实验三 傅里叶变换

  • 格式: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),验证傅立叶变换的相关 特性