累加!怎么会出错呢？

cumtsafe303

clear all
syms n k1 k2 k3
t=0:2:140
k1=1.5e-3
k2=3e-3
k3=6.58e-3
r1=symsum(n^(-2)*exp(-k1*t*(n^2)),n,1,inf)
r2=symsum(n^(-2)*exp(-k2*t*(n^2)),n,1,inf)
r3=symsum(n^(-2)*exp(-k3*t*(n^2)),n,1,inf)
r2=double(r2)
r3=double(r3)
r1=double(r1)
y2=1-6*pi^(-2)*r2
y1=1-6*pi^(-2)*r1
y3=1-6*pi^(-2)*r3
plot(t,y1,t,y2,t,y3)

feynmand

你的程序在我的电脑上可以正常运行，给出了结果，不知道你哪里出错了？
我用的版本是2009a

cumtsafe303

file:///d:/

cumtsafe303

??? Error: Assignment statements do not produce results. (Use == to test
for equality.)

Error in ==> D:\matlab\Matlab6p5FULL\toolbox\symbolic\@sym\double.m
On line 45  ==> D = reshape(eval(X),m,n);

Error in ==> D:\matlab\Matlab6p5FULL\work\fa.m
On line 12  ==> r1=double(r1)

cumtsafe303

我的是6.5的。

cumtsafe303

楼主,你解答的不对,我用09的版本试,还是出错.你能解释清吗?

cumtsafe303

??? Error using ==> mupadinit at 35
Unable to initialize kernel

Error in ==> mupadengine.mupadengine>mupadengine.evalin at 117
                [res,status] = mupadmex(statement,8);

Error in ==> C:\Program Files\MATLAB\R2009a\toolbox\symbolic\symengine.p>symengine at 29


Error in ==> sym.sym>sym.sym at 94
            eng = symengine;

Error in ==> syms at 77
      assignin('caller',x,sym(x));

cumtsafe303

上面是2009a的结果?

bainhome

楼主,你解答的不对,我用09的版本试,还是出错.你能解释清吗?

问题没解决之前，不要轻易下结论说别人对错。这应该是MATLAB新弄进来的“超级高档”符号计算引擎mupad造成的，如果用2008a以下的所有版本(maple引擎)就不会出现这种问题，目前此问题原因不明，但很可能是由于防火墙或者杀毒软件等block造成的。
可以关闭防火墙、安全卫士的保险箱、杀毒软件的所有监控，然后再试试，我也用2009a，可以出结果，不过我的方法过于狠毒——将MATLAB装在d盘，然后重新做系统，现在MATLAB算“绿色版”，好像可以避开无法开启mupad引擎的问题。

jinzi2006

算出的结果！

jinzi2006

我的是matlab R2009a啊!

zheng720

我觉得直接把r1=symsum(n^(-2)*exp(-k1*t*(n^2)),n,1,inf)
r2=symsum(n^(-2)*exp(-k2*t*(n^2)),n,1,inf)
r3=symsum(n^(-2)*exp(-k3*t*(n^2)),n,1,inf)
r2=double(r2)
r3=double(r3)
r1=double(r1)
该为
r1=subs(n^(-2)*exp(-k1*t*(n^2)),n,1,inf)
r2=subs(n^(-2)*exp(-k2*t*(n^2)),n,1,inf)
r3=subs(n^(-2)*exp(-k3*t*(n^2)),n,1,inf)
下面的直接画图就可以了，应该很容易做的，不知道我解的对不对？
画出来结果和10楼的一样。

cumtsafe303

谢谢，大家。算出来了。

ljelly

谢谢，大家。算出来了。
cumtsafe303 发表于 2009-5-8 14:03

分享一下吧，谢谢！

changshanyan

看来还需好好努力

changshanyan

看来还需好好努力