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

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

[讨论] 关于MATLAB的fft变换

  [复制链接]

该用户从未签到

尚未签到

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

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

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

×
ATP仿真得到的VFTO波形见附件,变换程序如下。
6 g  p. }1 A4 g: S目的是求取VFTO波形在50MHz范围内的幅频曲线5 w. a) ^1 ^( m1 ~
怎么感觉结果怪怪的?
/ Y# _3 F- z5 M% I7 }3 a- ?" |- A1 Q( h4 v0 ^. `
clear all;
$ `- ^! e: Z1 c. D/ ?6 Eclc;
, `" R8 B3 I# I! }- m* m6 v: u" D- V* ^
%Reading the inout waveform8 ^7 U8 C; b" [, }- Z$ {
fid=fopen('wave.dat','r+');
5 W0 A7 q$ D3 V2 SD=fscanf(fid,'%f');
! n5 f- z0 ]4 ?% J+ C' Z+ Q7 q. t7 R* T0 o% b- u6 R4 ~' B# ^
for i=1:numel(D)/2. L* E: S7 v3 W) D. W8 X
    t(i)=D(2*i-1);
5 x8 M/ b  w& ~! {    x(i)=D(2*i);0 ~/ u8 @9 u" q4 ~- @% W. l
end   
& ]# i" n9 }0 C2 \0 X: ?2 l8 B1 d: y$ W6 X! ~  ^. g2 I  `: ~1 l
2 @/ P$ _: `* [# G3 X. f& B
fs=500e+6;
& y5 Y3 y, C8 a  s1 fp=log2(numel(D)/2);
5 G6 T' _& w$ v# V; ynumber=2^p; %65536;  K5 Q" k) Z# P3 S; \
, \& T' n) P2 E6 Z. o& \+ d* d
y=fft(x,number);" R& L6 ?  ~8 z- }! _+ o# I
n=length(y);, g% }( B9 e8 S
m=abs(y)*2/n;
( e' ]$ ^7 Y, i5 _1 pf=(0:n/2-1)*fs/n;
! d: r+ o- P$ pplot(f, m(1:n/2),'k');

wave.rar

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

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

该用户从未签到

尚未签到

 楼主| 发表于 2009-6-25 08:57:17 | 显示全部楼层
hehe,大家都忙,没人答理呢
, B! o+ K" [3 o. O! R欢迎讨论
"真诚赞赏,手留余香"
还没有人打赏,支持一下
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
  • 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, 2025-6-7 03:14

    Powered by Discuz! X3.5 Licensed

    © 2001-2025 Discuz! Team.

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