xuanze.m

来自「小波变换用于半导体激光器可靠性分析的全部程序。我只完成了其中一部分」· M 代码 · 共 32 行

%装载实际测量信号
load e:\data\vfa043.dat;   
s=vfa043(512:1000);
ls=length(s);
%画出实际测量信号波形
figure(1);
plot(s,'LineWidth',2); 
%用db4小波对信号进行6层分解
[c,l]=wavedec(s,6,'db4');
%对分解结构[c,l]的各层低频部分进行重构
figure(2);
subplot(7,1,1);
plot(s);
ylabel('s');
for i=1:6
    decmp=wrcoef('a',c,l,'db4',7-i);
subplot(7,1,i+1);
plot(decmp);
    ylabel(['a',num2str(7-i)]);
end

%对分解结构[c,l]的各层高频部分进行重构
figure(3);
subplot(7,1,1);
plot(s);
ylabel('s');
for i=1:6
    decmp=wrcoef('d',c,l,'db4',7-i);
subplot(7,1,i+1);
plot(decmp);
    ylabel(['d',num2str(7-i)]);
end