声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1512|回复: 5

[FFT] 包络图中看不到 载波频率和边频

[复制链接]
发表于 2011-1-4 14:35 | 显示全部楼层 |阅读模式

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

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

x
F:\未命名.BMPsampleFreq   = 1000;%采样频率
sampleLength =1000;%采样点数
t=0:1/sampleFreq:(sampleLength-1)/sampleFreq;
Data=sin(2*pi*300*t).*sin(2*pi*50*t)+sin(2*pi*300*t);
Data1=Data+randn(1,length(Data));
time_plot_s = 0:1/sampleFreq:(sampleLength-1) / sampleFreq; %建立离散时间列向量
fft_plot_Hz = sampleFreq*(1:sampleLength/2)/sampleLength;
figure;
subplot(311)
plot(time_plot_s,Data1);
title('原始信号时域波形');
ylabel('振幅/m/s^2');
xlabel('时间/s');

subplot(312)
Y=fft(Data1);
Pyy=abs(Y)*2/sampleLength;%求幅值 乘上后面的2/N得到正确幅值
title('傅里叶变换的频谱图');
ylabel('振幅/m/s^2');
xlabel('频率/Hz');
f=(0:length(Pyy)-1)'*sampleFreq/length(Pyy);%进行对应的频率转换
plot(f(1:sampleLength/2),Pyy(1:sampleLength/2));%做频谱图

subplot(313)
y=hilbert(Data1) ;%da是滤波后的信号
y=abs(y);
y=y-mean(y);

p=abs(fft(y,sampleLength))*2/sampleLength;
title('包络图');
plot((0:sampleLength/2-1)/sampleLength*sampleFreq,p(1:sampleLength/2)); %fs是采样频率

回复
分享到:

使用道具 举报

 楼主| 发表于 2011-1-4 14:44 | 显示全部楼层
忘了写问题了 ,我的问题是在第二幅频谱图中可以看到300的载波频率与250和350的边频,
而在包络当中只看到了50赫兹的调制频率,看不到载波频率和边频,为什么呢?
怎么样才够让调制频率 、载波 边频都出现在包络谱中
发表于 2011-1-5 15:21 | 显示全部楼层
做了Hilbert解调肯定只有调制频率呀,想要看到载波就不要解调
 楼主| 发表于 2011-1-5 16:23 | 显示全部楼层
F:\QQ
这是西南交大硕士论文的硕士论文得到的图并且附上了解释
我是照着他做的,一样的信号,却得不到他这样的效果图
 楼主| 发表于 2011-1-5 17:01 | 显示全部楼层
本帖最后由 dongyun.wang 于 2011-1-5 17:06 编辑

这是西南交大硕士论文中的保包络图文中说到
“从仿真信号的频谱图但是从图中找不到调制频率,但是在包络谱中我们可以清晰的看到调制频率及其倍
,载波频率和边频被抑制,因此我们可以得出结论,用包络解调方法处理故障轴承信号提取其中的调制频率的方案是可行的”。
我怎样得到他这样的包络谱图?谢谢

QQ12.bmp
发表于 2011-1-6 10:49 | 显示全部楼层
怎么看不见大图?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-5-8 15:24 , Processed in 0.143237 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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