MATLAB计算二重数值积分(dblquad)？
第一，使用MATLAB计算下图中的二重积分。
第二，启动MATLAB，新建脚本（Ctrl+N），输入如下代码：close all; clear all; clcfigure(‘Position’,[50,50,600,500],’Name’,’double integration’)f=@(x,y) x*sin(y)-cos(x)+y-3;Q1=dblquad(f,0,2*pi,-pi,pi,1.0e-3)x=linspace(0,2*pi);y=linspace(-pi,pi);ff=x.*sin(y)-cos(x)+y-3;plot3(x,y,ff,’r’,’LineWidth’,3)grid on; view(35,35)xlabel(‘x’,’FontSize’,15); ylabel(‘y’,’FontSize’,15)其中Q1=dblquad(f,0,2*pi,-pi,pi,1.0e-3)采用默认方法quad计算二重积分，绝对计算精度设为1.0e-3。plot3(x,y,ff,’r’,’LineWidth’,3)是绘制被积分函数ff=x.*sin(y)-cos(x)+y-3的图像。
第三，保存和运行上述脚本，在命令行窗口（Command Window）得到如下结果：Q1 = -118.4351也就是说，该二重数值积分的结果为-118.4351.
第四，同时得到被积分函数x.*sin(y)-cos(x)+y-3的图像。
第五，在Q1=dblquad(f,0,2*pi,-pi,pi,1.0e-3)后面再添加一行命令，Q2=dblquad(f,0,2*pi,-pi,pi,1.0e-3,’quadl’)。也就是在绝对计算精度都是1.0e-3的情况下，使用quadl方法计算二重数值积分。
第六，保存和运行上述脚本，在命令行窗口（Command Window）得到采用quad方法（默认）和quadl方法计算的二重数值积分结果。Q1 = -118.4351Q2 = -118.4355也就是说，采用采用quad方法（默认）计算的二重数值积分结果为-118.4351，而采用quadl方法计算的二重数值积分结果为-118.4355，两者在小数点第四位略有差别。
延伸阅读
matlab怎么对矩阵数据积分？
1.
第一步,双击matlab软件图标,打开matlab软件,可以看到matlab软件的界面。
2.
第二步,使用syms命令,创建七个符号变量a、b、c、d、x、y、z。
3.
第三步,使用符号变量d、x,创建2×2符号矩阵A,其中A=[d^2,11*x;13*d*sin(x),21*x*sin(d)]。
4.
第四步,使用函数int(A,’x’),求解符号矩阵A中每个元素关于变量x的积分。
matlab中trapz(A,3)是什么意思？
trapz是梯形法数值积分 trapz(A,3)是对A的第三维数值积分 A原来是一个三维矩阵 大小为mxnxk B=trapz(A,3)返回一个mxn的二维矩阵 每个元素B(i,j)相当于原来的A(i,j,:)的积分 trapz(A,3)类似的于sum(A,3),都是将A的第三维的数据累加积分 只是trapz(A,3)用的是梯形面积，sum(A,3)用的是矩形面积
