声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2130|回复: 4

[FFT] 求助:绘制频谱图 有点急

[复制链接]
发表于 2006-10-29 16:26 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

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

x
如何绘制一周期为T,振幅为F0的矩形波的幅值频谱图和相位频谱图.
表达式为                 F0            0<t<1/2T
                     F(t)=   -F0          -1/2T<t<T  
希望哪位大哥给出程序 自己好琢磨
回复
分享到:

使用道具 举报

发表于 2006-10-29 19:21 | 显示全部楼层
照理说应该先滤波才能做FFT,但只要谱线数足够大的话频率反折现象也可忽略了.

clear all;
fs=256;
N=256;
T=fs*N;

F0=1;
s=[F0*ones(1,N/2) -F0*ones(1,N/2)];
sf=fft(s);
sf=sf(1:N/2)/N*2;

df=fs/N;
f=(0:N/2-1)*df;

figure();
subplot(211);
stem(f,abs(sf),'marker','none');

subplot(212);
stem(f,angle(sf),'marker','none');

[ 本帖最后由 yangzj 于 2006-10-30 12:46 编辑 ]
 楼主| 发表于 2006-10-30 09:09 | 显示全部楼层
不是很清楚每句语句的意思,能给出解释么?麻烦了
发表于 2010-9-17 11:37 | 显示全部楼层
还没有找到答案,继续
发表于 2010-9-27 11:14 | 显示全部楼层
clear all;
fs=256; %采样频率
N=256; %采样点数
T=fs*N;%信号周期

F0=1;
s=[F0*ones(1,N/2) -F0*ones(1,N/2)]%前办个周期个数据为F0,后半个周期数据位-F0
sf=fft(s);%傅里叶变换
sf=sf(1:N/2)/N*2;%转换为实际的幅值

df=fs/N;%FFT频率间隔
f=(0:N/2-1)*df;%对应的实际频率

figure();%作图
subplot(211);
stem(f,abs(sf),'marker','none');%幅度谱

subplot(212);
stem(f,angle(sf),'marker','none');%相位谱

都是基本的,楼主多读些程序就好了~
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-14 15:49 , Processed in 0.062132 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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