实验六 离散时间信号与系统.docVIP

实验六? 离散时间信号与系统 [实验目的] 1． 2． [实验原理与步骤] 我们所接触的信号大多为连续信号，为使之便于处理，往往要对其进行采样，对信号抽样并保证其能完全恢复，对抽样频率有一定的限制。 一. 离散时间正弦序列的MATLAB实现 解：由MATLAB产生的信号实际上是离散的。可以使用命令stem(n,x)作出离散序列x(n)各离散时间点n的取值。 A=2;N=12;pha=pi/6; omega=2*pi/12; n=10:10; x=A*sin(omega*n+pha); stem(n,x); ylabel(‘x(n)’);xlabel(‘Time index n’); title(‘离散时间正弦序列’)； 连续时间正弦信号如图1（a）所示，离散时间正弦序列如图1（b）所示。  图1（a）  图1（b） 若要表示具有特定采样频率的信号，需定义时间轴向量。例如： 用MATLAB实现如下： ％离散信号波形 f1=50;f2=100;fs=1000; t=0:1/fs:1;n=t*fs; y=sin(2*pi*f1*t)+2*sin(2*pi*f2*t); subplot(211); plot(t(1:50),y(1:50)); ％画y(t)的前0.05s的值 subplot(212); stem(n(1:50),y(1:50)); ％画y(k)的前50各样点值 图2 解：MATLAB程序实现如下： %采样定理例题 Dt=0.00005; t=-0.005:Dt:0.005; xa=exp(-1000*abs(t)); ％模拟信号 Wmax=2*pi*2000; K=500;k=0:1:K; W=k*Wmax/K; Xa=xa*exp(-j*t*W)*Dt; Xa=real(Xa); ％连续时间傅立叶变换 W=[-fliplr(W),W(2:501)]; ％频率从-Wmax to Wmax Xa=[fliplr(Xa), Xa(2:501)]; ％Xa介于-Wmax to Wmax subplot(421) plot(t*1000,xa) xlabel(时间（毫秒）) ylabel(xa(t)) title(模拟信号) subplot(422) plot(W/(2*pi*1000),Xa*1000) xlabel(频率(kHz)) ylabel(Xa(jw)) title(连续时间傅立叶变换) Ts=0.0002;n=-25:1:25; x=exp(-1000*abs(n*Ts)); ％离散信号 K=500; k=0:1:K; w=pi*k/K; X=x*exp(-j*n*w); X=real(X); ％离散时间傅立叶变换 w=[-fliplr(w), w(2:K+1)]; X=[fliplr(X), X(2:K+1)]; subplot(423) stem(n*Ts*1000,x) xlabel(时间（毫秒）) ylabel(xl(n)) title(离散信号Ts=0.2毫秒) subplot(424) plot(w/pi, X) xlabel(频率(弧度)) ylabel(Xl(w)) title(离散时间傅立叶变换) Ts=0.0005;n=-10:1:10; x=exp(-1000*abs(n*Ts)); ％离散信号 K=500; k=0:1:K; w=pi*k/K; X=x*exp(-j*n*w); X=real(X); ％离散时间傅立叶变换 w=[-fliplr(w), w(2:K+1)]; X=[fliplr(X), X(2:K+1)]; subplot(425) st