第十三章 数字信号设计.docVIP

第一章 已知LTI离散系统，x(n)=[1 1 1],h(n)=[0 1 2 3],求y(n)= x(n)* h(n) x=[1,1,1]; h=[0,1,2,3]; y=conv(x,h); subplot(2,2,1);stem([0:length(x)-1],x); ylabel(x(n));xlabel(Time index n); subplot(2,2,2);stem([0:length(h)-1],h); ylabel(h(n));xlabel(Time index n); subplot(2,1,2);stem([0:length(y)-1],y); ylabel(y(n)=x(n)*h(n));xlabel(Time index n); 第二章 用fourier函数,理论上求下列连续时间信号的频谱。 (1). 三角脉冲信号x1(t)= t=-2:0.1:2; x=tripuls(t,2,0); plot(t,x); syms t w xt=sym((t+1)*Heaviside(t+1)-2*t*Heaviside(t)+(t-1)*Heaviside(t-1)); Fw=fourier(xt,t,w); FFw=maple(convert,Fw,piecewise); FFp=abs(FFw); ezplot(FFp,[-10*pi 10*pi]) axis([-10*pi 10*pi 0 1]) (2). 单边指数信号x2(t)=eu(t) N=256; fs=6; Ts=1/fs; t=(0:N-1)*Ts; x=exp(-t); y=fft(x); mag=Ts*abs(y); Ws=2*pi*fs; w=(0:length(y)-1)*Ws/length(y); X=1./sqrt(w.^2+1); w1=w(1:length(y)/2); plot(w1,mag(1:length(y)/2)); xlable(频率 （弧度/秒）); ylable(幅度谱); 第四章 分别使用矩形窗、哈明窗、汉宁窗设计一个N=10的FIR低通和高通滤波器， 截频为。 1) 作出各滤波器的单位脉冲响应 2) 作出各滤波器的幅频响应并比较各滤波器的通带纹波和阻带纹波。 若当输入为，计算各滤波器的输出并作出响应波形。 N=10;M=N-1;wc=pi/3;n=0:M;N1=2048;n1=0:N1-1; b1=fir1(M,wc/pi,boxcar(N)); [H1,w]=freqz(b1,wc,512); H1_db=20*log10(abs(H1)); h1=ifft(H1,N1); b2=fir1(M,wc/pi,hamming(N)); [H2,w]=freqz(b2,wc,512); H2_db=20*log10(abs(H2)); h2=ifft(H2,N1); b3=fir1(M,wc/pi,hanning(N)); [H3,w]=freqz(b3,wc,512); H3_db=20*log10(abs(H3)); h3=ifft(H3,N1); figure(1) c=plot(w,H1_db,w,H2_db,r,w,H3_db,g); figure(2) subplot(3,1,1);stem(n1,real(h1));axis([0 25 -0.1 0.2]); subplot(3,1,2);stem(n1,real(h2),r);axis([0 25 -0.1 0.2]); subplot(3,1,3);stem(n1,real(h3),g);axis([0 25 -0.1 0.2]); b4=fir1(M,wc/pi,high,boxcar(N+1)); [H4,w]=freqz(b4,wc,512); H4_db=20*log10(abs(H4)); h4=ifft(H4,N1); b5=fir1(M,wc/pi,high,hamming(N+1)); [H5,w]=freqz(b5,wc,512); H5_db=20*log10(abs(H5)); h5=ifft(H5,N1); b6=fir1(M,wc/pi,high,hanning(N+1)); [H6,w]=freqz(b6,wc,512); H6_db=20*log10(abs(H6)); h6=ifft(H6,N1); figure(3) c=plot(w,H4_db,w,H5_db,r,w,H6_db,g); figure(4) subplot(3,1,1);stem(n1,real(h4));axis([0 25 -0.2 0.2]); subplot(3,1,2);stem(n1,real