信号处理与系统课程设计指导书分析.docVIP

《信号处理与系统分析》 课程设计指导书 南通大学 电子信息学院　 信息工程系 2013年5月 前言 《信号处理与系统》是南通大学杏林学院通信信息类专业的一门专业基础课程，其理论性强，是其它后续专业课程的基础。 开设该课程设计的重要意义在于：首先，从帮助学习《信号处理与系统》课程的角度讲，学生Matlab的应用方法，为后续专业课的学习打下了坚实的基础；另外实践环节使学生在综合使用现代电子信息技术和手段进行设计、制作创新能力。Matlab语言的编程方法； 二、熟悉用于《信号处理与系统分析》的Matlab主要函数的应用； 三、记录实验结果（包括波形和数据），撰写课程设计报告。 主要内容及步骤 一、连续系统的时域分析 信号的产生 （1）编写生成连续阶跃信号u(t-t0)及冲激信号δ(t-t0)的函数 ： function [x,y]=jieyue(t1,t2,t0) 和 function [x,y] =chongji(t1,t2,t0)，信号的时间变量取值区间为t1～t2，t0为阶跃点或冲激点处的时间，x为信号的时间向量，y为相应的信号值向量。（提示：冲激信号可以用时间宽度为dt、高度为1/dt的矩形脉冲来近似表示。当dt很小时，矩形脉冲信号可近似认为是冲激信号。在对该矩形脉冲信号采点取样后（设取样间隔为dt），信号值y的第1+(t0- t1)/dt个元素的值为1/dt，其余元素的值为0。） 参考程序: function [t,y1]=jieyue(t1,t2,t0) dt=0.01; ttt=t1:dt:t0-dt; tt=t0:dt:t2; t=t1:dt:t2; n=length(ttt); nn=length(tt); u=zeros(1,n); uu=ones(1,nn); y1=[u,uu]; return function [t,y2]=chongji(t1,t2,t0) dt=0.01; t=t1:dt:t2; n=length(t); y2(1:n)=0; y2(1,(t0-t1)/dt+1)=1/dt; 说明：上述定义的函数中， zeros(M,N)和ones(M,N)的功能是分别产生行为M以及列为N的矩阵。 以上两函数不可以直接运行，被调用后方可运行。 （2）调用上述函数，产生信号u(t-2),δ(t-4),e-tu(t),-6s≤t≤6s,并画出波形。 (观察用stairs和plot函数画图的区别) 举例：[t1,y1]=jieyue(-6,6,2); stairs(t1,y1) 信号的卷积 （1）编写求f(t)=f1(t)*f2(t)的函数： function [t,f]=myconv(f1,f2,t1,t2,dt), 其中t1、t2和t分别表示信号 f1(t)、f2(t)和f(t)的时间向量,dt代表信号取样点的时间间隔。 提示：将f1(t)和f2(t)离散化，卷积积分公式可转化为如下的形式：,其中表达式中的Δ代表信号离散化后取样点间的时间间隔，Δ越小，则误差越小。 而式中的可以通过调用Matlab现有的求序列卷积和的函数conv来实现。 另外，函数的返回值向量t由向量t1和t2确定。为了使t包含了f(t)的非零样点值的时间区间，则t的第一个元素t(1)= t1(1)+t2(1)，t的长度为t1和t2的长度之和减1。 参考程序: function [k,f]=myconv(f1,f2,k1,k2,p) %The function of compute f=f1*f2 % f: 卷积和序列f(k)对应的非零样值向量 % k： 序列f(k)的对应序号向量 % f1: 序列f1(k)非零样值向量 % f2: 序列f2(k)的非零样值向量 % k1: 序列f1(k)的对应序号向量 % k2: 序列f2(k)的对应序号向量 f=conv(f1,f2)*p; %计算序列f1与f2的卷积和f k0=k1(1)+k2(1); %计算序列f非零样值的起点位置 k3=length(f1)+length(f2)-2;%计算卷积和f的非零样值的宽度 k=k0:p:k0+k3*p; %确定卷积和f非零样值的序号向量 （2）调用上述已编函数myconv，求x1(t)=e-t[u(t)-u(t-2)]和x2(t)=δ