matlab的复频域分析

这个代码怎么修改啊，实在不会啊，哪位朋友帮帮我吧！😭
syms t s %求零输入响应
Yzis=(3s+13)/(s^2+3s+2);
yzi=ilaplace(Yzis);
xt=4exp(-2t)Heaviside(t);%求零状态响应
Xs=laplace(xt);
Yzss=Xs/(s^2+3s+2);
yxs=ilaplace(Yzss);
y=yzi+yxs;%求全响应
b=[0,1];a=[1,3,2];%求系统函数及零极点
sys=tf(b,a);
subplot(221);pzmap(sys);axis([-3,3,-2,2]);
subplot(222);impluse(b,a);%系统函数极点对时域波形的影响
w=-8pi:0.01:8pi;%画出系统的幅频和相频特性曲线
H=freqs(b,a,w);
subplot(223);plot(w,abs(H)),grid on
xlabel('/omega(rad/s)');
ylabel('|H(\omega)|');
title('幅频特性')
subplot(224);plot(w,angle(H)),grid on
xlabel('/omega(rad/s)');
ylabel('|phi(\omega)|');
title('相频特性')