ex3_4.asv

我们上课的资料

% 关闭所有图象窗口并清除内存变量
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,der3));
% 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=-10:1:10;
% K=2时系统的开环Nyquist曲线
figure(1);
nyquist(numo1, deno1,wt);
title('System Nyquist Charts with K=2')
grid;
figure(2);
step(numc1,denc1,t); 
title('System TimeResponsewithK：2／)；


    ％K：20时系统的开环Nycluist曲线
    nyquist(num02，den02，wt)；
    title('SystemNyquistChartswithK；20／)
    grid；
    ％系统的闭环阶跃响应
    y1：step(numcl，dencl，t)；
    y2：step(numc2，denc2，t)；    ，
    ％绘图
    subplot(211)，，
    Plot(t，y1)；
    title('System TimeResponsewithK：2／)；
    xlabel(，Time—8xr)；
    ylabel('Response—value')；
    grid；
    subplot(212)，
    少t(t，y2)；
    title('System TimeResponsewithK：20，)；
    ?xlabel('Time—sect)；
    ylabel(，Response—valuet)；
    grid；