实验一：用MATLAB产生时域离散信号

实验一：用MATLAB产生时域离散信号

一、实验目的

1、了解常用时域离散信号及其特点

2、掌握用MATLAB产生时域离散信号的方法

二、实验内容及步骤

1、编写程序，产生以下离散序列：

（1）f(n)=δ(n) (-3

n1=-3;n2=4;n0=0;

n=n1:n2;

x=[n==n0];

stem(n,x,'filled');

axis([n1,n2,0,1.1*max(x)]);

xlabel('时间(n)');ylabel('幅度x(n)');

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

（2）f(n)=u(n) (-5

n1=-5;n2=5;n0=0;

n=n1:n2;

x=[n>=n0];

stem(n,x,'filled');

axis([n1,n2,0,1.1*max(x)]);

xlabel('时间(n)');ylabel('幅度x(n)');

title('单位阶跃序列');

box

（3）f(n)= e(0.1+j1.6∏)n (0

n1=16;a=0.1;w=1.6*pi;

n=0:n1;

x=exp((a+j*w)*n);

subplot(2,2,1);plot(n,real(x));

title('复指数信号的实部');

subplot(2,2,3);stem(n,real(x),'filled');

title('复指数序列的实部');

subplot(2,2,2);plot(n,imag(x));

title('复指数信号的虚部');

subplot(2,2,4);stem(n,imag(x),'filled');

title('复指数序列的虚部');

box

（4）f(n)=3sin(nП/4) (0

f=1/8;Um=3;nt=3;

N=20;T=1/f;

dt=T/N;

n=0:nt*N-1;

tn=n*dt;

x=Um*sin(2*f*pi*tn);

subplot(2,1,1);plot(tn,x);

axis([0,nt*T,1.1*min(x),1.1*max(x)]);

ylabel('x(t)');

subplot(2,1,2);stem(tn,x);

axis([0,nt*T,1.1*min(x),1.1*max(x)]);

ylabel('x(n)');

box

2、一个连续的周期性方波信号频率为200Hz，信号幅度在-1～+1V之间，要求在图形窗口上显示其两个周期的波形。以4kHz的频率对连续信

号进行采样，编写程序生成连续信号和其采样获得的离散信号波形。

解 f=200;nt=2;

N=20;T=1/f;

dt=T/N;

n=0:nt*N-1;

tn=n*dt;

x=square(2*f*pi*tn,50);

subplot(2,1,1);plot(tn,x);

axis([0,nt*T,1.1*min(x),1.1*max(x)]);

ylabel('x(t)');

subplot(2,1,2);stem(tn,x);

axis([0,nt*T,1.1*min(x),1.1*max(x)]);

ylabel('x(n)');

Box

f=4000;nt=2;

N=20;T=1/f;

dt=T/N;

n=0:nt*N-1;

tn=n*dt;

x=square(2*f*pi*tn,50);

subplot(2,1,1);plot(tn,x);

axis([0,nt*T,1.1*min(x),1.1*max(x)]);

ylabel('x(t)');

subplot(2,1,2);stem(tn,x);

axis([0,nt*T,1.1*min(x),1.1*max(x)]);

ylabel('x(n)');

Box

2、思考题：通过例题程序，你发现采样频率Fs、采样点数N、采样时间间隔dt在程序编写中有怎样的联系？使用时需注意什么问题？

采样频率Fs与T的关系为Fs=1/T，而采样时间dt=T/N

所以dt=1/N*Fs

在使用时需要注意Fs与T的关系，不要将其混淆。

三、实验小结

通过本次的实验有以下两点收获：

（1） 学会了MATLAB软件的使用

（2） 了解常用时域离散信号及其特点