今天上数字信号处理的实验。。。做滤波器仿真,放上来玩玩:
- %模拟巴特沃斯低通滤波器设计
- clear;close all
- fp=500;fs=750;Ap=3;As=15; %输入滤波器指标
- [N,fc]=buttord(fp,fs,Ap,As,'s') %计算阶数N和3dB截止频率fc
- [B,A]=butter(N,fc,'s'); %设计低通巴特沃斯模拟滤波器
- [hf,f]=freqs(B,A,1024); %计算模拟滤波器频率响应,freqs为工具函数
- subplot(1,3,1);
- plot(f,20*log10(abs(hf)/abs(hf(1))))
- title('模拟巴特沃斯滤波器');
- grid;
- xlabel('f/Hz');
- ylabel('幅度(dB)');