matlab第五章 函数、导数和积分等.pptVIP

第五章 函数、导数与积分等 的数值计算;1.函数值的计算;解： x=1:2:5;y=2:2:6; z=log(x.^2+y.^2)+2.*x.*y.*tan(y./x); x,y,z x = 1 3 5 y = 2 4 6 z = -7.1307 102.3804 158.4400;2.多项式的计算;多项式相乘：conv p1=[5,-4,3,-2,1,-6];p2=[1,0,-2,3]; p3=conv(p1,p2) p3 = 5 -4 -7 21 -17 7 -8 15 -18 多项式相除：[q,r]=deconv(p1,p2) p1=[5,-4,-7,21,-17,7,-8,15,-18];p2=[1,0,-2,3]; [q,r]=deconv(p1,p2) q = 5 -4 3 -2 1 -6 r = 0 0 0 0 0 0 0 0 0 p2=[1,0,-2,5];[q,r]=deconv(p1,p2) q = 5 -4 3 -12 9 -32 r = 0 0 0 0 0 0 70 -94 142 qq=poly2sym(q) qq =5*x^5-4*x^4+3*x^3-12*x^2+9*x-32 rr=poly2sym(r) rr =70*x^2-94*x+142;;解： p=[1,1];q=[1,-6,11,-6]; [a,b,r]=residue(p,q) a = 2.0000 -3.0000 1.0000 b = 3.0000 2.0000 1.0000 r = [] ;解： p=1;q=[1,0,0,0,-1]; [a,b,r]=residue(p,q) a = -0.2500 0.2500 -0.0000 + 0.2500 -0.0000 - 0.2500i b = -1.0000 1.0000 0.0000 + 1.0000 0.0000 - 1.0000i r = [];将简单分式之和并为有理分式;3.导数的数值计算;例8 [x,y]=meshgrid(0:0.2:1,0:0.3:1.5); f=x.^2+x.*y; [fx,fy]=gradient(f,0.2,0.3);fx fx = 0.2000 0.4000 0.8000 1.2000 1.6000 1.8000 0.5000 0.7000 1.1000 1.5000 1.9000 2.1000 0.8000 1.0000 1.4000 1.8000 2.2000 2.4000 1.1000 1.3000 1.7000 2.1000 2.5000 2.7000 1.4000 1.6000 2.0000 2.4000 2.8000 3.0000 1.7000 1.9000 2.3000 2.7000 3.1000 3.3000;例9 [x,y]=meshgrid(-2:0.2:2); f=x.*exp(-x.^2-y.^2); [fx,fy]=gradient(f,0.2,0.2); contour(f) hold on quiver(fx,fy) hold off;4.定积分的数值计算;例11 计算;Newton Cotes法：[S,nf]=quad8(‘F’,a,b,tol) 例12：同11。 解： function y=f(x) y=exp(-x.^2/2); format long; [S3,nf]=quad8(Newton,-8,8) S3 = 2.50662827463083 nf = 168;例13 分别用quad函数和quad8函数求定积分的近似值，并在相同的积分精