Question

Matlab程序注释

（1）num=[2547];
den=[1,2005];
g=tf(num,den);
phai=feedback(g,1);
num=[571.2];
den=[1,0];
gc=tf(num,den);
ggc=gc*g;
phaic=feedback(ggc,1);
bode(ggc);
grid on;
（2）num=conv([800],[1/40 1]);
den=conv([1 0],conv([1/3/pi 1],[1/1000 1]));
a=tf(num,den);
phaic=feedback(a,1);
t=0:0.01:3;
u=pi/6*sin(3*pi*t);
lsim(phaic,u,t);
grid on;
希望注释的详细一些，谢谢啦

Answer 1

（1）num=[2547];
den=[1,2005];
g=tf(num,den); %构造传递函数g=2547/(s+2005)
phai=feedback(g,1); %引入反馈phai=2547/(s+2005+2547)
num=[571.2];
den=[1,0];
gc=tf(num,den); %gc=571.2/s,积分环节
ggc=gc*g;
phaic=feedback(ggc,1); %phaic=1455000/(s^2 + 2005 s + 1455000)
bode(ggc); 画bode图
grid on;
（2）num=conv([800],[1/40 1]); %conv做卷积，这里实际上就相当于800*[1/40 1]
den=conv([1 0],conv([1/3/pi 1],[1/1000 1]));
a=tf(num,den);
phaic=feedback(a,1);
t=0:0.01:3; %时间
u=pi/6*sin(3*pi*t); %输入u(t)
lsim(phaic,u,t); %系统响应图
grid on;
应该没什么难理解的阿，很简单的东西