test.m

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

%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(311)
plot(pcsym10)
subplot(312)
plot(gaopints)
subplot(313)
plot(gaopints1)