已知系统的微分方程，求系统的零极点和单位冲激响应和频率响应曲线的MATLAB代码

好的，以下是求解系统零极点、单位冲激响应和频率响应曲线的MATLAB代码：

假设系统的微分方程为：

dy/dt + 2y = 3u(t)

其中u(t)为单位阶跃函数。

1. 求解系统的零极点：

% 定义系统的传递函数
sys = tf([3], [1 2]);

% 求解系统的零极点
[z, p] = zero(sys);
disp('系统的零点为：');
disp(z);
disp('系统的极点为：');
disp(p);

2. 求解系统的单位冲激响应：

% 定义单位冲激函数
t = 0:0.01:10;
u = zeros(size(t));
u(1) = 1;

% 求解系统的单位冲激响应
h = impulse(sys, t);

% 绘制单位冲激响应曲线
figure;
plot(t, h);
xlabel('时间');
ylabel('系统响应');
title('系统的单位冲激响应');

3. 求解系统的频率响应曲线：

% 定义系统的频率范围
w = logspace(-2, 2, 1000);

% 求解系统的频率响应曲线
[mag, phase] = bode(sys, w);

% 绘制频率响应曲线
figure;
subplot(2,1,1);
semilogx(w, 20*log10(mag));
grid on;
xlabel('频率(rad/s)');
ylabel('幅值(dB)');
title('系统的频率响应曲线');
subplot(2,1,2);
semilogx(w, phase);
grid on;
xlabel('频率(rad/s)');
ylabel('相位(deg)');

运行以上代码后，即可求得系统的零极点、单位冲激响应和频率响应曲线。