求助，积分运算

ibeing

1. u1=1.136;  sigmal1=log(1.066); fN1=0.025;
   
2. u2=1.461; sigmal2=log(1.056);  fN2=0.283;
   
3. nf1=fN1/(fN1+fN2); nf2=fN2/(fN1+fN2);
   
4. dx=0.01;x=(0.6:dx:2);
   
5. k=(x.^3-1).*(exp(0.0022./(200.*x))./0.85-1);
   
6. cl1(1,:)=exp(-(log(x./u1)).^2./(2*sigmal1^2))./(x.*sigmal1*sqrt(2*pi)).*nf1;
   
7. cl2(1,:)=exp(-(log(x./u2)).^2./(2*sigmal2^2))./(x.*sigmal2*sqrt(2*pi)).*nf2;
   
8. gfpdfl=cl1+cl2;
   
9. hold on ;
   
10. plot(k,gfpdfl,'b');
    
11. hold off
    
12. axis ([-0.2 1.2  0 6]);
    
13. xlabel('Kappa')
    
14. syms x
    
15. gm=int(gfpdfl,x,0,1)
    
16. a=inf;b=-inf;
    
17. syms k;
    
18. z=int(gfpdfl,k,a,b);

复制代码

上面程序是计算并画了一个对数正态分布谱图，然后希望用计算的对数正态分布函数分别计算积分值gm和z，gm表示的是该对数正态分布曲线的积分平均，理论上应该是一个数值，为什么结果出来是一个很大的数组？z也没有结果。像这种积分程序到底该怎么写？请大家指教。

迷鹿

不能给楼主帮助 小小支持一下

ibeing

迷鹿 发表于 2014-1-6 16:02

                               
登录/注册后可看大图

不能给楼主帮助 小小支持一下

谢谢，问题解决了