带波纹和阻带衰减,单位为dB。函数返回值n为模拟滤波器的最小阶数;wc为模拟滤波器的截止频率,单位为rad/s)。因本设计中提供的滤波器参数较少,故忽略了滤波器最小阶数的选择,在结果分析中再对典型阶数进行讨论。Р 其次是模拟滤波器性能测试的问题。可从频域、时域两方面分析滤波器特性,其中频域分析可用滤波器的传递函数描述,绘出其幅度和相位曲线,可采用MATLAB信号处理工具箱提供的函数abs给出系统的幅频曲线和angle给出系统的相频曲线;时域分析可用脉冲(冲激)响应来描述,绘出其脉冲响应和阶跃响应曲线,可分别采用MATLAB信号处理工具箱提供的函数[[y,t]=]impulse(H)给出该系统的模拟脉冲响应和[[y,t]=]step(H)来得到该系统的阶跃响应。Р 6Р 武汉理工大学《数字信号处理》课程设计说明书Р 3程序设计Р 3.1程序编写Р 结合理论分析阶段的研究,得到该模拟高通滤波器的MATLAB程序如下: %Butterworth 模拟低通滤波器原型设计Р n=0:0.01:2; %设定频率点Р for ii=1:4 %定义循环,产生不同阶数的曲线Р switch iiР case 1,N=2;Р case 2,N=5;Р case 3,N=10;Р case 4,N=30;Р endР [z,p,k]=buttap(N); %调用Butterworth模拟低通滤波器原型函数Р [b,a]=zp2tf(z,p,k); %将零点极点增益形式转换为传递函数形式Р [H,w]=freqs(b,a,n); %按n指定的频率点给出频率响应Р magH2=(abs(H)).^2;Р subplot(211);Р hold on;Р plot(w,magH2);Р endР xlabel('w/wc');Р ylabel('|H(jw)|^2');