实验二：多项式与线性方程组求解.docVIP

实 验 报 告 课程名称： MATLAB语言与应用技术 实验名称： 多项式与线性方程求解 院 （系）： 机电学院 专业班级： 工程机械1101 姓 名： 甘超 学 号： 110730123 指导教师： 郑建校 2013 年11 月4日 实验二：多项式与线性方程求解 实验地点：大楼五楼八号机房 试验时间：2013年11月4日 实验目的：1.掌握MATLAB中的多项式创建及其运算 2. 掌握MATLAB中线性方程组的求解 3.熟悉MATLAB的基本特性以及软件操作 实验内容： 多项式的表示和创建 例：多项式5*x^4-4*x^3-2*x+的行矢量表示 p=[5 -4 0 -2 1]; poly2sym(p,x) ans = 5*x^4-4*x^3-2*x+1 多项式的基本运算 例：求多项式p1和p2的加减乘除 P1=3*x^4+2x^2-8x+5 P2=3x^3+x^2+6 p=[5 -4 0 -2 1]; poly2sym(p,x) ans = 5*x^4-4*x^3-2*x+1 p1=[3 0 2 -8 5]; p2=[0 3 1 0 6]; p1+p2 ans = 3 3 3 -8 11 p1-p2 ans = 3 -3 1 -8 -1 w=conv(p1,p2) w = 0 9 3 6 -4 7 17 -48 30 3.求多项式的根 p2=[0 3 1 0 6]; roots(p2) ans = -1.3814 0.5240 + 1.0831i 0.5240 - 1.0831i 4.多项式的微分·积分和求值 p=[3 3 3 -8 11] p = 3 3 3 -8 11 polyder(p) ans = 12 9 6 -8 pp=[1 2 3 4] pp = 1 2 3 4 polyint(pp) ans = 0.2500 0.6667 1.5000 4.0000 0 polyval(p,2) ans = 79 二．线性方程组求解 1.恰定方程组求解 A=[10 -2 -1;-2 10 -1;-1 -2 5]; b=[3 15 10]; x=inv(A)*b x = 1 2 3 x=A\b x = 1 2 3 2线性超定方程组求解 A=[1 2 -3;3 4 2;-6 7 0;4 -2 3]; b=[5 2 8 6]; x=pinv(A)*b x = 0.3115 1.0189 0.2910 x=A\b x = 0.3115 1.0189 0.2910 3.线性欠定方程组求解 A=[1 -2 3;0 1 -1;-1 0 -1;1 -3 4]; b=[4 -3 -4 1]; x=pinv(A)*b x = 2.2549 1.2157 1.0392 y=A\b Warning: Rank deficient, rank = 2 tol = 4.6151e-015. y = 3.4706 0 -0.1765 norm(x) ans = 2.7645 norm(y) ans = 3.4751 4.齐次线性方程组求解 A=[1 -2 3 -4;0 1 -1 1;-1 0 -1 2 ;1 -3 4 -5]; a=null(A) a = -0.4565 -0.6770 0.5667 0.1102 0.6770 -0.4565 0.1102 -0.5667 a*a ans = 1.0000 -0.0000 -0.0000 1.0000 实验总结：多项式在代数中占有重要地位，它是最简单最容易计算的函数，但同时也是其他绝大多数复杂函数的构成基础。MATLAB提供了多项式的表达以及创建和各种运算方法，处理起来非常简单。同时在有关线性方程组的问题处理中，MATLAB也