声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2252|回复: 5

[小波] 利用小波变换进行时频分析时,如何进行选取采样频率?

[复制链接]
发表于 2013-2-25 22:55 | 显示全部楼层 |阅读模式

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

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

x

fs=512

fs=512
12.jpg 源程序:

clear;
clc;
fs=1024;  %采样频率,可根据需要设置采样频率,一般为2的整数次幂
% %f1=100;                        
% f2=200;
% t=0:1/fs:1;
R=xlsread('Gravity');%读取数据
t=R(:,1);
s=R(:,2);
%s=10*sin(2*pi*f1*t)+15*sin(2*pi*f2*t);  %两个不同频率正弦信号合成的仿真信号
%s=sin(2*pi*f1*t)/2;  %两个不同频率正弦信号合成的仿真信号
%%%%%%%%%%%%%%%%%小波时频图绘制%%%%%%%%%%%%%%%%%%
wavename='cmor3-3';
totalscal=256;                    %尺度序列的长度,即scal的长度
wcf=centfrq(wavename);            %小波的中心频率
cparam=2*wcf*totalscal;           %为得到合适的尺度所求出的参数
a=totalscal:-1:1;  
scal=cparam./a;                   %得到各个尺度,以使转换得到频率序列为等差序列
coefs=cwt(s,scal,wavename);       %得到小波系数
f=scal2frq(scal,wavename,1/fs);   %将尺度转换为频率
imagesc(t,f,abs(coefs));          %绘制色谱图
colorbar;
%colormap(flipud(gray));% 黑白显示,打印比较清晰
xlabel('时间 t/s');
ylabel('频率 f/Hz');
title('小波变换时频图');

所读取的数据中就是一些随时间变化的离散值。当采样频率fs=1024时,得到图

当采样频率fs=1024时,得到图
两张图的纵坐标的变化范围随fs的变化进行变化 这是怎么回事?

回复
分享到:

使用道具 举报

发表于 2013-2-26 09:38 | 显示全部楼层
纵坐标就是频率轴,你采样频率fs设为1024,那么频率轴的最高频率就是奈奎斯特频率是fs/2,这个和fft是一样的。
发表于 2013-2-26 10:06 | 显示全部楼层
不明白,你写的程序是一个图fs1024的那个,下面的图就应该是对的,但是第一个图你是怎么来的呢,好像是没说清楚吧?
 楼主| 发表于 2013-2-26 11:47 | 显示全部楼层

第一个图是fs=512时所绘制的图,第二个是fs=1024时所绘制的图
发表于 2013-3-12 21:55 | 显示全部楼层
本帖最后由 寂寞的部落 于 2013-3-12 22:03 编辑

我也遇到相同的问题,这个与你数据的采样频率有关,采样频率设成你数据信号的采样频率
发表于 2013-9-26 10:12 | 显示全部楼层
请问,此处的采样频率是不是收集数据时的采样频率呢,我采集数据时的采样频率是1250Hz,这里是不是要把fs设为2500呢

补充内容 (2013-10-10 09:03):
“totalscal=256;                    %尺度序列的长度,即scal的长度”
请问256如何计算得到
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-11 00:56 , Processed in 0.092947 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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