连续信号的谱分析(FFT)

linux_li · 发表于 2005-12-25 18:47:01

对于频谱很宽的信号，为了防止采样后产生的频谱混叠失真，一般采用滤波器去除幅度较小的高频部分。因为这里主要进行带宽有限的谱分析，为了分析的方便，采用了标准正弦信号，其他信号可以做进一步的推广。
相关参数（下面的几个概念一定要清楚，否则采集系统是不可控制下工作的，所做的分析没有任何意义）
要分析信号的最高频率采样频率对信号的采样点数采样持续时间 DFT变换区间长度N
等间隔采样时间T

[此贴子已经被作者于2005-12-25 10:56:24编辑过]

linux_li · 发表于 2005-12-25 18:48:19

<

>谱分析范围的分析与仿真
<

0cm 0cm 0pt"> <o:p></o:p>
<

0cm 0cm 0pt 18pt; TEXT-INDENT: -18pt; mso-list: l0 level1 lfo1; tab-stops: list 18.0pt">1．谱分析范围的分析与仿真<o:p></o:p>
谱分析的范围受采样频率<v:shapetype> <v:stroke joinstyle="miter"></v:stroke><v:formulas><v:f eqn="if lineDrawn pixelLineWidth 0"></v:f><v:f eqn="sum @0 1 0"></v:f><v:f eqn="sum 0 0 @1"></v:f><v:f eqn="prod @2 1 2"></v:f><v:f eqn="prod @3 21600 pixelWidth"></v:f><v:f eqn="prod @3 21600 pixelHeight"></v:f><v:f eqn="sum @0 0 1"></v:f><v:f eqn="prod @6 1 2"></v:f><v:f eqn="prod @7 21600 pixelWidth"></v:f><v:f eqn="sum @8 21600 0"></v:f><v:f eqn="prod @7 21600 pixelHeight"></v:f><v:f eqn="sum @10 21600 0"></v:f></v:formulas><v:path o:connecttype="rect" gradientshapeok="t" o:extrusionok="f"></v:path><o:lock aspectratio="t" v:ext="edit"></o:lock></v:shapetype><v:shape><v:imagedata></v:imagedata></v:shape>的影响，其影响如下图所示<o:p></o:p>
在<v:shape> <v:imagedata></v:imagedata></v:shape>是<v:shape> <v:imagedata></v:imagedata></v:shape>20倍的情况下，信号频谱的情况，基本表现良好,采样获取信号的时域基本没有什么失真<o:p></o:p>
<v:shape><v:imagedata></v:imagedata></v:shape><o:p></o:p>
在<v:shape> <v:imagedata></v:imagedata></v:shape>是<v:shape> <v:imagedata></v:imagedata></v:shape>5倍的情况下，信号频谱的情况，基本表现良好,采样获取信号的时域显示有一定程度的失真<o:p></o:p>
<v:shape><v:imagedata></v:imagedata></v:shape> 在<v:shape> <v:imagedata></v:imagedata></v:shape>是<v:shape> <v:imagedata></v:imagedata></v:shape>2倍的情况下，信号频谱混叠,采样获取信号的时域显示完全失真<o:p></o:p>
<v:shape><v:imagedata></v:imagedata></v:shape><o:p></o:p>
结论：如上图所列出的，谱分析的范围受采样频率的限制，因此在进行谱分析中，必须满足采样定理，为了获取更好的效果，一般应该采取采样频率是信号最高频率的五倍以上。

linux_li · 发表于 2005-12-25 18:50:08

<

0cm 0cm 0pt 18pt; TEXT-INDENT: -18pt; mso-list: l0 level1 lfo1; tab-stops: list 18.0pt">1．频率分辨率的分析与仿真<o:p></o:p>
<

0cm 0cm 0pt 42pt; TEXT-INDENT: -21pt; mso-list: l0 level2 lfo1; tab-stops: list 42.0pt">a) 定义：谱分析中能够分辨的两个频谱分量的最小间隔，用频率采样间隔F来表示。换句话说，就是能够分析的信号中的最低频率，如果信号中含有最低频率为<v:shapetype> <v:stroke joinstyle="miter"></v:stroke><v:formulas><v:f eqn="if lineDrawn pixelLineWidth 0"></v:f><v:f eqn="sum @0 1 0"></v:f><v:f eqn="sum 0 0 @1"></v:f><v:f eqn="prod @2 1 2"></v:f><v:f eqn="prod @3 21600 pixelWidth"></v:f><v:f eqn="prod @3 21600 pixelHeight"></v:f><v:f eqn="sum @0 0 1"></v:f><v:f eqn="prod @6 1 2"></v:f><v:f eqn="prod @7 21600 pixelWidth"></v:f><v:f eqn="sum @8 21600 0"></v:f><v:f eqn="prod @7 21600 pixelHeight"></v:f><v:f eqn="sum @10 21600 0"></v:f></v:formulas><v:path o:connecttype="rect" gradientshapeok="t" o:extrusionok="f"></v:path><o:lock aspectratio="t" v:ext="edit"></o:lock></v:shapetype><v:shape><v:imagedata></v:imagedata></v:shape>，采样后要保持该频率成分，则采样长度<v:shape> <v:imagedata></v:imagedata></v:shape><v:shape><v:imagedata></v:imagedata></v:shape><v:shape><v:imagedata></v:imagedata></v:shape>。<o:p></o:p>
<

0cm 0cm 0pt 21pt"><v:shape><v:imagedata></v:imagedata></v:shape> 因为<v:shape> <v:imagedata></v:imagedata></v:shape>=<v:shape> <v:imagedata></v:imagedata></v:shape>,所以<v:shape> <v:imagedata></v:imagedata></v:shape><o:p></o:p>
b) 从上面公式可以看出提高频率分辨率的一种方法，就是在不改变采样频率的情况下，增加采样时间<v:shape> <v:imagedata></v:imagedata></v:shape>,也就是增加<v:shape> <v:imagedata></v:imagedata></v:shape>。<o:p></o:p>
下面给出<v:shape> <v:imagedata></v:imagedata></v:shape>不同情况下仿真结果:<o:p></o:p>
<v:shape><v:imagedata></v:imagedata></v:shape>=100,也就是采样时间为0.1s，F=10Hz,主瓣很宽，分辨率很低<o:p></o:p>
<v:shape><v:imagedata></v:imagedata></v:shape><v:shape><v:imagedata></v:imagedata></v:shape>=500,也就是采样时间为0.5s，F=2Hz,主瓣明显变窄，分辨率得到很大的提高<o:p></o:p>
<v:shape><v:imagedata></v:imagedata></v:shape><v:shape><v:imagedata></v:imagedata></v:shape>=1000,也就是采样时间为1s，F=1Hz,主瓣变窄，分辨率提高<o:p></o:p>
<v:shape><v:imagedata></v:imagedata></v:shape><o:p></o:p>
c) 上面的方法实现了分辨率的提高，但是考虑到如果采样长度<v:shape> <v:imagedata></v:imagedata></v:shape>取得较长，但在T不变的情况下，采样点数<v:shape> <v:imagedata></v:imagedata></v:shape>增多，使计算机的工作量增大。这里给出提高频率分辨率的另外一种方法，就是在保持<v:shape> <v:imagedata></v:imagedata></v:shape>不变的情况下，降低采样频率，这样虽然提高了频率分辨率，但是引起了谱分析范围的减少。由于这个是直观的描述就不再做仿真了。<o:p></o:p>综上所述，采集信号的时间间隔和长度，对得出的结论正确与否有直接影响，在一般情况下，分析中为了处理的方便，一般采用固定<v:shape> <v:imagedata></v:imagedata></v:shape>的方式，因此它的频率分析范围取决于采样间隔T（保证满足采样定理）,最低的分析频率取决于采样的长度<v:shape> <v:imagedata></v:imagedata></v:shape>,对即要求高的分析范围又要求频率分辨率的系统，一般采取高的采样率和合适的采样时间来完成。

linux_li · 发表于 2005-12-25 18:51:21

<

0cm 0cm 0pt 18pt; TEXT-INDENT: -18pt; mso-list: l0 level1 lfo1; tab-stops: list 18.0pt">1．关于DFT变换区间长度N的仿真<o:p></o:p>
<

0cm 0cm 0pt; TEXT-INDENT: 18pt"> <o:p></o:p>
<

0cm 0cm 0pt; TEXT-INDENT: 18pt">DFT变换区间长度N=1024下的频谱情况<o:p></o:p>
<v:shapetype><v:stroke joinstyle="miter"></v:stroke><v:formulas><v:f eqn="if lineDrawn pixelLineWidth 0"></v:f><v:f eqn="sum @0 1 0"></v:f><v:f eqn="sum 0 0 @1"></v:f><v:f eqn="prod @2 1 2"></v:f><v:f eqn="prod @3 21600 pixelWidth"></v:f><v:f eqn="prod @3 21600 pixelHeight"></v:f><v:f eqn="sum @0 0 1"></v:f><v:f eqn="prod @6 1 2"></v:f><v:f eqn="prod @7 21600 pixelWidth"></v:f><v:f eqn="sum @8 21600 0"></v:f><v:f eqn="prod @7 21600 pixelHeight"></v:f><v:f eqn="sum @10 21600 0"></v:f></v:formulas><v:path o:connecttype="rect" gradientshapeok="t" o:extrusionok="f"></v:path><o:lock aspectratio="t" v:ext="edit"></o:lock></v:shapetype><v:shape><v:imagedata></v:imagedata></v:shape><o:p></o:p>
DFT变换区间长度N=65536下的频谱情况<o:p></o:p>
<v:shape><v:imagedata></v:imagedata></v:shape> 
 <o:p></o:p>

linux_li · 发表于 2005-12-25 18:51:56

<

0cm 0cm 0pt 18pt; TEXT-INDENT: -18pt; mso-list: l0 level1 lfo1; tab-stops: list 18.0pt">1．谱分析的误差问题<o:p></o:p>
<

0cm 0cm 0pt">①混叠现象:<o:p></o:p>
<

0cm 0cm 0pt">在对联系信号进行处理的时间，只有先对其采样后才可以利用DFT进行谱分析。采样频率f必须满足采样定理，这个问题在上面 谱分析范围的分析与仿真 中已经做了完整的说明。<o:p></o:p>
②栅栏效应：<o:p></o:p>
由于栅栏效应，有可能漏掉大的频谱分量，为了解决这个问题，可以采用改变序列长度N，从图上可以清楚看出，在N=1024漏掉的某些频谱分量在N=65536的情况下被检测出来了。<o:p></o:p>
③截断效应：<o:p></o:p> 由FT理论可以知道，若信号的频谱有限宽，那么持续的时间无限长，由于DFT不可能无限长，所以截断后的序列与原来序列必然有差别，这个差别表现在谱分析上主要有两点：

linux_li · 发表于 2005-12-25 18:52:30

<

0cm 0cm 0pt">ⅰ泄露：<o:p></o:p>
<

0cm 0cm 0pt"> 根据FT理论，正弦信号在频谱分析上表现为离散谱线，从上图可以明显看出，离散谱线向附近展宽，显然这种展宽使谱的分辨率降低<o:p></o:p>
<

0cm 0cm 0pt">ⅱ 谱间干扰：<o:p></o:p>
 在主谱线的两边形成很多的旁瓣，在图中可以看出，旁瓣信号也是比较强的，由于频率接近，把旁瓣误以为另一信号的谱线，从而造成假信号，会对谱分析产生较大的偏差。<o:p></o:p>
 <o:p></o:p>

linux_li · 发表于 2005-12-25 18:58:10

没有办法，无法实现图片传输，只有这样了，需要的可以和我联系，我这里有DOC格式的，呵呵，希望大家都交流，共同进步

jm2009 · 发表于 2007-3-6 16:21:59

给我发一份吧！

Email：jm2009@21cn.com

谢谢！

zxh0294 · 发表于 2007-3-17 14:13:02

也发一份给我吧！
Email：zxh0294@yahoo.com.cn
xiexie

wh_hejia · 发表于 2007-3-22 18:58:21

看了你的讲解大开眼界！
也发一份给我吧！
Email：wh_hejia@yahoo.com.cn
谢谢大侠了！

mahairui · 发表于 2007-3-23 20:11:30

也发一份给我吧！
Email：mahairui@sohu.com

desertwolf · 发表于 2007-4-2 08:43:22

谢谢了！！！！！！！！！！！！！！！！

ykn0310331 · 发表于 2007-4-13 18:00:49

楼主，想向你请教一下信号频域分析中的时移特性和频移特性怎么编VI,拜托了～～～～～～～～～～～～～

pjpj7777 · 发表于 2007-5-31 19:36:39

请也给我一份吧
pjpj7777@163.com
谢谢!!!!

kuansso · 发表于 2007-12-21 10:03:19

请给我发一份，谢谢!!email]kuansso@yahoo.com.cn[/email]

账号		自动登录	找回密码
密码			注册

LabVIEW+单片机学习套件全套教程资料下载[免费]	LabVIEW论坛精华列表贴	USB0816数据采集卡	《LabVIEW宝典》
LabWindows/CVI论坛精华贴	NET0816以太网数据采集卡	RC0210远程设备授权系统关闭关停	锁定打开设备户外分布式数据采集
NET1624低速高精度以太网数据采集卡	WIFI0824SD无线WIFI网络数据采集卡	脱机运行 SD存储小尺寸微型串口采集	远程采集安卓手持移动采集纪录仪