实验一序列的产生及绘图.docVIP

实验一 序列的产生及绘图 一、实验目的 1．熟悉信号处理软件MATLAB的使用。 2．离散信号的基本运算实现。 3．了解基本序列及复杂序列的产生方法。 4．运用卷积方法观察系统的时域特性。 5．掌握线性时不变系统的频域表示方法。 二、实验内容 1．熟悉扩展函数 2．运行例题程序 3．编程实现下列内容 （1）利用扩展函数产生序列并画图 (a) -5=n=5 (b) 和 0=n=50 w(n)为白噪声 函数为 w=randn(size(n)) （2）设线性移不变系统的抽样响应为 输入序列为 求系统输出 y(n)并画图 提示: 输出为输入和抽样响应的卷积 三、实验报告要求 1．记录例题程序的实验结果、图形。 2．写出自己编写的程序并记录结果、图形 %如果要了解哪个函数的应用方法请用help命令 如help zreos %本软件中 * 表示乘法, 卷积用函数 conv 或修改后的卷积 conv_m %以下是7个扩展函数 %扩展函数1~7的用法和该软件自带函数用法一致,即在调用时要将实参代入 %例：应用扩展函数3需要输入x1(n),x2(n)的值。 %在Command Window（命令窗口）中输入 % n1=1:5; % n2=2:6; % x1=[1 3 5 7 9]; % x2=[2 4 6 8 10]; % [y,n]=sigadd(x1,n1,x2,n2) 即可得两序列相加的结果 %7个扩展函数要分别存到不同的文件中,并且文件名要和该扩展函数的函数名一致 %如产生单位取样序列的函数所存文件的文件名必须为 impseq %1.单位取样序列 x(n)=delta(n-n0) 要求n1=n0=n2 function[x,n]=impseq(n0,n1,n2) n=[n1:n2]; x=[(n-n0)==0]; %2.单位阶跃序列 x(n)=u(n-n0) 要求n1=n0=n2 function[x,n]=stepseq(n0,n1,n2) n=[n1:n2]; x=[(n-n0)=0]; %3.信号加 y(n)=x1(n)+x2(n) %find函数：找出非零元素的索引号 %x1:第一个序列的值，n1：x1的索引号，任意整数，可以不连续以及不按顺序排列 %x2:第二个序列的值，n2：x2的索引号，任意整数，可以不连续以及不按顺序排列 function[y,n]=sigadd(x1,n1,x2,n2) n=min(min(n1),min(n2)):max(max(n1),max(n2)); y1=zeros(1,length(n)); y2=y1; y1(find((n=min(n1))(n=max(n1))))=x1; y2(find((n=min(n2))(n=max(n2))))=x2; y=y1+y2; %4.信号乘 y(n)=x1(n)*x2(n) %find函数：找出非零元素的索引号 %x1:第一个序列的值，n1：x1的索引号，任意整数，可以不连续以及不按顺序排列 %x2:第二个序列的值，n2：x2的索引号，任意整数，可以不连续以及不按顺序排列 function[y,n]=sigmult(x1,n1,x2,n2) n=min(min(n1),min(n2)):max(max(n1),max(n2)); y1=zeros(1,length(n)); y2=y1; y1(find((n=min(n1))(n=max(n1))))=x1; y2(find((n=min(n2))(n=max(n2))))=x2; y=y1.*y2; %5.移位 y(n)=x(n-n0) %x:移位前序列的值，m：x的索引号，任意整数，可以不连续以及不按顺序排列 %y:移位后序列的值，n：y的索引号，任意整数，可以不连续以及不按顺序排列 function[y,n]=sigshift(x,m,n0) n=m+n0; y=x; %6.翻褶 y(n)=x(-n) %x: 翻褶前序列的值，x的索引号n：任意整数，可以不连续以及不按顺序排列 %y: 翻褶后序列的值，y的索引号n：任意整数，可以不连续以及不按顺序排列 function[y,n]=sigfold(x,n) y=fliplr(x); n=-fliplr(n); %7.修改后的卷积(系统自带卷积函数(conv)只能得到输出值,而无法表示取值范围(n),修改后的卷积既给出了卷积值,也给出了它的取值范围.) function[y,ny]=conv_m(x,nx,h,nh) n