北京科技大学电工实验报告.docxVIP

北京科技大学电工实验报告 　　实验四线性系统的频域分析　　一.实验目的　　熟练掌握运用MATLAB命令绘制控制系统伯德图的方法；　　了解系统伯德图的一般规律及其频域指标的获取方法；　　熟练掌握运用伯德图分析控制系统稳定性的方法；　　设计超前校正环节并绘制Bode图；　　设计滞后校正环节并绘制Bode图。　　二.实验原理及内容　　1、频率特性函数G(j?(来自:写论文网:北京科技大学电工实验报告))。　　设线性系统传递函数为：　　b0sm?b1sm?1?????bm?1s?bmG(s)?nna0s?a1s?????an?1s?an　　则频率特性函数为：　　b0(j?)m?b1(j?)m?1?????bm?1(j?)?bmG(jw)?nna0(j?)?a1(j?)?????an?1(j?)?an　　由下面的MATLAB语句可直接求出G(jw)。　　i=sqrt(-1)%求取-1的平方根　　GW=polyval(num，i*w)./polyval(den，i*w)　　其中为系统的传递函数模型。而w为频率点构成的向量，点右除运算符表示操作元素点对点的运算。从数值运算的角度来看，上述算法在系统的极点附近精度不会很理想，甚至出现无穷大值，运算结果是一系列复数返回到变量GW中。　　2、用MATLAB作伯德图　　控制系统工具箱里提供的bode()函数可以直接求取、绘制给定线性系统的伯德图。当命令不包含左端返回变量时，函数运行后会在屏幕上直接画出伯德图。如果命令表达式的左端含有返回变量，bode()函数计算出的幅值和相角将返回到相应的矩阵中，这时屏幕上不显示频率响应图。命令的调用格式为：　　[mag,phase,w]=bode(num,den)　　[mag,phase,w]=bode(num,den,w)　　或　　[mag,phase,w]=bode(G)　　[mag,phase,w]=bode(G,w)　　矩阵mag、phase包含系统频率响应的幅值和相角，这些幅值和相角是在用户指定的频率点上计算得到的。用户如果不指定频率w,MATLAB会自动产生w向量，并根据w向量上各点计算幅值和相角。这时的相角是以度来表示的，幅值为增益值，在画伯德图时要转换成分贝值，因为分贝是作幅频图时常用单位。可以由以下命令把幅值转变成分贝：magdb=20﹡log10(mag)　　绘图时的横坐标是以对数分度的。为了指定频率的范围，可采用以下命令格式：logspace(d1，d2)　　或　　logspace(d1，d2，n)　　第一种格式是在指定频率范围内按对数距离分成50等分的，即在两个十进制数　　向量中的点数50是一个默认?1?10d和?2?10d之间产生一个由50个点组成的分量，12　　值。例如要在?1?弧度/秒与?2?100弧度/秒之间的频区画伯德图，则输入命令时，d1?log10(?1),d2?log10(?2)在此频区自动按对数距离等分成50个频率点，返回到工作空间中，即　　w=logspace(-1，2)　　要对计算点数进行人工设定，则采用第二种格式。例如，要在?1?1与?2?1000之间产生100个对数等分点，可输入以下命令：　　w=logspace(0，3，100)　　在画伯德图时，利用以上各式产生的频率向量w，可以很方便地画出希望频率的伯德图。由于伯德图是半对数坐标图且幅频图和相频图要同时在一个绘图窗口中绘制，因此，要用到半对数坐标绘图函数和子图命令。　　对数坐标绘图函数　　利用工作空间中的向量x，y绘图，要调用plot函数，若要绘制对数或半对数坐标图，只需要用相应函数名取代plot即可，其余参数应用与plot完全一致。命令公式有：　　semilogx(x，y，s);只对x轴进行对数变换，y轴仍为线性坐标。　　semilogy(x，y，s);y轴取对数变换的半对数坐标图。　　Loglog(x，y，s);全对数坐标图，即x轴和y轴均取对数变换。　　子图命令　　MATLAB允许将一个图形窗口分成多个子窗口，分别显示多个图形，这就要用到subplot()函数，其调用格式为：　　subplot(m，n，k)　　该函数将把一个图形窗口分割成m×n个子绘图区域，m为行数，n为列数，用户可以通过参数k调用各子绘图区域进行操作，子图区域编号为按行从左至右编号。对一个子图进行的图形设置不会影响到其它子图，而且允许各子图具有不同的坐标系。例如，subplot(4，3，6)则表示将窗口分割成4×3个部分。在第6部分上绘制图形。MATLAB最多允许9×9的分割。　　【范例6-3】给定单位负反馈系统的开环传递函数为：G(s)?10(s?1)，试画出伯德图。s(s?7)　　利用以下M