北邮DSP第一次实验作业.docVIP

实验一： 单位冲击序列 程序：n1=input(n1=); n2=input(n2=); n=n1:n2; k=0; nt=length(n);nk=abs(k-n1)+1;x=zeros(1,nt);x(nk)=1;subplot(2,1,1);stem(n,x,filled );title(delta(n));m1=input(m1=);m2=input(m2=)m=m1:m2; k=input(k=);mt=length(m); mk=abs(k-m1)+1; y=zeros(1,mt); y(mk)=1; subplot(2,1,2);stem(m,y,filled ); title(delta(m-k)); 单位阶跃响应 程序：n1=input(n1=); n2=input(n2=);n=n1:n2; k=0; nt=length(n); nk=abs(k-n1)+1;x=[zeros(1,nk-1),ones(1,nt-nk+1)]; stem(n,x,filled ); 正弦序列 程序：nt=input(nt=) ;N=input(N=);n=0:nt*N-1;f=input(f=); T=1/f;dt=T/N;A=input(A=);fai=input(fai=);tn=n*dt;x=A*sin(2*pi*f*tn+fai);stem(tn,x,fille 复指数序列 程序：n=0:N-1;N=input(N=);w=input(w=);x=exp(j*w*n); subplot(3,1,1);stem(n,x,filled);title(x=exp(jwn));subplot(3,1,2);stem(n,real(x),filled);title(real(x));subplot(3,1,3);stem(n,imag(x),filled);title(imag(x)); 指数序列 程序：N=input(N=);n=1:N-1;a=input(a=);x=a.^n; stem(n,x,filled );title(exponent); 复指数序列的实部是余弦序列，虚部是正弦序列 实验二、离散系统的冲击响应与阶跃响应 1、y[n]+ 0.75y[n ?1]+ 0.125y[n ? 2] = x[n]? x[n ?1] 程序：a=[1,0.75,0.125];b=[1,-1,0];N=20;n=0:N-1;hn=impz(b,a,n); gn=dstep(b,a,n)subplot(1,2,1);stem(n,hn,filled);title(h(n));subplot(1,2,2);stem(n,gn,filled);title(g(n)); 左边为冲击响应 右边阶跃响应 2、y[n] = 0.25{x[n ?1]+ x[n ? 2]+ x[n ?3]+ x[n ? 4]} 程序：a=[1,0.75,0.125];b=[1,-1,0];N=20;n=0:N-1;hn=impz(b,a,n); gn=dstep(b,a,n)subplot(1,2,1);stem(n,hn,filled);title(h(n));subplot(1,2,2);stem(n,gn,filled);title(g(n)); 左边为冲击响应 右边阶跃响应 实验三：求系统零极点及幅频响应 程序：b=input(b=);a=input(a=);rz=roots(b);rp=roots(a); subplot(2,1,2);zplane(b,a);n=(0:500)*pi/500;[h,w]=freqz(b,a,n);subplot(2,1,1),plot(n/pi,abs(h));title(Magnitude);