科学工程计算与matlab编程3.pptVIP

3.5.2.2 第二类曲面积分 又称对坐标的曲面积分 可转化成第一类曲面积分 若曲面由参数方程给出 例： 的上半部，且积分沿椭球面的上面。 ％引入参数方程 x=a*sin(u)*cos(v); y=b*sin(u)*sin(v); z=c*cos(u)， u[0,pi/2], v[0,2*pi]. syms u v; syms a b c positive; x=a*sin(u)*cos(v); y=b*sin(u)*sin(v); z=c*cos(u); A=diff(y,u)*diff(z,v)-diff(z,u)*diff(y,v); I=int(int(x^3*A,u,0,pi/2),v,0,2*pi) I = 2/5*pi*a^3*c*b Sum((2*y(1:end-1,:)+diff(y)).*diff(x))/2 格式： S=trapz(x,y) 例： x1=[0:pi/30:pi]; y=[sin(x1) cos(x1) sin(x1/2)]; x=[x1 x1 x1]; S=sum((2*y(1:end-1,:)+diff(y)).*diff(x))/2 S = 1.9982 0.0000 1.9995 S1=trapz(x1,y) % 得出和上述完全一致的结果 S1 = 1.9982 0.0000 1.9995 例： 画图 x=[0:0.01:3*pi/2, 3*pi/2]; % 这样赋值能确保 3*pi/2点被包含在内 y=cos(15*x); plot(x,y) % 求取理论值 syms x, A=int(cos… (15*x),0,3*pi/2) A = 1/15 随着步距h的减小，计算精度逐渐增加： h0=[0.1,0.01,0.001,0.0001,0.00001,0.000001]; v=[]; %v定义为一个数组 for h=h0, x=[0:h:3*pi/2, 3*pi/2]; y=cos(15*x); I=trapz(x,y); v=[v; h, I, 1/15-I ];end %1/15-I表示用解析法与用trapz数值法的误差 format long; v v = 0.10000000000000 0.05389175150076 0.01277491516591 0.01000000000000 0.06654169546584 0.00012497120083 0.00100000000000 0.06666541668004 0.00000124998663 0.00010000000000 0.06666665416667 0.00000001250000 0.00001000000000 0.06666666654167 0.00000000012500 0.00000100000000 0.06666666666542 0.00000000000125 3.4.2 单变量数值积分问题求解 梯形公式 格式：(变步长）(Fun:函数的字符串变量） y=quad(Fun,a,b) y=quadl(Fun,a,b) % 求定积分 y=quad(Fun,a,b, ) y=quadl(Fun,a,b, ) %限定精度的定积分求解，默认精度为10－6。后面函数算法更精确，精度更高。 例： 第三种：匿名函数(MATLAB 7.0) 第二种：inline 函数 第一种，一般函数方法 函数定义被积函数： y=quad(‘c3ffun’,0,1.5)%将函数名用‘’括起来。 y = 0.9661 用 inline 函数定义被积函数： f=inline(2/sqrt(pi)*exp(-x.^2),x); y=quad(f,0,1.5) y = 0.9661 运用符号工具箱： syms x, y0=vpa(int(2/sqrt(pi)*exp(-x^2),0,1.5),60) y0 = .966105146475310713936933729949905794996224943257461473285749 y=quad(f,0,1.5,1e-20) % 设置高精度，但该方法失效，用quad不能设置精度，但是用quadl的话可以设置精度。 提高求解精度： y=quadl(f,0,1.5,1e-20) y = 0.9661 abs(y-y0) ans = .6402