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

zhangliuchun

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

先谢过。
; a3 p* U* x2 ^5 T2 ?" B2 P
8 n1 p* E! ?+ [7 b8 G9 ]& G我的邮箱是zhangliuchun@yahoo.com.cn

zhangliuchun

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

zhangliuchun

欢迎指点，讨论，谢谢。

2 l7 ^* p* g  g: E( H  e8 H共同进步

zhangliuchun

已自行解决。

zhangliuchun

4# zhangliuchun

& {/ x+ m: R/ v, C7 p& x; F8 c
自己解决该问题的代码：

3 W$ A5 d4 Y' O7 s/ g: U2 i( kclear
- k+ ^9 I5 T% y" O8 _load qs_切空变_vfto520090507
) X( N! R& F9 [" w5 Btime=t;
6 z. _- d) Q4 Q0 v4 G2 ?1 kdata=vMoa4;
1 B* H( ]2 Z5 l9 V%FFT
% N   = data_num;
N=length(time);
Ts=time(2) - time(1);
) n+ P; x1 r+ K% R* {$ Cfs =1/Ts;
7 J$ j2 t4 Y! w9 s; C' o3 Bfp=fs*(0 : N/2)/N;
Y=fft(data);
. T; C4 O5 c- I2 _& W$ g+ W
" p7 L- I3 K' D  U! E8 ?8 U
* e6 P2 R$ v2 ~  V6 T) @' U8 jPyy=abs(Y(1 : N/2+1))*2/N;%求幅谱,乘上后面的2/N得到正确幅值；
+ X! C! C. D& Y$ K% i
f= fs*(0:N/2)/N;



figure;hold on;
plot(f(2:end),Pyy(2:end),'k','LineWidth',1.5)
( Z% n; j- \9 z: C$ m  Pxlabel(' 频  率 (Hz) ');
5 x, O! [7 }6 dylabel('幅值 ');
title('电压信号vMoa4的幅频谱图');
9 T* j3 z6 w: C# N* E
请各位指正。

maisigete

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

maisigete

回复 5# zhangliuchun


    请问画图的时候为什么要从第二个点开始画？

Ericwjc

回复 7# maisigete
# P9 F" R, ^. M" C4 [
" Q4 R- u" x0 |2 G  g4 K& w7 ~: q除以N是由于DFT是N个点相加求和，乘以2是由于出现了正、负频率且正负频率的幅值均为原来的一半，这个在频谱校正里有讲到。

maisigete

回复 8# Ericwjc

. J+ d' l5 {3 U7 Z& j
; w1 a  ]3 ]$ y! z1 a/ J5 L    请问，在哪本书里什么内容有讲这个？我对这个不太了解
    还有，如果我进行傅里叶分解是为了选择频点在频域内计算，计算完了以后把频域的结果再反变换到时域，那么在傅里叶分解的时候是否还要乘以2/N？如果乘了，在反变换之前是否要对频域内各频点下的结果除以2/N再反变换？
    非常感谢

Ericwjc

回复 9# maisigete

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