q710.m

<MATLAB及在电子信息工程中的应用>原程序,

%《MATLAB在电子信息课程中的应用》第七章例7.10程序q710
% 时域采样及其频谱
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月

clear;close all;
fs=10000; fs1=1000; fs2=400; fs3=200;	% 设置四种采样频率
t=0:1/fs:0.1;		% 采集信号长度为0.1秒
A=444.128; a=50*sqrt(2)*pi; b=a;
xa=exp(-a*t).*sin(b*t); 
k=0:511; f=fs*k/512;		%由wk=2πk/512=2πfT求得模拟频率f
Xa=dtft(xa,2*pi*k/512);   % 近似模拟信号频谱
T1=1/fs1;t1=0:T1:0.1;	% 采集信号长度为0.1秒                             
x1=A*exp(-a.*t1).*sin(b*t1);		%1kHz采样序列x1(n)
X1=dtft(x1,2*pi*k/512);      % x1(n)的512点dtft
T2=1/fs2;t2=0:T2:0.1;		% 采集信号长度为0.1秒
x2=A*exp(-a.*t2).*sin(b.*t2);	%400Hz采样序列x2(n)   
X2=dtft(x2,2*pi*k/512);      % x2(n)的512点dtft
T3=1/fs3;t3=0:T3:0.1;		% 采集信号长度为0.1秒
x3=A*exp(-a.*t3).*sin(b.*t3); 	%200Hz采样序列x3(n)                
X3=dtft(x3,2*pi*k/512);      % x3(n)的512点dtft
figure(1);
subplot(3,2,1);plot(t,xa);
axis([0,max(t),min(xa),max(xa)]);title('模拟信号');
xlabel('t(s)');ylabel('Xa(t)');line([0,max(t)],[0,0])
subplot(3,2,2);plot(f,abs(Xa)/max(abs(Xa)));
title('模拟信号的幅度频谱');axis([0,500,0,1])
xlabel('f (Hz)');ylabel('|Xa(jf)|');
subplot(3,2,5);stem(t1,x1,'.');
line([0,max(t1)],[0,0]);axis([0,max(t1),min(x1),max(x1)])
title('采样序列x1(n)（fs1=1kHz）');xlabel('n');ylabel('X1(n)');
f1=fs1*k/512;
subplot(3,2,6);plot(f1,abs(X1)/max(abs(X1)));
title('x1(n)的幅度谱');xlabel('f (Hz)');ylabel('|X1(jf)|');
figure(2);
subplot(3,2,1);stem(t2,x2,'.');
line([0,max(t2)],[0,0]);axis([0,max(t2),min(x2),max(x2)]);
title('采样序列x2(n)（fs2=400Hz）');xlabel('n');ylabel('X2(n)');
f=fs2*k/512;
subplot(3,2,2);plot(f,abs(X2)/max(abs(X2)));
title('x2(n)的幅度谱');xlabel('f (Hz)');ylabel('|X2(jf)|');
subplot(3,2,5);stem(t3,x3,'.');
line([0,max(t3)],[0,0]);axis([0,max(t3),min(x3),max(x3)]);
title('采样序列x3(n)（fs3=200Hz）');xlabel('n');ylabel('X3(n)');
f=fs3*k/512;
subplot(3,2,6);plot(f,abs(X3)/max(abs(X3)));
title('x3(n)的幅度谱');xlabel('f (Hz)');ylabel('|X3(jf)|')