matlab信号处理小波变换


% 定义信号信息
fs=2^6;    %采样频率
dt=1/fs;    %采样间隔
timestart=-8;
timeend=8;
t=(0:(timeend-timestart)/dt-1)*dt+timestart;
L=length(t);
z=4*sin(2*pi*linspace(6,12,L).*t);%为待变换信号
%matlab自带的小波变换
figure(1)
[wt,f,coi] = cwt(z,'amor',fs);
pcolor(t,f,abs(wt));shading interp

clc
clear all
close alldata12_20=csvread('E:\experiment\5-5\12-1.csv',22,0);
t12=data12_20(:,1);
y12_20=data12_20(:,2);
z12_20=data12_20(:,3);dt=data12_20(2,1)-data12_20(1,1);    %采样间隔
fs=1/dt;   %采样频率
L=length(t12);
figure
plot(t12, z12_20)
%方法一
wavename='cmor3-3';
totalscal=256;
Fc=centfrq(wavename); % 小波的中心频率
c=2*Fc*totalscal;
scals=c./(1:totalscal);
f=scal2frq(scals,wavename,1/fs); % 将尺度转换为频率
coefs=cwt(z12_20,scals,wavename); % 求连续小波系数
figure
imagesc(t12,f,abs(coefs));
set(gca,'YDir','normal')
colorbar;
xlabel('时间 t/s');
ylabel('频率 f/Hz');
title('方法一');
ylim([0,1e6]);
%方法二
figure(1)
[wt,f,coi] = cwt(z12_20,'amor',fs);
pcolor(t12,f,abs(wt));shading interp
ylim([0,1e6]);
title('方法二');%方法三
[coefs_healthy,frequencies] = cwt(z12_20,1:256,'morl',dt);
figure
surf(t12,frequencies,abs(coefs_healthy))
shading flat
colormap('jet')
view([0,0,1])
ylim([0,1e6]);
title('方法三');

fls的代码

clc
clear all
% close alldata12_1=csvread('.\5-6-2\12-1.csv',22,0);
t12_1=data12_1(:,1);
y12_1=data12_1(:,2);data12_2=csvread('.\5-6-2\12-2.csv',22,0);
t12_2=data12_2(:,1);
y12_2=data12_2(:,2);data34_1=csvread('.\5-6-2\34-1.csv',22,0);
y34_1=data12_1(:,2);data12_2=csvread('.\5-6-2\34-2.csv',22,0);
y34_2=data12_2(:,2);dt=data12_1(2,1)-data12_1(1,1);    %閲囨牱闂撮殧
fs=1/dt;   %閲囨牱棰戠巼
% figure
% plot(t12_1, y12_1)[coefs_1,frequencies] = cwt(y12_1,4:4:1024,'morl',dt);
[coefs_2,~] = cwt(y12_2,4:4:1024,'morl',dt);
[coefs_3,~] = cwt(y34_1,4:4:1024,'morl',dt);
[coefs_4,~] = cwt(y34_2,4:4:1024,'morl',dt);coefs_1h =zeros(256,9999);
coefs_2h =zeros(256,9999);
coefs_3h =zeros(256,9999);
coefs_4h =zeros(256,9999);
for i=1:256coefs_1h(i,:) = abs(hilbert(coefs_1(i,:)));coefs_2h(i,:) = abs(hilbert(coefs_2(i,:)));coefs_3h(i,:) = abs(hilbert(coefs_3(i,:)));coefs_4h(i,:) = abs(hilbert(coefs_4(i,:)));
endaa=abs(coefs_1h-coefs_2h)./max(max(coefs_1h));
bb=abs(coefs_3h-coefs_4h)./max(max(coefs_3h));
figure
% subplot(2,1,1)
% surf(t12_1,frequencies,abs(coefs_1h-coefs_2h))
surf(t12_1,frequencies,coefs_1h)
shading interp;
view([0,0,1])
ylim([0,1e6]);
%
% subplot(2,1,2)
% % surf(t12_1,frequencies,abs(coefs_3h-coefs_4h))
% surf(t12_1,frequencies,bb)
% shading interp;
% view([0,0,1])
% ylim([0,1e6]);% figure
% subplot(3,1,1)
% surf(t12_1,frequencies,abs(coefs_1))
% shading interp;
% %shading flat
% %colormap('jet')
% view([0,0,1])
% ylim([0,1e6]);
% % title('12缁勫皬娉㈠彉鎹 );
% %set(gca,'CLim',[a,b])
% % colorbar;caxis(temp1);
% subplot(3,1,2)
% surf(t12_1,frequencies,abs(coefs_2))
% shading interp;
% %shading flat
% %colormap('jet')
% view([0,0,1])
% ylim([0,1e6]);
% subplot(3,1,3)
% surf(t12_1,frequencies,abs(coefs_1-coefs_2))
% shading interp;
% %shading flat
% %colormap('jet')
% view([0,0,1])
% ylim([0,1e6]);

参考文章
参考文章

matlab信号处理小波变换相关推荐

MATLAB信号处理工具箱函数列表分类
**现将MATLAB信号处理工具箱函数进行分组,便于记忆查询和长期回顾.(只解释基本用途,具体用法请在help目录下查询)** Waveform Generation(波形产生) chairp: 产生 ...
matlab中udt函数,《MATLAB信号处理超级学习手册》——2.5 离散时间信号中的运算...
本节书摘来自异步社区<MATLAB信号处理超级学习手册>一书中的第2章,第2.5节,作者:MATLAB技术联盟 , 史洁玉著,更多章节内容可以访问云栖社区"异步社区"公 ...
matlab 信号处理教程,MATLAB信号处理仿真实验_教程-学习文件.pdf
MATLAB信号处理仿真实验从观察正弦波开始用绳量给我的地界坐落在佳美之处我的产业实在美好杜伟韬 duweitao@ 广播电视数字化教育部工程研究中心 2013年 4月23 日于北京定福庄 ...
《MATLAB信号处理超级学习手册》——2.7 本章小结
本节书摘来自异步社区<MATLAB信号处理超级学习手册>一书中的第2章,第2.7节,作者:MATLAB技术联盟 , 史洁玉著,更多章节内容可以访问云栖社区"异步社区"公 ...
MATLAB信号处理---学习小案例(10)---Z反变换
定义X(z)X(z)X(z)的Z反变换(IZT)为x(n)=12πj∮CX(z)zn−1dzx(n)=\frac{1}{2\pi j} \oint_{C}X(z)z^{n-1}dzx(n)=2πj1 ...
MATLAB信号处理——信号与系统的分析基础(1)
1-1 信号的概念信号是传递信息的函数,可以表示为一个或几个独立变量的函数.按照时间.幅值的连续与离散可以分为: (1)连续时间信号:时间连续,幅值可以连续也可以离散 (2)模拟信号:时间连续,幅值 ...
MATLAB信号处理——信号与系统的分析基础（2）
z1-3 离散时间序列 1)单位采样序列单位采样序列(也叫单位脉冲序列)δ(n)定义为它与单位冲激函数(狄拉克函数)很像,不同的是,狄拉克函数当n=0时,取值为无穷大. 脉冲序列可以通过下面命令来 ...
MATLAB基于小波变换的语音信号去噪算法改进
MATLAB基于小波变换的语音信号去噪算法改进概述 0. 需要调用的子函数 0.1 Gnoisegen函数 0.2 snrr函数 1. 语音信号输入和加噪 1.1 语音信号输入 1.2 语音信号加噪 ...
MATLAB信号处理——信号与系统的分析基础（3）
1-4信号的基本运算前面说过,序列是一种函数,函数我们就可以对它进行运算,而且很多运算法则都是一些基础的函数也有的. 1)序列相加与相乘相加:在对应的点处把对应的序列值相加即可 ...
MATLAB信号处理——信号与系统的分析基础（5）
1-6 LTI(线性时不变系统)的时域分析 LTI具有4个基本性质:1.线性:2.时不变性:3.微分性:4.积分性,它可以用下面的线性常微分方程来描述.https://baike.baidu.com/ ...

matlab信号处理小波变换

matlab信号处理小波变换相关推荐

最新文章

热门文章