实验一信号数字信号处理实验报告
七七 程序及图像程序及图像 1.1.分析理想采样信号序列的特性。分析理想采样信号序列的特性。 1.1.理想采样信号序列理想采样信号序列>> n=0:50;A=444.128;a=50*sqrt(2)*pi; >> w0=50*sqrt(2)*pi;T=1/1000; >> x=A*exp(-a*n*T).*sin(w0*n*T); >> k=-25:25; W=(pi/12.5)*k; f=(1/25)*k*1000; >> X=x*(exp(-j*pi/12.5).(n'*k); >> magX=abs(X);angX=angle(X); >> subplot(3,1,1);stem(x);title('x 的时域') >> subplot(3,1,2);stem(f,magX);title('x 的幅度谱') >> subplot(3,1,3);stem(f,angX);title('x 的相位谱')>> T=1/300; >> x=A*exp(-a*n*T).*sin(w0*n*T); >> magX=abs(X);angX=angle(X); >> subplot(3,1,1);stem(x);title('x 的时域') >> subplot(3,1,2);stem(f,magX);title('x 的幅度谱') >> subplot(3,1,3);stem(f,angX);title('x 的相位谱')>> T=1/200; >> x=A*exp(-a*n*T).*sin(w0*n*T); >> X=x*(exp(-j*pi/12.5).(n'*k); >> magX=abs(X);angX=angle(X); >> subplot(3,1,1);stem(x);title('x 的时域') >> subplot(3,1,2);stem(f,magX);title('x 的幅度谱') >> subplot(3,1,3);stem(f,angX);title('x 的相位谱')没有出现“混淆”现象,采样频率均大于信号频率的二倍。2.2.离散信号,系统和系统响应的分析离散信号,系统和系统响应的分析(1)>> n=0:50;xb=1 zeros(1,50);k=-25:25; >> Xb=xb*(exp(-j*pi/25).(n'*k); >> magX=abs(Xb);angX=angle(Xb); >> subplot(3,1,1);stem(xb);title('xb 的时域') >> subplot(3,1,2);stem(magX);title('xb 的幅度谱') >> subplot(3,1,3);stem(angX);title('xb 的相位谱')>> n=1:50; hb=zeros(1,50); >> hb(1)=1;hb(2)=2.5;hb(3)=2.5;hb(4)=1;k=-25:25; >> X=hb*(exp(-j*pi/25).(n'*k); >> magX=abs(X);angX=angle(X); >> subplot(3,1,1);stem(hb);title('hb 的时域') >> subplot(3,1,2);stem(magX);title('hb 的幅度谱') >> subplot(3,1,3);stem(angX);title('hb 的相位谱')>> y=conv(xb,hb); >> subplot(1,1,1);stem(y);title('y 的时域')>> n=0:99; >> k=-50:50; >> X=y*(exp(-j*pi/25).(n'*k); >> magX=abs(X); >> angX=angle(X); >> subplot(3,1,2);stem(magX);title('y 的幅度谱') subplot(3,1,3);stem(angX);title('y 的相位谱')(2)>> n=0:50;xc=ones(1,10) zeros(1,41);>> subplot(3,1,1);stem(n,xc);title('xc 的时域')>> axis(0 50 0 1.2); >> k=-25:25; >> X=xc*(exp(-j*pi/25).(n'*k); >> magX=abs(X);angX=angle(X); >> subplot(3,1,2);stem(magX);title('xc 的幅度谱') >>subplot(3,1,3);stem(angX);title('xc 的相位谱')>> ha=ones(1,10) zeros(1,41); >> y=conv(xc,ha);subplot(3,1,1); >>stem(y);title('线性卷积 y 的时域')>> n=0:100;k=-50:50; >> X=y*(exp(-j*pi/25).(n'*k); >> magX=abs(X);angX=angle(X); >> subplot(3,1,2);stem(magX);title('y 的幅度谱') >>subplot(3,1,3);stem(angX);title('y 的相位谱')(3)(3)>> n=0:50;A=1;a=0.4;w=2.0734;T=1; >> x=A*exp(-a*n*T).*sin(w*n*T); >> k=-25:25;W=(pi/12.5)*k; >> f=(1/25)*k*1000; >> X=x*(exp(-j*pi/12.5).(n'*k); >> magX=abs(X);angX=angle(X); >> subplot(3,1,1);stem(x);title('x 的时域') >> subplot(3,1,2);stem(f,magX);title('x 的幅度谱') >> subplot(3,1,3);stem(f,angX);title('x 的相位谱')>> ha=ones(1,10) zeros(1,41); >> y1=conv(x,ha); >> subplot(3,1,1);stem(y1);title('线性卷积 y1 的时域') >>n=0:100;k=-50:50; >>X=y1*(exp(-j*pi/25).(n'*k); >>magX=abs(X);angX=angle(X); >>subplot(3,1,2);stem(magX);title('y1 的幅度谱') >>subplot(3,1,3);stem(angX);title('y1 的相位谱')>> n=0:50;A=1;a=0.1;w=2.0734;T=1; >>x=A*exp(-a*n*T).*sin(w*n*T); >> ha=ones(1,10) zeros(1,41); >>y2=conv(x,ha); >>subplot(3,1,1);stem(y2);title('线性卷积 y2 的时域') >>n=0:100;k=-50:50;X=y2*(exp(-j*pi/25).(n'*k); >>magX=abs(X);angX=angle(X); >>subplot(3,1,2);stem(magX);title('y2 的幅度谱') >>subplot(3,1,3);stem(angX);title('y2 的相位谱')>> n=0:50;A=1;a=0.1;w=1.2516;T=1; >>x=A*exp(-a*n*T).*sin(w*n*T); >> ha=ones(1,10) zeros(1,41); >>y3=conv(x,ha); >>subplot(3,1,1);stem(y3);title('线性卷积 y3 的时域') >>n=0:100;k=-50:50;X=y3*(exp(-j*pi/25).(n'*k); >>magX=abs(X);angX=angle(X); >>subplot(3,1,2);stem(magX);title('y3 的幅度谱') >>subplot(3,1,3);stem(angX);title('y3 的相位谱')3>> n=0:50;>> A=444.128; >> a=50*sqrt(2)*pi; >> w0=50*sqrt(2)*pi; >> T=1/1000; >> x=A*exp(-a*n*T).*sin(w0*n*T); >> ha=ones(1,10) zeros(1,41);>> k=-25:25;X=x*(exp(-j*pi/12.5).(n'*k); >> magX=abs(X); >> subplot(3,2,1);stem(magX);title('输入信号的幅度谱');>> angX=angle(X); >> subplot(3,2,2);stem(angX);title('输入信号的相位谱')>> Hb=hb*(exp(-j*pi/12.5).(n'*k); >> magHb=abs(Hb); >> subplot(3,2,3);stem(magHb);title('系统响应的幅度谱');>> angHb=angle(Hb); >> subplot(3,2,4);stem(angHb);title('系统响应的相位谱')>> y=conv(x,ha);>> n=1:101;k=1:99; >> Y=y*(exp(-j*pi/12.5).(n'*k); >> magY=abs(Y); >> subplot(3,2,5);stem(magY);title('输出信号的幅度谱');>> angY=angle(Y); >> subplot(3,2,6);stem(angY);title('输出信号的相位谱')