分析理想采樣信號序列的特性

分析理想采樣信號序列的特性T=1/1000n=0:50;A=444.128;a=50*sqrt(2.0)*pi;T=(1/1000);w0=50*sqrt(2.0)*pi;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);subplot(3,1,1);stem(n,x);title('理想采樣信號序列')magX=abs(X);subplot(3,1,2);stem(f,magX);title('理想采樣信號序列的幅度譜')angX=angle(X);subplot(3,1,3);stem(f,angX);title('理想采樣信號序列的相位譜')T=1/300n=0:50;A=444.128;a=50*sqrt(2.0)*pi;T=(1/300);w0=50*sqrt(2.0)*pi;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);subplot(3,1,1);stem(n,x);title('理想采樣信號序列')magX=abs(X);subplot(3,1,2);stem(f,magX);title('理想采樣信號序列的幅度譜')angX=angle(X);subplot(3,1,3);stem(f,angX);title('理想采樣信號序列的相位譜')（3)n=0:50;A=444.128;a=50*sqrt(2.0)*pi;T=(1/200);w0=50*sqrt(2.0)*pi;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);subplot(3,1,1);stem(n,x);title('理想采樣信號序列')magX=abs(X);subplot(3,1,2);stem(f,magX);title('理想采樣信號序列的幅度譜')angX=angle(X);subplot(3,1,3);stem(f,angX);title('理想采樣信號序列的相位譜')離散信號、系統(tǒng)和系統(tǒng)響應(yīng)的分析n=0:50;x=[1zeros(1,50)];subplot(3,1,1);stem(n,x);title('單位脈沖信號序列')k=-25:25;X=x*(exp(-j*pi/12.5)).^(n'*k);magX=abs(X);subplot(3,1,2);stem(magX);title('單位脈沖信號的幅度譜')angX=angle(X);subplot(3,1,3);stem(angX);title('單位脈沖信號的相位譜')特定沖激串hb(n)N=1:50;x=zeros(1,50);x(1)=1;x(2)=2.5;x(3)=2.5;x(4)=1;closeall;subplot(3,1,1);stem(x);title('特定沖擊串')k=-25:25;X=x*(exp(-j*pi/12.5)).^(n'*k);magX=abs(X);subplot(3,1,2);stem(magX);title('特定沖擊串的幅度譜')angX=angle(X);subplot(3,1,3);stem(angX);title('特定沖擊串的相位譜')系統(tǒng)響應(yīng)n=1:50;hb=zeros(1,50);hb(1)=1;hb(2)=2.5;hb(3)=2.5;hb(4)=1;closeall;subplot(3,1,1);stem(hb);title('系統(tǒng)hb[n]')n=0:50;x=[1zeros(1,50)];subplot(3,1,2);stem(x);title('輸入信號x[n]')y=conv(x,hb);subplot(3,1,3);stem(y);title('輸出信號y[n]')矩形序列xc(n)n=0:50;x=[ones(1,10),zeros(1,41)];subplot(3,1,1);stem(n,x);title('矩形序列')axis([05001.2]);k=-25:25;X=x*(exp(-j*pi/12.5)).^(n'*k);magX=abs(X);subplot(3,1,2);stem(magX);title('矩形序列的幅度譜')angX=angle(X);subplot(3,1,3);stem(angX);title('矩形序列的相位譜')ha(n)n=0:50;x=[ones(1,10),zeros(1,41)];subplot(3,1,1);stem(n,x);title('ha(n)序列')axis([05001.2]);k=-25:25;X=x*(exp(-j*pi/12.5)).^(n'*k);magX=abs(X);subplot(3,1,2);stem(magX);title('ha(n)序列的幅度譜')angX=angle(X);subplot(3,1,3);stem(angX);title('ha(n)序列的相位譜')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('系統(tǒng)響應(yīng)的幅度譜')；angHb=angle(Hb);subplot(3,2,4);stem(angHb);title('系統(tǒng)響應(yīng)的相位譜')；n=1:99;k=1:99;Y=y*(exp