登录论坛

查看完整版本 : 求助!两种方法算出来不一致


未注册
2008-07-29, 13:50
第一种方法:
global D1 D2 T;
t=1;
spaceT=0.05;
for T=0.001:spaceT:1.2
....(中间是关于T的函数的构造:Z(T))
t=t+1;
a(t)=T;
y(t)=Z;
end
plot(a,y,'k');%画出Z(T)与T的图。

第二种方法:
global D1 D2 Tt;
sym T;%把T当所符号
t=1;
....(中间是关于T的函数的构造:Z(T)) 和第一种方法中间的一模一样
spaceT=0.05;
for Tt=0.001:spaceT:1.2
t=t+1;
a(t)=Tt;
y(t)=subs(Z,T,Tt); %用Tt来代符号T
end
plot(a,y,'r');

算出来的结果竟然不一样,第一种方法的数据是对的,Z的值从0开始,从小到大,最后为0.3647;但为什么第二种就不对了呢,是从1开始逐渐减小,到最后一个点的值一样为0.3647。