信号、系统及系统响应实验。

figure(1) xb=unitimp(n); figure(2) hb=tedingxulie(n); y=conv(xb,hb); k1=1:99; N=length(y); n1=1:N;

[magX,angX]=fpres(y,n1,k1); figure(3) subplot(3,1,1); stem(y);

title('系统响应'); xlabel('时间'); ylabel('幅值');

subplot(3,1,2); stem(magX)

title('系统响应幅度谱') xlabel('时间'); ylabel('幅值');

subplot(3,1,3); stem(angX)

title('系统响应相位谱') xlabel('频率'); ylabel('幅值');

%%%%%%%%%卷积定理验证 clc; clear all; n=0:50; A=1; fs=1; T=1/fs; a=0.1; w0=1.2516;

xa=A*exp(-a*n*T).*sin(w0*n*T); 的采样信号

k=-25:25;

X=xa*(exp(-j*pi/12.5)).^(n'*k); magX=abs(X); figure(1) subplot(3,2,1);

stem(magX);title('输入信号xa的幅度谱'); angX=angle(X); subplot(3,2,2);

stem(angX);title('输入信号xa的相位谱'); n1=1:10;

ha=sign(sign(10-n1)+1); k1=-25:25;

Hb=ha*(exp(-j*pi/12.5)).^(n1'*k1); magHb=abs(Hb); subplot(3,2,3);

stem(magHb);title('系统响应ha的幅度谱');

angHb=angle(Hb); subplot(3,2,4);

stem(angX);title('系统响应ha的相位谱'); y=conv(xa,ha); n2=1:60;k2=1:60;

Y=y*(exp(-j*pi/12.5)).^(n2'*k2); magY=abs(Y); subplot(3,2,5);

stem(magY);title('输出信号y的幅度谱'); angY=angle(Y); subplot(3,2,6);

stem(angY);title('输出信号y的相位谱'); XHb=X.*Hb; figure(2) subplot(2,1,1);

stem(magY);title('输出信号y的幅度谱'); subplot(2,1,2);

stem(abs(XHb));title('xa的幅度谱与hb的% 理想幅度谱相乘');