计算化学公社

标题: 一个matlab程序运行不出来，希望高手看看 [打印本页]

作者
Author: zyj19831206 时间: 2017-8-31 13:52
标题: 一个matlab程序运行不出来，希望高手看看
function csegout
%微观离析反应器计算
cmean1 = quad(@fun1,0,10)
cmean2 = quad(@fun2,10,14)
cmean3 = quad(@fun3,14,24)
cmean = cmean1+cmean2+cmean3

function y=fun1(t)
syms t
k=2.4e-3;
cA0=1600;
y=cA0.*(-0.000412308.*t.^3+0.00538945.*t.^2-0.00246543.*t-...
0.000100973)/(1+k*cA0.*t);

function y=fun2(t)
syms t
k=2.4e-3;
cA0=1600;
y=cA0.*(-0.019084.*t+0.292621)/(1+k*cA0.*t);

      function y=fun3(t)
      syms t
      k=2.4e-3;
      cA0=1600;
      y=cA0.*(0.0000103077.*t.^3-0.000479624.*t.^2+0.00425519.*t+...
      0.0315744)/(1+k*cA0.*t);

作者
Author: kyuu 时间: 2017-8-31 19:47
本帖最后由 kyuu 于 2017-8-31 22:09 编辑

cmean1 = quad(@fun1,0,10)
cmean2 = quad(@fun2,10,14)
cmean3 = quad(@fun3,14,24)
cmean = cmean1+cmean2+cmean3

function y=fun1(t)
k=2.4e-3;
cA0=1600;
y=cA0.*(-0.000412308.*t.^3+0.00538945.*t.^2-0.00246543.*t-...
0.000100973)./(1+k*cA0.*t);
end

function y=fun2(t)
k=2.4e-3;
cA0=1600;
y=cA0.*(-0.019084.*t+0.292621)./(1+k*cA0.*t);
end

function y=fun3(t)
      k=2.4e-3;
      cA0=1600;
      y=cA0.*(0.0000103077.*t.^3-0.000479624.*t.^2+0.00425519.*t+...
      0.0315744)./(1+k*cA0.*t);
end

作者
Author: zyj19831206 时间: 2017-8-31 21:07

kyuu 发表于 2017-8-31 19:47
cmean1 = quad(@fun1,0,10)
cmean2 = quad(@fun2,10,14)
cmean3 = quad(@fun3,14,24)

还是不行啊，为什么？

作者
Author: kyuu 时间: 2017-8-31 21:44
本帖最后由 kyuu 于 2017-9-3 16:10 编辑

zyj19831206 发表于 2017-8-31 21:07
还是不行啊，为什么？

复制我的代码，别在你自己的上边改，你点乘写的不对，乘出来的向量维度输入与输出不一致
cmean1 =
42.4139

cmean2 =
9.0121

cmean3 =
2.6709

cmean =

54.0970

欢迎光临计算化学公社 (http://bbs.keinsci.com/)