实验三 matlab数值运算,实验三matlab数值计算.doc

实验三matlab数值计算.doc

实验三MATLAB数值计算一、实验目的熟悉MATLAB多项式的运用。(1)多项式的求值、求根和部分分式展开(2)多项式的乘除法和微积分(3)多项式拟合和插值二、实验内容和步骤1多项式求值函数POLYVAL可以用来计算多项式在给定变量时的值，是按数组运算规则进行计算的。语法POLYVALP,S说明P为多项式,S为给定矩阵。【例1】计算PX3X22X1多项式的值。P321POLYVALP,2计算X2时多项式的值ANS17X0053POLYVALP,X计算X为向量时多项式的值ANS1000027500600001075001700002475003400002多项式求根ROOTS用来计算多项式的根。语法RROOTSP说明P为多项式；R为计算的多项式的根，以列向量的形式保存。与函数ROOTS相反，根据多项式的根来计算多项式的系数可以用POLY函数来实现。语法PPOLYR【例2】计算多项式PXX36X272X27的根以及由多项式的根得出系数。P167227ROOTSP计算多项式的根ANS1212295734503884POLY1212295734503884计算多项式的系数ANS10000600007200002700113特征多项式对于一个方阵S，可以用函数POLY来计算矩阵的特征多项式的系数。特征多项式的根即为特征值，用ROOTS函数来计算。语法PPOLYS说明S必须为方阵；P为特征多项式。【例3】根据矩阵来计算的特征多项式系数。A123456780PPOLYAA123456780P1000060000720000270000RROOTSPR1212295734503884程序分析PX36X272X27为矩阵A的特征多项式，121229,57345和03884为矩阵S的特征根。4部分分式展开用RESIDUE函数来实现将分式表达式进行多项式的部分分式展开。KSNPSR2PSR1ASB语法R,P,KRESIDUEB,A说明B和A分别是分子和分母多项式系数行向量；R是R1R2RN留数行向量；P为P1P2PN极点行向量；K为直项行向量。【例4】将表达式进行部分分式展开。B5327A4083R,P,KRESIDUEB,AB5327A4083R141670665313320P157371164404093K12500程序分析表达式展开结果为。2501493S16S05573S42多项式的乘除法和微积分1多项式的乘法和除法多项式的乘法语法PCONVPL,P2说明P是多项式P1和P2的乘积多项式。多项式的除法语法Q,RDECONVPL,P2说明除法不一定会除尽，会有余子式。多项式P1被P2除的商为多项式Q，而余子式是R。【例5】计算表达式X32X23X410X220X30。U1234V102030CCONVU,VC1040100160170120Q,RDECONVC,UQ102030R0000002多项式的微分和积分多项式的微分由POLYDER函数实现。MATLAB没有专门的多项式积分函数，但可以用P/LENGTHP11,K的方法来完成积分，K为常数。【例6】求多项式3X26X9X22X的微分和积分A369B120KPOLYDERA,BK12364218SLENGTHK11S4321P1K/S,0多项式积分,常数K0P1312211803多项式拟合和插值1多项式拟合多项式曲线拟合是用一个多项式来逼近一组给定的数据，使用POLYFIT函数来实现。拟合的准则是最小二乘法，即找出使最小的FX。2N1IIIYFX语法PPOLYFITX,Y,N说明X、Y向量分别为N个数据点的横、纵坐标；N是用来拟合的多项式阶次；P为拟合的多项式，P为N1个系数构成的行向量。【例7】对多项式Y2X3X25X10曲线拟合。经过一阶、二阶和三阶拟合的曲线如图2所示。图一阶、阶和三阶拟合曲线X110P21510Y1POLYVALP,XY1COLUMNS1THROUGH7163270142260436682COLUMNS8THROUGH10101014321960P1POLYFITX,Y1,1一阶拟合P120480005224000P2POLYFITX,Y1,2二阶拟合P232000014720001816000P3POLYFITX,Y1,3三阶拟合P32000010000500001000002插值运算插值运算是根据数据点的规律，找到一个多项式表达式可以连接两个点，插值得出相邻数据点之间的数值。1一维插值一维插值是指对一个自变量的插值，INTERP1函数是用来进行一维插值的。语法YIINTERP1X,Y,XI,’’说明X、Y为行向量；XI是插值范围内任意点的X坐标，YI则是插值运算后的对应Y坐标；是插值函数的类型，“LINEAR”为线性插值默认，“NEAREST”为用最接近的相邻点插值，“SPLINE”为三次样条插值，“CUBIC”为三次插值。【例8】经过线性插值、三次样条插值计算出横坐标为95的对应纵坐标，如图所示。图线性插值和三次样条插值X110P21510Y1POLYVALP,XY1COLUMNS1THROUGH7163270142260436682COLUMNS8THROUGH10101014321960Y11INTERP1X,Y1,95线性插值Y111696Y12INTERP1X,Y1,95, SPLINE 三次样条插值Y1216822二维插值二维插值是指对两个自变量的插值，INTERP2函数是用来进行二维插值的。语法ZIINTERP2X,Y,Z,XI,YI,’’说明是插值函数的类型有，“LINEAR”为双线性插值默认，“NEAREST”为用最接近点插值，“CUBIC”为三次插值。AVERAGE{2}AVERAGE_SCORES思考题1计算PXX410X32X24X3多项式的在X5值。2表达式进行部分分式展开。10S3计算表达式5X34X27X19X5的根以及由多项式的根得出系数。4求多项式5X24X93X26X的微分和积分