hanzaofenjie.m

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

%产生fbm信号，H=0.4
%初始化randn发生器
randn('state',1);
%设置参数H和采样长度lg
H=0.4;
lg=1000;
x1=wfbm(H,lg);
m1=mean(x1);
x=x1-m1;
t=1:lg;
%生成白噪声信号
N=1000;
t=1:N;
y1=5*(rand(1,N)-0.5);
m2=mean(y1);
y=y1-m2;
s=x+y;
t=1:1000;
%用db4小波进行五层分解
[c,l]=wavedec(s,5,'db4');
%对分解结构[c,l]的各层低频部分进行重构
a5=wrcoef('a',c,l,'db4',5);
subplot(7,1,1);
plot(s);
subplot(7,1,2);
plot(a5);
%对分解结构[c,l]的各层高频部分进行重构
for i=1:5
    decmp=wrcoef('d',c,l,'db4',6-i);
subplot(7,1,i+2);
plot(decmp);
    ylabel(['d',num2str(6-i)]);
end