数字信号处理_2012_11_3.doc

1、实验二%编程练习：（1 ） （2）nux21nux2（3 ） （4）n（5） （6）12nux uxn12分别对应 x1,x2,x3,x4,x5,x6.时间：2012-19-18%在 NOTEBOOK 中的程序代码为：n=-1:14;x1=(1./2).n.*uDT(n);subplot(3,2,1);stem(n,x1,fill);xlabel(n);title(x1);axis(-1,10,-0.1,1.1); grid on;x2=2.n.*uDT(n);subplot(3,2,2);stem(n,x2,fill);xlabel(n);title(x2);axis(-1,10,-0.1,

2、500);grid on;x3=(-1./2).n.*uDT(n);subplot(3,2,3);stem(n,x3,fill);xlabel(n);title(x3);axis(-1,10,-1.1,1.1) ;grid on;x4=(-2).n.*uDT(n);subplot(3,2,4);stem(n,x4,fill),xlabel(n);title(x4);axis(-1,10,-100,100);grid on;x5=2.(n-1).*uDT(n-1);subplot(3,2,5); stem(n,x5,fill);xlabel(n);title(x5);axis(-1,10,-0.

3、1,100);grid on;x6=(1./2).(n-1).*uDT(n);subplot(3,2,6);stem(n,x6,fill);xlabel(n);title(x6);axis(-1,10,-0.1,1.1) grid on ; 实验二各序列波形图0 2 4 6 8 1000.51nx10 2 4 6 8 100200400nx20 2 4 6 8 10-101nx30 2 4 6 8 10-1000100nx40 2 4 6 8 10050100nx50 2 4 6 8 1000.51nx6实验四b=2,-1.6,-0.9;a=1,-2.5,1.96,-0.48;subplot(

4、121)zplane(b,a)legend(零点, 极点)title(极点与单位圆的位置关系)grid onsubplot(122)impz(b,a,20)title(单位取样响应)grid on -1 -0.5 0 0.5 1-1-0.500.51Real PartImaginaryPart与与与与与与与与与与与0 5 10 1501234567n (samples)Amplitude与与与与与与与与与与b1=1,-1;a1=1,-0.9,-0.65,0.873,0;subplot(121)zplane(b1,a1)legend(零点, 极点)title(零极点和单位圆的位置关系)subpl

5、ot(122)impz(b1,a1,30)title(单位取样响应)grid on -1 -0.5 0 0.5 1-1-0.500.512Real PartImaginaryPart与与与与与与与与与与与与0 5 10 15 20 25-0.8-0.6-0.4-0.200.20.40.60.81n (samples)Amplitude与与与与与与与与与与b=1,0,0;a=1,-3./4,1./8;H,w=freqz(b,a,400,whole);Hm=abs(H);Hp=angle(H);subplot(211)plot(w,Hm)grid onxlabel(omega(rad/s)ylab

6、el(Magnitude)title(离散系统幅频特性曲线)subplot(212)plot(w,Hp)grid onxlabel(omga(rad/s)ylabel(Phase)title(离散系统相频特性曲线) 0 1 2 3 4 5 6 70123(rad/s)Magnitude与与与与与与与与与与0 1 2 3 4 5 6 7-101mga(rad/s)Phase与与与与与与与与与与实验五编程练习中的程序代码是：%1、模拟信号 ，以 进行抽样，求 N 分别)8cos(5)4sin(*2)(tttxnt01.1N为 (1) N=45 (2) N=50 (3) N=55 (4) N=60

7、点 DFT 的幅值谱。%t=0:0.01:1; x=2*sin(4*pi*t)+5*cos(8*pi*t); Y1=fft(x,45); f1=100*(0:44)/45;subplot(4,1,1)plot(f1,Y1(1:45);title(N=45)Y2=fft(x,50); f2=100*(0:49)/50;subplot(4,1,2)plot(f2,Y2(1:50);title( N=50)Y3=fft(x,55); f3=100*(0:54)/55;subplot(4,1,3)plot(f3,Y3(1:55);title(N=55 )Y4=fft(x,60); f4=100*(0:

8、59)/60; subplot(4,1,4)plot(f4,Y4(1:60);title( N=60) 其波形如下：0 10 20 30 40 50 60 70 80 90 100-1000100 N=450 10 20 30 40 50 60 70 80 90 100-2000200 N=500 10 20 30 40 50 60 70 80 90 100-2000200 N=55 0 10 20 30 40 50 60 70 80 90 100-1000100 N=60%2、模拟信号 ，以 进行抽样， N=64，)8cos(5)4sin(*2)(tttxnt01.1N并在信号中加入噪声 w

9、(t),比较有无噪声时的信号谱。%t=0:0.01:1; x1=2*sin(4*pi*t)+5*cos(8*pi*t); Y1=fft(x1,64); f1=100*(0:63)/64; subplot(2,1,1)plot(f1,Y1(1:64);title(无噪声时的信号谱)x2=2*sin(4*pi*t)+5*cos(8*pi*t)+rand(size(t); Y2=fft(x2,64); f2=100*(0:63)/64; subplot(2,1,2)plot(f2,Y2(1:64);title(加入噪声时的信号谱) 其波形如下：0 10 20 30 40 50 60 70 80 90 100-50050 与 与 与 与 与 与 与 与0 10 20 30 40 50 60 70 80 90 100-50050100 与 与 与 与 与 与 与 与 与