数字信号实验二.docVIP

. . 自定义一个长度为8点的信号，信号幅度值也由自己任意指定，对该信号作DTFT，分别画出幅度谱和相位谱. 解：代码如下： n=0:7; x=(0.9*exp(j*pi/3)).^n; w=-pi:pi/200:pi; X= x* exp(-j).^(n*w); magX=abs(X); angX=angle(X); subplot(2,1,1); plot(w/pi,magX); xlabel(w/pi); ylabel(幅度|X|); title(幅频特性); subplot(2,1,2); plot(w/pi,angX); xlabel(w/pi); ylabel(相位(rad/π)); title(相频特性); 实现P138的题2.27,2.28,2.30,2.34 2.27 解：归一化因子为A=1/(1+2N)； （1）当N = 5时： clc;clear all N = 5;M = 2*N+1;n = -20:20; x = [zeros(1,15),ones(1,M),zeros(1,15)]; omega = -pi:2*pi/1500:pi; X = sin(0.5*M*omega)./sin(0.5*omega)/M; subplot(211); plot(omega/pi,X); xlabel(\omega/\pi); title(频谱); ylabel(X(e^j^\omega));axis([-1,1,min(X),1]);grid; （2）当N = 15时： clc;clear all N = 15;M = 2*N+1;n = -20:20; x = [zeros(1,15),ones(1,M),zeros(1,15)]; omega = -pi:2*pi/1500:pi; X = sin(0.5*M*omega)./sin(0.5*omega)/M; subplot(211); plot(omega/pi,X); xlabel(\omega/\pi); title(频谱); ylabel(X(e^j^\omega));axis([-1,1,min(X),1]);grid; （3）当N = 25时： clc;clear all N = 25;M = 2*N+1;n = -20:20; x = [zeros(1,15),ones(1,M),zeros(1,15)]; omega = -pi:2*pi/1500:pi; X = sin(0.5*M*omega)./sin(0.5*omega)/M; subplot(211); plot(omega/pi,X); xlabel(\omega/\pi); title(频谱); ylabel(X(e^j^\omega));axis([-1,1,min(X),1]);grid; （4）当N = 100时： clc;clear all N = 100;M = 2*N+1;n = -20:20; x = [zeros(1,15),ones(1,M),zeros(1,15)]; omega = -pi:2*pi/1500:pi; X = sin(0.5*M*omega)./sin(0.5*omega)/M; subplot(211); plot(omega/pi,X); xlabel(\omega/\pi); title(频谱); ylabel(X(e^j^\omega));axis([-1,1,min(X),1]);grid; 2.28 （1）h(n)=0.9^| n |u(n) 解：代码如下： clc;clear all b = 1; a = [1,-0.9]; [H,w] = freqz(b,a,1024,whole); magX = abs(H);angX = angle(H); subplot(2,1,1); plot(w/pi,magX);grid; xlabel(w/pi); ylabel(幅值); title(幅频特性); subplot(2,1,2); plot(w/pi,angX);grid; xlabel(以\pi为单位的频率);title(相频特性);ylabel(弧度)； (2) h(n)=[0.3^n+0.5^n]u(n) 解：代码如下： clc;clear all b = [2,-0.8]; a = [1,-0.8,0.15]; [H,w] = freqz(b,a,1024,whole); magX = abs(H);angX = angle(H); subplot(2,1,1); plot(w/pi,magX);grid; xlabel(w/pi); ylabel(幅值); title(幅频特性); sub