我这个方程的那里有问题？

hitlgq 发布于 2007-04-18 16:51， 1513 次点击

帮我看看吧！高手些，谢谢
function F=fcht(t)
Ka=0.0028;
Q=47595;
p=0.1;
peq=2.5;
T=250;
R=8.314;
t=0:100;
F=1-exp(-4*1*(Ka*exp(-Q/R/T)*(p-peq))^(1/2)*t^(1/2));
plot(t,F,'r-')
xlabel('t');
ylabel('F');
grid on ;

10 回复

尘埃1222007-04-18 18:35

F=1-exp(-4*1*(Ka*exp(-Q/R/T)*(p-peq))^(1/2)*t^(1/2));
改为：
F=1-exp(-4*1*(Ka*exp(-Q/R/T)*(p-peq))^(1/2).*t.^(1/2));

程序能通过，但结果中有复数，你可以看看方程是否描述正确

当有矩阵运算时，最好用点乘，点除的运算符号，这样可以避免出错

hitlgq2007-04-18 20:18

但是我哪个方程没有数组和矩阵啊，点乘不行啊。最后的结果应该是一条曲线F=1-exp（-0.1976*t^0.5)

hitlgq2007-04-18 20:21

2楼能不能把你的邮箱给我，我把方程的推导过程给，你帮我看看好吗？

hitlgq2007-04-18 20:21

回复：（hitlgq）我这个方程的那里有问题？

2楼能不能把你的邮箱给我，我把方程的推导过程给，你帮我看看好吗？

尘埃1222007-04-19 08:38

呵呵，你把推导过程传上来，将程序意思大概描述下，问题应该不大的

邮箱个人信息里面有

hitlgq2007-04-19 18:52

回复：（尘埃122）呵呵，你把推导过程传上来，将程序...

只有本站会员才能查看附件，请登录

非常地感谢你，这上面有向您这样热心的确实不多，方程的推导在
附件里，里面有方程的说明，谢谢。。。

hitlgq2007-04-19 18:53

回复：（尘埃122）呵呵，你把推导过程传上来，将程序...

其中n=1

hitlgq2007-04-19 18:55

回复：（尘埃122）呵呵，你把推导过程传上来，将程序...

边界值F（0）=0

#10

尘埃1222007-04-20 09:13

呵呵，你确定方程推导对了？

我有些看迷糊了

#11

hitlgq2007-04-20 13:40

回复：（尘埃122）呵呵，你确定方程推导对了？我有些看...

应该没有问题，我做毕业设计，这个方程是导师给我的