符号运算太复杂导致matlab出错

beyondhxf

我用了一个50*50次的循环，每次循环计算一个很复杂的积分，到中途出错了
以前做的18*18时候，没有出错，说明不是程序的错误。

1. 
   
2. ??? Error using ==> reshape
   
3. To RESHAPE the number of elements must not change.
   
4. Error in ==> sym.maple at 94
   
5.       result = reshape(result,size(varargin{3}));
   
6. Error in ==> sym.int at 51
   
7.    r = reshape(maple('map','int',f(:),[x.s '=(' a.s ')..(' b.s ')']),size(f));
   

复制代码

请问如何改进我的程序，以便让它能运行通？

补充？难道说在调用maple做符号运算的时候，会积累什么错误之类的？调用次数太多就出错？maple就崩溃？

[ 本帖最后由 beyondhxf 于 2009-10-26 13:15 编辑 ]

ChaChing

符号运算来此学的, 没这方面经验!
仅有出错提示无出错代码, 真没办法, 待高人路过

friendchj

result = reshape(result,size(varargin{3}));根据提示，应该是result和size(varargin{3})的元素大小不一样。最好能贴出程序代码，可以简化贴出，能说明问题即可。

beyondhxf

原帖由 beyondhxf 于 2009-10-26 13:11 发表

                               
登录/注册后可看大图

我用了一个50*50次的循环，每次循环计算一个很复杂的积分，到中途出错了
以前做的18*18时候，没有出错，说明不是程序的错误。
??? Error using ==> reshape
To RESHAPE the number of elements must not change. ...

贴出相关的程序，仅是说明性的，不全面。

1. 
   
2. syms sita
   
3. for i=1:50
   
4.      A(1,i)=B(sita)
   
5.      。。。
   
6.      A(4,i)=B(sita)
   
7. end
   
8. %这里，B函数是一个关于sita的函数，返回一个关于sita的表达式(比较复杂，恐怕是致错的根源)，是sym型的。
   
9. for i=1:50
   
10. for j=1:50
    
11. C(i,j)=int(A(1,i)*A(2,j),-1,1);
    
12. end
    
13. end
    
14. %由于循环次数是50*50，比较多，在算到一半的时候会出现如上提示的错误。而如果循环次数在18*18次，则不出错。
    
15. 本人猜测，这是由于调用maple内核过多而导致的错误，并不是程序本身的语法错误。
    

复制代码

friendchj

参考一下：http://forum.vibunion.com/forum/ ... mp;page=1#pid441572