作者johnny4003 ()
看板MATLAB
标题[问题]在计算中有关多次数列叠代的问题!!!!
时间Fri Aug 24 01:02:48 2012
目前在计算有关扩散的问题,由於时间系数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: 180.177.10.89