傅里叶变换和非周期信号的频谱
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第10章傅里叶变换和非周期信号的频谱
10.1利用fourier 函数求下列信号的傅里叶变换F (j ω),并利用ezplot 函数绘出其幅度频谱|F (j ω)|和相位频谱φ(ω)。观察比较三个信号的幅频特性和相频特性,并利用傅里叶变换的性质加以解释。
(1)t
t t f ππ2)2sin()(1= (2))
2(2)]2(2sin[)(2--=t t t f ππ
(3)23]
2)
2sin([)(t t t f ππ=
(1)syms t im re phase;
f = sin(2*pi*t)/(2*pi*t);
Fw = fourier(f)
subplot(311);
ezplot(f);
axis([-2.5 2.5 -pi pi]);
xlabel('时域波形');
subplot(312);
ezplot(abs(Fw));
axis([-2.5 2.5 -pi pi]);
xlabel('幅度谱');
im = imag(Fw);
re = real(Fw);
phase = atan(im/re)
subplot(313);
ezplot(phase);
axis([-2.5 2.5 -pi pi]);
xlabel('相位谱');
(2)
syms t im re phase;
f =
sin((2*pi*(t-2))/(2*pi*(t
-2));
Fw = fourier(f)
subplot(311);
ezplot(f);
axis([-2.5 2.5 -pi pi]);
xlabel('时域波形');
subplot(312);
ezplot(abs(Fw));
axis([-2.5 2.5 -pi pi]);
xlabel('幅度谱'); im = imag(Fw);
re = real(Fw);
phase = atan(im/re) subplot(313);
ezplot(phase);
axis([-2.5 2.5 -pi pi]); xlabel('相位谱');
(3)
syms t im re phase;
f =
[sin(2*pi*t)/(2*pi*t)]^2; Fw = fourier(f)
subplot(311);
ezplot(f);
axis([-2.5 2.5 -pi pi]); xlabel('时域波形'); subplot(312);
ezplot(abs(Fw));
axis([-2.5 2.5 -pi pi]); xlabel('幅度谱');
im = imag(Fw);
re = real(Fw);
phase = atan(im/re) subplot(313);
ezplot(phase);
axis([-2.5 2.5 -pi pi]); xlabel('相位谱');