matlab数字带阻滤波器设计以及DSP设计(精)(word文档良心出品).docxVIP

数字滤波器设计是用硬件或者软件实现的一种算法，这个算法是为了达到滤波的 目的而对数字输入信号进行运算产生数字输出信号。数字滤波器这个词是指执行 滤波算法的特定硬件或者软件程序，数字滤波器经常作用的对象是数字化的模拟 信号，或者刚好是存储在计算机存储器里代表某些变量的数据。 下图给出了一个具有模拟输入信号和输出信号的实时数字滤波器的简化框图。带 限模拟信号x(t)被周期地抽样，然后通过量化转化成一数字序列 x(n)(n=0, 1, )。 输人 亠 ADC 1 1 ■ 1 OAC Mt yio ■ 数字处理器依据滤波器的计算算法执行滤波运算，把输入系列 x(n)映射到输出 系列y(n)。DAC把数字滤波后的输出转化成模拟值，这些模拟值接着被模拟滤波 器平滑，并且消去不想要的高频分量? 利用MATLAB的强大运算功能，基于 MATLAB信号处理工具的数字滤波器设计 法可以快速有效的设计由软件组成的常规数字滤波器，设计方便、快捷，极大的 减轻了工作量。在设计过程中可以对比滤波器特性，随时更改参数，以达到滤波 器设计的最优化。 1、MATLB仿真源程序 fs=15000;T=1/fs; rp=1;rs=40; wp1=0.11*pi;wp2=0.81*pi;ws1=0.31*pi;ws2=0.61*pi; % 数字带阻滤波器技术指 标 wc1=(2/T)*tan(wp1/2); % 频率预畸变 wc2=(2/T)*ta n(wp2/2);wr1=(2/T)*ta n(ws1/2);wr2=(2/T)*ta n(ws2/2); w0=sqrt(wc1*wc2);B=wc2-wc1; wp=1; %归一化通带截止频率 ws=wp*(wr1*B)/(w0A2-wr1A2); % 归一化阻带 截止频率 [N,wc]=buttord(wp,ws,rp,rs,s); %求滤波器阶数和3dB截止频率 [Z,P,K]=buttap(N); [Md,Nd]=zp2tf(Z,P,K); %将零极点形式转换为传输形式 [M,N]=lp2bs(Md,Nd,w0,B); %对低通滤波器进行频率变换，转换为带阻滤波器 [h,w]=freqs(M,N); %模拟带阻滤波器的幅频响应 xlabel(频率/Hz);ylabel(幅度 );title(模拟带阻滤波器); [b,a]=bilinear(M,N,15000); % 对模拟滤波器双线性变换 figure(1); freqz(b,a);[H,W]=freqz(b,a); % 绘出频率响应 axis([0,1,-100,20]); figure(2); plot(W*fs/(2*pi),abs(H));grid on; T I T I回 xlabel(频率 /Hz);ylabel(幅值); title(数字滤波器幅频响应 ■ 口-200■JQO-600File Edit View [ntert TooU Desktop Window Help DQtl昌k題窃?富□「因吧 ■ 口 -200 ■JQO -600 File Edit View [ntert TooU Desktop Window Help DQtl昌k題窃?富□「因 吧01 FJ Figure 1 -??」6誓p)岀5希￡ -HH^WUS1于 o i a s Normalized Frequency (xi rwV^ample) |H(ejOmega)| ); 指标参数： sl=0.5062 rad, sl=0.5536 rad, pl=0.3750 rad, pu=0.6750 rad, p=1dB,s=40 dB 由MATLAB可得出系统函数系数： Bz=[ 0.3521 0.1819 1.0416 0.3576 1.0416 0.1819 0.3521] Az=[ 1.0000 0.3531 1.0547 0.3038 0.7232 0.0645 0.0095] 系统函数为： H(Z)= (0.3521+0.1819 *z-1+1.0416*z-2+ 0.3576*z-3+ 1.0416*z-4+ 0.1819*z +0.3521*z”(1.0000+ 0.3531*z+ 1.0547*z+ 0.3038*z+ -4-5-60.7232*z+ 0.0645*z+ 0.0095*z) DSP带阻源程序： #i ncludemath.h #defi ne IIRNUMBER 30 #defi ne SIGNAL1F 2200 #defi ne SAMPLEF 8000 #define PI 3.1415926 -5-6-1-2-3 float In putWave(); float IIR(); float