补充：频域特征值提取的MATLAB代码实现（小波分析）

之前的文章信号频域分析方法的理解（频谱、能量谱、功率谱、倒频谱、小波分析）中提到了离散小波分解的例子，其参考代码如下：

t_s = 0.005; %采样周期
t_start = 0.001; %起始时间
t_end = 10;     %结束时间
t = t_start : t_s : t_end;
y = 10*sin(2*pi*0.5*t)+3*sin(2*pi*10*t)+1*sin(2*pi*20*t)+3*randn(1,length(t));  %生成信号
len = length(y);
%生成突变信号
y2 = 50*sin(2*pi*50*t);
for i = 1: lenif i>=601&&i<=604y(i) = y(i)+y2(i);elsey(i) = y(i);end
end
figure
plot(y) %绘制原始信号
[c,l] = wavedec(y,5,'db5');
%重构1~5层细节函数
d5 = wrcoef('d',c,l,'db5',5);
d4 = wrcoef('d',c,l,'db5',4);
d3 = wrcoef('d',c,l,'db5',3);
d2 = wrcoef('d',c,l,'db5',2);
d1 = wrcoef('d',c,l,'db5',1);
%重构1~5层近似函数
a5 = wrcoef('a',c,l,'db5',5);
a4 = wrcoef('a',c,l,'db5',4);
a3 = wrcoef('a',c,l,'db5',3);
a2 = wrcoef('a',c,l,'db5',2);
a1 = wrcoef('a',c,l,'db5',1);
figure
subplot(4,2,1);
plot(a1)
subplot(4,2,2);
plot(d1)
subplot(4,2,3);
plot(a2)
subplot(4,2,4);
plot(d2)
subplot(4,2,5);
plot(a3)
subplot(4,2,6);
plot(d3)
subplot(4,2,7);
plot(a4)
subplot(4,2,8);
plot(d4)

结果如下：

原始信号

离散小波分析结果

左侧四行是1~4阶的近似信号，右侧四行是1~4阶的细节信号。

欢迎关注我的专栏与信号处理有关的那些东东

欢迎关注我的公众号“括号的城堡”，微信号为“khscience”，会有更多有趣的东西分享。

补充：频域特征值提取的MATLAB代码实现（小波分析）相关推荐

matlab小波分析工具箱原理与应用_补充：频域特征值提取的MATLAB代码实现（小波分析）...
之前的文章信号频域分析方法的理解(频谱.能量谱.功率谱.倒频谱.小波分析)中提到了离散小波分解的例子,其参考代码如下: t_s 结果如下: 原始信号离散小波分析结果左侧四行是1~4阶的近似信号,右 ...
matlab频域怎么提取特征值,补充：频域特征值提取的MATLAB代码实现（小波分析）...
之前的文章信号频域分析方法的理解(频谱.能量谱.功率谱.倒频谱.小波分析)中提到了离散小波分解的例子,其参考代码如下: t_s = 0.005; %采样周期 t_start = 0.001; %起始时 ...
时域特征值提取的MATLAB代码实现（均方根、峰值因子、脉冲因子、裕度因子、峭度因子、波形因子和偏度等）
时域计算代码相对简单,下列代码主要来自求助matlab时域分析程序-峰值.均方根.裕度.峭度_百度知道.不过原文中求"峭度因子"的方法不对,即"Kr = sum(y.^4 ...
频域特征指标及其MATLAB代码实现（重心频率、均方频率、均方根频率、频率方差、频率标准差）
以前我们梳理过常用的有量纲时域特征指标和无量纲时域特征指标,其实特征提取方法中,一些频谱特征指标也是十分有用的. 下述频域特征指标是根据信号的功率谱求得的.在故障诊断或者其他应用场景下,可以与时域指标 ...
【图像分割】基于计算机视觉实现胸部CT肺质提取附matlab代码
1 内容介绍在现代医学领域中,医学影像处理技术随着计算机科学和影像技术的进步,已经成为医学领域重要的一个分支.室外光照度不均.CT自身空间分辨率和层厚参数.人体组织器官蠕动等诸多外界因素造成了医学X ...
（论文加源码）基于时频域特征分析和SVM分类器的DEAP脑电信号情感状态识别（matlab代码）（四分类）
论文和源码见个人主页:https://download.csdn.net/download/qq_45874683/85069822 (论文加源码)(matlab代码)基于时频域特征分析和SVM分类器 ...
matlab 提取最小值,MATLAB中取最小值的代码
matlab 最小二乘法平面拟合的代码 p=[xyones(length(x),1)]\z这样得到的就是平面z=p(1)*x+p(2)*y+p(3)祝你学习快乐!再问:xyz轴中如何显示所拟合的平面 ...
fortran转换 matlab代码,将Fortran77代码转换为Matlab代码以查找特征值/向量
我将Fortran 77中的书面代码转换为Matlab代码.该函数使用QL算法计算矩阵的特征值和特征向量.由于某些原因,我不能在matlab中使用eig函数的结果.这种方法得到的特征值与eig函数得到 ...
频域平滑滤波 matlab,用MATLAB实现频域平滑滤波以及图像去噪代码
用MATLAB实现频域平滑滤波以及图像去噪代码提问时间:2018-01-18 19:06:50 是数字图象处理的实验,麻烦高人给个写好的代码,希望能在重要语句后面附上一定的说明,只要能在MATLAB ...

补充：频域特征值提取的MATLAB代码实现（小波分析）

补充：频域特征值提取的MATLAB代码实现（小波分析）相关推荐

最新文章

热门文章