声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1251|回复: 4

[FFT] 矩形波的幅频图怎么画啊

[复制链接]
发表于 2009-9-22 19:06 | 显示全部楼层 |阅读模式

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

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

x
t=0:0.001:1;
x = square(2*pi*10*t,50);
plot(t,x)

大家帮帮忙啊,由这几句程序生成的矩形波,它的幅频图怎么画啊,程序怎么编写啊?怎么也做不出来了:@L ,请大侠们帮帮忙吧,谢谢了:handshake !
回复
分享到:

使用道具 举报

发表于 2009-9-23 08:47 | 显示全部楼层
幅频图的绘制与波形是什么有关吗?
建议动手搜下
 楼主| 发表于 2009-9-23 11:21 | 显示全部楼层

回复 沙发 ChaChing 的帖子

tau_T=1/4;                                            % 占空比1/4
omega_0=2*pi;                                         % 基波频率
n0=0;n1=15;
n=n0:n1;
F_n=[(sin(tau_T*pi*n))/n*pi]*exp(-j*tau_T*n.*pi);
F_n=2*F_n;
Fn_max=max(abs(F_n));
Fn_min=min(abs(F_n));
subplot(2,1,1),stem(n,abs(F_n),'.');            
axis([n0 n1 Fn_min-0.1 Fn_max+0.1]);
line([n0 n1],[0 0],'color','r');                 
title('单边幅度频谱');   
请问一下,我这个程序怎么画出来幅度是一个值啊,哪里错了呢?
 楼主| 发表于 2009-9-23 11:25 | 显示全部楼层

回复 沙发 ChaChing 的帖子

t = 0:1/1000:0.05;
y = square(2*pi*100*t,20);%f=100hz
Y=fft(y,512);%1024为采样点数
df = 1000/1024;%频率分辨率
f=df*(0:511)
n=1:512
plot(f,abs(Y))
ylable('magnitude')
xlabel('frequency (Hz)')
周期函数的频谱应该是离散的,我这个画出来的怎么不是离散的呢,还有幅值好像不对啊,到底哪里错了啊,初学,请帮帮忙吧!我上网查了可是还是解决不了。
发表于 2009-9-23 11:42 | 显示全部楼层
个人水平专业有限, 建议楼主看看这里有没有你要问的问题--信号处理专栏话题索引
http://forum.vibunion.com/forum/thread-49413-1-1.html
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-13 19:23 , Processed in 0.064843 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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