test.asv

来自「这是一个高频滤波程序」· ASV 代码 · 共 18 行

%high frequency enhancement 1-0.95(Z-1)
clear;
load pcsym10.dat
%函数调用的方法实现滤波
 B=[1 -0.95];    
 A=1;
 gaopints1=filter(B,A,pcsym10);
%循环的方法实现滤波
%gaopints(1,1)=0;
%n=length(pcsym10)
%for k=2:n
  % gaopints(1,k)=pcsym10(1,k)-0.95*pcsym10(1,k-1);
%end

subplot(211)
plot(pcsym10)
subplot(212)
plot(gaopints)