在matlab中如何编程实现干扰信号的fft转换？

zhangliuchun

各位好，请教一下，如何把我附件中的电压信号实现FFT变换，画出其幅频曲线？

5 O3 |- c2 E! [  o% ]! O先谢过。

我的邮箱是zhangliuchun@yahoo.com.cn

zhangliuchun

欢迎讨论，欢迎指点，谢谢。

zhangliuchun

欢迎指点，讨论，谢谢。

/ q" p. d! k% h) K共同进步

zhangliuchun

已自行解决。

zhangliuchun

4# zhangliuchun


自己解决该问题的代码：
) \4 v& _- b7 \! |0 |
" i% ], k1 H; X4 J/ eclear
load qs_切空变_vfto520090507
6 _2 y. f' j5 Stime=t;
data=vMoa4;
3 q  e* k( a; _  D+ p1 ]3 ^%FFT
% g' Z3 r$ w, n6 o5 d4 c% N   = data_num;
N=length(time);
5 r; o7 d# ^: s) F: DTs=time(2) - time(1);
fs =1/Ts;
! E1 o; Q* X. U6 }! ]; l4 R' Ffp=fs*(0 : N/2)/N;
8 @0 j& ]; U( v1 W6 i1 T) mY=fft(data);


Pyy=abs(Y(1 : N/2+1))*2/N;%求幅谱,乘上后面的2/N得到正确幅值；
/ U6 |  I# E2 \3 v! h7 }
f= fs*(0:N/2)/N;
6 Q. F! H: w+ T8 }& o
2 c8 K( r$ ~3 y& e% `
2 Y& ^* r" Z" ]7 g
3 N7 T) e. E! bfigure;hold on;
plot(f(2:end),Pyy(2:end),'k','LineWidth',1.5)
4 h# n+ {7 }2 l4 L: J- Q5 Ixlabel(' 频  率 (Hz) ');
" Y9 L! l4 ~# {9 @ylabel('幅值 ');
title('电压信号vMoa4的幅频谱图');
" Y: A3 ?# {. ~' w9 X
请各位指正。

maisigete

"Pyy=abs(Y(1 : N/2+1))*2/N;%求幅谱,乘上后面的2/N得到正确幅值；"
5 J4 P$ h; h0 E7 G6 R7 J) R请问为什么要乘2/N？我在网上看到别的程序中也有这样的，不过也看到过不乘的

maisigete

回复 5# zhangliuchun

) a* d* ]5 b7 X" }. X+ y7 j$ g
! x% ~1 M9 X! {) H8 K, a- f    请问画图的时候为什么要从第二个点开始画？

Ericwjc

回复 7# maisigete
) c7 w* W: E# I6 ?) f1 j& `
除以N是由于DFT是N个点相加求和，乘以2是由于出现了正、负频率且正负频率的幅值均为原来的一半，这个在频谱校正里有讲到。

maisigete

回复 8# Ericwjc


6 V: T3 S4 n* }2 Z    请问，在哪本书里什么内容有讲这个？我对这个不太了解
    还有，如果我进行傅里叶分解是为了选择频点在频域内计算，计算完了以后把频域的结果再反变换到时域，那么在傅里叶分解的时候是否还要乘以2/N？如果乘了，在反变换之前是否要对频域内各频点下的结果除以2/N再反变换？
    非常感谢

Ericwjc

回复 9# maisigete
8 m* h* D( p/ x: ?% x3 ]+ ?/ n
为了得到准确的幅值才乘以2/N的。实际的DFT的定义有两种，一种是在正变换里除以N，另一种是在逆变换里除以N。应该说DFT只是一种数学的定义，它可以分析频率成分（实际分析中存在频谱混叠，所以比较精确的场合需要进行校正），但得出的相应频率成分的幅值时需要修正，如乘以2/N。至于这个系数放在正变换里还是逆变换里，关系不大，只要对应上就可以了。
" i5 E/ E8 j2 R4 o' Y你可以看下《离散频谱分析校正理论与技术》这本书，个人感觉不错。