泰克频谱仪FSK解调后滤波处理

1、%解调+中值滤波

%程序中红色标注是需要根据具体情况修改的参数

PhaseDiff=diff(unwrap(angle(ValueData)));

DouPhaDiff=double(PhaseDiff);

medFiltData=medfilt1(DouPhaDiff,10);

%整形成非标准方波

save=zeros(1,31800);

for i=10:32000

for index=1:19

posNum=0;

negNum=0;

if(medFiltData(i-10+index)>=0)

posNum=posNum+1;

else negNum=negNum+1;

end

if(posNum>negNum)

save(i)=0.2;

else save(i)=-0.2;

end

end

end

%方波整形

save1=repmat(save,1,1);%复制矩阵

for k=2:32000

if(save(k-1)==-0.2 && save(k)==0.2)

para=0;

for l=1:10

if(save(k+l)==0.2)

para=para+1;

end

end

if(para==10)

save1(k)=0.2;

else

for l=0:10

save1(k+l)=-0.2;

end

end

elseif (save(k-1)==0.2 && save(k)==-0.2)

para=0;

for l=1:10

if(save(k+l)==-0.2)

para=para+1;

end

end

if(para==10)

save1(k)=-0.2;

else

for l=0:10

save1(k+l)=0.2;

end

end

end

end

medData=medfilt1(save1,10);%中值滤波