对于固定位置的干扰(比如带内的单音干扰),可以采用当前的值减去统计出来的平均值。
clear
close all
a=randn(10^4,1024);
a_rmv_dc=zeros(10^4,1024);
for cntx=1:10^4
for cnty=1:1024
if(mod(cnty,128) == 0)
a(cntx,cnty)=10+randn;
end
end
end
a_s0=sum(a);
a_avg=a_s0./10^4;
for cntx=1:10^4
for cnty=1:1024
a_rmv_dc(cntx,cnty)=a(cntx,cnty)-a_avg(cnty);
end
end
plot(a_rmv_dc(1,1:1024),'r*')
hold on
plot(a(1,1:1024),'b')