Matlab---基本

求零状态响应

系统的微分方程为
y”(t) + 2y’(t) + 77y(t) = f(t)
在t≥0接入激励f(t)=10sin(2πt)，求零状态响应。

sys = tf([1],[1 2 77]);
t = 0:0.01:5;
f = 10*sin(2*pi*t);
y = lsim(sys, f, t);
plot(t,y);
xlabel('Time(sec)')
ylabel('y(t)')

系统的冲激响应和阶跃响应

求以下系统的冲激响应和阶跃响应。
7y”(t) + 4y’(t) + 6y(t) = f’(t)+ f(t)

a=[7 4 6]; 
b=[1 1]; 
subplot(2,1,1)
impulse(b,a); 
subplot(2,1,2)
step(b,a);