Matlab编程与应用第五讲.PDF

Matlab编程与应用 第五讲 中国科技大学信息学院 陆伟 luwei@ustc.edu.cn 本讲内容  part1: Matlab多项式运算  part2: 微分方程数值解法  part3: 符号计算 Part1: Matlab多项式运算 Matlab中多项式表示方法：  n 阶多项式用长度为 n+1的向量表示， 缺少的幂次项系数为 0 ! p x a x a x ax (a)  n n1  n n1 1 0 [ ,a a, , , a ]a  1 n n 1 0 Matlab 中表示为向量： 例 3 2 ： [2 -1 0 3 ] 2 x 3x  系数中的零不能省！  多项式显示: poly2sym(p,’x’) 多项式加减运算 多项式加减就是其所对应的系数向量的加减运算  对于次数相同的多项式，可以直接对其系数向量进行 加减运算；  如果两个多项式次数不同，则应该把低次多项式中系 数不足的高次项用 0 补足，然后进行加减运算。 3 2 [2, 1, 0, 3]  p  2x x  3 例： 1 p  x22  1 [ 0, 0, 2,[ 2, 1] 1] p p  x x2  x  2 3 4 2 [2, 1, 2, 4]  1 2 多项式加减运算  练习：编一个自动补零的多项式加减函数 p = mypolyadd(a,b) 输入：a,b为待相加多项式对应的向量（可以是列向量） 输出：相加后多项式对应的向量。 多项式乘法、除法运算 k =conv (p,q)  多项式乘法运算： 3 2 2 1x 2 x 3x  例：计算多项式 和 的乘积 p=[2,-1,0,3]; q=[2,1]; k=conv(p,q);  多项式除法运算：[k,r]=deconv (p,q) k p q r 其中 返回的是多项式 除以 的商， 是余式。 [k,r