MATLAB求一元函数的数值积分(quad,quadl,quadv)

时间:2024-10-16 14:12:37

1、第一,计算下图两个一元函数的数值积分。

MATLAB求一元函数的数值积分(quad,quadl,quadv)

2、第二,启动MATLAB,新建脚本(Ctrl+N),输入如下代码:close all; clear all; 艘早祓胂clcformat compactfigure('Position',[50,50,600,500],'Name','sin(0.5x)')f1=@(x) sin(0.5*x);[Q1a,Fcnt1a]=quad(f1,0,2)[Q1b,Fcnt1b]=quadl(f1,0,2)fplot(f1,[0,2]);%-------------------%figure('Position',[100,100,600,500],'Name','e^x+1')f2=@(x) exp(x)+1;[Q2a,Fcnt2a]=quad(f2,-1,2)[Q2b,Fcnt2b]=quadl(f2,-1,2)fplot(f2,[-1,2]);其中f1=@(x) sin(0.5*x),f2=@(x) exp(x)+1是建立被积分的函数(形式为函数句柄/匿名函数)。[Q,Fcnt]=quad(function,a,b),quadl(function,a,b)是分别采用quad和quadl两种方法求数值积分。其中被积函数function要为上述建立的函数句柄/匿名函数,a,b分别积分上限,[Q,Fcnt]分别返回数值积分的结果和函数计算的次数。

MATLAB求一元函数的数值积分(quad,quadl,quadv)

3、第三,保存和运行上述脚本,在命令行窗口(Command Window)得到如下结果:Q1a = 0.9194Fcnt1a = 13Q1b = 0.9194Fcnt1b = 18Q2a = 10.0212Fcnt2a = 37Q2b = 10.0212Fcnt2b = 18可以看出,采用quad( )和quadl两种方法计算函数(1)的数值积分都是0.9194,函数(2)的数值积分都是10.0212,只是计算次数不同。例如采用quad( )计算函数(1)的数值积分,计算次数为13,采用quadl( )计算次数为18。

MATLAB求一元函数的数值积分(quad,quadl,quadv)

4、第四,下图为被积分函数(1)的图像。

MATLAB求一元函数的数值积分(quad,quadl,quadv)

5、第五,下图为被积分函数(2)的图像,其中fplot( )绘图命令可以绘制函数句柄/匿名函数的图像。

MATLAB求一元函数的数值积分(quad,quadl,quadv)
© 手抄报圈