2010-2011-01_Matlab复习.docVIP

2010-2011-01 Matlab复习 专题一 启动与退出Matlab 命令窗口及使用（默认窗口default） 简单计算：+、-、*、/、^、sin、pi、sqrt、exp 命令行编辑：调回已输入的内容 命令历史窗口 工作路径和搜索路径 帮助 专题2 多项式 多项式的根 求多项式的根（即使多项式等于0的解）是经常遇到的问题。在Matlab中，一个多项式是用多项式系数行向量表示，向量中的系数按照其所对应的自变量的阶次的降序进行排列。 例如：多项式在Matlab中的表示 p=[1 -12 0 25 116] p = 1 -12 0 25 116 注意：不存在，所以系数为0 然后利用函数roots来求多项式的根。 r=roots(p) r = 11.7473 2.7028 -1.2251 + 1.4672i -1.2251 - 1.4672i 如果知道一个多项式的根，能不能构建相应的多项式呢？ 用poly函数从一个根向量中构建一个多项式向量。 pp=poly(r) pp = 1.0000 -12.0000 -0.0000 25.0000 116.0000 多项式的乘法 多项式的乘法实际上是多项式系数向量之间的卷积（convolution）Matlab中使用conv函数。 例如：求和的乘积 a=[1 2 3 4];b=[1 4 9 16];c=conv(a,b) c = 1 6 20 50 75 84 64 结果 多项式的加法 Matlab没有提供专门的函数执行多项式加法。 如果两个多项式向量长度相同，则多项式加法就是将两个多项式直接相加。 例子同上 d=a+b d = 2 6 12 20 如果两个多项式阶次不同时，补0 例如： e=c+[0 0 0 d] e = 1 6 20 52 81 96 84 多项式的除法 Matlab中用deconv函数 例如： [q,r]=deconv(c,b) q = 1 2 3 4 r = 0 0 0 0 0 0 0 q:quotient(商) ；； r:remainder(余数) 多项式的微分和积分 Matlab中微分操作函数polyder，积分操作函数polyint 例如： g=[1 6 20 48] g = 1 6 20 48 h=polyder(g) h = 3 12 20 f=polyint(h) f = 1 6 20 0 多项式求值 使用函数polyval 例如：对多项式在[-1,3]之间求值 p=[1 4 -7 -10];x=linspace(-1,3); v=polyval(p,x); plot(x,v) 曲线拟合 曲线拟合是进行数据分析时经常遇到的问题，它是根据一组或多组测量数据找出一条数学上可描述的曲线的过程。这条曲线有时候将穿过测量的数据点，而有时候将会非常接近于但不会穿过测量的数据点。评价一条曲线是否准确地描述了测量数据的最通用的方法，是看测量数据点与该曲线上对应点之间的平方误差是否达到最小，这种曲线拟合的方法称为最小二乘曲线拟合。 Matlab提供了函数polyfit用于实现最小二乘多项式曲线拟合。 例如： p1、p3的结果是多项式，然后计算多项式的值作图 专题3 矩阵分析和操作 矩阵的创建 如果是列向量，[]中间用； 创建二维数组 元素的访问 多维数组的创建 矩阵的逆 inv(A) 求解线性方程组 也可以用左除，A\b 计算行列式值 det(A) 元素的排序 [B,p]=sort(A,dim) dim=1,按列排序，dim=2,按行排序，B是排序后的矩阵，p是记录B中的元素在A中的位置。 例题参见实验指导书 数据分析 专题4 数据插值 有时候，我们并不仅仅对给定点处的函数值感兴趣，还对这些数据点中间的某些点处的函数值感兴趣。当我们无法快速对这些中间点执行函数运算时，就需要用到数据插值的概念。 其实数据插值的概念并不难理解，在我们用Matlab进行绘图时，就用到了数据插值。由于Matlab在默认情况下需要用连续的直线来连接各个数据点，那么这些数据点之间的各个数据就需要通过线性插值的方法获得。线性插值就是数据插值的一种。该方法建设两个数据点之间的中间值都落在这两个数据点连成的直