資源簡介
復小波變換提取小波脊線,可較好識別頻率
代碼片段和文件信息
clc;clear;close?all
t1=0:1/40:7.5-1/40;
t2=7.5:1/40:25.6-1/40;
t=[t1t2];
x1=sin(1.6*2*pi*t1)+sin(3*2*pi*t1+0.04*2*pi*t1.^2+4*sin(0.1*2*pi*t1));
x2=sin(1.6*2*pi*(t2-7.5)+0.055*2*pi*(t2-7.5).^2)+sin(3*2*pi*t2+0.04*2*pi*t2.^2+4*sin(0.1*2*pi*t2));
x=[x1x2];
N=length(x);
y=awgn(x11);
fs=40;fc=3;totalscal=200;
FreqBins=linspace(1/fs6/fstotalscal);
fw=FreqBins*fs;
scales=fc./FreqBins;
coes_y=cwt(yscales‘cmor?4-3‘);
imagesc(tscalesabs(coes_y));
figure(2)
imagesc(tfwabs(coes_y));
a=zeros(1N);
b=zeros(1N);
b_p=zeros(1N);
[a(1)b(1)]=max(abs(coes_y(81:1211)));
temp1=80;b_p(1)=b(1)+temp1;
for?i=2:N
????[a(i)b(i)]=max(abs(coes_y(temp1+b(i-1)-10:temp1+b(i-1)+10i)));
????temp1=b_p(i-1)-10-1;
????b_p(i)=b(i)+temp1;
end
figure(3)
plot(tfw(b_p))
c=zeros(1
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件??????84917??2010-08-12?11:04??原始信號小波量圖.png
?????文件??????38786??2010-08-12?11:04??提取出的小波脊線.png
?????文件???????1846??2010-08-10?09:19??wtsimulation.m
-----------?---------??----------?-----??----
???????????????125549????????????????????3
- 上一篇:matlab程序仿真信號SG平滑濾波
- 下一篇:gps基本原理及其matlab仿真
評論
共有 條評論