小波变换在谐波检测中的应用【求助】

hitdely

%%参考书====时频分析技术及其应用

clear;

clc;

close all;

fs=1024;
%采样频率

f1=100;

f2=200;

t=0:1/fs:0.2;

s=sin(2*pi*f1*t)+sin(2*pi*50*t)+sin(2*pi*f2*t);
%不同频率正弦信号合成的仿真信号

figure;

plot(t,s);

ylabel('幅值/V');

xlabel('时间/s');

[c,l]=wavedec(s,3,'dmey');%db系列 dmey

%对信号做层数为3的多尺度分解

a3=wrcoef('a',c,l,'dmey',3);

a2=wrcoef('a',c,l,'dmey',2);

a1=wrcoef('a',c,l,'dmey',1);

%显示逼近系数

figure;

subplot(3,2,1);

plot(t,a3,'LineWidth',2);

ylabel('a3');

subplot(3,2,3);

plot(t,a2,'LineWidth',2);

ylabel('a2');

subplot(3,2,5);

plot(t,a1,'LineWidth',2);

ylabel('a1');

xlabel('时间/t');

%重构第1-3层细节系数

% d1=a1-a2;

% d2=a2-a3;

% d3=a3;

d3=wrcoef('d',c,l,'coif3',3);

d2=wrcoef('d',c,l,'coif3',2);

d1=wrcoef('d',c,l,'coif3',1);

%显示细节系数

subplot(3,2,2);

plot(t,d3,'LineWidth',2);

ylabel('d3');

subplot(3,2,4);

plot(t,d2,'LineWidth',2);

ylabel('d2');

subplot(3,2,6);

plot(t,d1,'LineWidth',2);

ylabel('d1');

xlabel('时间/t');

请大家帮忙分析一下 我不清楚‘逼近系数’ 与‘细节系数’是什么意思？如何对应50Hz 100Hz和200Hz 第一个应该是50Hz 而第二个第三个分量是什么意思 谢谢各位

chenlu1986

帮你顶一下，期待小波高手回答

jiyae16241626

初学振动分析

看不懂，希望高手多多指教

wangyanfang678

逼近系数’是低频信号 ；‘细节系数’是高频信号；
a1d1分别是第一层（分解原信号s）分解后的低频和高频信号
a2d2分别是第二层(第一层分解后的低频信号a1)分解后的低频和高频信号
a3d3分别是第三层(第二层分解后的低频信号a2)分解后的低频和高频信号

wangyanfang678

a1是0-200赫兹的信号
d1是200-400赫兹的信号
a2是0-100赫兹的信号
d2是100-200赫兹的信号
a3是0-50赫兹的信号
d3是50-100赫兹的信号


d3

hitdely

回复 6 # wangyanfang678 的帖子

哦 谢谢 有点明白了

hitdely

回复 6 # wangyanfang678 的帖子

你好 请问如何从小波变换中提取各分量的幅值 谢谢

yummyking

为什么我还没升级