澳大利亚(二)学生.ppt

机电工程系 机电工程系 第6章 数值计算 多项式、曲线拟合和插值、数值微积分、线性代数 6.1 多项式计算 6.1.1多项式的表达和创建 6.1.2多项式四则运算和微积分运算 6.1.3多项式的求值、求根和部分分式展开 6.1.4两个有限长序列的卷积 6.1.1多项式的表达和创建 高数中多项式的表达： Matlab约定以上的降幂多项式用一个长度为（n+1）的降幂系数行数组表示，如下所示 注：多项式中缺少的幂次项，用0补位 6.1.1多项式的表达和创建 例 创建多项式 6.1.2多项式四则运算和微积分运算 多项式的加法和减法 多项式加法/减法 对应项系数加/减 1.阶次相同可直接相加/减； 2.阶次不同补阶再相加/减； 3.也可定义函数：ployadd； 6.1.2多项式四则运算和微积分运算 例2.调用polyadd函数完成两个多项式的相加运算。 6.1.2多项式四则运算和微积分运算 建立polyadd.m文件名 并存储到搜索路径 function[poly]=polyadd(poly1,poly2) %polyadd(poly1,poly2)完成两个阶次可能不同的多项式相加 if length(poly1)length(poly2) short=poly1; long=poly2; else short=poly2; long=poly1; end mz=length(long)-length(short); if mz0 poly=[zeros(1,mz),short]+long; else poly=long+short; end 以上函数polyadd(文件名：polyadd.m), 由密西根大学J.Shriver编写 6.1.2多项式四则运算和微积分运算 编写M脚本文件exam6_2.m如下 clear a=[1,2,0,-5,5]; b=[1,2,3]; c=polyadd(a,b) 运行exam6_2.m，结果如下 结果即多项式 6.1.2多项式四则运算和微积分运算 例3 完成例2中两个多项式b(x)和a(x)相减运算。 6.1.2多项式四则运算和微积分运算 6.1.2多项式四则运算和微积分运算 多项式乘法和除法 多项式相乘是一个卷积过程，两式相乘可以通过两式系数的卷积来完成; Matlab中，用函数conv实现,调用格式为： c=conv(f1,f2)。 多项式除法是乘法的逆运算， Matlab中，用函数deconv实现,调用格式为： [q,r]=deconv(f1,f2)。 6.1.2多项式四则运算和微积分运算 例4 完成例2中两个多项式相乘运算 6.1.2多项式四则运算和微积分运算 编写M脚本文件exam6_4.m如下： clear a=[1,2,0,-5,5]; b=[1,2,3]; e=conv(a,b) 6.1.2多项式四则运算和微积分运算 例5 完成例2中两个多项式相除运算 6.1.2多项式四则运算和微积分运算 编写M脚本文件exam6_5.m如下： clear a=[1,2,0,-5,5]; b=[1,2,3]; [q,r]=deconv(a,b) 6.1.2多项式四则运算和微积分运算 多项式微分和积分 多项式的微分由polyder函数实现； 调用格式为 dp = polyder(P),它的功能是求多项式P的微分。 dp = polyder(P,Q)----------是求P*Q的微分。 [p, q] = polyder(P,Q)------是求P/Q的微分 微分的分子存入p，分母存入q。 matlab中没有专门的积分函数，但可以用指令[p./(length(p):-1:1),k]来实现，其中k为常数。 6.1.2多项式四则运算和微积分运算 例6 完成例2中多项式b(x)微分和积分运算 6.1.3多项式求值、求根和部分分式展开 多项式求值 多项式求值函数： 按数组规则计算 （polyval） 按运算规则计算（polyvalm） polyval调用格式: y= polyval(p,x)------求多项式p在x点的值， x 也可以为一数组 polyvalm调用格式 y= polyvalm(p,x)------求多项式p对于矩阵x的值 要求x必须是方阵，x如果