xiao4.m

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

load db8;
w=db8;
subplot(421);stem(w);grid;
title('原始尺度滤波器');
[Lo_D,Hi_D,Lo_R,Hi_R]=orthfilt(w);
subplot(423);stem(Lo_D);
title('分解的低通滤波器');grid;
subplot(424);stem(Hi_D);
title('分解的高通滤波器');grid;
subplot(425);stem(Lo_R);
title('重构的低通滤波器');grid;
subplot(426);stem(Hi_R);
title('重构的高通滤波器');grid;
df=[Lo_D;Hi_D];
rf=[Lo_R;Hi_R];
disp('检查双正交性的结果，df*rf=:');
id=df*df'
fftld=fft(Lo_D);ffthd=fft(Hi_D);
freq=[1:length(Lo_D)]/length(Lo_D);
subplot(427);plot(freq,abs(fftld));
title('转换模量：低通');grid;
subplot(428);plot(freq,abs(ffthd));
title('转换模量：高通');grid;