201332000110吕丹实验一.docVIP

实验一 离散时间信号和系统 一、实验目的 1、熟悉MATLAB的软件环境，学习MATLAB的编程方法 2、学习离散时间典型信号及其运算的MATLAB实现 3、学习离散时间系统的MATLAB实现并掌握求解系统对输入的响应的方法 二、实验原理 1、典型序列 （1）单位抽样序列impseq(n0,n1,n2) （2）单位阶跃序列stepseq(n0,n1,n2) （3）实指数序列a.^n （4）复指数序列exp((a+bj)*n) （5）正弦或余弦序列sin or cos 序列的运算 序列加sigadd(x1,n1,x2,n2) 序列乘sigmult(x1,n1,x2,n2) 数乘a*x 序列移位sigshift(x,m,n0) 序列反折sigfold(x,n) 序列卷积 系统的实现及响应的求解 差分方程： 在MATLAB中y(n)=filter(b,a,x) 系统的频率响应 三、实验内容 1、画出下列序列在给定区间的波形图 解：a. n=[-5:5]; %给定n的取值范围是-5《n《5 x=2*impseq(-2,-5,5)-impseq(4,-5,5); %求两个单位抽样序列的差 subplot(2,2,1); %将多个图画到一个平面表示是图排成行表示图排成列，表示图所在的位置从左到右从上到下的第一个位置。 stem(n,x); %绘制二维离散数据的火柴杆图 xlabel(n); %定义x轴为n ylabel(x(n)); %定义y轴为x(n) 图形 b. x(n)=n[u(n)-u(n-10)]+10*exp(-0.3(n-10))[u(n-10)-u(n-20)],0《n《20 n=[0:20]; %给定n的取值范围是0《n《20 x1=n.*(stepseq(0,0,20)-stepseq(10,0,20)); %求两个单位阶跃序列的差和矩阵序列的积 x2=10*exp(-0.3*(n-10)).*(stepseq(10,0,20)-stepseq(20,0,20)); %求两个单位阶跃序列的差和指数序列的积 x=x1+x2; %求两个定义序列的和 subplot(2,2,2); %将多个图画到一个平面表示是图排成行表示图排成列，表示图所在的位置从左到右从上到下的第个位置。stem(n,x); %输出离散序列，(plot连续)? xlabel(n); ylabel(x(n)); 图形 单位阶跃序列的波形在t=0处发生跃变，所以在t=0此点常常不予定义或定义为t=0.5 n=[0:50]; %给定n的取值范围是0《n《50 x=cos(0.04*pi*n)+0.2*randn(size(n)); %randn(size(A))，返回一个和A有同样维数大小的随机数组 subplot(2,2,1); stem(n,x); xlabel(n); ylabel(x(n)); 图形 给出如下差分方程y(n)-y(n-1)+0.9y(n-2)=x(n); a.计算并画出脉冲响应h(n),n=-20,......,100 b.计算并画出阶跃响应s(n),n=-20,......,100 c.由此h(n)规定的系统是否稳定？ 解：由给定的差分方程可知，其系数数组为： a.MATLAB程序 b=[1]; a=[1,-1,0.9]; x=impseq(0,-20,120); %单位抽样序列的起点为-20，终点为120，在n=0点处生成一个单位脉冲 n=[-20:120]; %给定n的取值范围是-20《n《120 h=filter(b,a,x); %输入x为滤波前序列，h为滤波结果序列，b/a提供滤波器系数，b为分子，a为分母