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

 找回密码
 立即加入
搜索
查看: 1844|回复: 2

关于MATLAB的fft变换

 荐 [复制链接]

该用户从未签到

尚未签到

发表于 2009-6-24 17:52:25 | 显示全部楼层 |阅读模式

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

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

×
ATP仿真得到的VFTO波形见附件,变换程序如下。
# c; P1 Y* N' v) |* ]! \目的是求取VFTO波形在50MHz范围内的幅频曲线
" T. @2 ^5 M, `怎么感觉结果怪怪的?
# m7 m- o* k6 Z. h1 e* h, x5 f6 Z% U, Q% V) C
clear all;$ W7 b; S/ R, [2 l0 Y5 N! {
clc;
8 A; ], U. g$ G% {$ ~8 {! u4 A
8 v5 X! f. p" N9 A%Reading the inout waveform
% o: m% `' [4 w$ zfid=fopen('wave.dat','r+');
/ W5 D- _1 P- s% \D=fscanf(fid,'%f');8 o$ Q% a; g: K, p0 `. L+ q4 A

7 R, E( ?; s# p; ~% W' Cfor i=1:numel(D)/2" h. K; }" J' ~0 g2 z, C5 n6 E
    t(i)=D(2*i-1);
/ L5 Z" H8 N! W. t    x(i)=D(2*i);+ G& b3 C! v* [' H( x/ K& `
end    & ], g4 O: J9 _2 P/ ^; |2 U

. p8 X  [$ ^- c3 Z  C3 l  K' q( u" |
fs=500e+6;6 E/ g  g+ x3 B* J) ~
p=log2(numel(D)/2);8 v4 ~& Y( e  L, t
number=2^p; %65536;8 P+ T) C$ Y, @* A2 Z: C5 l
/ r9 x; y' R( R3 m
y=fft(x,number);
( Y6 X1 l. ~* g( in=length(y);" E2 O' Z8 @7 S+ x) ~
m=abs(y)*2/n;0 d: M% T* L3 X( `' d* L! ^
f=(0:n/2-1)*fs/n;
  I1 a, C/ e0 \* F; |' Jplot(f, m(1:n/2),'k');

wave.rar

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

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

该用户从未签到

尚未签到

 楼主| 发表于 2009-6-25 08:57:17 | 显示全部楼层
hehe,大家都忙,没人答理呢  ^. ~  n' I7 }% G2 x* j& |
欢迎讨论
"真诚赞赏,手留余香"
还没有人打赏,支持一下
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • TA的每日心情
    开心
    2018-12-28 19:16
  • 签到天数: 1 天

    连续签到: 1 天

    [LV.1]初来乍到

    累计签到:1 天
    连续签到:1 天
    发表于 2009-6-28 19:53:13 | 显示全部楼层
    我感觉程序没有错,只是直流分量应该除以n而不是n/2
    "真诚赞赏,手留余香"
    还没有人打赏,支持一下
    帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
    您需要登录后才可以回帖 登录 | 立即加入

    本版积分规则

    招聘斑竹

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

    GMT+8, 2026-10-9 07:35

    Powered by Discuz! X3.5 Licensed

    © 2001-2026 Discuz! Team.

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