实验五Matlab频域分析2013汇编.doc

实验 MATLAB频域特性分析系统的频率响应是在正弦信号作用下系统的稳态输出响应。对于线定常系统，在正弦信号作用下，稳态输出是与输入同频率的正弦信号，仅是幅值和相位不同。设系统传函为，其频率特性 例对系统，在输入信号和下可由求系统的输出信号》num=2;den=[1 2 3]; 》G=tf(num,den); 》t=0:0.1:6*pi; 》=sin(t);/ u=sin(3*t); 》y=lsi(G,u,t); 》plot(t,,t,y) 运行程序显示系统响应如图5-所示。 a) sin t的响应 b) sin (3t)的响应 图5-1 正弦信号输入系统的稳态响应 5.2用绘制极坐标图 频率特性中的奈奎斯特图是奈奎斯特（Nyquist）稳定性判据的基础。反馈控制系统稳定的充分必要条件为：奈奎斯特曲线逆时针包围点的次数等于系统开环右极点个数。 调用Matlab中nyquist() 函数可绘出奈奎斯特图，其调用格式为： 或 式中，；为用户提供的频率范围；re为极坐标的实部；im为极坐标的虚部。若不指定频率范围，则为。在输入指令中，如果缺省了左边的参数说明，奈奎斯特函数将直接生成奈奎斯特图；当命令包含左端变量时，即时，则奈奎斯特函数将只计算频率响应的实部和虚部，并将计算结果放在数据向量和中。在此情况下，只有调用plot函数和向量、，才能生成奈奎斯特图。 例5-2 设系统的传递函数为，绘制其奈奎斯特图。 解 程序如下： 》奈奎斯特，则绘制其奈奎斯特图的程序如下： 》num=[1000];den=[1,8,17,10]; 》nyquist (num,den);grid 或 num=[1000];den=[1,8,17,10]; sys=tf(num,den); nyquist (sys); grid 运行程序，显示奈奎斯特奈奎斯特奈奎斯特》v=[-10,0,-1.5,1.5]; 》axis(v) a) b) 图5-3 奈奎斯特局部图 例5-4 设某系统的开环传递函数为 则绘制其奈奎斯特图的程序如下： 》num=10*[1, 4, 4]; 》den=conv([1, 1],[1, -2,9]); 》nyquist (num,den); 》grid 或 num=10*[1,4,4]; 》den=conv([1,1],[1,-2,9]); sys=tf(num,den); nyquist (sys); grid 运行程序，显示得图5-4 a)。 若规定实轴、虚轴范围（10，10），（-10，10），则绘制其奈奎斯特图的程序如下： 》num=10*[1 4 4]; den=conv([1 1],[1 -2 9]); 》nyquist (num,den); 》axis([-10,10,-10,10]) 运行程序，显示得图5-4 b)。 a） b） 图 5-4 奈奎斯特曲线 5.3 用bode(sys)画对数坐标图 伯德图由对数幅频和对数相频两张图构成轴采用对数分度，幅值为对数增益即分贝()，相位为线性分度。Matlab中绘制伯德图的函数为bode( )，其调用格式为 [mag,phase,]=bode(num,den,w) 或 sys=tf(num,den); bode(sys) 式中，G(s)=num/den频率自动选择范围从，若选择频率范围，可应用logspace()函数，其格式为 式中，表示最小频率10ab表示最大频率10b，n表示10a～10b之间频率点数。 ，则绘制其伯德图的Matlab程序如下: num=5*[0.1 1]; den=conv([1 0],conv([0.5 1],[1/2500 0.6/50 1])); bode(num,den) 函数 的应用。程序如下： w=logspace(-1,4,300); %确定频率范围及点数 [mag,phase,w]=bode(n