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

实验一 离散信号的产生与运算实验任务：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） 已知，试显示、和在区间的波形。n=0:20;x=n.*sin(n);x1=fliplr(x);n1=-fliplr(n);x2=n/2.*sin(n/2);subplot(4,1,1),plot(n,x);title( X(n)ylabel(x(n)subplot(4,1,2),plot(n1,x1);title(X(-n)ylabel(x(n)subplot(4,1,3),plot(n,-x);title(-X(n)ylabel(x(n)subplot(4,1,4),plot(n,x2);title(X(n/2)ylabel(x(n)（6） 编写Matlab程序，描绘下列信号序列的卷积波形： n = 0:10x1 = n=0; x2 = exp(0.2*n).*x1;subplot(2,1,1), plot(n,x1); title(x1(n)subplot(2,1,2), plot(n,x2); title(x2(n) （7）已知一个II