小波变换在谐波检测中的应用【求助】
本帖最后由 hitdely 于 2011-4-18 16:34 编辑%%参考书====时频分析技术及其应用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');=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 而第二个第三个分量是什么意思 谢谢各位 帮你顶一下,期待小波高手回答 {:{40}:} 看不懂,希望高手多多指教 逼近系数’是低频信号 ;‘细节系数’是高频信号;
a1d1分别是第一层(分解原信号s)分解后的低频和高频信号
a2d2分别是第二层(第一层分解后的低频信号a1)分解后的低频和高频信号
a3d3分别是第三层(第二层分解后的低频信号a2)分解后的低频和高频信号
a1是0-200赫兹的信号
d1是200-400赫兹的信号
a2是0-100赫兹的信号
d2是100-200赫兹的信号
a3是0-50赫兹的信号
d3是50-100赫兹的信号
d3 回复 6 # wangyanfang678 的帖子
哦 谢谢 有点明白了 回复 6 # wangyanfang678 的帖子
你好 请问如何从小波变换中提取各分量的幅值 谢谢 为什么我还没升级
页:
[1]