23t.m

包含50个数字信号处理方面的matlab程序

clc;
close all;
clear all;
Rp=3;
Rs=18;
fs=2000;
ws1=0.2*pi;
wp1=0.3*pi;
wp2=0.4*pi;
ws2=0.5*pi;
wp3=(2*fs)*tan(wp1/2);
wp4=(2*fs)*tan(wp2/2);
ws3=(2*fs)*tan(ws1/2);
ws4=(2*fs)*tan(ws2/2);
wp=[wp3,wp4];
ws=[ws3,ws4];
[n,wn]=buttord(wp,ws,Rp,Rs,'s');
[b,a]=butter(n,wn,'bandpass','s');
[bz,az]=bilinear(b,a,fs);
[db,mag,pha,grd,w]=freqz_m(bz,az);
plot((w*2000)/(2*pi),db);
grid;