matlab 用fft作谱分析,matlab关于用fft作谱分析

matlab关于用fft作谱分析

Matlab 关于用 FFT 作 谱分析中文摘要： matlab是一种用于算法开发、数据可视化、数据分析以及数值 计算的高级技术计算语言和交互式环境。本文利用matlab的强大计 算功能和信号数据处理功能，用FFT对连续信号和时域离散信号进 行谱分析。 英文摘要：matlab is used in algorithm development, data visualization, data analysis and numerical calculation of high-level technical computing language and interactive environment.。In this paper, matlab powerful computing capabilities and signal data processing functions, using FFT of a continuous signal and the time-domain discrete signal spectrum analysis. 关键词：FFT、连续信号、时域离散信号、matlab 一、引言本文研究在matelab中用FFT对连续信号和时域离散信号进行 谱分析的方法，了解可能出现的分析误差及其原因。除了数字信号 变换和滤波的功能外，matlab还具有很强大的图像处理功能，数学 计算功能等。而且不光是在物理学方面，在数学，化学，生物学， 经济学，甚至是社会学也有很多应用。2、傅立叶变换 1、离散傅立叶变换(DFT)离散傅立叶级数变换是周期序列，有N个独立的数值，所以他的 许多特性可以通过有限长序列沿拓来得到。对于一个长度为N的有限 长序列x(n)，也即小(n)只在0～(N-1)各点上有非零值，即        其他 , 0 1 0 ), ( ) ( N n n x n x 把序列x(n)以N为周期进行周期沿拓得到周期序列 ，则有 ) ( ~ n x        其他 , 0 1 0 ), ( ~ ) ( N n n x n x 所以，有限长序列x(n)的离散傅立叶变换(DFT)为          1 0 1 0 , ) ( )] ( [ ) ( N n kn N N k W n x n x DFT k X 逆变换为          1 0 1 0 , ) ( 1 )] ( [ ) ( N n kn N N k W k X N k X IDFT n x 编制信号产生子程序，产生以下典型信号供谱分析：    n R n x 4 1  ， 1  n 3 0   n，    n x 2 n  8 7 4   n， 其它 n 0， n  4 3 0   n，    n x 3 n  3 7 4   n 0 ， 其它 n   n n x 4 cos 4    n n x 8 sin 5    t t t t x    20 cos 16 cos 8 cos 6    应当注意，如果给出的是连续信号 ，则首先要根据其最高频率   t x a 确定采样速率 以及由频率分辨率选择采样点数 N，然后对其进行 s f 软件采样(即计算 ， )，产生对应序列 。对     nT x n x a  1 0    N n   n x 信号 ，频率分辨率的选择要以能分辨开其中的三个频率对应的   t x 6 谱线为准则。对周期序列，最好截取周期的整数倍进行谱分析，否 则有可能产生较大的分析误差编写主程序 下图给出了主程序框图，供参考。本实验提供 FFT 子程序和通用绘 图子程序。 开始 读入长度 N 调用信号产生子程序产生实验信号 调用绘图子程序(函数)绘制时间序列波形图 调用 FFT 子程序(函数)计算信号的 DFT 主程序框图 %x6(n)=cos(8пt)+cos(16пt)+cos(20пt) fs=64kHz, N=16,32,64 m=( FFT 点数= ); n=0:(m-1); subplot(2,2,1); x6=cos(pi*n/8)+cos(pi*n/4)+cos(pi*n*5/16); stem(n,x6, . ); % stem:二维火柴棍图 xlabel( n ); ylabel( x6(n) ); title( x6(n)=cos(8пt)+cos(16пt)+cos(20пt)的函数 ); subplot(2,2,2); xa=fft(x6,16); i=0:15; stem(i,abs(xa), . ); % abs:求模，绝对值 调用绘图子程序(函数)绘制 曲线   k X 结束xlabel( k ); ylabel( x6(n) ); title( x6(n)的 16点 FFT ); subplot(2,2,3); % m 行 n列个图像，第 p个区域 xb=fft(x6,32); i=0:31; stem(i,abs(xb), . ); xlabel( k ); ylabel( x6(n) ); title( x6(n)的 32点 FFT ); subplot(2,2,4); xc=fft(x6,64); i=0:63; stem(i,abs(xc), . ); xlabel( k ); ylabel( x6(n) ); title( x6(n)的 64点 FFT );直接运行程序，按照实验内容及程序提示键入 1~8，分别对 ~ 及 、 进行谱分析。输出   n x 1   n x 6       n x n x n x 5 4 7         n jx n x n x 5 4 8   的波形及其 8点 DFT 和 16点 DFT， 的 16点、32点     n x n x 5 1 ~   n x 6 和 64点采样序列及其 DFT。1、 及其 8点和 16点 DFT   n x 1 2、 及其 8点和 16点 DFT   n x 23、 及其 8点和 16点 DFT   n x 3 4、 的 8点和 16点波形及其 DFT   n x 45、 的 8点和 16点波形及其 DFT   n x 5 6、 的 16点、32点和 64点采样序列波形及其 DFT   n x 6