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

zhangliuchun

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

. Y. l  H  w3 n; {先谢过。
0 H7 J1 ~$ r9 I1 Y8 C
9 E* ?' _9 P. i6 J$ C我的邮箱是zhangliuchun@yahoo.com.cn

zhangliuchun

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

zhangliuchun

欢迎指点，讨论，谢谢。
# B) l6 ?& F* ~5 A  V
) M5 V6 B( J  A) k; M& l7 v4 |共同进步

zhangliuchun

已自行解决。

zhangliuchun

4# zhangliuchun

/ p" Z) N% P2 m
1 q: b9 h! }1 U2 o; Z+ ~2 O8 C自己解决该问题的代码：
* H1 k' f# ]$ @$ j9 ?0 _
clear
& T. \- ?' W8 L; Z1 Kload qs_切空变_vfto520090507
time=t;
3 M5 n) S# ?% d( gdata=vMoa4;
9 A. d4 T* _" v, O6 V%FFT
" z5 ~2 h1 |; W' B2 f% N   = data_num;
N=length(time);
Ts=time(2) - time(1);
fs =1/Ts;
# p+ I3 y0 ]) K& Pfp=fs*(0 : N/2)/N;
Y=fft(data);


Pyy=abs(Y(1 : N/2+1))*2/N;%求幅谱,乘上后面的2/N得到正确幅值；

f= fs*(0:N/2)/N;

% x2 S1 W  c) o  {0 T9 {8 o

( N1 [! r) d7 R# D+ r1 S5 Nfigure;hold on;
plot(f(2:end),Pyy(2:end),'k','LineWidth',1.5)
7 d4 T! b# [- u, k6 G7 xxlabel(' 频  率 (Hz) ');
ylabel('幅值 ');
' C% P4 u. N, Utitle('电压信号vMoa4的幅频谱图');

请各位指正。

maisigete

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

maisigete

回复 5# zhangliuchun
  {& [3 J4 f2 ~5 Z) G) d, L
! P. n1 F' K. A" o
    请问画图的时候为什么要从第二个点开始画？

Ericwjc

回复 7# maisigete
) w& G1 V& f) }- {
) K8 V6 k+ f8 @! m  T$ c除以N是由于DFT是N个点相加求和，乘以2是由于出现了正、负频率且正负频率的幅值均为原来的一半，这个在频谱校正里有讲到。

maisigete

回复 8# Ericwjc
! f. l! _7 E+ [$ `( ^/ m

; K. S- [8 D& O2 D4 [1 H    请问，在哪本书里什么内容有讲这个？我对这个不太了解
8 ^2 x8 h4 j; F% V% w, u- S    还有，如果我进行傅里叶分解是为了选择频点在频域内计算，计算完了以后把频域的结果再反变换到时域，那么在傅里叶分解的时候是否还要乘以2/N？如果乘了，在反变换之前是否要对频域内各频点下的结果除以2/N再反变换？
/ @9 }! J- E& G2 J5 i/ v9 F# c3 U    非常感谢

Ericwjc

回复 9# maisigete
. K- Z9 N) C7 w, K
# W7 }+ B' s9 w7 @为了得到准确的幅值才乘以2/N的。实际的DFT的定义有两种，一种是在正变换里除以N，另一种是在逆变换里除以N。应该说DFT只是一种数学的定义，它可以分析频率成分（实际分析中存在频谱混叠，所以比较精确的场合需要进行校正），但得出的相应频率成分的幅值时需要修正，如乘以2/N。至于这个系数放在正变换里还是逆变换里，关系不大，只要对应上就可以了。
1 G9 _8 i+ i# R) Y# {6 i5 L  _你可以看下《离散频谱分析校正理论与技术》这本书，个人感觉不错。