Matlab系列之离散信号的产生与运算的实验报告.docVIP

实验一 离散信号的产生与运算 实验任务： 1、在Matlab中运行各例题程序，理解时域离散信号的产生方法和离散序列的运算性质，了解各条语句的意义。 2、完成以下习题： （1）编写程序产生下列离散序列： ① n1 = -3; n2 = 4; n0 = 0; n = n1:n2; x = [n==n0]; % x在n=n0时为1，其余为0 stem(n,x,filled); %filled:序列圆心处用实心圆表示 axis([n1,n2,0,1.1*max(x)]) title(单位抽样序列) xlabel(time(n)) ylabel(Amplitude:x(n)) ② n1 = -5; n2 = 5; n0 = 0; n = n1:n2; x = [n=n0]; % x在n=n0时为1，其余为0 stem(n,x,filled); axis([n1,n2,0,1.1*max(x)]) title(单位阶跃序列) xlabel(time(n)) ylabel(Amplitude:x(n)) ③ n1 = 16; sigma = 0.1; omega = 1.6; n = 0:n1; x = exp((sigma+i*omega)*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(复指数序列的虚部) ④ f = 1/8; Um = 3; nt = 2; 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); title(连续时间正弦信号和离散正弦序列) ylabel(x(t)) subplot(2,1,2),stem(tn,x); ylabel(x(n)) ⑤ n = -20:20; f = sinc(n/5); subplot(2,1,1),plot(n,f); title(sinc信号) subplot(2,1,2),stem(n,f); title(sinc序列) （2）一个连续的周期性三角波信号频率为50Hz，信号幅度在0～+2V之间，在窗口上显示2个周期的信号波形，对信号的一个周期进行16点采样来获得离散信号。试显示连续信号和采样获得的离散信号波形。（用sawtooth产生三角波） f = 50; Um = 1; nt = 2; N = 16; T = 1/f; dt = T/N; n = 0:nt*N-1; tn = n*dt; x = (Um+1)*sawtooth (2*f*pi*tn); subplot(2,1,1),plot(tn,x); title(连续时间三角信号和离散三角序列) ylabel(x(t)) subplot(2,1,2),stem(tn,x); ylabel(x(n)) （3） 用Matlab实现下列信号序列： n1 = -5; n2 = 5; n = n1:n2; x1= [n==-3]; % x在n=-3时为1，其余为0 x2= [n==4]; % x在n=4时为1，其余为0 x=x1+2*(x2) stem(n,x,filled); %filled:序列圆心处用实心圆表示 axis([n1,n2,0,1.1*max(x)]) title(单位抽样序列) xlabel(time(n)) ylabel(Amplitude:x(n)) （4） 已知，显示、和在区间的波形。 f = 1/10; Um = 3; nt = 2; N = 20; T = 1/f; dt = T/N; n = 0:nt*N-1; tn = n*dt; x = Um*cos(2*f*pi*tn); subplot(3,1,1),plot(tn,x); title(连续时间余弦信号X(n)) ylabel(x(t)) n=n+3; subplot(3,1,2),plot(tn-3,x); title(连续时间余弦信号X(n-3)) ylabel(x(t)) n=n-6; subplot(3,1,3),plot(tn+3,x); title(连续时间余弦信号X(n+3)) ylabel(x(t)) （5） 已知，试显示、、和在区间的波形