zuoye.m

局部放电的n-q-phi和dltU处理程序

load sinsin
%画出原始图像
figure(1)
subplot(111);image(X);colormap(map);
title('原始图像');
axis square
%尺度为2，小波函数为db1时的尺度分解
[c,s]=wavedec2(X,2,'db1');
sizex=size(X);
sizec=size(c)
val_s=s
%从不同的方位从分解结构[c,s]中提取尺度2的系数
chd2=detcoef2('h',c,s,2);
cvd2=detcoef2('v',c,s,2);
cdd2=detcoef2('d',c,s,2);
sizecd2=size(chd2)
ca2=appcoef2(c,s,'db1',2);
figure(2)
%画出尺度2的低频信号图像
sizeca1=size(ca2)
subplot(221);image(ca2);colormap(map);
title('图像的低频部分 ');
axis square
%画出尺度2的高频信号图像
subplot(222);image(chd2);
title('高频的水平部分 ');
axis square
subplot(223);image(cvd2);
title('高频的垂直部分 ');
axis square
subplot(224);image(cdd2);
title('高频的斜线部分 ');
axis square