傅里叶变换及应用Word格式.docx

1、定理设齐（”是以HOCvs ）为周期的实函数、且在I 2d丿上满足狄利克雷条件、即齐（-在一个周期上满足：（1）连续或只有有限个第一类间断点；（2）只有有限个极值点.则在连续点处J有cos n cot + b/j1 -%二不住方CCOS讪d= 12）r單（用吨*）在间断点4处，（1）式右端级数收敛于2/inCXf bi CiXcoso S s2y于島 ex T ur,ee+ 0 C - n f “心 * + “人-斶2 2-处）晋+Ee +e厶W=1旦厂_ 5 叽_ +也0叩 ”T,2,3,则XaW二工c严“X 42CtX+丿+（J0-inCX, +c”+Q6XdlCX=Co+（qg +C2W

2、 +代JnCX9+%+丿（2）式称为傅里叶级数的复指数形式，具有明显的物理意义.容易证明5可以合写成一个式子，即q 冷直齐（Wgai,2,.）傅里叶积分任何一个非周期函数/”）都可看成是由某个周期函数齐（）当T-+8时转化而来的.即imx/T（/） =/（r）由公式（2）、（3）得T R 了耐吨严9可知/（r）=71iml 心（少叫卜/ w=x T令马二gg = _（DI则T或A于是f（t）= lim 八）T卄TJArYdt/?=-xuf w=-x 0 4r t令0 （）=2 吐 CH dre（ ,H故（4）03” ）T女）二土匚/ （少叫（1）严注意到当即TTS 时从而按照积分的定义，（4）

3、可以写为:或者/（/）=2匚匚/ （少叫咔叫0公式（5）称为函数/（ ）的傅氏积分公式.定理若/ （*）在（揺，+8）上满足条件：仃）/ （”在任一有限区间上满足狄氏条件；（2） / C）在无限区间（-8, +8）上绝对可 积，即收敛，则在几）的连续点成里；而在/ （“的间断点匚处/ （ /。+0） +/亿一0）应以 来代替.上述定理称为傅氏积分定理可以证明，当/”）满足傅氏积分定理条件时，公式（5）可以写为三角形式，即一x/（r）cos） ； ylabel（* y/m ）；Subplot（212） W=fft（w）；W=fftshift（W）；plot（t,abs（W） /br, 1 Lin

4、eWidth,;title （冲激函数的傅里叶变换）；xlabel （； ylabel（A/m*）;其时域图像和频域图像如图1所示倔1 4C-图1冲激函数的时域和频谱图像分析:从图中可以看出，冲激信号的频率为0处的分量最大，然后向两端快速衰减，表明脉冲信号 中实际占主导地位的其实是直流分量。余弦信号我们已经知道，任何信号都可以分解成为不同频率的正或余弦信号的叠加，那么现在研究余弦 信号的时域和频域特性（3）。用Mat lab可以产生余弦信号并分析其频谱的特性。Mat lab 程序：M=10 ；N=2 ；t=linspace（-10,10,N） xcos=cos（3*t） Jsubplot（21

5、1）plot（t,xcos）；title C余弦信号的时域图像）；xlabel（ t/sr）;ylabel（ y/m*） subplot（212）plot（ t ,abs（ fftshi ft（ fft（xcos） ； t i 11 eC 余弦信号的频域图像）xlabel（w/（rad/s） 余弦信号的时域图像与频域图像如图2所示 M 鼻 j 2 0 2A &9 10图2余弦函数的时域和频谱频率突变信号频率突变信号在现实生活总很常见，下面用Mat lab来产生频率突变信号和分析其傅里叶变 换。%Mat lab 程序:t=linspace（-10,10,N）:sl=find（t= 0）；x（s2

6、） =cos（2*pi*3*t（s2） subplot （211）；plot（l.x）；titleC频率突变信号）； xlabel（rt/s*）；ylabel（? y/m1）subplot（212） X= fft（x）；plot（trabs（X） Ititled频率突变信号的傅里叶变换图像*）； xlabel（r f/hz*）;ylabeK y/m*）其图像如图3所示图3频率突变信号的时域和频谱图象分析：频率突变信号的频率在3和5的位置对应的幅值特别高。因此标记出这两 个频 谱峰值对应的频率分量，正好可以验证信号的频率成份。高斯信号在信号中，常会伴随着来声，而高斯噪声是常见的棗声，研究它的特性

7、对于消除噪声有很大 的意义。Mat lab程序如下:M 二 10 ；N=2AM ；t=linspace（-10,10, N）；a=l/4 ；g=exp（-a*t. *2）;plot（t,g）titleC高斯信号的时域图像）；xlabelCt/s1）;y/m ）；subplot（212）G=fft （g）:G=fftshift（G）:plot （t.abs（G） titleC高斯信号的频域图像J xlabel（* f/Hz J ； ylabel（* y/m*）；高斯信号的时域和频域图像如图4所示-5 4 2 0 2 4 6 S 10 高廟馆号隔嶽tEBd 44a4 6dx图4高斯信号的时域和频域

8、图像 ，WX图像分析：这是一个正态分布函数，具有单峰性，归一性。其傅立叶变换函数的图象中，只有频率为0 的地方有极大的峰值，说明小概率时间发生的机会是极小的，越向原点，时间发生的可能性越 大。随机序列研究随机序列住有很大的意义，在数字信号的传输过程中，往往会产生噪声，而噪声并是随 机序列研究其特性对消除棗声有很大的意义利用MATLAB很容易产生两类随机信号:Rand（hN）在区间0.1上产生N点均匀分布的随机序列Randn（hN）r生均值为0,方差为1的高斯随机序列，也就是白噪声序列例如下图表示点数为32点的均匀分布的随机序列与高斯随机序列，其血t lab仿真结果如图下 所示，其中图和图分别表

9、示序列一和序列二的时域和频域图像。用Mat 1 ab产生的随即序列和其傅里叶变换的程序如下图所示clear al 1 ；N二32 ；x rand=rand（1,N）:x randn=randn（1,N）；xn=0 :figured）subplot （2,1,1）; stem（xn, x rand） ； titleC系列 1 的时域图像subplot （2.1,2）； stem（xniabs（fftshift （fft （x rand）: t i tle（!系列 1 的频域 图像 Jfigure （2） stem（xn,x_randn）: t i t le（1 系列 2 的时域图像）subplo

10、t （2.1,2） ； stem（xn（abs（fftshift （fft （x randn）: title（* 系列 2 的频域图像J系列1的时域羽惊10$0.60 102CjJ（Vc?7:T） 5 10 15 20 25 30 35系列1的腹域图像151 1 11055 10202530图序列一的时域和频域图像系列啲时域图像系列2的顿域囱像$厂64oTIG ?ITQ O0 6 10 15 20 26 30 36图序列二的时域和频域图像男女声音的辨别男女声音有什么区别，怎么样区别男女声音，下面用傅里叶变换分析一段男女芦音，观察结 果，并分析。首先用录音器分别录两段男女声音音频，分别保存为和。

11、在Mat lab中用Y, Fs=wavread（ f i 1 ename）即可直接读取 Mat 1 ab程序N=5000 ；aAwavreadC ,N）；bAwavreacK*.N）；A=fftshift（fft（a）；B 二 fftshift（fft（b）；subplot（223） plot（abs（A） ； titleC女生清唱傅里叶变换）subplot（224）:plot （abs （B） ； ti tie （男生清唱傅里叶变换） subplot（221）；plot （a） ； ti11 eC 女声subplot（222） plot（b） * title（男声）运行结果如图6所示由于男女声音音高不同，其频谱属于不同频段范围。因此可以用上述方法才辨别男女声音。如图所示，将直流分量至于中间，明显可以看出女声的高频分量多，男声的低频分量少。由此可以得出，此分析结果正确。女声6040A女生清唱傅里叶变挽 U20004000 6000图6男女声频谱分析0 2000 4000 6000参考文献:1 李红著复变函数与积分变换北京：高等教育出版杜，19992张元林积分变换M 北京：高等教育出版社,2003.3胡广书数字信号处理一一理论.算法与实现M 清华大学出版社，1997.4飞思科技产品研发中心.Matlab7辅助信号处理技术与应用M. 2005.