Re: [問題]在計算中有關多次數列疊代的問題!!!!

[TOTOROBOY (小胖)]

把D1 改成D1(i)
或者是直接改成下面程式碼


clear
  clc
  Ea=0.45;
  D0=7*10^(-6);
  k=8.617*10^(-5);
  i=0:1:5400 ;
  y1=873-i*(873/5400)
  D=D0*exp(-Ea./(k*y1)); % Diffusion coefficient
  D1=sqrt(2*D.*i); % Diffusion length

  plot(i,D1,'color',[0 0 1]);


※ 引述《johnny4003 ()》之銘言：
:   目前在計算有關擴散的問題，由於時間係數i為0~5400，
:   代入y1(i)的值再代入D(y1)，而所得的D(y1)需要再乘上
:   i，才能得到結果的D1。
:   例如i=10
:   先計算出y1(10)，將y1(10)代入D方程式得D(y1(10))，然
:   後再將其D(y1(10))*10開根號求結果。
:   請問下列方程式哪邊有問題?謝謝~~
:   clear
:   clc
:   Ea=0.45;
:   D0=7*10^(-6);
:   k=8.617*10^(-5);
:   for i=0:1:5400 ; //時間
:   y1=873-i*(873/5400)
:   D=D0*exp(-Ea/(k*y1)); // Diffusion coefficient
:   D1=sqrt(2*D*i); // Diffusion length
:   end
:   plot(i,D1,'color',[001]);

--
※ 發信站: 批踢踢實業坊(ptt.cc) 
◆ From: 1.162.16.86
推 johnny4003:感謝感謝~~