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

zhangliuchun

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

) }, c- E, u( I先谢过。

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

zhangliuchun

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

zhangliuchun

欢迎指点，讨论，谢谢。
$ {+ x2 L3 Z0 `$ R; y( M
; o1 F; }& y0 q- O" J& |8 O8 W& ~共同进步

zhangliuchun

已自行解决。

zhangliuchun

4# zhangliuchun
$ L# H# {$ `% C' `: l. e
3 y/ A7 _! P* p3 E# M& _1 f; \
自己解决该问题的代码：
* w. f7 Y" U; u) h8 s  m- o3 f7 _
clear
; E. m/ E! K9 U$ m2 @load qs_切空变_vfto520090507
9 I3 j# R7 E, `. K& btime=t;
& X- K9 V1 H! @7 y) G3 Z* Adata=vMoa4;
4 O3 z# N- I  p. X7 _9 `9 g. J; C3 {%FFT
% N   = data_num;
N=length(time);
( m* }9 Z) `0 z/ ^/ ~# ?Ts=time(2) - time(1);
fs =1/Ts;
" @! G( u; A' U/ _fp=fs*(0 : N/2)/N;
2 x4 O- p0 A2 Z( a$ SY=fft(data);
( j$ e! m9 L! j8 y! {9 j) x
% Y) R" i& C2 I4 R+ o0 \
Pyy=abs(Y(1 : N/2+1))*2/N;%求幅谱,乘上后面的2/N得到正确幅值；

5 _, n* u2 J- cf= fs*(0:N/2)/N;
3 _) k, V- d( t% X0 d8 G* R7 I, _0 }

  A6 P: {+ C7 T1 z/ `  W
1 V% m& B% S  Ifigure;hold on;
plot(f(2:end),Pyy(2:end),'k','LineWidth',1.5)
+ c: F) z$ g: q  uxlabel(' 频  率 (Hz) ');
ylabel('幅值 ');
title('电压信号vMoa4的幅频谱图');
, B- U- h5 G( ?0 z% B8 K. L
请各位指正。

maisigete

"Pyy=abs(Y(1 : N/2+1))*2/N;%求幅谱,乘上后面的2/N得到正确幅值；"
, x/ S" B8 q& X9 @$ Y# E. o6 B请问为什么要乘2/N？我在网上看到别的程序中也有这样的，不过也看到过不乘的

maisigete

回复 5# zhangliuchun
/ c/ A9 z+ {0 e; ^7 [5 K$ y

1 E/ Q: b) |0 }" X3 a7 `    请问画图的时候为什么要从第二个点开始画？

Ericwjc

回复 7# maisigete
3 [/ P) g7 Z9 L0 k/ w
除以N是由于DFT是N个点相加求和，乘以2是由于出现了正、负频率且正负频率的幅值均为原来的一半，这个在频谱校正里有讲到。

maisigete

回复 8# Ericwjc


    请问，在哪本书里什么内容有讲这个？我对这个不太了解
0 k+ @9 T% X* e3 h( B    还有，如果我进行傅里叶分解是为了选择频点在频域内计算，计算完了以后把频域的结果再反变换到时域，那么在傅里叶分解的时候是否还要乘以2/N？如果乘了，在反变换之前是否要对频域内各频点下的结果除以2/N再反变换？
    非常感谢

Ericwjc

回复 9# maisigete

为了得到准确的幅值才乘以2/N的。实际的DFT的定义有两种，一种是在正变换里除以N，另一种是在逆变换里除以N。应该说DFT只是一种数学的定义，它可以分析频率成分（实际分析中存在频谱混叠，所以比较精确的场合需要进行校正），但得出的相应频率成分的幅值时需要修正，如乘以2/N。至于这个系数放在正变换里还是逆变换里，关系不大，只要对应上就可以了。
, q6 l4 d, ]# b) S* {1 v你可以看下《离散频谱分析校正理论与技术》这本书，个人感觉不错。