请教大家 matlab 积分问题

bettyshen

此式为要积分的公式。其中f是表达式很长的函数。我用int 积分，结果总是不对。
如下的错误。不知怎么改呢？
??? Error using ==> reshape
To RESHAPE the number of elements must not change.

请大家指教。谢谢

ChaChing

Ref : 常见的程序出错问题整理 http://forum.vibunion.com/thread-46001-1-1.html, 14F
From http://forum.vibunion.com/home-space-uid-63979-do-blog-id-18250.html

qibbxxt

建议你用quadl等函数进行数值积分

bettyshen

谢谢大家的回答，我用quadl积分 出现下面的提示。Warning: System is inconsistent. Solution does not exist  请问怎么解决呢

bettyshen

。    。

bettyshen

回复 3 # qibbxxt 的帖子

1. clear
   
2. clc
   
3. 
   
4. w0=1.0;w1=2.0;w2=2.0;alpha=1.0;beta1=0.01;beta2=0.02;D1=0.01;D2=0.01;x1=0.3;x2=0.3;
   
5. 
   
6. g=w0^2.*a+alpha.*a.^3;
   
7. n=(alpha.*a.^2/4)./(w0^2+3*alpha*a.^2/4);
   
8. b0=(w0^2+3*alpha*a.^2/4).^0.5;
   
9. b2=(w0^2+3*alpha*a.^2/4).^0.5*n/2;
   
10. 
    
11. S12=D1/pi()./(((2*b0).^2-w1^2).^2+4*x1^2*w1^2.*(2*b0).^2); %%S1(2*b0)
    
12. S14=D1/pi()./(((4*b0).^2-w1^2).^2+4*x1^2*w1^2.*(4*b0).^2); %%S1(4*b0)
    
13. 
    
14. S21=D2/pi()./((b0.^2-w2^2).^2+4*x2^2*w1^2*b0.^2);          %%S2(b0)
    
15. S23=D2/pi()./(((3*b0).^2-w2^2).^2+4*x2^2*w2^2*(3*b0).^2);  %%S2(3*b0)
    
16. 
    
17. g1=w0^2.*u+alpha*u.^3;
    
18. n1=(alpha*u.^2/4)/(w0^2.*u+3*alpha*u.^2/4);
    
19. b01=(w0^2+3*alpha*u.^2/4).^0.5;
    
20. b21=(w0^2+3*alpha*u.^2/4).^0.5.*n1/2;
    
21. 
    
22. SS12=D1/pi()./(((2*b01).^2-w1^2).^2+4*x1^2*w1^2*(2*b01).^2); %%S1(2*b0)
    
23. SS14=D1/pi()./(((4*b01).^2-w1^2).^2+4*x1^2*w1^2*(4*b01).^2); %%S1(4*b0)
    
24. SS21=D2/pi()./((b01.^2-w2^2).^2+4*x2^2*w1^2*b01.^2);          %%S2(b0)
    
25. SS23=D2/pi()./(((3*b01).^2-w2^2).^2+4*x2^2*w2^2*(3*b01).^2);  %%S2(3*b0)
    
26. 
    
27. 
    
28. m1=-u.^2./8./g1.*(-beta1.*(4*w0^2+5*alpha.*u.^2./2)+beta2*u.^2.*(w0^2+3*alpha*u.^2/4))+...
    
29.    pi()*u.^2./32./g1.*((2*b01).*SS12.*(1./(u+3/4.*u.^2).^(1/2).*u.^2./(u+u.^3).*(1+3/2*u)+4*(u+3/4*u.^2).^(1/2).*u./(u+u.^3)-2*(u+3/4*u.^2).^(1/2).*u.^2./(u+u.^3).^2.*(1+3*u.^2))+...
    
30.      b21.*SS14.*(-1/16./(u+3/4*u.^2).^(3/2).*u.^2.*(1+3/2*u)+1/4./(u+3/4*u.^2).^(1/2).*u))+...
    
31.    pi()*u.^3./16./g1.^2.*((2*b01).*(2*b01+2*b21).*SS12+2*(b21).*b21.*SS14)+...
    
32.     pi()*u./8./g1.*((2*b01-b21).*SS21.*((1/2./(u+3/4*u.^2).^(1/2).*(1+3/2*u)+1/16./(u+3/4*u.^2).^(3/2).*u.^2.*(1+3/2*u)-1/4./(u+3/4*u.^2).^(1/2).*u).*u./(u+u.^3)+...
    
33.         ((u+3/4*u.^2).^(1/2)-1/8./(u+3/4.*u.^2).^(1/2).*u.^2)./(u+u.^3)-((u+3/4*u.^2).^(1/2)-1/8./(u+3/4.*u.^2).^(1/2).*u.^2).*u./(u+u.^3).^2.*(1+3*u.^2))+...
    
34.         (2*b21).*SS23.*(-1/8./(u+3/4*u.^2).^(3/2).*u.^3./(u+u.^3).*(1+3/2.*u)+3/4./(u+3/4.*u.^2).^(1/2).*u.^2./(u+u.^3)-1/4./(u+3/4*u.^2).^(1/2).*u.^3./(u+u.^3).^2.*(1+3*u.^2)))+...
    
35.     pi()*u./8./(g1).^2.*((2*b01-b21).*(2*b01+b21).*SS21+3*((b21).^2).*SS23);
    
36. 
    
37. B=@(a)pi()/16*a.^4./g.^2.*(4*b0.^2.*S12+b2.^2.*S14)+pi()/14*a.^4./g.^2.*(4*b0.^2.*S21+b2.^2.*S23);
    
38. BB=@(a)1./B;  %被积函数1的推导
    
39. 
    
40. B1=pi()*u.^4/16./g1.^2.*(4*b01.^2.*SS12+b21.^2.*SS14)+pi()*u.^4/14./g1.^2.*(4*b01.^2.*SS21+b21.^2.*SS23);
    
41. 
    
42. y=@(u)2*m1./B1;         %被积函数2的推导
    
43. 
    
44. q=quadl(BB.*exp(quadl(y,u,0,a)),a,0,inf);

复制代码

程序不能计算，不知道问题出在哪里，请帮忙指教，谢谢