canshuguji.m

小波变换用于半导体激光器可靠性分析的全部程序。我只完成了其中一部分

%装载实际测量信号
load e:\data\vfa020.dat;   
s=vfa020(512:1000); 
%用db10小波对合成信号进行5层分解
[c,l]=wavedec(s,5,'db10');;
%从小波分解结构[c,l]中提取低频概貌系数
ca5=appcoef(c,l,'db10',5);
%从小波分解结构[c,l]中提取高频细节系数
cd5=detcoef(c,l,5);
cd4=detcoef(c,l,4);
cd3=detcoef(c,l,3);
cd2=detcoef(c,l,2);
cd1=detcoef(c,l,1);
%求合成信号小波系数的样本方差
v5=sum((cd5).^2)/length(cd5);
v4=sum((cd4).^2)/length(cd4);
v3=sum((cd3).^2)/length(cd3);
%估计谱参数r的值
r=log2(abs((v5-v4)/(v4-v3)));
H=(r-1)/2;H