第五章数值计算.docVIP

第五章 数值计算 5.1多项式的计算 在全部初等函数里，多项式是结构最为简单，应用最为广泛的一类函数，有必要单独列出进行讨论。 1．多项式的表示法 在MATLAB中将多项式采用降幂排列，并可表示为向量的形式。 例如，多项式p(x)=x3-12x+5可用向量表示为 p(x)=[1 0 –12 5] 2．多项式在给定点的值 利用MATLAB中的polyval函数，可以直接求出多项式在给定点的值。 y=polyval(p,x0),其中p为多项式，x0为给定点。 【例5－1】求多项式p=[4 -3 2 1]在点x=[1 2 3 4 5]处的值 输入 p=[4 -3 2 1]; x=[1 2 3 4 5] y=polyval(p,x); %求多项式在x处的值。 x,y 运行后可得 x = 1 2 3 4 5 y = 4 25 88 217 436 3．多项式方程求根 利用MATLAB中的roots函数，可以直接求得多项式方程的根。 p(xi)=0的根xi用函数xi=roots(p). 【例5－2】求多项式方程的根 输入 p1=[1 0 -12 5];p2=[1 0 -12 20]; x1=roots(p1); x2=roots(p2); x1,x2 运行后可得 x1 = -3.6562 3.2332 0.4230 x2 = -4.1072 2.0536 + 0.8075i 2.0536 - 0.8075i 4．多项式相乘 两多项式相乘，可以利用MATLAB中的conv(convolution)函数来求得. p3=conv(p1,p2), 5、多项式相除： 两多项式与相除，可以利用MATLAB中的deconv 函数来求得，其调用格式为[q,r]=deconv(p1,p2),其中p1为被除式，p2为除式，q为商式，r为余式。 【例5－3】多项式相乘与相除 输入 p1=[5 -4 3 -2 1 -6]; p2=[1 0 -2 3]; p3=conv(p1,p2); [q,r]=deconv(p1,p2); p3,q,r 运行后可得 p3 = 5 -4 -7 21 -17 7 -8 15 -18 q = 5 -4 13 r = 0 0 0 -25 39 -45 6有理分式的分解与合并： 有理分式同多项式有着紧密的联系，依照代数定理，任何有理分式都可分解为部分分式（或称最简分式）之和 式中pm(x)与qn(x)为m次与n次多项式，当m≥n时，r(x)为商式，当mn时，r(x)为空，bi为方程qn(x)=0的根。ai为对应的实数或复数。 将有理分式进行分解，是一件十分复杂和繁琐的工作。如果我们调用MATLAB里的residule函数，将会使问题变得简单容易，其调用格式为 [a,b,r]=residule(p,q); 有理分式分解. [p,q]=residule(a,b,r). 将简单分式之和合并为有理分式. 【例5－4】将有理分式分解为最简分式之和。 输入 p=[1 0 0 0];q=[1 1 -2]; [a,b,r]=residue(p,q); a,b,r 运行后可得 a = 2.6667 0.3333 b = -2 1 r = 1 -1 即： 【例5－5】将有理分式分解为最简分式之和。 输入 p=[1 1]; q=[1 -6 11 -6]; [a,b,r]=residue(p,q); a,b,r 运行后可得 a = 2.0000 -3.0000 1.0000 b = 3.0000 2.0000 1.0000 r = [] 故求得 【例5－6】将简单分式之和合并为有理分式。 输入 a=[2.667 0.333]; b=[-2 1]; r=[1 -1]; [p,q]=residue(a,b,r); p,q 运行后可得 p = 1.0000 0 0 -0.0010 q = 1 1 -2 5．2导数的数值计算 导数的数值计算或称数值求导，是指导数值的一类近似计算方法，它不必去求出导函数（需要符号求导运算），而是利用函数值的差分来构造出近似导数值。 下面为求导的几个常见命令： 1、yx=diff(y)./diff(x)（即y′=△y/△x） 其中diff=differential，差别的，微分。 2、[fx,f