控制系统仿真第2章 控制系统分析.ppt

2．多项式的运算 MATLAB语言提供了大量关于多项式运算的内部函数。表2-10列出了常用的多项式运算函数。 【例2-27】多项式的运算。 A=[1 2 3 4]; 　 B=[3 -5 8 23 -7]; 　 roots(A) % 求多项式A的根。 　　 ans = -1.6506 -0.1747 + 1.5469i -0.1747 - 1.5469i polyder(B) % 多项式B求导。 　　 ans = 12 -15 16 23 　 conv(A,B) % 多项式A与B相乘。 　　 ans = 3 1 7 36 43 87 71 -28 3. 多项式的曲线拟合 在数值分析中，曲线拟合就是用解析表达式逼近离散数据，即离散数据的公式化。 【例2-28】 　　%输入自变量x数组 　　x = [0 3.3 4 5 6 6.4 7 7.4 8 8.6 9 10]; 　　％输入对应自变量x数组的多项式值域y数组 　　y = [6.016 5.608 5.359 5.360 5.470 5.420 5.350 5.338 5.260 5.205 5.110 5.000]; 　　％进行多项式拟合，定义多项式最高次幂为３ 　　n = 3; 　　p = polyfit(x, y, n); 　　xi = linspace (0, 10, 1000); 　　z = polyval (p, xi); plot (x, y, +b , x, y, r, xi, z, :g )　％绘制多项式拟合曲线，其中红色实线为实验数据绘制的曲线，绿色为采用多项式拟合方法绘制的曲线。 　 运算结果如图2-6所示。 　　 图2-6 曲线拟合curve fitting 多项式插值是指根据给定的有限个样本点，产生另外的估计点以达到数据更为平滑的效果。该技巧在信号处理与图像处理上应用广泛。所用指令有一维的interp1、二维的interp2、三维的interp3 【例2-29】一维插值 　　 x=[-2 1 5 10 20]; 　　 y=[1 9 11 20 24]; 　　 xi=3; 　　 yi=interp1(x,y,xi, linear) 　　 yi= 　　 10 2.7 MATLAB常用绘图命令 2.7.1 二维图形的绘制 二维图形的绘制是MATLAB语言图形处理的基础，常用的绘图命令是plot。 　1. 单曲线 　　语句格式plot(x,y) ：其中x和y为长度相同的向量，分别用于存储x坐标和y坐标数据。 　　 y=[0, 0.48,0.84,1,0.91,0.6,0.14]; 　　 plot(y) 生成的图形见图2-7，是以序号1、2、…、6为横坐标、数组y的数值为纵坐标画出的折线。 　 x=linspace(0,2*pi,30); % 生成一组线性等距的数值 y=sin(x); plot(x,y) 生成的图形见图2-8，是[0 2]上30个点连成的光滑的正弦曲线 图2-7 单曲线绘图实例1 图2-8单曲线绘图实例2 2．多重线 在同一个画面上可以画许多条曲线，只需多给出几个数组，语句格式为plot(x1,y1,x2,y2，…，xn,yn)。　　 多重线的另一种画法是利用hold命令。在已经画好的图形上，若设置hold on，MATLA将把新的plot命令产生的图形画在原来的图形上。而命令hold off 将结束这个过程。 【例2-30】 　　 t=0:pi/100:6; y=sin(t); plot(t,y) 先画好图2-9，然后用下述命令增加cos(x)的图形，也可得到图2-10。 　　 hold on 　　 y1=sin(t+0.25); plot(t,y1) 　　 y2=sin(t+0.5); plot(t,y2) 　　 hold off 图2-9多重线绘图实例1 图2-10多重线绘