设为首页收藏本站|繁體中文 快速切换版块

 找回密码
 立即加入
搜索
查看: 2991|回复: 11

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

[复制链接]

该用户从未签到

尚未签到

发表于 2009-5-6 13:54:19 | 显示全部楼层 |阅读模式

马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!

您需要 登录 才可以下载或查看,没有账号?立即加入

×
各位好,请教一下,如何把我附件中的电压信号实现FFT变换,画出其幅频曲线?
  f5 {( g5 t9 [" T3 G" p. Z! K# Z+ I
先谢过。
2 z" o' M5 U5 G# O* ]
' N  U5 b6 `+ M我的邮箱是zhangliuchun@yahoo.com.cn

qs_切空变_vfto5.rar

448.31 KB, 下载次数: 8, 下载积分: 威望 -2 点, 学分 -5 点

"真诚赞赏,手留余香"
还没有人打赏,支持一下
楼主热帖
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

该用户从未签到

尚未签到

 楼主| 发表于 2009-5-6 14:39:13 | 显示全部楼层
欢迎讨论,欢迎指点,谢谢。
"真诚赞赏,手留余香"
还没有人打赏,支持一下
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

该用户从未签到

尚未签到

 楼主| 发表于 2009-5-6 14:40:02 | 显示全部楼层
欢迎指点,讨论,谢谢。
9 v% M3 w- s( v# j; \* k' Y6 k  e
共同进步
"真诚赞赏,手留余香"
还没有人打赏,支持一下
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

该用户从未签到

尚未签到

 楼主| 发表于 2009-5-11 16:31:23 | 显示全部楼层
已自行解决。
"真诚赞赏,手留余香"
还没有人打赏,支持一下
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

该用户从未签到

尚未签到

 楼主| 发表于 2009-5-11 17:13:18 | 显示全部楼层
4# zhangliuchun 5 N' j2 |3 w$ Q1 p5 Z, d) `3 z
; J0 l& G, q: x- P+ y/ S, h- }* }
" F& r$ Y7 m# U* s6 U# g$ @
自己解决该问题的代码:
8 `% l" ]( ^3 d3 e. [2 I; N1 K/ R
clear; {) N4 I- r3 i. f$ |& Q
load qs_切空变_vfto520090507; O+ `& R' u" F  ^/ n8 V4 |
time=t;  \9 A2 @! z! }2 Z$ w
data=vMoa4;; o* Q7 Q" e2 t. s8 ~* d! `- N( `
%FFT/ h& `1 u5 D) h4 N6 k
% N   = data_num;) R, L& ~6 u8 {& w2 ^
N=length(time);
2 t( B( z9 ]7 M- u+ W! r- w. aTs=time(2) - time(1);
9 E$ K& g# L6 Z! [3 Ifs =1/Ts;$ d% ?9 I: A# I0 m4 d
fp=fs*(0 : N/2)/N;2 M- G/ i* B+ E* z
Y=fft(data);! ^6 l- O" ^# z
" s( Z% e) D2 ~
0 E3 X) N; o3 W. e
Pyy=abs(Y(1 : N/2+1))*2/N;%求幅谱,乘上后面的2/N得到正确幅值;
( k/ d+ x, f7 Y+ r  X
) Q& S) o" J8 p! j1 V3 Kf= fs*(0:N/2)/N;
$ |2 z; U) \6 O* ~0 ^) S# v. s- o, P

, D- g% h1 V2 J5 J8 U
0 l0 K0 N1 \/ H; M+ H- x. h. t1 Efigure;hold on;. S. B2 ?6 e7 j% ^  r* Y4 M3 G: j, E
plot(f(2:end),Pyy(2:end),'k','LineWidth',1.5)
2 a. s8 P6 b9 k! b& R  yxlabel(' 频  率 (Hz) ');
9 |+ m$ u' g+ w2 r7 k' jylabel('幅值 ');, O. Z$ x3 q6 n8 S; r- ~' i
title('电压信号vMoa4的幅频谱图'); 3 Q+ X- @) t" H  H# k9 e: {, ]

' W# i) G+ `) {0 @. j请各位指正。
"真诚赞赏,手留余香"
还没有人打赏,支持一下
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

该用户从未签到

尚未签到

发表于 2012-5-16 16:00:05 | 显示全部楼层
"Pyy=abs(Y(1 : N/2+1))*2/N;%求幅谱,乘上后面的2/N得到正确幅值;"
! @; _! t+ l+ h" F) W7 q" A. u请问为什么要乘2/N?我在网上看到别的程序中也有这样的,不过也看到过不乘的
"真诚赞赏,手留余香"
还没有人打赏,支持一下
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

该用户从未签到

尚未签到

发表于 2012-5-16 16:30:02 | 显示全部楼层
回复 5# zhangliuchun
7 v4 P' ~" }. B0 Z2 C  D( @" s
$ g2 f& d5 e8 L- `1 x) z( R' F8 ]3 _1 f; z, P
    请问画图的时候为什么要从第二个点开始画?
"真诚赞赏,手留余香"
还没有人打赏,支持一下
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    无聊
    2020-12-23 16:05
  • 签到天数: 34 天

    连续签到: 1 天

    [LV.5]常住居民I

    累计签到:34 天
    连续签到:1 天
    发表于 2012-5-16 18:20:00 | 显示全部楼层
    回复 7# maisigete
    6 @# N+ r  g4 R# m$ V" u$ x5 h1 m, f. e8 e
    除以N是由于DFT是N个点相加求和,乘以2是由于出现了正、负频率且正负频率的幅值均为原来的一半,这个在频谱校正里有讲到。

    评分

    参与人数 1威望 +2 学分 +2 收起 理由
    maisigete + 2 + 2 多谢帮助 很专业啊

    查看全部评分

    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

    该用户从未签到

    尚未签到

    发表于 2012-5-17 09:56:37 | 显示全部楼层
    回复 8# Ericwjc " y' X' g0 f/ S% S7 T( G
    5 F: a3 z- A+ c! ^
    2 i6 I: @* k/ b# ]
        请问,在哪本书里什么内容有讲这个?我对这个不太了解
    : n* o9 {* ]3 b0 Z1 i# n    还有,如果我进行傅里叶分解是为了选择频点在频域内计算,计算完了以后把频域的结果再反变换到时域,那么在傅里叶分解的时候是否还要乘以2/N?如果乘了,在反变换之前是否要对频域内各频点下的结果除以2/N再反变换?
    ( h6 r! I* W( L& L    非常感谢
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    无聊
    2020-12-23 16:05
  • 签到天数: 34 天

    连续签到: 1 天

    [LV.5]常住居民I

    累计签到:34 天
    连续签到:1 天
    发表于 2012-5-17 10:57:52 | 显示全部楼层
    回复 9# maisigete ! W- Z( D: K& `
    3 T, b* t, E. y' x3 Q9 `" [
    为了得到准确的幅值才乘以2/N的。实际的DFT的定义有两种,一种是在正变换里除以N,另一种是在逆变换里除以N。应该说DFT只是一种数学的定义,它可以分析频率成分(实际分析中存在频谱混叠,所以比较精确的场合需要进行校正),但得出的相应频率成分的幅值时需要修正,如乘以2/N。至于这个系数放在正变换里还是逆变换里,关系不大,只要对应上就可以了。# [8 e* U  @, F7 a& V
    你可以看下《离散频谱分析校正理论与技术》这本书,个人感觉不错。
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
    您需要登录后才可以回帖 登录 | 立即加入

    本版积分规则

    招聘斑竹

    小黑屋|手机版|APP下载(beta)|Archiver|电力研学网 ( 赣ICP备12000811号-1|赣公网安备36040302000210号 )|网站地图

    GMT+8, 2025-4-5 21:28

    Powered by Discuz! X3.5 Licensed

    © 2001-2025 Discuz! Team.

    快速回复 返回顶部 返回列表