第7讲Matlab数值计算.pptVIP

Matlab Language Matlab数值计算 第7讲 * Matlab Language * 6.1 多项式(polynomial) 多项式的MATLAB表达 多项式由一个行向量表示 该向量元素是该多项式的系数 且按降幂次序排列 如：多项式x4－12x3＋25x＋116由行向量： p=[1 -12 0 25 116]表示。 注意，必须包括具有零系数的项。 求解多项式的根？ roots指令 * Matlab Language * 6.1 多项式(polynomial) 举例：求解多项式x4－12x3＋25x＋116的根 p=[1 -12 0 25 116] p = 1 -12 0 25 116 r=roots(p) r = 11.7473 2.7028 -1.2251 + 1.4672i -1.2251 - 1.4672i MATLAB按惯例规定，多项式是行向量，根是列向量 * Matlab Language * 6.1 多项式(polynomial) 已知多项式的根，求对应的多项式 ？ 能！ 使用poly指令 举例：由上例所得的根求其多项式 pp=poly(r) pp = 1.0000 -12.0000 -0.0000 25.0000 116.0000 即：x4－12x3＋25x＋116 * Matlab Language * 6.1 多项式(polynomial) 多项式的乘法(conv指令) 举例：多项式a(x)=x3＋2x2＋3x＋4 和b(x)= x3＋4x2＋9x＋16的乘积。 a=[1 2 3 4] ; b=[1 4 9 16]; c=conv(a, b) c = 1 6 20 50 75 84 64 两个以上的多项式的乘法需要重复使用conv. * Matlab Language * 6.1 多项式(polynomial) 多项式加法: MATLAB没有提供进行加法运算的函数。 如果两个多项式向量大小相同，标准的数组加法有效。 把多项式a(x)与上面给出的b(x)相加。 d=a+b d = 2 6 12 20 结果: d(x)= 2x3＋6x2＋12x＋20 当两个多项式阶次不同，低阶的多项式必须用首零填补，使其与高阶多项式有同样的阶次。 考虑上面多项式c和d相加： e=c+[0 0 0 d] e = 1 6 20 52 81 96 84 结果: d(x)= x6+6x5+20x4+52x3+81x2+96x+84 * Matlab Language * 6.1 多项式(polynomial) 问题：编写一个多项式加法运算的函数文件 明确需求 定义一个函数，需要有两个输入参数（比如：p1、p2），一个输出参数p_out。 p1、p2表示两个待计算的多项式。 p_out表示两个多项式的求和结果 函数的内部处理 如p1、p2两参数大小相等，则直接相加: p_out=p1+p2 如p1、p2两参数大小不等： If length(p1)length(p2) P2前面要补0元素，使p1、p2两参数大小相等 否则 P1前面要补0元素，使p1、p2两参数大小相等 * Matlab Language * 6.1 多项式(polynomial) 函数文件的编写 function p_out=poly_sum(p1, p2) % calculate the sum of two polynomials if length(p1)==length(p2) p_out=p1+p2; elseif length(p1)length(p2) p2=[zeros(1,length(p1)-length(p2)), p2]; p_out=p1+p2; else p1=[zeros(1,length(p2)-length(p1)), p1]; p_out=p1+p2; end * Matlab Language * 6.1 多项式(polynomial) 多项式的除法(deconv) 举例说明： c(x)=x6＋6x5＋20x4＋50x3＋75x2＋84x＋64 除以b(x)= x3＋4x2＋9x＋16 c=[1 6 20 50 75 84 64]; b=[1 4 9 16