信号系统实验1

电气学科大类级《信号与控制综合实验》课程实验报告(基本实验一:信号与系统基本实验)姓名学号专业班号同组者1 学号专业班号同组者2 学号专业班号指导教师日期实验成绩评阅人综合实验和实验报告要求信号与控制综合实验，是集多门技术基础课程以及其它延伸课程理论于一体的综合性实验课程，需要综合多门学科理论知识和实验方法来体现，因此，实验目的不是简单的课程理论验证和练习，而是综合应用、研究开发、设计创新。

应采用尽可能好的设计，使所设计的电路和系统达到要实现的功能，步骤和方案自行拟定，实现对设计思路的实验验证。

完成多个实验项目的，应将实验内容整理综合后写成一份总报告，以利于锻炼整理归纳和总结能力，在总报告中以第二级标题形式依次写下所完成的实验项目、内容及实验设计过程。

实验报告按“题目、目录、正文（分所完成的各实验项目）、结论、心得与自我评价、参考文献”6个部分撰写；正文主要包括以下几个内容：任务和目标、总体方案设计（原理分析与方案设计特点，选择依据和确定）、方案实现和具体设计（过程）、实验设计与实验结果、结果分析与讨论。

（格式方面请注意：每个图应该有图号和图名，位于图的下方，同一图号的分图应在同一页，不要跨页；每个表应该有表号和表名，位于表的上方，表号表名与表（数据）也应在同一页，不要跨页；建议各部分题目采用四号黑体、设计报告内容文字采用小四号宋体）注：报告中涉及实验指导书或教材内容，只需注明引用位置，不必在报告中再加以阐述。

不得不加引用标记地抄袭任何资料。

每一基本实验部分按计划学时100分成绩计算（100％），需要完成60分的实验项目；实验报告、设计部分和创新研究内容另外计分（分别为10％、20％和10％）。

再按照学时比例与本课程其它部分实验综合成为总实验成绩。

每一部分实验均为：基本实验：0～60分，考核基本理论的掌握和基本操作技能、实验室道德规范；实验报告：0～10分，考核思考和总结表述能力；完成设计性实验：0～20分，评价设计能力；完成创新性实验：0～10分，鼓励创新。

实验总结班级：10电子班学号：1039035 姓名：田金龙这学期的实验都有：信号的时域分析、线性时不变系统的时域分析、连续时间信号系统的频域分析、连续时间在连续时间信号的频域LTI系统的复频域分析、连续时间LTI系统的频域分析。

在这学期的学习中学习了解到很多关于信号方面的处理方法加上硬件动手的实践能力，让我对课堂上所学到的知识有了更深层次的理解也加深了所学知识的印象。

下面则是对每次实验的分析和总结：实验一：信号的时域分析在第一次试验中进行信号的时域分析还有的就是学会使用MATLAB软件来利用它实现一些相关的运算并且绘制出相关的信号图。

在时域分析中掌握连续时间信号和离散时间信号的描述方法，并能够实现各种信号的时域变化和运算。

了解单位阶跃信号和单位冲激信号的拓展函数，以便于熟悉这两种函数在之后的程序中的应用。

在能够对简单信号的描述的前提下，通过一些简单的程序，实现信号的分析，时域反相，时域尺度变换和周期信号的描述。

clear,close alldt=0.01;t=-2:dt:2;x=u(t);plot(t,x)title('u signal u(t)')grid on连续时间信号的时域分析后，则是离散时间信号的仿真。

通过对连续时间信号的描述和对离散时间信号的描述，发现它们的不同之处在于对时间的定义和对函数的图形描述。

在离散时间信号的图形窗口描述时，使用的是stem（n，x）函数。

在硬件实验中，使用一些信号运算单元，加法器，减法器，倍乘器，反相器，积分器和微分器。

输入相应的简单信号，观察通过不同运算单元输出的信号。

实验二：线性时不变系统的时域分析在线性时不变系统的时域分析中主要研究的就是信号的卷积运算，学会进行信号的卷积运算和MATLAB对卷积运算的实现。

而系统则通常是由若干部件或单元组成的一个整体，根据系统所处理的信号不同，系统又有多种不同的分类。

而在学习总最常研究的则是线性时不变系统，而线性时不变系统则是形同同时满足齐次性和叠加性。

matlab软件仿真实验（信号与系统）（1）《信号与系统实验报告》学院：信息科学与⼯程学院专业：物联⽹⼯程姓名：学号：⽬录实验⼀、MATLAB 基本应⽤实验⼆信号的时域表⽰实验三、连续信号卷积实验四、典型周期信号的频谱表⽰实验五、傅⽴叶变换性质研究实验六、抽样定理与信号恢复实验⼀MATLAB 基本应⽤⼀、实验⽬的：学习MATLAB的基本⽤法，了解 MATLAB 的⽬录结构和基本功能以及MATLAB在信号与系统中的应⽤。

⼆、实验内容：例⼀已知x的取值范围，画出y=sin(x)的图型。

x=0:0.05:4*pi;y=sin(x);plot(y)例⼆计算y=sin(π/5)+4cos(π/4)例三已知z 取值范围，x=sin(z);y=cos(z);画三维图形。

z=0:pi/50:10*pi;x=sin(z);y=cos(z);plot3(x,y,z)xlabel('x')ylabel('y')zlabel('z')例四已知x的取值范围，⽤subplot函数绘图。

参考程序：x=0:0.05:7;y1=sin(x);y2=1.5*cos(x);y3=sin(2*x);y4=5*cos(2*x);subplot(2,2,1),plot(x,y1),title('sin(x)')subplot(2,2,2),plot(x,y2),title('1.5*cos(x)')subplot(2,2,3),plot(x,y3),title('sin(2*x)')subplot(2,2,4),plot(x,y4),title('5*cos(2*x)')连续信号的MATLAB表⽰1、指数信号：指数信号Ae at在MATLAB中可⽤exp函数表⽰，其调⽤形式为：y=A*exp(a*t) (例取 A=1,a=-0.4)参考程序：A=1;a=-0.4;t=0:0.01:10;ft=A*exp(a*t);plot(t,ft);grid on;2、正弦信号：正弦信号Acos(w0t+?)和Asin(w0t+?)分别由函数cos和sin表⽰，其调⽤形式为：A*cos(w0t+phi) ;A*sin(w0t+phi) (例取A=1，w0=2π，?=π/6) 参考程序：A=1;w0=2*pi; phi=pi/6; t=0:0.001:8;ft=A*sin(w0*t+phi);plot(t,ft);grid on ;3、抽样函数：抽样函数Sa(t)在MATLAB中⽤sinc函数表⽰，其定义为：sinc(t)=sin(πt)/( πt)其调⽤形式为：y=sinc(t)参考程序：t=-3*pi:pi/100:3*pi;ft=sinc(t/pi);plot(t,ft);grid on;4、矩形脉冲信号：在MATLAB中⽤rectpuls函数来表⽰，其调⽤形式为：y=rectpuls(t,width),⽤以产⽣⼀个幅值为1，宽度为width,相对于t=0点左右对称的矩形波信号，该函数的横坐标范围由向量t决定，是以t=0为中⼼向左右各展开width/2的范围，width的默认值为1。

实验一 连续时间信号分析一、实验目的（一）掌握使用Matlab 表示连续时间信号1、学会运用Matlab 表示常用连续时间信号的方法2、观察并熟悉常用信号的波形和特性（二）掌握使用Matlab 进行连续时间信号的相关运算1、学会运用Matlab 进行连续时间信号的时移、反褶和尺度变换2、学会运用Matlab 进行连续时间信号微分、积分运算3、学会运用Matlab 进行连续时间信号相加、相乘运算4、学会运用Matlab 进行连续时间信号卷积运算二、实验条件装用Matlab R2015a 的电脑。

三、实验内容1、利用Matlab 命令画出下列连续信号的波形图。

（1）)4/3t (2cos π+ 程序：t=-3:0.01:3; ft=2*cos(3*t+pi/4); plot(t,ft)图像：（2）)t (u )e 2(t--程序：t=-6:0.01:6; ut=(t>=0);ft=(2-1*exp(-t)).*ut; plot(t,ft)图像：（3）)]2()(u )][t (cos 1[--+t u t π 程序：t=-6:0.01:6; ut=(t>=0); ut2=(t>=2);ft=(1+cos(pi*t)).*(ut-ut2); plot(t,ft)图像：2、利用Matlab 命令画出复信号)4/t (j 2e )t (f π+=的实部、虚部、模和辐角。

程序：t=0:0.01:20;ft=2*exp(1j*(t+pi/4));subplot(2,2,1);plot(t,real(ft));title('ʵ²¿');axis([-0.5,20,-2.5,2.5]); subplot(2,2,2);plot(t,imag(ft));title('Ð鲿');axis([-0.5,20,-2.5,2.5]); subplot(2,2,3);plot(t,abs(ft));title('Ä£');axis([-0.5,20,-0.5,2.5]); subplot(2,2,4);plot(t,angle(ft));title('·ø½Ç');axis([-0.5,20,-3.5,3.5]);图像：3、已知信号的波形如下图所示：试用Matlab 命令画出()()()()2332----t f t f t f t f ，，，的波形图。

信号与系统实验一连续时间系统卷积的数值计算实验目的1 加深对卷积概念及原理的理解;2 掌握借助计算机计算任意信号卷积的方法.实验原理()()()tototftoftf d21⎰∞∞--=卷积实验流程图源程序＃include 〈stdio。

h〉float u（float t）｛while(t〉=0）return（1）;while(t<0）return(0）;｝float f1(float t)｛return（u(t+2)—u(t-2));｝float f2（float t)｛return（t*(u（t）-u(t-2))+（4-t)＊（u(t-2）—u（t-4））)；｝main（)｛float t,i，j，result=0;for(i=—2;i<=6;i=i+0.1）｛result=0；for(j=0;j<=4;j=j+0.1）result+=f2（j）*f1(i—j)*0。

1;printf("%.1f\t%.2f\t",i，result);}printf （"\n”）；}实验数据—2。

0 0。

00 —1.9 0.01 —1.8 0。

03 —1。

7 0.06 -1。

6 0.10—1。

5 0.15 -1。

4 0。

21 —1.3 0。

28 -1。

2 0。

36 —1.1 0。

45-1.0 0.55 -0.9 0.66 -0。

8 0。

78 -0.7 0。

91 -0。

6 1。

05-0.5 1.20 -0.4 1。

36 -0.3 1.53 -0。

2 1.71 —0。

1 1.90 0.0 2.10 0。

1 2。

29 0。

2 2.47 0。

3 2。

64 0。

4 2。

800。

5 2。

95 0。

6 3.09 0.7 3.22 0。

8 3。

34 0.9 3。

451。

0 3。

55 1。

1 3。

64 1.2 3。

72 1.3 3.79 1。

4 3.85 1。

5 3。

90 1.6 3.94 1.7 3。

1.已知系统的差分方程如下式：y1(n)=0.9y1(n-1)+x(n)程序编写如下：(1)输入信号x(n)=R10 (n)，初始条件y1(-1)=1，试用递推法求解输出y1(n)；a=0.9; ys=1; %设差分方程系数a=0.9,初始状态: y(-1)=1xn=ones(1,10); %矩型序列R10(n)=u(n)-u(n-10)，定义其宽度为0~9n=1:35; %设差分方程系数a=0.9,初始状态: y(-1)=1xn=sign(sign(10-n)+1);B=1;A=[1,-a]; %差分方程系数xi=filtic(B,A,ys); %由初始条件计算等效初始条件输入序列xiyn=filter(B,A,xn,xi); %调用filter解差分方程,求系统输出y(n)n=0:length(yn)-1;subplot(2,1,1);stem(n,yn,'linewidth',2); axis([-5,15,0,8]); grid ontitle('图(a) y1(n)=0.9y1(n-1)+x(n) 初始条件y1(-1)=1 ');xlabel('n');ylabel('y(n)')(2) 输入信号x(n)=R10 (n)，初始条件y1(-1)=0，试用递推法求解输出y1(n)。

a=0.9; ys=0; %设差分方程系数a=0.9,初始状态: y(-1)=1xn=ones(1,10); %矩型序列R10(n)=u(n)-u(n-10)B=1;A=[1,-a]; %差分方程系数xi=filtic(B,A,ys); %由初始条件计算等效初始条件输入序列xiyn=filter(B,A,xn,xi); %调用filter解差分方程,求系统输出y(n)n=0:length(yn)-1;subplot(2,1,2);stem(n,yn, 'linewidth',2); axis([-5,15,0,8]); grid ontitle('图(b) y1(n)=0.9y1(n-1)+x(n) 初始条件y1(-1)=0 ');xlabel('n');ylabel('y(n)') 图形输出如下：-505101502468图(a) y1(n)=0.9y1(n-1)+x(n) 初始条件y1(-1)=1ny (n )-55101502468图(b) y1(n)=0.9y1(n-1)+x(n) 初始条件y1(-1)=0ny (n )2. 已知系统差分方程为： y 1(n )=0.9y 1(n -1)+x (n ) 用递推法求解系统的单位脉冲响应h (n )，要求写出h (n )的封闭公式，并打印h (n )~n 曲线。

1.已知系统的差分方程如下式：y1(n)=0.9y1(n-1)+x(n)程序编写如下：(1)输入信号x(n)=R10 (n)，初始条件y1(-1)=1，试用递推法求解输出y1(n)；a=0.9; ys=1; %设差分方程系数a=0.9,初始状态: y(-1)=1xn=ones(1,10); %矩型序列R10(n)=u(n)-u(n-10)，定义其宽度为0~9n=1:35; %设差分方程系数a=0.9,初始状态: y(-1)=1xn=sign(sign(10-n)+1);B=1;A=[1,-a]; %差分方程系数xi=filtic(B,A,ys); %由初始条件计算等效初始条件输入序列xiyn=filter(B,A,xn,xi); %调用filter解差分方程,求系统输出y(n)n=0:length(yn)-1;subplot(2,1,1);stem(n,yn,'linewidth',2); axis([-5,15,0,8]); grid ontitle('图(a) y1(n)=0.9y1(n-1)+x(n) 初始条件y1(-1)=1 ');xlabel('n');ylabel('y(n)')(2) 输入信号x(n)=R10 (n)，初始条件y1(-1)=0，试用递推法求解输出y1(n)。

a=0.9; ys=0; %设差分方程系数a=0.9,初始状态: y(-1)=1xn=ones(1,10); %矩型序列R10(n)=u(n)-u(n-10)B=1;A=[1,-a]; %差分方程系数xi=filtic(B,A,ys); %由初始条件计算等效初始条件输入序列xiyn=filter(B,A,xn,xi); %调用filter解差分方程,求系统输出y(n)n=0:length(yn)-1;subplot(2,1,2);stem(n,yn, 'linewidth',2); axis([-5,15,0,8]); grid ontitle('图(b) y1(n)=0.9y1(n-1)+x(n) 初始条件y1(-1)=0 ');xlabel('n');ylabel('y(n)') 图形输出如下：-505101502468图(a) y1(n)=0.9y1(n-1)+x(n) 初始条件y1(-1)=1ny (n )-55101502468图(b) y1(n)=0.9y1(n-1)+x(n) 初始条件y1(-1)=0ny (n )2. 已知系统差分方程为： y 1(n )=0.9y 1(n -1)+x (n ) 用递推法求解系统的单位脉冲响应h (n )，要求写出h (n )的封闭公式，并打印h (n )~n 曲线。

实验一连续时间信号的时域和频域分析一. 实验目的：1. 熟悉MATLAB 软件平台。

2. 掌握MATLAB 编程方法、常用语句和可视化绘图技术。

3. 编程实现常用信号及其运算MATLAB 实现方法。

4. 编程实现常用信号的频域分析。

二. 实验原理：1、连续时间信号的描述：（1）向量表示法连续信号是指自变量的取值范围是连续的，且对于一切自变量的取值，除了有若干个不连续点之外，信号都有确定的值与之对应。

严格来说，MATLAB 并不能处理连续信号，而是用等时间间隔点的样值来近似表示连续信号。

当取样时间间隔足够小时，这些离散的样值就能较好地近似连续信号。

矩阵是MATLAB 进行数据处理的基本单元，矩阵运算是MATLAB 最重要的运算。

通常意义上的数量（也称为标量）在MATLAB 系统中是作为1×1 的矩阵来处理的，而向量实际上是仅有一行或者一列的矩阵。

通常用向量表示信号的时间取值范围，如t = -5:5，但信号x(t)、向量t 本身的下标都是从1 开始的，因此必须用一个与向量x 等长的定位时间变量t，以及向量x，才能完整地表示序列x(t)。

在MATLAB 可视化绘图中，对于以t 为自变量的连续信号，在绘图时统一用plot 函数；而对n 为自变量的离散序列，在绘图时统一用stem 函数。

（2）符号运算表示法符号对象(Symbolic Objects 不同于普通的数值计算)是Matlab 中的一种特殊数据类型，它可以用来表示符号变量、表达式以及矩阵，利用符号对象能够在不考虑符号所对应的具体数值的情况下能够进行代数分析和符号计算(symbolic math operations)，例如解代数方程、微分方程、进行矩阵运算等。

符号对象需要通过sym 或syms 函数来指定, 普通的数字转换成符号类型后也可以被作为符号对象来处理.我们可以用一个简单的例子来表明数值计算和符号计算的区别: 2/5+1/3 的结果为0.7333(double 类型数值运算), 而sym(2)/sym(5)+sym(1)/sym(3)的结果为11/15, 且这里11/15 仍然是属于sym 类型, 是符号数。