《信号与系統》仿真实验 实验指示书要点 实验三 信号的频谱分析
1 《信号与系统》仿真实验 实验指示书要点 实验三 信号的频谱分析
系统仿真的要点 关于快速傅立叶变换的t函数的 调用 频谱分析举例
2 系统仿真的要点 • 关于快速傅立叶变换的fft函数的 调用 • 频谱分析举例
、数的调用 ft(xDFT的快速算法 (1)数据长度为N; 若N为2的整数次幂,则是基2的 DFT,运算速度较快; 若N不为2的整数次幂,则运算速 度较慢。 (2)若数据实际长度小于N,则自动补 些零,使之长度等于N
3 一、fft函数的调用 • fft(x)—DFT的快速算法 （1） 数据长度为N； 若N为2的整数次幂，则是基-2的 DFT，运算速度较快； 若N不为2的整数次幂，则运算速 度较慢。 （2）若数据实际长度小于N，则自动补 一些零，使之长度等于N
复指数座标 正变换 X(k+1)=∑x(m+1W减 反变换 x(m+1)=∑X(k+1) 其中
4 • 复指数座标 – 正变换 – 反变换 其中 N n nk WN X k x n 1 ( 1) ( 1) N n nk WN x n X k 1 ( 1) ( 1) N j N W e 2
ff(x或f(x,N)的调用编程 t=(0:1/255:1); x=exp(-0.9)*sin(2 pi 60*t) y=fft(x) ry-real(y a=abs(y) plot(t, ry, r,t, a, b-") grI d 这里是归一化频率
5 • fft(x) 或fft(x,N)的调用编程 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 -20 -10 0 10 20 30 40 50 t=(0:1/255:1); x=exp(-0.9)*sin(2*pi*60*t); y=fft(x); ry=real(y); a=abs(y); plot(t,ry, 'r-' ,t,a, 'b--'); grid 这里是归一化频率