一个matlab程序运行不出来，希望高手看看

zyj19831206

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);

kyuu

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

zyj19831206

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

还是不行啊，为什么？

kyuu

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

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

cmean2 =
9.0121


cmean3 =
2.6709

cmean =

54.0970