数字信号处理_2012_11_3.doc

实验二%编程练习：（1） （2） （3） （4） （5） （6）分别对应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,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.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 ; 实验二各序列波形图实验四b=2,-1.6,-0.9;a=1,-2.5,1.96,-0.48;subplot(121)zplane(b,a)legend(零点,极点)title(极点与单位圆的位置关系)grid onsubplot(122)impz(b,a,20)title(单位取样响应)grid on b1=1,-1;a1=1,-0.9,-0.65,0.873,0;subplot(121)zplane(b1,a1)legend(零点,极点)title(零极点和单位圆的位置关系)subplot(122)impz(b1,a1,30)title(单位取样响应)grid on 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)ylabel(Magnitude)title(离散系统幅频特性曲线)subplot(212)plot(w,Hp)grid onxlabel(omga(rad/s)ylabel(Phase)title(离散系统相频特性曲线) 实验五编程练习中的程序代码是：%1、模拟信号，以 进行抽样，求 N分别为 (1) N=45 (2) N=50 (3) N=55 (4) N=60点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:59)/60; subplot(4,1,4)plot(f4,Y4(1:60);title( N=60) 其波形如下：%2、模拟信号，以 进行抽样， N=64，并在信号中加入噪声w(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