% ZV 两脉冲输入整形 仿真
wn = 30; % 无阻尼固有频率
b = 0.2; % 阻尼比
G1 = tf(wn*wn,[1 2*wn*b wn*wn]); % 原系统
step(G1)

wn = 30; % 模型失配参数
b = 0.18;
K = exp(-b*pi/(sqrt(1-b*b)));
T = pi/((wn*sqrt(1-b*b)));
A1 = 1/(1+K);
A2 = K/(1+K);

Index = 15; % 时间因子
TN = 500; % 线条点数
TNm = Index*TN; % 线条总点数
u1(1:TN,1) = A1; % 第一个脉冲
u1(TN+1:TNm,1) = A2+A1; % 第二个脉冲
Ts = linspace(0,Index*T,TNm)';
y = lsim(G1,u1,Ts);
y1 = step(G1,Ts);

plot(Ts,u1,Ts,y,Ts,y1);
title('阶跃响应');legend('输入整形脉冲','未加输入整形','输入整形');
axis([0 1 0 1.6]);