波器时只需要代入相应的实现函数即可。下面我们将给出 FIR 数字滤波器的主要程序。通信系统综合设计与实践第 9页%========================= F IR低通滤波器======================= Ft=2000; % 采样频率 fpts=[100 120]; % 通带边界频率 fp=100Hz ,阻带截止频率 fs=120Hz mag=[1 0]; dev=[0.01 0.05]; % 通带波动 1% ,阻带波动 5% [n21,wn21,beta,ftype]=kaiserord(fpts,mag,dev,Ft);%kaiserord 估计采用凯塞窗设计的 FIR 滤波器的参数 b21=fir1(n21,wn21,Kaiser(n21+1,beta)); %由 fir1 设计滤波器[h,w]=freqz(b21,1); % 得到频率响应 plot(w/pi,abs(h)); grid on title('FIR 低通滤波器'); 0 0. 1 0. 2 0. 3 0. 4 0. 5 0. 6 0. 7 0. 8 0. 91 0 0. 2 0. 4 0. 6 0. 8 1 1. 2 1. 4 F IR 低通滤波器如图 6.6 FIR 低通滤波器三 AM解调 t=-1:0.00001:1; A0=10; % 载波信号振幅 A1=5; % 调制信号振幅 A2=3; % 已调信号振幅 f=3000; % 载波信号频率 w0=2*f*pi; m=0.15; % 调制度 k=0.5 ; %DSB 前面的系数 mes=A1*cos(0.001*w0*t); % 调制信号 Uam=A2*(1+m*mes).*cos((w0).*t); %AM 已调信号 Dam=Uam.*cos(w0*t); %对 AM 调制信号进行解调 subplot(2,1,1);
全文预览
