gaojian 发表于 2007-5-29 22:02

求助:各位同仁能否提供一个小波分析输出时频图的相关例子!??

如题!
感激不尽!

[ 本帖最后由 zhangnan3509 于 2007-6-1 17:06 编辑 ]

hector1982 发表于 2007-5-30 09:17

原帖由 gaojian 于 2007-5-29 22:02 发表 http://www.chinavib.com/forum/images/common/back.gif
如题!感激不尽!
典型的例子在很多相关的小波分析书中都有介绍:
我这有个小程序也是从书中得到的,可以供参考
%含噪的三角波与正弦波的组合
%生成正弦信号
clc;
clear;
close all;
N=1000;
t=0:0.001:0.6;
x=sin(2*pi*50*t)+sin(2*pi*120*t);
y=x+2*randn(size(t));
figure;
subplot(211)
plot(t,y,'LineWidth',2);
subplot(212)
plot(1000*t(1:50),y(1:50));
%xlabel('sample number n');
%ylabel('amplitude A');
%wavelet 1D decompose
=wavedec(y,7,'db10');
% reconstruct 1-7 layer approach coefficient.
a7=wrcoef('a',c,l,'db10',7);
a6=wrcoef('a',c,l,'db10',6);
a5=wrcoef('a',c,l,'db10',5);
a4=wrcoef('a',c,l,'db10',4);
a3=wrcoef('a',c,l,'db10',3);
a2=wrcoef('a',c,l,'db10',2);
a1=wrcoef('a',c,l,'db10',1);

% show approach coefficient
figure
subplot(7,1,1);
plot(a7,'LineWidth',2);
ylabel('a7');
subplot(7,1,2);
plot(a6,'LineWidth',2);
ylabel('a6');
subplot(7,1,3);
plot(a5,'LineWidth',2);
ylabel('a5');
subplot(7,1,4);
plot(a4,'LineWidth',2);
ylabel('a4');
subplot(7,1,5);
plot(a3,'LineWidth',2);
ylabel('a3');
subplot(7,1,6);
plot(a2,'LineWidth',2);
ylabel('a2');
subplot(7,1,7);
plot(a1,'LineWidth',2);
ylabel('a1');
xlabel('sample sequence n');
% reconstruction 1-7 layer details coefficient
d7=wrcoef('d',c,l,'db10',7);
d6=wrcoef('d',c,l,'db10',6);
d5=wrcoef('d',c,l,'db10',5);
d4=wrcoef('d',c,l,'db10',4);
d3=wrcoef('d',c,l,'db10',3);
d2=wrcoef('d',c,l,'db10',2);
d1=wrcoef('d',c,l,'db10',1);
% show detail coefficient
figure
subplot(7,1,1);
plot(d7,'LineWidth',2);
ylabel('d7');
subplot(7,1,2);
plot(d7,'LineWidth',2);
ylabel('d6');
subplot(7,1,3);
plot(d5,'LineWidth',2);
ylabel('db5');
subplot(7,1,4);
plot(d5,'LineWidth',2);
ylabel('d4');
subplot(7,1,5);
plot(d3,'LineWidth',2);
ylabel('d3');
subplot(7,1,6);
plot(d2,'LineWidth',2);
ylabel('d2');
subplot(7,1,7);
plot(d1,'LineWidth',2);
ylabel('d1');
xlabel('sample squence n');

%%%%%%%%%%%%%%%hilbert
fs=1000;
Y2=fft(a1,512);
Pyy2=Y2.*conj(Y2)/512;
f2=1000*(0:256)/512;
y=hilbert(a1);
ydata=abs(y);
y=y-mean(y);
nfft=1000;
p=abs(fft(ydata,nfft));

y1=hilbert(d1);
ydata1=abs(y1);
y1=y1-mean(y1);
nfft=1000;
p1=abs(fft(ydata1,nfft));


figure
subplot(311)
plot((0:nfft/2-1)/nfft*fs,p(1:nfft/2));
subplot(312)
plot((0:nfft/2-1)/nfft*fs,p1(1:nfft/2));
subplot(313)
plot(f2,Pyy2(1:257));

[ 本帖最后由 hector1982 于 2007-5-30 09:36 编辑 ]

form 发表于 2007-5-30 15:26

本帖最后由 wdhd 于 2016-9-19 13:03 编辑

原帖由 hector1982 于 2007-5-30 09:17 发表
%含噪的三角波与正弦波的组合
%生成正弦信号
clc;
clear;
close all;
N=1000;
t=0:0.001:0.6;
x=sin(2*pi*50*t)+sin(2*pi*120*t);
y=x+2*randn(size(t));
figure;
subplot(211)
plot(t,y,'LineWidth ...
这个好像不是搂住需要的,粗略翻了下书,help wpviewcf。找本书看看

zhlong 发表于 2007-5-30 15:37

楼主要得到时频图可以用 tftb时频工具箱中的 tfrscalo()函数

hector1982 发表于 2007-5-30 19:51

本帖最后由 wdhd 于 2016-9-19 13:03 编辑

原帖由 zhlong 于 2007-5-30 15:37 发表
楼主要得到时频图可以用 tftb时频工具箱中的 tfrscalo()函数
tfrscalo() 是matlab中自带的函数吗?怎么我的matlab7.0没有啊?

zhlong 发表于 2007-5-30 19:56

tftb时频工具箱 可以到E盘下载 http://free.ys168.com/?zhlong

hector1982 发表于 2007-5-30 20:01

xiexie

zhlong 发表于 2007-5-30 20:06

下面是我做的关于时频分析方法比较的一个例子,信号由3个线性调频信号组成,其中一个还调幅,另外两个频率很接近。图1(a)上面是信号的波形,下面3个是它的组成成份的时域波形。做时频分析用到的函数都是tftb时频工具箱中的,里面有连续小波变换(CWT)的时频图。

gaojian 发表于 2007-5-31 00:03

回复 #8 zhlong 的帖子

感谢各位的帮助!
谢谢!

sky163 发表于 2007-5-31 16:44

小波的时频图也是用 tfrscalo()函数实现的吗?
那个函数不是用在画STFT的时频图吗?

STFT做时频分析时用specgram()函数可以吗?

zhlong 发表于 2007-5-31 16:53

STFT做时频分析时用specgram()函数可以吗?
可以
小波的时频图也是用 tfrscalo()函数实现的吗
是的
%        =TFRSCALO(X,T,WAVE,FMIN,FMAX,N,TRACE) computes
%        the scalogram (squared magnitude of a continuous wavelet
%        transform).

sky163 发表于 2007-5-31 21:31

再请教下:用tfrscalo()函数做实现小波时频分析,编程的时候步骤是怎样的?

zhlong 发表于 2007-5-31 21:45

回复 #12 sky163 的帖子

下载tftb时频工具箱,help tfrscalo

jing123 发表于 2007-5-31 22:26

好像会一直提示错误

jing123 发表于 2007-5-31 22:42

图里的STFT和小波的时频图是怎么弄出来的,

可不可以把程序贴出来让大家看一下
页: [1] 2
查看完整版本: 求助:各位同仁能否提供一个小波分析输出时频图的相关例子!??