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

第五章 函数、导数与积分等 的数值计算 函数值、多项式、导数、定积分、二重积分、非线性方程（组）、级数的计算 1.函数值的计算 解： x=pi/6; y=sin(x)+sin(2.*x)+sin(3.*x); x,y x = 0.5236 y = 2.3660 2.多项式的计算 多项式的表示法：p=[1,0,-12,5] 多项式在给定点的值：p=[4,-3,2,1]; y=polyval(p,5) y = 436 多项式方程求根:roots p1=[1,0,-12,5];p2=[1,0,-12,20];x1i=roots(p1) x1i = -3.6562 3.2332 0.4230 x2i=roots(p2) x2i = -4.1072 2.0536 + 0.8075i 2.0536 - 0.8075i 多项式相乘：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,0,0,0];q=[1,1,-2]; [a,b,r]=residue(p,q) a = 2.6667 0.3333 b = -2 1 r = 1 -1 解： 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 = [] 将简单分式之和并为有理分式 例6 a=[2.6667,0.3333]; b=[-2,1];r=[1,-1]; [p,q]=residue(a,b,r) p = 1.0000 0 0 -0.0001 q = 1 1 -2 3.导数的数值计算 [fx,fy]=gradient(f, x, y) yx=diff(y)./diff(x) 例7 x=0:0.5:3;f=x.^2;fx=gradient(f,0.5) fx = 0.5000 1.0000 2.0000 3.0000 4.0000 5.0000 5.5000 4.定积分的数值计算 Simpson法：[S,nf]=quad(‘F’,a,b,tol) 例11 计算 解： function y=f(x) y=exp(-x.^2/2); [S2,nf]=quad(simp12,-8,8) S2 = 2.5066 nf = 105 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.5066282746