电气系数字信号实验程序和图

1、 1.5 1 0.5 0 0 50 100 150 200 250 300 350 400 450 500 1.5 1 0.5 0 100 110 120 130 140 150 160 170 180 190 200 实验：6-1.1 window=chebwin(35,30; b=fir1(34,0.48,'high',window; freqz(b,1,512; 50 Magnitude (dB 0 -50 -100 -150 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 Normalized Frequency (´p rad/sampl

2、e 0.9 1 1000 Phase (degrees 0 -1000 -2000 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 Normalized Frequency (´p rad/sample 0.9 1 实验6-1.2: f=0 0.6 0.6 1; m=1 1 0 0; b=fir2(30,f,m; h,w=freqz(b,1,128; plot(f,m,w/pi,abs(h; 1.4 1.2 1 0.8 0.6 0.4 0.2 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 实验6-2： wn=0.4 0.7;

3、n=38; b1=fir1(n,wn,'stop' win=chebwin(n+1,30; b2=fir1(n,wn,'stop',win; h1,w1=freqz(b1,1,512,2; h2,w2=freqz(b2,1,512,2; subplot(2,1,1; plot(w1,20*log(abs(h1; subplot(2,1,2; plot(w2,20*log(abs(h2; 50 0 -50 -100 -150 -200 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 100 0 -100 -200 -300 0 0.

4、1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 实验 6-3： f=0 0.2 0.2 0.4 0.4 0.6 0.6 1; m=1 1 0 0 1 1 0 0; b=fir2(40,f,m; h,w=freqz(b,1,256; plot(f,m,':',w/pi,abs(h; legend('hope','real' 50 0 -50 -100 -150 -200 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.5 hope real 1 0.5 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 实验6-4： fs=8000; f=1500 2000; a=1 0; dev=0.01 0.1; n,f0,a0,w=remezord(f,a,dev,fs; b=remez(n,f0,a0,w; h,w1=freqz(b,1,1024,fs; plot(w1,20*log(abs(h; 20 0 -20 -40 -60 -80 -100 -120 -140 -160 0 500 1000 1500 2000 2500 3000 3500 4000 【注】 ：感谢龙