lx1_2.m

matlab与控制系统仿真的PPt多媒体教程和源码

%lx1_2.m
clear
clc
close all
K=10;num=[1];den=[1 1 0];
[nums dens]=series(K,1,num,den);
w=logspace(-1,2,200);
[mag,phase,w]=bode(nums,dens,w);
%第一步：计算K=10时，未校正开环系统的相角裕度pm
[gm,pm,wcg,wcp]=margin(mag,phase,w);
%第二步：计算所需要的附加超前相角
phi_hope=55;
phi=(phi_hope-pm)*pi/180;%弧度单位
%第三步：根据sin(Qm)=(a-1)/(a+1),计算校正网络参数a
alpha=(1+sin(phi))/(1-sin(phi))
%第四步：计算10loga，在未校正系统的Bode图上，
%        确定与幅值增益－10loga对应的频率wm
%        通过绘制-10loga线来确定wm
M=-10*log10(alpha)*ones(length(w),1);
%生成一个随w变化的M直线
semilogx(w,20*log10(mag),w,M)
grid on
xlabel('frequency[rad/s]')
ylabel('mag[dB]')
hold on