|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
fs=512
源程序:
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的变化进行变化 这是怎么回事?
|
|