请教关于积分后求和的问题

江山美人

编了下面一个小程序，当Ns=7之内的时候都可以运行，当N=8以上时就出现如下错误：错误使用 mupadmex
Error in MuPAD command: The number of arguments is incorrect. [has]
  Evaluating: sum


出错 sym/symsum (line 114)
   rSym = mupadmex('symobj::map',fsym.s,'symobj::symsum',x.s,a.s,b.s);
请问各位大神是什么原因，谢谢。

1. syms gamma k
   
2. Ns=4;
   
3. Bs=0.05;
   
4. R=0.2;
   
5. H=0.002;
   
6. n=3;
   
7. alpha=2*pi/Ns;
   
8. beta=Bs/(R+H/2);
   
9. fa=(cos(n*gamma))^2;
   
10. f11=int(fa,gamma,k*alpha,k*alpha+beta);
    
11. f11=vpa(f11);
    
12. T11=symsum(f11,k,0,Ns-1);
    
13. z=vpa(T11)

复制代码

chybeyond

我把你问题简化发到中文论坛，可能是一个bug。参考帖子：http://www.ilovematlab.cn/thread-297092-1-2.html

江山美人

chybeyond 发表于 2014-6-21 09:44
我把你问题简化发到中文论坛，可能是一个bug。参考帖子：http://www.ilovematlab.cn/thread-297092-1-2.htm ...

谢谢，中文论坛我看了

江山美人

chybeyond 发表于 2014-6-21 09:44
我把你问题简化发到中文论坛，可能是一个bug。参考帖子：http://www.ilovematlab.cn/thread-297092-1-2.htm ...

这个问题怎么报错啊，没有做过

江山美人

chybeyond 发表于 2014-6-21 09:44
我把你问题简化发到中文论坛，可能是一个bug。参考帖子：http://www.ilovematlab.cn/thread-297092-1-2.htm ...

打错了， 是怎么报告啊

chybeyond

江山美人 发表于 2014-6-22 09:19
打错了， 是怎么报告啊

不太清楚，去官网看看

江山美人

chybeyond 发表于 2014-6-22 12:09
不太清楚，去官网看看

我在7.4.0版本里面，算出来不会报错，可能高版本的还不够完善。

江山美人

这个问题突然想到了一个方法，共享一下，把求和分成两部分

1. syms gamma k
   
2. Ns=8;
   
3. Bs=0.05;
   
4. R=0.2;
   
5. H=0.002;
   
6. n=3;
   
7. alpha=2*pi/Ns;
   
8. beta=Bs/(R+H/2);
   
9. fa=(cos(n*gamma))^2;
   
10. f11=int(fa,gamma,k*alpha,k*alpha+beta);
    
11. f11=vpa(f11);
    
12. T1=symsum(f11,k,0,Ns-5);
    
13. T2=symsum(f11,k,Ns-4,Ns-1);
    
14. T11=T1+T2;
    
15. z=vpa(T11)

复制代码

就有结果了