1、七七 程序及图像程序及图像 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,
2、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(
3、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)
4、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.0
5、734;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;
6、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=ang
7、le(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(输出信号的相位谱)
《实验一信号数字信号处理实验报告》由会员j****9分享,可在线阅读,更多相关《实验一信号数字信号处理实验报告》请在金锄头文库上搜索。