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

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

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

[复制链接]

该用户从未签到

尚未签到

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

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

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

×
各位好,请教一下,如何把我附件中的电压信号实现FFT变换,画出其幅频曲线?) t9 E7 U$ A' |

! A$ n% y3 J6 Q# @- r先谢过。: y8 m  ^0 W& v

7 ?" d1 Z- O4 n' y$ E: T我的邮箱是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 o( _! ?- ~: j* t- y. C- J
% s. m7 T% u3 U+ f
共同进步
"真诚赞赏,手留余香"
还没有人打赏,支持一下
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

该用户从未签到

尚未签到

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

该用户从未签到

尚未签到

 楼主| 发表于 2009-5-11 17:13:18 | 显示全部楼层
4# zhangliuchun ) v5 \/ J$ |/ z# o9 l, I/ D$ P
3 R  `! `0 S' y, G

0 _4 U% }  Z# j7 n- T' [自己解决该问题的代码:# T7 K4 o. r3 [
; z4 E+ _* [$ G5 i& g
clear4 c, Z2 s3 l- F3 g- n$ |1 W: e: k
load qs_切空变_vfto520090507
  s- E! e7 L% J: Ntime=t;* F8 Y$ B* L/ R3 p
data=vMoa4;7 i8 O$ A( ?! w. a' G6 g
%FFT  \2 S& ?5 x- D, l  d8 ^9 g
% N   = data_num;
  ]8 \  X! _7 uN=length(time);
3 N. r+ v7 P' K* s& W) gTs=time(2) - time(1);! Z1 M: n4 P" D- b$ [) y! U- B
fs =1/Ts;; N2 O5 U- q- u& n- ^' X1 |
fp=fs*(0 : N/2)/N;
1 s, h1 C3 t6 |+ t' [Y=fft(data);6 m1 J# m, ]' ~( R6 f

6 z7 ~7 m% \# i3 Y, m3 }. d
' n5 F. C* r+ fPyy=abs(Y(1 : N/2+1))*2/N;%求幅谱,乘上后面的2/N得到正确幅值;# t2 Z5 H2 r+ x; A

$ _4 d  S( }4 uf= fs*(0:N/2)/N;# R1 N& p" g2 ^& D: _) E- W0 n
4 e4 d' q$ g6 K/ R
# B" o  n' @, `- g  r
6 w' R( C0 g6 o
figure;hold on;
/ w# ]6 n7 M* R" J7 H1 `4 C; d+ w' Oplot(f(2:end),Pyy(2:end),'k','LineWidth',1.5)
; V3 a6 L& d/ G: N( d$ O$ z+ L9 Ixlabel(' 频  率 (Hz) ');
8 |" s( b+ T( _1 E# lylabel('幅值 ');3 b. K* U! {# o" z
title('电压信号vMoa4的幅频谱图');
4 s) l+ T5 d4 `& h4 C$ t
" U5 U- t: Z7 p1 i3 y; s; G请各位指正。
"真诚赞赏,手留余香"
还没有人打赏,支持一下
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】

该用户从未签到

尚未签到

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

该用户从未签到

尚未签到

发表于 2012-5-16 16:30:02 | 显示全部楼层
回复 5# zhangliuchun # X6 d' E2 N3 f9 }$ g3 P0 S
6 ]9 J; u6 L" q# z! s: G

& G9 Y" ^( x9 J) T- H5 t, L6 F9 c    请问画图的时候为什么要从第二个点开始画?
"真诚赞赏,手留余香"
还没有人打赏,支持一下
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    无聊
    2020-12-23 16:05
  • 签到天数: 34 天

    连续签到: 1 天

    [LV.5]常住居民I

    累计签到:34 天
    连续签到:1 天
    发表于 2012-5-16 18:20:00 | 显示全部楼层
    回复 7# maisigete 5 k) t# [6 o) H/ ^8 v, i2 R+ N
    . @! i6 m. P8 b% f; @6 f( r0 y
    除以N是由于DFT是N个点相加求和,乘以2是由于出现了正、负频率且正负频率的幅值均为原来的一半,这个在频谱校正里有讲到。

    评分

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

    查看全部评分

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

    该用户从未签到

    尚未签到

    发表于 2012-5-17 09:56:37 | 显示全部楼层
    回复 8# Ericwjc 9 ^" |3 x7 M+ m+ V
    ' r( f, {1 a' N4 q- ^/ ]5 h+ J

    ( A: O" e! h6 r) F' c8 n    请问,在哪本书里什么内容有讲这个?我对这个不太了解
    ! n( S( j0 i% {    还有,如果我进行傅里叶分解是为了选择频点在频域内计算,计算完了以后把频域的结果再反变换到时域,那么在傅里叶分解的时候是否还要乘以2/N?如果乘了,在反变换之前是否要对频域内各频点下的结果除以2/N再反变换?
    - \+ Z$ S" d% \# C1 _' t2 a5 o    非常感谢
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    无聊
    2020-12-23 16:05
  • 签到天数: 34 天

    连续签到: 1 天

    [LV.5]常住居民I

    累计签到:34 天
    连续签到:1 天
    发表于 2012-5-17 10:57:52 | 显示全部楼层
    回复 9# maisigete - ]$ ?& `+ u2 r0 o+ D$ b

    " f) w' m1 A) Y+ i( n为了得到准确的幅值才乘以2/N的。实际的DFT的定义有两种,一种是在正变换里除以N,另一种是在逆变换里除以N。应该说DFT只是一种数学的定义,它可以分析频率成分(实际分析中存在频谱混叠,所以比较精确的场合需要进行校正),但得出的相应频率成分的幅值时需要修正,如乘以2/N。至于这个系数放在正变换里还是逆变换里,关系不大,只要对应上就可以了。
    0 v: x+ A& \4 e) ^% n1 \你可以看下《离散频谱分析校正理论与技术》这本书,个人感觉不错。
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
    您需要登录后才可以回帖 登录 | 立即加入

    本版积分规则

    招聘斑竹

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

    GMT+8, 2025-4-5 12:35

    Powered by Discuz! X3.5 Licensed

    © 2001-2025 Discuz! Team.

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