lpf.m

你也matlab实现抑制载波单边带调幅的仿真

function [t,st]=lpf(f,sf,B)
% 本函数为通带为B的低通滤波器
% 入口参数：f:频率样本向量
%           sf:频谱向量
%           B:低通滤波器带宽
% 出口参数：t:时间向量
%           st:输出信号向量
df=f(2)-f(1);
T=1/df;
hf=zeros(1,length(f));
bf=[-floor(B/df):floor(B/df)]+floor(length(f)/2);
hf(bf)=1;
yf=hf.*sf;
[t,st]=F2T(f,yf);
st=real(st);