实验一离散信号的MATLAB实现..docVIP

实验一 离散信号的MATLAB实现 一、目的要求 1)学习典型的离散时间信号的MATLAB实现方法； 2)学习离散时间序列的基本运算：相加、相乘、移位等； 3)学习噪声的产生； 4)掌握两个序列的卷积和相关运算 5)掌握离散系统单位脉冲响应的求解。 二、实验内容 1)典型的离散信号的表示方法：用matlab产生单位抽样信号δ(n) 、单位阶跃序列u(n) 、矩形序列RN(n) 、 实指数序列 anu(n)、正弦序列sin(ωn)。编写程序、并画出图形。 n=-20:20;n0=0;n1=10; w0=pi/4; x=[(n-n0)==0]; x1=[(n-n0)=0]; x2=[(n-n0)=0 (n-n1)=0]; x3=0.9.^n.*x1; x4=sin(w0*n); subplot(511);stem(n,x);axis([ -20 20 0 2]);ylabel(\sigma(n)); subplot(512);stem(n,x1);axis([ -20 20 0 2]);ylabel(u(n)); subplot(513);stem(n,x2);axis([ -20 20 0 2]);ylabel(BN(n)); subplot(514);stem(n,x3);axis([ -20 20 0 2]);ylabel(anu(n)); subplot(515);stem(n,x4);axis([ -20 20 -2 2]);ylabel(sin(w0n)); xlabel(n); 2)离散信号的基本运算：对序列x(n)={2,3,4,1,2,5} ,n=0,1,2,3,4,5,的移位、乘法、加法、翻转及尺度变换。编写程序、并画出图形。 n=-10:10;k=2;N=length(n); x=zeros(size(n));x1=zeros(size(n)); y=zeros(size(n));y1=zeros(size(n)); y2=zeros(size(n));y3=zeros(size(n)); y4=zeros(size(n)); x(11:16)=[2 3 4 1 2 5];x1(11:16)=[1 2 3 4 5 6]; y(1+k:N)=x(1:N-k); y1(1:N-k)=x(1+k:N); y2=x.*x1; y3=x+x1; y4=fliplr(x); y5=x(1:2:N);n1=fix(n(1:2:N)/2); subplot(421);stem(n,x);xlabel(n);ylabel(x(n)); subplot(422);stem(n,x1);xlabel(n);ylabel(x1(n)); subplot(423);stem(n,y);xlabel(n);ylabel(x(n-2)); subplot(424);stem(n,y1);xlabel(n);ylabel(x(n+2)); subplot(425);stem(n,y2);xlabel(n);ylabel(x(n)*x1(n)); subplot(426);stem(n,y3);xlabel(n);ylabel(x(n)+x1(n)); subplot(427);stem(n,y4);xlabel(n);ylabel(x(-n)); subplot(428);stem(n1,y5);xlabel(n);ylabel(x(2n)); 3)噪声的产生：产生方差为1，2，0.5的白噪声。 a=[0.5 1 2];N=1000; x=a*randn(1,N); subplot(311);plot(1:N,x(1,:)); subplot(312);plot(1:N,x(2,:)); subplot(313);plot(1:N,x(3,:)); 4)用matlab求系统的单位脉冲响应 、频率响应和零极点分布。编写程序、并画出图形。 二阶离散系统的差分方程 : y(n)-2.5y(n-1)+y(n-2)=f(n)-2f(n-1) Y(z)=(1-2z-1)/(1-2.5z-1+z-2) b=[1 -2]; a=[1 -2.5 1]; [h,t]=impz(b,a,40); [H,W]=freqz(b,a,256,whole,1); Hr=abs(H); Hphase=angle(H);Hphase=unwrap(Hphase); subplot(221);stem(t,h);xlabel(n);ylabel(h(n)); subplot(222);plot(W,Hr);xlabel(w);ylabel(H(ejw)); subplot(22