基础篇2,矩阵与数组.pptVIP

向量点积 函数 dot 格式 C = dot(A,B) %若A、B为向量，则返回向量A与B的点积，A与B长度相同；若为矩阵，则A与B有相同的维数。 C = dot(A,B,dim) %在dim维数中给出A与B的点积 在数学上，两向量的叉乘是一个过两相交向量的交点且垂直于两向量所在平面的向量。在Matlab中，用函数cross实现。 函数 cross 格式 C = cross(A,B) %若A、B为向量，则返回A与B的叉乘，即C=A×B，A、B必须是3个元素的向量；若A、B为矩阵，则返回一个3×n矩阵，其中的列是A与B对应列的叉积，A、B都是3×n矩阵。 C = cross(A,B,dim) %在dim维数中给出向量A与B的叉积。A和B必须具有相同的维数，size(A,dim)和size(B,dim)必须是3。 matlab语言把多项式表达成一个行向量，该向量中的元素是按多项式降幂排列的。 N+1个元素构成的行向量表示n阶多项式，按降幂排列 P(x)=anxn+a n-1 x n-1+…+a1x+a0 则用P=[an,a n-1,…a0]来表示 例:a=[1 2 3;4 5 6;7 8 0]; p=poly(a) p =1.00 -6.00 -72.00 -27.00 p是多项式p(x)=x3-6x2-72x-27的matlab描述方法，我们可用： p1=poly2str(p,‘x’) — 函数文件，显示 数学多项式的形式 p1 =x^3 - 6 x^2 - 72 x – 27 matlab规定多项式系数向量用行向量表示，一组根用列向量表示。 （3）*多项式的拟合 多项式拟合又称为曲线拟合，其目的就是在众多的样本点中进行拟合，找出满足样本点分布的多项式。这在分析实验数据，将实验数据做解析描述时非常有用。 命令格式：p=polyfit(x,y,n)，其中x和y为样本点向量，n为所求多项式的阶数，p为求出的多项式。 例如：卫星真实轨道，并非二体问题 操作： x0=0:0.1:1； y0=[-.447 1.978 3.11 5.25 5.02 4.66 4.01 4.58 3.45 5.35 9.22]; p=polyfit(x0,y0,3) 利用M文件编程： Clc; clear； x=0:0.1:2*pi; %生成样本点x y=sin(x)+0.5*rand(size(x)); %生成样本点y，通过随机矩阵 p=polyfit(x,y,3) %拟合出多项式（3阶） y1=polyval(p,x); %求多项式的值 plot(x,y,+,x,y1,-r) %绘制多项式曲线，以验证结果 了解 多项式插值是指根据给定的有限个样本点，产生另外的估计点以达到数据更为平滑的效果。该技巧在信号处理与图像处理上应用广泛。 插值的定义——是对某些集合给定的数据点之间函数的估值方法。 当不能很快地求出所需中间点的函数时，插值是一个非常有价值的工具。 Matlab提供了一维、二维、 三次样条等许多插值选择 所用指令有一维的interp1、二维的interp2、三维的interp3。这些指令分别有不同的方法（method），设计者可以根据需要选择适当的方法，以满足系统属性的要求。Help polyfun可以得到更详细的内容。 练习了解 y=interp1(xs,ys,x,’method’) nearest：执行速度最快，输出结果为直角转折 linear：默认值，在样本点上斜率变化很大 spline：最花时间，但输出结果也最平滑 cubic：最占内存，输出结果与spline差不多 %curve interpolation ys=[0 0.9 0.6 1 0 0.1 -0.3 -0.7 -0.9 -0.2]; %已有的样本点ys xs=0:length(ys)-1; %已有的样本点xs x=0:0.1:length(ys)-1;%新的样本点x y1=interp1(xs,ys,x,nearest); %插值产生新的样本点y1 y2=interp1(xs,ys,x,linear); %插值产生新的样本点y2 plot(xs,ys,‘+k’,x,y1,‘:r’,x,y2,‘-m’) %分别绘制不同方法产生的曲线； legend(sampled point,nearest,linear,spline,cubic)  * * 第3章 矩阵与数组 哈尔滨工业大学 2010年春季学期 《计算机辅助设计》 《Matlab/simulink教程及其在飞行器动态仿真中的应用》 航 天 学 院 航天工程系 高长生：综合楼414cortu