ex3_4.m

我们上课的资料

% 关闭所有图象窗口并清除内存变量
close all;
clear all;
% K=2时系统的开环传递函数
numo1=[0 0 0 2];
den1=[1 0];
den2=[1 1];
den3=[0.1 1];
deno1=conv(den1,conv(den2,den3));
% K=20时系统的开环传递函数
numo2=[0 0 0 20];
deno2=deno1;
% 系统的闭环传递函数
numc1=numo1;
denc1=deno1+numo1;
numc2=numo2;
denc2=deno2+numo2;
% 仿真时间及角频率初始化
t=0:0.1:15;
wt=0:10;
% K=2时系统的开环Nyquist曲线和闭环阶跃响应
figure(1);
nyquist(numo1, deno1,wt);
title('System Nyquist Charts with K=2')
figure(2);
step(numc1,denc1,t); 
title('System Time Response with K=2');
grid;
% K=20时系统的开环Nyquist曲线和闭环阶跃响应
figure(3);
nyquist(numo2, deno2,wt);
title('System Nyquist Charts with K=20')
figure(4);
step(numc2,denc2,t); 
title('System Time Response with K=20');
grid;