example6.m

来自「控制系统的matlab程序六个。含有波特图」· M 代码 · 共 12 行

i=0;
y=zeros(1,400);
t=[0:0.1:39.9]';
den=conv([1 0],conv([0.2 1],[3 1]));
for k=1:0.1:2.9
    g0=tf(k,den);
    f=feedback(g0,1);
    i=i+1;
    y(i,:)=step(g,t);
end
surf(y)
shading flat