标签：数值积分 end 积分 Isum nargin Matlab dx 无穷

Matlab无穷上限数值积分的近似方法

经常遇到从0积分到无穷的数值积分问题

首先要求被积函数有无穷限积分

数值积分可直接采用matlab 的 integral

有时候积分过程计算很复杂，可以取0到M（M是较大的数）的积分值代替无穷限积分

代码如下

function Isum = quadToInf(fun,a,dx0,tol,method)
if nargin < 2 ,a=0 ;end
if nargin < 3 ,dx0=0.5 ;end
if nargin < 4 ,tol = 5e-4 ;end
if nargin < 5 ,method = 1 ;end
j=0;dx = dx0;Isum = 0;x1 = a; maxint = 35;
%fprintf('\n j dx x2 I_j Isum\n');
while j<maxint
    x2 = x1 + dx;
    switch method
    case 1 
        I = integral(fun,x1,x2);
    case 2 
        I = quad(fun,x1,x2);
    otherwise 
        fprintf('method = %d not allowed',method);
    end
    Isum = Isum + I;
    %fprintf('%4d %8.1f %8.1f %12.8f %12.8f\n',j,dx,x2,I,Isum);
    if j>5 && abs(I/Isum) < tol
        break; 
    end
j = j+1;x1 = x2;dx = 2*dx;
end
end

标签：数值积分,end,积分,Isum,nargin,Matlab,dx,无穷
来源： https://www.cnblogs.com/elapsetor/p/15642510.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。