精度和速度要比quad函数高，而trapz函数用的是梯形积分法，精度低，适用于数值函数 和光滑性不好的函数的积分计算. 例7用数值积分的方法计算∫xd 解先用trapz(仪，Y)来求x'dr,相应的Matlab程序为 >》x=-1:0.1:2:y=x.^3:%取积分步长为0.1 >》trapz(k,y) 结果为ans=3.7575 现在把步长取到0.01，则有 >x=-1:0.01:2:y=x.^3: >trapz(x,y) 结果为ans=3.7501 通过上面的运算结果可以看到，利用不同的步长进行计算，可以得到不同精度的近似值， 如果取步长为0.01，则输出结果为3.7500，与步长取到0.01时的计算值接近. 可以利用符号积分命令计算，相应的Mat1ab程序为 >>clear all >syms x: >>int(x3.X.-1.2) 结果为ans=15/4 下面用quad和quadl米计算x'dr，首先建立名为jifen.m的m文件。 function y=jifen(x) Y=x>3: 然后在命令窗口输入以下Mat1ab程序 >>quad('jifen',-1,2) 结果为ans=3.7500. 利用quadl函数同样可以得到上面结果。 例8计第颈值积分血，并用符号积分命令求解，及案输出结果 解利用rapz函数计算数值积分 r,相应的Matlab程序为 1+x >x=0:0.1:1:y=sin(x.2)./(1+x): >》trapz(x,y) 结果为an5=0.1811. 当步长变为0.01时运算结果为0.1808， 利用quad函数求解，先建立名为jifen1.m的m文件，程序为 function y=jifen(x) y=sin(x.^2）./(1+x): 在命令窗口进行调用，程序为 5 精度和速度要比 quad 函数高，而 trapz 函数用的是梯形积分法，精度低，适用于数值函数 和光滑性不好的函数的积分计算. 例 7 用数值积分的方法计算 2 3 1 x xd − . 解 先用 trapz(X,Y)来求 2 3 1 x xd − ，相应的 Matlab 程序为 >> x=-1:0.1:2；y=x.^3； %取积分步长为 0.1 >> trapz(x,y) 结果为 ans =3.7575 现在把步长取到 0.01，则有 >> x=-1:0.01:2；y=x.^3； >> trapz(x,y) 结果为 ans =3.7501 通过上面的运算结果可以看到，利用不同的步长进行计算，可以得到不同精度的近似值. 如果取步长为 0.01，则输出结果为 3.7500，与步长取到 0.01 时的计算值接近． 可以利用符号积分命令计算，相应的 Matlab 程序为 >> clear all >> syms x； >> int(x^3,x,-1,2) 结果为 ans =15/4 下面用 quad 和 quadl 来计算 2 3 1 x xd − ，首先建立名为 jifen.m 的 m 文件. function y=jifen(x) Y=x>^3； 然后在命令窗口输入以下 Matlab 程序 >> quad('jifen',-1,2) 结果为 ans =3.7500. 利用 quadl 函数同样可以得到上面结果. 例 8 计算数值积分 2 1 0 sin d 1 x x + x  ，并用符号积分命令求解，观察输出结果. 解 利用 trapz 函数计算数值积分 2 1 0 sin d 1 x x + x  ，相应的 Matlab 程序为 >> x=0:0.1:1； y=sin(x.^2)./(1+x)； >> trapz(x,y) 结果为 ans =0.1811. 当步长变为 0.01 时运算结果为 0.1808. 利用 quad 函数求解，先建立名为 jifen1.m 的 m 文件，程序为 function y=jifen(x) y=sin(x.^2)./(1+x)； 在命令窗口进行调用，程序为