信号与系统实验 01 - | 连续系统频域分析的MATLAB实现

Posted 2021-07-14 Neutionwei

一、实验目的

1. 掌握连续时间信号频谱特性的MATLAB分析方法；

2. 掌握连续系统的频率响应MATLAB分析方法方法。

二、实验原理

1. 连续时间信号的频谱---傅里叶变换

syms t w f
f=exp(-2*t)*sym('heaviside(t)');
F=fourier(f)
subplot(3,1,1);ezplot(f,[0:2,0:1.2]);
subplot(3,1,2);ezplot(abs(F),[-10:10]);
subplot(3,1,3);ezplot(angle(F),[-10:10]);

syms t w
F=1/(1+w^2);
f=ifourier(F,w,t)
ezplot(f)

2. 连续系统的频域分析

三、实验内容

syms t w f F y Y s S
f=exp(-t)*sym('heaviside(t)');
F=fourier(f)
subplot(3,3,1);ezplot(f);title('波形');
subplot(3,3,2);ezplot(abs(F),[-20:20]);title('f:幅频特性');
subplot(3,3,3);ezplot(angle(F),[-20:20]);title('f:相频特性');
s=cos(10*t);
x=f*s;
X=fourier(x,t,w)
subplot(3,3,4);ezplot(x);title('x:波形');
subplot(3,3,5);ezplot(abs(X),[-20:20]);title('x:幅频特性');
subplot(3,3,6);ezplot(angle(X),[-20:20]);title('x: 相频特性');
H=sym('heaviside(w+10)-heaviside(w-10)');
Y=X*H;y=ifourier(Y,w,t)
subplot(3,3,8);ezplot(abs(Y),[-20:20]);title('y: 幅频特性');
subplot(3,3,9);ezplot(angle(Y),[-20:20]);title('y: 相频特性');

b=[1,0];
a=[1,5,6];
freqs(b,a,0:0.01:1000);