实验三 傅里叶变换
- 格式:doc
- 大小:30.50 KB
- 文档页数:2
实验三 傅里叶变换
一 实验目的
(1) 掌握连续时间信号傅立叶变换的数值计算方法;
(2) 熟悉基本信号的频域转换;
(3) 熟悉傅立叶变换的性质。
二 实验原理
1. 连续时间信号傅立叶变换的数值计算:
信号 f(t)的傅立叶变换定义为:
F(jw)= ∫ ∞ ?∞ f (t )e ? jwt dt = lim τ →0 n = ?∞ ∑ f ( nτ ) e
∞ ? jwn τ τ (3.1)
当取 τ 足够小时,上式的近似情况可以满足实际需要。若信号 f(t)是时限的,或当 |t|大于某个给定值时,f(t)的值已衰减得很厉害,可以近似地看成时限信号时,则(3.1) 式中的 n 可取有限值,设为 N,并对(3.1)式中的频率 w 进行取样,设
wk = 2π k Nτ N ?1 n=0 (3.2) ? jwk nτ
F (k ) = τ ? ∑ f (nτ )e , 0≤k≤N (3.3)
用 MATLAB 实现时,其要点是要正确生成 f(t)的 N 个样本 f(nτ)的向量 f,
及向量 e ? jwk nτ ,两向量的内积(即两矩阵的乘积)结果即完成式(3.3)的计算。
2. MATLAB 的 Symbolic Math Toolbox 提供了能直接求解傅立叶变换及逆变换的 函数 fourier(f)和 ifourier(F)。
在调用函数 fourier()和 ifourier()之前,要用 syms 命令对所用到的变量进行说明, 即要将这些变量说明成符号变量。对 fourier()中的函数 f 及
ifourier()中的函数 F,也要 用符号定义符 sym 将 f 或 F 说明为符号表达式; f 或 F 是 MATLAB 中的通用表达式, 若 则不必用 sym 加以说明。
三 实验内容
1..已知门信号 f(t)=g2(t)= 0 |t|>1 |t|<1 1, 利用数值计算的方法,求其傅立叶变换F(jw)。
实现该过程的 MATLAB 命令程序如下:
R=0.01;t=-2:R:2;
f=Heaviside(t+1)-Heaviside(t-1);
W1=2*pi*5; %频率带宽
N=500;k=-N:N;W=k*W1/N;
F1=f*exp(-j*t'*W)*R;
F=real(F1);
subplot(2,1,1);plot(t,f)
xlabel('t');ylabel('f(t)'); title('f(t)=u(t+1)-u(t-1)');
axis([-2 2 0 1.2])
subplot(2,1,2);plot(W,F)
xlabel('W');ylabel('F(w)');
title('f(t)的傅氏变换 F(w)');
axis([-40 40 -0.5 2])
2. 试用 fourier()函数求下列信号的傅立叶变换 F(jw), 并画出|f(jw)|。
(1). f(t)=te-3tu(t);
(2) f(t)=sgn(t)= 1, t>0 -1, t<0
3. 试用 ifourier()函数求 F ( jw) = ? 4. f1(t)的波形如下图所示。 f(t) 1
2ω 的逆傅立叶变换并画出波形。
4 +ω2 2 -3 0 3 t f(t)=f1(t-2)cos(100t), 试用连续信号的傅立叶变换数值算法,求:
(1) f1(t)的傅立叶变换 F1(jw)的|F1(jw)|及 φ1(w);
(2) f1(t-2)的傅立叶变换 F2(jw)的| F2(jw)|及 φ2(w);
(3) f(t)的傅立叶变换 F(jw) 的|F(jw)|及 φ(w);
分别比较|F1(jw)|、|F2(jw)|、|F(jw)|及 φ1(w)、φ2(w)、φ(w),验证傅立叶变换的相关 特性