[信息与通信]数字信号实验代码.docVIP

%时域采样与频域采样实验 figure(1) A=444.18; a=50*pi*sqrt(2);omega=50*pi*sqrt(2);%模拟信号 Tp=64/1000; Fs=1000;%抽样频率 n=0:Tp*Fs-1; y=A*exp(-a*n*1/Fs).*sin(omega*n*1/Fs); subplot(2,2,1); stem(n,y,.); xlabel(抽样频率为1000HZ时信号Y的抽样图形); y1=fft(y,64);%频域抽样用fft计算DFT。。。。。。。 subplot(2,2,2); n=0:length(y1)-1; stem(n,y1,.); xlabel(抽样频率为1000HZ时信号y的64点得FFT); subplot(2,2,3); n1=0:length(angle(y))-1; plot(n1,angle(y)); xlabel(抽样频率为1000HZ时信号Y的相频特性图); figure(2) A=444.18; a=50*pi*sqrt(2);omega=50*pi*sqrt(2);%模拟信号 Tp=64/1000; Fs=300;%抽样频率 n=0:Tp*Fs-1; y=A*exp(-a*n*1/Fs).*sin(omega*n*1/Fs);//两向量相乘用点乘 y2=[y,zeros(1,46)]; subplot(2,2,