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

magX=abs(X); angX=angle(X);

%%%%%%理想的采样信号 clc; clear all; n=0:50; A=444.128; fs=200; T=1/fs;

a=50*sqrt(2.0)*pi; w0=50*sqrt(2.0)*pi;

x=A*exp(-a*n*T).*sin(w0*n*T); k=-25:25;

[magX,angX]=fpres(x,n,k); figure(1) subplot(3,1,1); stem(x);

title('理想采样信号序列（fs=1000hz）');xlabel('时间'); ylabel('幅值');

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

title('理想采样信号序列幅度谱') xlabel('时间'); ylabel('幅值');

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

title('理想采样信号序列相位谱') xlabel('频率'); ylabel('幅值');

%%%%%%%单位脉冲序列 function xb=unitimp(n) % n=1:50; N=length(n); xb=zeros(1,N); xb(1)=1; k=-25:25;

[magX,angX]=fpres(xb,n,k);

subplot(3,1,1);

stem(xb)

title('单位脉冲序列');

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

title('单位脉冲序列幅度谱') xlabel('时间'); ylabel('幅值');

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

title('单位脉冲序列相位谱') xlabel('频率'); ylabel('幅值');

%%%%%%%%%指定序列

function hb=tedingxulie(n) % n=1:50; N=length(n); hb=zeros(1,N);

hb(1)=1;hb(2)=2.5;hb(3)=2.5;hb(4)=1;k=-25:25;

[magX,angX]=fpres(hb,n,k); subplot(3,1,1); stem(hb)

title('指定序列');

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

title('指定序列幅度谱') xlabel('时间'); ylabel('幅值');

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

title('指定序列相位谱') xlabel('频率'); ylabel('幅值');

%%%%%%%%%系统响应 clc; n=1:50; k=-25:25;