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

zhangliuchun

各位好，请教一下，如何把我附件中的电压信号实现FFT变换，画出其幅频曲线？
# i& {  [' [$ I
先谢过。

$ q# O$ M2 l- }0 g9 u: X我的邮箱是zhangliuchun@yahoo.com.cn

zhangliuchun

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

zhangliuchun

欢迎指点，讨论，谢谢。
: y& h, L: R' q8 J7 i2 s1 u* F4 }
  R6 o9 p$ t" [/ I$ m( ~8 j共同进步

zhangliuchun

已自行解决。

zhangliuchun

4# zhangliuchun

- f# {# i2 C" b. \7 V
自己解决该问题的代码：
  s" @& ^9 q' m- T2 m( ~
clear
load qs_切空变_vfto520090507
0 I& U/ _! J9 n5 Xtime=t;
data=vMoa4;
- D8 W( o5 p! i, y# {# L+ F5 ]%FFT
7 x  h  e/ u  S+ s- ~: \% N   = data_num;
: G: _9 M, S$ n( QN=length(time);
Ts=time(2) - time(1);
4 d0 f1 A: B5 X0 k, Q7 Lfs =1/Ts;
fp=fs*(0 : N/2)/N;
Y=fft(data);


& s/ l! E# R% }. @: c: P& a8 i& }Pyy=abs(Y(1 : N/2+1))*2/N;%求幅谱,乘上后面的2/N得到正确幅值；
& b* q. i% K# P6 a; _: i# s& V
* A7 d' }4 d: X$ w0 h/ of= fs*(0:N/2)/N;

6 e  ~4 r' a* O# H
3 g. h! E8 u8 J9 H1 s+ W
& i( T9 L: f3 l- S% G6 Cfigure;hold on;
plot(f(2:end),Pyy(2:end),'k','LineWidth',1.5)
6 m/ l4 A0 J% nxlabel(' 频  率 (Hz) ');
- `! b; W1 U  u7 p  @$ s5 aylabel('幅值 ');
% j" a- D9 g. n1 w3 Rtitle('电压信号vMoa4的幅频谱图');
+ t% @6 X* x7 T/ }9 f9 _" h
) l8 V5 _( ~& U* d请各位指正。

maisigete

"Pyy=abs(Y(1 : N/2+1))*2/N;%求幅谱,乘上后面的2/N得到正确幅值；"
. K: N. w$ l, I& }- p( L2 m* ^请问为什么要乘2/N？我在网上看到别的程序中也有这样的，不过也看到过不乘的

maisigete

回复 5# zhangliuchun


: o; P- W1 T3 S# E" U# u# M9 c    请问画图的时候为什么要从第二个点开始画？

Ericwjc

回复 7# maisigete
# S' P6 g4 _2 W+ n3 n+ V
* x+ ^( {: F' `( o: \除以N是由于DFT是N个点相加求和，乘以2是由于出现了正、负频率且正负频率的幅值均为原来的一半，这个在频谱校正里有讲到。

maisigete

回复 8# Ericwjc

. g1 x( r% h: ~, E" [2 q) `
    请问，在哪本书里什么内容有讲这个？我对这个不太了解
    还有，如果我进行傅里叶分解是为了选择频点在频域内计算，计算完了以后把频域的结果再反变换到时域，那么在傅里叶分解的时候是否还要乘以2/N？如果乘了，在反变换之前是否要对频域内各频点下的结果除以2/N再反变换？
    非常感谢

Ericwjc

回复 9# maisigete
- `8 Q7 {, H6 K$ [: p2 w
为了得到准确的幅值才乘以2/N的。实际的DFT的定义有两种，一种是在正变换里除以N，另一种是在逆变换里除以N。应该说DFT只是一种数学的定义，它可以分析频率成分（实际分析中存在频谱混叠，所以比较精确的场合需要进行校正），但得出的相应频率成分的幅值时需要修正，如乘以2/N。至于这个系数放在正变换里还是逆变换里，关系不大，只要对应上就可以了。
5 d' x6 M! T+ P你可以看下《离散频谱分析校正理论与技术》这本书，个人感觉不错。