MATLAB -1 数值计算.ppt

第五章 MATLAB的数值计算 1.多项式计算 2. 方程组求解 3. 函数及其数值分析 4.数据处理 5.数据插值 6.曲线拟合 7.傅里叶变换 1.多项式 运算 多项式表达方式的约定 创建多项式的方法 多项式运算函数 多项式表达方式的约定 MATLAB约定降幂多项式P(x)=a0xn+ a1xn-1 +…+an-1x+an用以下系数矢量（系数行向量）表示：p= [a0, a1 ,…,an-1,an],即把多项式的各项系数依降幂次序排放在行向量的元素位置上。 注意I：假如多项式中缺某幂次项，则应认为该幂次项的系数为零。 注意II：多项式加减时低阶多项式必须首零填补使得与高阶多项式阶次相同。 创建多项式的方法 系数矢量的直接输入法 在命令窗直接输入多项式的系数矢量， 由根矢量创建多项式 由给定的根矢量创建多项式，由函数poly实现。 特征多项式输入法 由矩阵的特征多项式取得，由函数poly实现。 多项式的字符标示：poly2str 多项式求根 n次多项式具有n个根，当然这些根可能是实根，也可能含有若干对共轭复根。MATLAB提供的roots函数用于求多项式的全部根，其调用格式为： x=roots(P) 其中P为多项式的系数向量，求得的根赋给向量x，即x(1),x(2),…,x(n)分别代表多项式的n个根。 例: 求多项式x4+8x3-10的根。 命令如下： A=[1,8,0,0,-10]; x=roots(A) 若已知多项式的全部根，则可以用poly函数建立起该多项式，其调用格式为： P=poly(x) 若x为具有n个元素的向量，则poly(x)建立以x为其根的多项式，且将该多项式的系数赋给向量P。 多项式运算：乘运算 例:a(x)=x2+2x+3; b(x)=4x2+5x+6; c = (x2+2x+3)(4x2+5x+6) a=[1 2 3];b=[4 5 6]; c=conv(a,b)=conv([1 2 3],[4 5 6]) c = 4.00 13.00 28.00 27.00 18.00 p=poly2str(c,x) p = 4 x^4 + 13 x^3 + 28 x^2 + 27 x + 18 多项式除运算deconv 多项式微分 matlab提供了polyder函数多项式的微分。 命令格式： polyder(p): 求p的微分 polyder(a,b): 求多项式a,b乘积的微分 [p,q]=polyder(a,b): 求多项式a,b商的微分 例：a=[1 2 3 4 5]; poly2str(a,x) ans = x^4 + 2 x^3 + 3 x^2 + 4 x + 5 b=polyder(a) b = 4 6 6 4 poly2str(b,x) ans =4 x^3 + 6 x^2 + 6 x + 4 多项式积分 多项式积分命令函数polyint 1 polyint(p) 多项式p的积分，常数项为0 2 polyint(p,k)多项式p的积分，常数项为k 有理多项式的运算函数 多项式之比的分式展开 [r,p,k]=residue(a,b) 从展开分式得到多项式 [a,b]=residue(r,p,k) 注意residue函数的可逆性 多项式的求值 MATLAB提供了两种求多项式值的函数：polyval与polyvalm，它们的输入参数均为多项式系数向量P和自变量x。两者的区别在于前者是代数多项式求值，而后者是矩阵多项式求值。 代数多项式求值 polyval函数用来求代数多项式的值，其调用格式为： Y=polyval(P,x) 若x为一数值，则求多项式在该点的值；若x为向量或矩阵，则对向量或矩阵中的每个元素求其多项式的值。 例 已知多项式x4+8x3-10，分别取x=1.2和一个2×3矩阵为自变量计算该多项式的值。 例: 画出上例x在-5到5间的曲线 矩阵多项式求值【了解】 polyvalm函数用来求矩阵多项式的值，其调用格式与polyval相同，但含义不同。polyvalm函数要求x为方阵，它以方阵为自变量求多项式的值。 设A为方阵，P代表多项式x3-5x2+8，那么polyvalm(P,A)的含义是： A*A*A-5*A*A+8*eye(size(A)) 而polyval(P,A)的含义是： A.*A.*A-5*A.*A+8*ones(size(A)) 例 仍以多项式x4+8x3-10为例，取一个2×2矩阵为自变量分别用polyval和polyvalm计算该多项式的值。 2. 方程组求解 代数方程组求解 方法1 矩阵的除运算 （教材） 左除 求逆 （了解：各种方程类型 恰定，超定，欠