第四章初等代数+第五章线性代数.pptVIP

例5.2.5、 三、初等变换法 机动 目录 上页 下页 返回 结束 将方程组化为阶梯形方程组，如果最后一个方程是零等于一个非零数，则方程组无解，否则有解。在有解的情况下，如方程的个数等于未知量个数，则方程组有唯一解；如果方程的个数小于未知量的个数，则方程组有无穷多个解。 命令: rref(A), 功能: 将A化为阶梯矩阵。 例5.2.6. 求矩阵 机动 目录 上页 下页 返回 结束 解：通过Matlab运算，分析知：原方程组等价于 * (参考L. P552 第四节,一. ) * 运行时, 点击按钮“性质7”, 可显示性质7. 第4章 初等代数运算命令与例题 第一节 一、多项式的表达 Matlab中用多项式的系数组成的行向量来表示多项式 机动 目录 上页 下页 返回 结束 多项式的表达与运算 第4章 多项式: 命令poly2str(p,’x’)，可以给出多项式的习惯形式。 二、 多项式的运算 机动 目录 上页 下页 返回 结束 （1）多项式的加法（减法） 例：将x^2+1与x^2+x+1作和，并与 多项式x^3-x+2作差。 命令形式：a+b或a-b （2）多项式乘法 命令形式：conv(a,b) 注：如果使用命令p= conv(a,b)，则可以得到乘积结果存放在变量p中，p可以选用任何其他字母。 例：求多项式x^2+1和x^3-x+3的乘积，并用习惯形式表达。 机动 目录 上页 下页 返回 结束 a=[1 0 1];b=[1 0 -1 3]; conv(a,b) 解： poly2str(ans,x) （3）多项式除法 机动 目录 上页 下页 返回 结束 命令形式：[q,r]=deconv(a,b) 功能：做a除以b的运算，计算结果得到的商是q,余子式是r。q,r都是向量且可以选用任何其他字母。 例：求多项式x^4+7x^3+16x^2+18x+8与x+3的商及余子式。 （4）求方阵的特征多项式 命令形式：poly(A) 机动 目录 上页 下页 返回 结束 例：求3阶方阵A=[1,2,3;4,15,60;7,8,9]的特征多项式。 （5）对多项式求导 命令形式：polyder(a) 例：求3多项式x^4+7x^3+16x^2+18x+8 的导数。 （6）对多项式乘积求导 命令形式：polyder(a,b) 机动 目录 上页 下页 返回 结束 例：求x^2+1和x^3-x+3的乘积的导数。 （7）多项式数据拟合 命令形式：polyfit(x,y,n) 功能：求由向量x,y对应的数据点的n次多项式拟合函数。 例4.1.1：求由以下数据的6次拟合多项式 x 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 y 2.1 2.3 2.5 2.9 3.2 3.3 3.8 4.1 4.9 5.4 5.8 (8)计算x为向量或矩阵点处多项式的值 命令形式：polyval(a,x) 功能：计算多项式a的变量在点阵x处的值。X 可以为向量或矩阵，计算结果是与x同维的向量或矩阵。 例：求多项式3*x^2+2*x+1在x=1,4,5,8时的值。 例：随机产生一个3阶方阵，并求出多项式4*x^3-3*x+12在此方阵处的值。 第二节 有理多项式的运算 有理函数的运算函数 例：求多项式x^5+x^4-6x^3+2x^2+4除以x^3+4x^2+x-7的导数。 机动 目录 上页 下页 返回 结束 函数名称 功能简介 [Num, Den]=polyder(p1, p2) [r, p, k]=residue(a, b) [a, b]=residue(r, p, k) 对有理分式(p1/p2)求导数 部分分式展开式 部分分式组合 第三节 代数式的符号运算 机动 目录 上页 下页 返回 结束 指令 含义 p=factor(s) p=expand(s) p=collect(s) p=collect(s,’v’) p=simplify(s) sn=subs(s,’old’,’new’) r=vpa(sn) p是对s定义的多项式进行因式分解的结果 p是对s定义的多项式进行展开的结果 把s中x的同幂项系数进行合并 把s中v的同幂项系数进行合并 对s进行简化 这两条实现代数式的求值问题。sn是变量替换后的符号表达式的变量名，s为替换前，old为被替换变量，new为替换变量，r为最终求得结果。 机动 目录 上页