实验内容
syms t
y1=cos(10*t)*heaviside(t)
y2=exp(-t)*cos(10*t)*heaviside(t)
Y1=laplace(y1) %拉普拉斯变化
Y2=laplace(y2) %拉普拉斯变化
figure(1); %第一个窗口
a1=[0 1 0];b1=[1 0 100];%a1为零点函数 b1为极点函数
zplane(a1,b1);
title('y1');
figure(2);%第一个窗口
a2=[0 1 1];b2=[1 2 101] %a1为零点函数 b1为极点函数
zplane(a2,b2);
title('y2');
%收敛半径不变,收敛域会改变
syms s w t j
h(s)=s^2/[(s+1)^2+1];
figure(1);%第一个窗口
a1=[1 0 0];b1=[1 2 2];%a1为零点函数 b1为极点函数
zplane(a1,b1);
%画h(t)
figure(2);
h(t)=ilaplace(h(s)) %拉普拉斯逆变化
ezplot(h(t));
%画h(jw)
figure(3);
w=0:0.01:20;
H=freqs(a1,b1,w);
plot(w,abs(H));
title('h(jw)');
%系统为高通