数值积分与matlab求解程序.pptVIP

Matlab多重积分 （二）用MATLAB 数值计算三重积分 调用格式一：Q3=triplequad(FUN,a,b,c,d,p,q) 调用格式二：Q3=triplequad(FUN,a,b,c,d,p,q,tol) 调用格式三：Q3=triplequad(FUN,a,b,c,d,p,q,tol,@QUADL) 调用格式四： Q3=triplequad(FUN,a,b,c,d,p,q,tol,@MYQUADF) 调用格式五： Q3=triplequad(FUN,a,b,c,d,p,q,tol,@QUADL,P1,P2,...) 调用格式六：Q3=triplequad(FUN,a,b,c,d,p,q,[ ]，[ ],P1,P2,...) 与 Q3=triplequad(FUN,a,b,c,d,p,q,1.e-6,@QUAD,P1,P2,...) 相同. 例11 分别用MATLAB 函数triplequad的调用格式二和三计算 ? ??? ? 的值，取误差限为tol= 10-4 ，并将计算结果与精确值比较.其中V是三维长方体区域2≤x≤2,2≤y≤2,0≤z≤4. 解： 建立并保存被积函数的M 文件 function u=integrnd1(x,y,z) u=x+exp(y)+sin(z); 在MATLAB 工作窗口输入程序 a= -2;b=2;c=-2;d=2;p=0;q=4; Q3=triplequad(@integrnd1,a,b,c,d,p,q,1.e-4) QL3=triplequad(@integrnd1,a,b,c,d,p,q,1.e-4,@quadl) syms x y z f=x+exp(y)+sin(z); jfz=int(f,z,p,q); jfy=int(jfz,y,c,d); jfx=int(jfy,x,a,b); I3=double(jfx), Juewu3=abs(I3-Q3) JuewuL3=abs(I3-QL3) 运行后屏幕显示如下 Q3 = 1.425178451284647e+002, QL3 = 1.425178451284647e+002 I3 = 1.425178309849224e+002, Juewu3 = 1.414354233020276e-005 JuewuL3 = 1.414354233020276e-005 Thank You！ MATLAB求解数值积分 MATLAB求解连续函数积分 引言 我们知道,若函数f(x)在区间[a,b]上连续且其原函数为F(x),则可用Newton-Leibnitz公式 求定积分的值 , Newton-Leibnitz公式无论在理论上还是在解决实际问题上都起了很大作用，但它并不能完全解决定积分的计算问题，因为积分学涉及的实际问题极为广泛，而且极其复杂，在实际计算中经常遇到以下三种情况： 1 Matlab求解连续函数积分 ?(1) 被积函数f(x)并不一定能够找到用初等函数的 有限形式表示的原函数F(x)，例如： Newton-Leibnitz公式就无能为力了 (2) 还有被积函数f(x)的原函数能用初等函数表示，但表达式太复杂，例如函数 并不复杂，但积分后其表达式却很复杂，积分后其原函数F(x)为： (3) 被积函数f(x)没有具体的解析表达式, 其函数 关系由表格或图形表示。 对于这些情况, 要计算积分的准确值都是十分困难的。由此可见, 通过原函数来计算积分有它的局限性, 因而研究一种新的积分方法来解决Newton-Leibniz公式所不能或很难解决的积分问题, 这时需要用数值解法来建立积分的近似计算方法。 将积分区间细分,在每一个小区间内用简单函数代替复杂函数进行积分，这就是数值积分的思想，用代数插值多项式去代替被积函数发f(x)进行积分是本章讨论数值积分的主要内容。 1.1 定积分的Matlab符号计算 例1 由 y =sin x， y =cos x , x=-1/2,x=3/2所围成的平面区域D.求平面区域D 的面积S. 解 输入作函数图形的程序 x=-1:0.001:2; F1= sin(x); F2=cos(x); plot(x ,F1,b-,x ,F2,g-), axis([-1,pi/4+1,-1.3,1.3]), xlabel(x), ylabel(y), title(y=sinx , y=cosx 和x=-0.5及x=1.5所围成的平面区域的图形) 运行后屏幕显示图形. 求平面区域D 的面积S.输入计算面积S 的程序 syms x f1= cos(x)-sin(x);