Matlab语言及其在电子信息工程中应用第6讲(V1.2).pptVIP

Matlab语言及其在电子信息工程中的应用 主讲：许钢 （V1.2） 第六讲 MATLAB多项式函数及符号计算 一、多项式 二、符号表达式和符号矩阵的创建 三、符号矩阵的基本运算 四、函数画图、求导、求零点、极值 第六讲 MATLAB多项式函数及符号计算 一、多项式 第六讲 MATLAB多项式函数及符号计算 一、多项式 第六讲 MATLAB多项式函数及符号计算 一、多项式 第六讲 MATLAB多项式函数及符号计算 一、多项式 第六讲 MATLAB多项式函数及符号计算 一、多项式 第六讲 MATLAB多项式函数及符号计算 一、多项式 第六讲 MATLAB多项式函数及符号计算 一、多项式 第六讲 MATLAB多项式函数及符号计算 一、多项式 第六讲 MATLAB多项式函数及符号计算 一、多项式 第六讲 MATLAB多项式函数及符号计算 一、多项式 第六讲 MATLAB多项式函数及符号计算 一、多项式 第六讲 MATLAB多项式函数及符号计算 二、符号表达式和符号矩阵的创建 第六讲 MATLAB多项式函数及符号计算 二、符号表达式和符号矩阵的创建 第六讲 MATLAB多项式函数及符号计算 二、符号表达式和符号矩阵的创建 第六讲 MATLAB多项式函数及符号计算 三、符号矩阵的基本运算 第六讲 MATLAB多项式函数及符号计算 四、函数画图、求导、求零点、极值 第六讲 MATLAB多项式函数及符号计算 四、函数画图、求导、求零点、极值 第六讲 MATLAB多项式函数及符号计算 四、函数画图、求导、求零点、极值 第六讲 MATLAB多项式函数及符号计算 四、函数画图、求导、求零点、极值 第六讲 MATLAB多项式函数及符号计算 四、函数画图、求导、求零点、极值 第六讲 MATLAB多项式函数及符号计算 （第六讲 完） *Matlab语言及其在电子信息工程中的应用 　　在MATLAB 里，多项式用其系数向量表示： 对应于向量： 1、多项式的创建 创建 生成的多项式的系数向量： poly(A) 　　%创建矩阵A的特征多项式。 2、多项式的常用函数 poly(q) %构造一个以q向量为根的多项式； poly(A) %得出方阵A的特征多项式； roots(p) %返回多项式的根向量 注1：多项式的零系数项要填上零。 polyval(p,x) %求多项式p在某点x处的值; 3、多项式的加减法 （1）依向量加法 例：a=[1,2,3,4]; b=[1,4,9,16]; d=a+b 系统回复 d= 2 6 12 20 当两个多项式次数不同时，要用0补足 例：c=[1,-5,0,1,-5] %四次多项式 e=c+[0,d] %将三次多项式d前面加一个0，补足成四次多项式，再相加。 系统回复 e= 1 -3 6 13 15 结果是： 3、多项式的加减法 （2）实现多项式加法的M-文件 例： function p=mmpadd(a,b) %MMPADD Polynomial addition if nargin2 error(Not enough input arguments) end a=a(:); %转变成行向量 b=b(:); na=length(a); nb=length(b); p=[zeros(1,nb-na) , a]+[zeros(1,na-nb) , b]; 4、卷积及其解卷积 卷积对应求多项式的乘积，解卷积则对应多项式的除法。 conv(a,b) %卷积：求多项式a,b的乘积 例： a=[1,2,3,4]; b=[1,4,9,16]; c=conv(a,b) 4、卷积及其解卷积 [q,r]=deconv(c,a) %解卷积 %用c(x)除以a(x)反求出上例 %中b(x),其中q是商，r是余式。 例： a=[1,2,3,4]; b=[1,4,9,16]; c=conv(a,b) [q,r]=deconv(c,a) c(3)=c(3)+2; [q,r]=deconv(c,a) 5、多项式曲线拟合 polyfit(x,y,n) 用最小二乘法，拟合出所给数据x,y的n次多项式系数。 例： x=[1,2,3,4,5]; y=[5.5,43.1,128,290.7,498.4]; plot(x,y,r*); hold on p=polyfit(x,y,4) x1=1:0.1:5; y1=polyval(p,x1); plot(x1,y1); hold off 6、多项式的展开 [r,p,k]=residue(b,a) %