如何在不修改原传递函数的前提下绘制零点多于极点系统的阶跃响应

问题原因 #

你定义的传递函数 G(s) = -RC·s 是纯微分环节，属于非因果系统（零点数1个 > 极点数0个），MATLAB内置的step()函数仅支持仿真因果系统的时间响应，因此触发报错。

不修改传递函数的解决方案 #

纯微分器的阶跃响应是冲激函数（y(t) = -RC·δ(t)），无法通过step()直接绘制，需手动构造响应并可视化：

方法1：使用符号工具箱绘制冲激示意图 #

% 定义参数
R = 47*10^3;
C = 1*10^-6;
RC = R*C;

% 符号化定义阶跃响应
syms t
y = -RC * dirac(t);

% 绘制冲激函数
fplot(y, [-0.1 0.1]);
title('微分器阶跃响应（冲激函数）');
xlabel('时间 t (s)');
ylabel('输出 y(t)');
grid on;

方法2：用窄脉冲近似冲激（工程可视化） #

如果不需要严格的符号冲激，可通过极窄的矩形脉冲近似：

% 定义参数
R = 47*10^3;
C = 1*10^-6;
RC = R*C;

% 构造时间轴
t = linspace(-0.1, 0.1, 10000);
y = zeros(size(t));
% 在t=0附近设置窄脉冲，近似冲激
delta_t = 1e-4;
idx = abs(t) < delta_t/2;
y(idx) = -RC / delta_t; % 保证脉冲面积等于冲激强度

% 绘制近似响应
plot(t, y);
title('微分器阶跃响应（冲激近似）');
xlabel('时间 t (s)');
ylabel('输出 y(t)');
ylim([-500, 100]); % 调整y轴范围以看清脉冲
grid on;

补充说明 #

纯微分器的阶跃响应本质是冲激函数，这是因为阶跃信号的导数就是冲激信号。实际工程中不存在理想微分器，通常会加入惯性环节使其成为因果系统，但你要求不修改原传递函数，因此只能通过上述方式可视化理想微分器的响应特性。

内容的提问来源于stack exchange，提问作者user10799430