信号与系统实验 上机报告
信号与系统实验报告总结
信号与系统实验实验一常用信号的观察方波:正弦波:三角波:在观测中,虚拟示波器完全充当实际示波器的作用,在工作台上连接AD1为示波器的输入,输入方波、正弦波、三角波信号时,可在电脑上利用软件观测到相应的波形,其纵轴为幅值可通过设置实现幅值自动调节以观测到最佳大小的波形,其横轴为时间,宜可通过设置实现时间自动调节以观测到最佳宽度的波形。
实验四非正弦周期信号的分解与合成方波DC信号:DC信号几乎没有,与理论相符合,原信号没有添加偏移。
方波基波信号:基波信号为与原方波50Hz信号相对应的频率为50Hz的正弦波信号,是方波分解的一次谐波信号。
方波二次谐波信号:二次谐波信号频率为100Hz为原方波信号频率的两倍,幅值较一次谐波较为减少。
方波三次谐波信号:三次谐波信号频率为150Hz为原方波信号的三倍。
幅值较一二次谐波大为减少。
方波四次谐波信号:四次谐波信号的频率为200Hz为原方波信号的四倍。
幅值较三次谐波再次减小。
方波五次谐波信号:五次谐波频率为250Hz为原方波信号的五倍。
幅值减少到0.3以内,几乎可以忽略。
综上可知:50Hz方波可以分解为DC信号、基波信号、二次、三次、四次、五次谐波信号…,无偏移时即无DC信号,DC信号幅值为0。
分解出来的基波信号即一次谐波信号频率与原方波信号频率相同,幅值接近方波信号的幅值。
二次谐波、三次谐波、四次谐波、五次谐波依次频率分别为原方波信号的二、三、四、五倍,且幅值依次衰减,直至五次谐波信号时几乎可以忽略。
可知,方波信号可分解为多个谐波。
方波基波加三次谐波信号:基波叠加上三次谐波信号时,幅值与方波信号接近,形状还有一定差异,但已基本可以看出叠加后逼近了方波信号。
方波基波加三次谐波信号加五次谐波信号:基波信号、三次谐波信号、五次谐波信号叠加以后,比基波信号、三次谐波信号叠加后的波形更加接近方波信号。
综上所述:方波分解出来的各次谐波以及DC信号,叠加起来以后会逼近方波信号,且叠加的信号越多,越是接近方波信号。
信号与系统实验四实验报告
实验四 时域抽样与频域抽样一、实验目的加深理解连续时间信号的离散化过程中的数学概念和物理概念,掌握时域抽样定理的基本内容。
掌握由抽样序列重建原连续信号的基本原理与实现方法,理解其工程概念。
加深理解频谱离散化过程中的数学概念和物理概念,掌握频域抽样定理的基本内容。
二、 实验原理时域抽样定理给出了连续信号抽样过程中信号不失真的约束条件:对于基带信号,信号抽样频率sam f 大于等于2倍的信号最高频率m f ,即m sam f f 2≥。
时域抽样是把连续信号x (t )变成适于数字系统处理的离散信号x [k ] ;信号重建是将离散信号x [k ]转换为连续时间信号x (t )。
非周期离散信号的频谱是连续的周期谱。
计算机在分析离散信号的频谱时,必须将其连续频谱离散化。
频域抽样定理给出了连续频谱抽样过程中信号不失真的约束条件。
三.实验内容1. 为了观察连续信号时域抽样时抽样频率对抽样过程的影响,在[0,0.1]区间上以50Hz 的抽样频率对下列3个信号分别进行抽样,试画出抽样后序列的波形,并分析产生不同波形的原因,提出改进措施。
)102cos()(1t t x ⨯=π答: 函数代码为: t0 = 0:0.001:0.1;x0 =cos(2*pi*10*t0);plot(t0,x0,'r')hold onFs =50;t=0:1/Fs:0.1;x=cos(2*pi*10*t); stem(t,x); hold offtitle('连续信号及其抽样信号')函数图像为:)502cos()(2t t x ⨯=π同理,函数图像为:)0102cos()(3t t x ⨯=π同理,函数图像为:由以上的三图可知,第一个图的离散序列,基本可以显示出原来信号,可以通过低通滤波恢复,因为信号的频率为20HZ,而采样频率为50>2*20,故可以恢复,但是第二个和第三个信号的评论分别为50和100HZ,因此理论上是不能够恢复的,需要增大采样频率,解决的方案为,第二个信号的采样频率改为400HZ,而第三个的采样频率改为1000HZ,这样可以很好的采样,如下图所示:2. 产生幅度调制信号)200cos()2cos()(t t t x ππ=,推导其频率特性,确定抽样频率,并绘制波形。
信号与系统测试报告
信号与系统测试报告在进行信号与系统测试时,我们主要关注信号的特性以及系统的响应。
通过测试,我们可以验证系统的性能是否符合设计要求,以及信号是否能够正确地传输和处理。
本次测试旨在评估系统的频率响应、时域响应和稳定性等方面的表现,以确保系统能够准确、稳定地工作。
我们对系统的频率响应进行了测试。
通过输入不同频率的信号,我们可以观察系统对不同频率信号的响应情况。
测试结果显示,系统在特定频率范围内表现良好,能够准确地传输信号并保持稳定。
然而,在高频率下系统的响应有所下降,需要进一步优化以提高高频响应能力。
我们对系统的时域响应进行了测试。
通过输入不同形状的信号,如方波、正弦波等,我们可以观察系统对信号的延迟、失真等情况。
测试结果显示,系统在时域上能够准确地响应输入信号,并且延迟较小,失真程度也较低。
这表明系统具有良好的时域特性,能够满足实际应用中的需求。
我们还对系统的稳定性进行了测试。
通过输入不同幅度的信号,我们可以观察系统的稳定性和抗干扰能力。
测试结果显示,系统在输入信号幅度较小的情况下表现稳定,但在输入信号幅度较大时出现了一定程度的失真。
这提示我们需要进一步优化系统的动态范围,以提高系统的稳定性和抗干扰能力。
综合以上测试结果,我们可以得出结论,系统在频率响应、时域响应和稳定性等方面表现良好,能够满足大多数实际应用的需求。
然而,仍有一些方面需要进一步优化,如提高高频响应能力、优化动态范围等。
通过持续的测试和优化,我们相信系统将能够更好地满足用户的需求,并在实际应用中发挥更大的作用。
总的来说,信号与系统测试是确保系统正常工作的重要环节。
通过不断测试和优化,我们可以提高系统的性能和稳定性,确保系统能够准确、稳定地传输和处理信号。
希望通过本次测试报告的分享,能够帮助更多的人了解信号与系统测试的重要性,促进系统技术的进步和发展。
《信号与系统》课程实验报告
《信号与系统》课程实验报告《信号与系统》课程实验报告一图1-1 向量表示法仿真图形2.符号运算表示法若一个连续时间信号可用一个符号表达式来表示,则可用ezplot命令来画出该信号的时域波形。
上例可用下面的命令来实现(在命令窗口中输入,每行结束按回车键)。
t=-10:0.5:10;f=sym('sin((pi/4)*t)');ezplot(f,[-16,16]);仿真图形如下:图1-2 符号运算表示法仿真图形三、实验内容利用MATLAB实现信号的时域表示。
三、实验步骤该仿真提供了7种典型连续时间信号。
用鼠标点击图0-3目录界面中的“仿真一”按钮,进入图1-3。
图1-3 “信号的时域表示”仿真界面图1-3所示的是“信号的时域表示”仿真界面。
界面的主体分为两部分:1) 两个轴组成的坐标平面(横轴是时间,纵轴是信号值);2) 界面右侧的控制框。
控制框里主要有波形选择按钮和“返回目录”按钮,点击各波形选择按钮可选择波形,点击“返回目录”按钮可直接回到目录界面。
图1-4 峰值为8V,频率为0.5Hz,相位为180°的正弦信号图1-4所示的是正弦波的参数设置及显示界面。
在这个界面内提供了三个滑动条,改变滑块的位置,滑块上方实时显示滑块位置代表的数值,对应正弦波的三个参数:幅度、频率、相位;坐标平面内实时地显示随参数变化后的波形。
在七种信号中,除抽样函数信号外,对其它六种波形均提供了参数设置。
矩形波信号、指数函数信号、斜坡信号、阶跃信号、锯齿波信号和抽样函数信号的波形分别如图1-5~图1-10所示。
图1-5 峰值为8V,频率为1Hz,占空比为50%的矩形波信号图1-6 衰减指数为2的指数函数信号图1-7 斜率=1的斜坡信号图1-8 幅度为5V,滞后时间为5秒的阶跃信号图1-9 峰值为8V,频率为0.5Hz的锯齿波信号图1-10 抽样函数信号仿真途中,通过对滑动块的控制修改信号的幅度、频率、相位,观察波形的变化。
《信号与系统》实验三
三:
源程序:
(1):τ/T=1/4时的周期矩形脉冲的幅度谱和相位谱:
n=-20:20;
F=zeros(size(n));
forii=-20:20
F(ii+21)= sin(ii*pi/4)/(ii*pi+eps);
end
F(21)=1/4;
实验
内容
1.求图1所示周期信号( , )的傅里叶级数,用Matlab做出其前3、9、21、45项谐波的合成波形与原信号作比较,并做出其单边幅度谱和相位谱。
图1 周期为2的三角脉冲信号
2. 求图2所示的单个三角脉冲( )的傅里叶变换,并做出其幅度谱和相位谱。
图2 单个三角脉冲
3. 求不同占空比下周期矩形脉冲的幅度谱和相位谱,例如 、 。
y=1/4;
forn=1:m
y=y+4/(n*n*pi*pi)*(1-cos(n*pi/2)).*cos(n*pi.*t);
end
源代码:
t=-6:0.01:6;
d=-6:2:6;
fxx=pulstran(t,d,'tripuls');
f1=fourierseries(3,t);
f2=fourierseries(9,t);
n=1:10;
a=zeros(size(n));
fori=1:10
a(i)=angle(4/(i*i*pi*pi)*(1-cos(i*pi/2)))
end
n=0:pi:9*pi
stem(n,a,'fill','linewidth',2);
axis([0,9*pi,-0.2,0.2])
信号与系统实验报告
信号与系统实验报告一、实验目的(1) 理解周期信号的傅里叶分解,掌握傅里叶系数的计算方法;(2)深刻理解和掌握非周期信号的傅里叶变换及其计算方法;(3) 熟悉傅里叶变换的性质,并能应用其性质实现信号的幅度调制;(4) 理解连续时间系统的频域分析原理和方法,掌握连续系统的频率响应求解方法,并画出相应的幅频、相频响应曲线。
二、实验原理、原理图及电路图(1) 周期信号的傅里叶分解设有连续时间周期信号()f t ,它的周期为T ,角频率22fT,且满足狄里赫利条件,则该周期信号可以展开成傅里叶级数,即可表示为一系列不同频率的正弦或复指数信号之和。
傅里叶级数有三角形式和指数形式两种。
1)三角形式的傅里叶级数:01212011()cos()cos(2)sin()sin(2)2cos()sin()2n n n n a f t a t a t b t b t a a n t b n t 式中系数n a ,n b 称为傅里叶系数,可由下式求得:222222()cos(),()sin()T T T T nna f t n t dtb f t n t dtTT2)指数形式的傅里叶级数:()jn tn nf t F e式中系数n F 称为傅里叶复系数,可由下式求得:221()T jn tT nF f t edtT周期信号的傅里叶分解用Matlab进行计算时,本质上是对信号进行数值积分运算。
Matlab中进行数值积分运算的函数有quad函数和int函数。
其中int函数主要用于符号运算,而quad函数(包括quad8,quadl)可以直接对信号进行积分运算。
因此利用Matlab进行周期信号的傅里叶分解可以直接对信号进行运算,也可以采用符号运算方法。
quadl函数(quad系)的调用形式为:y=quadl(‘func’,a,b)或y=quadl(@myfun,a,b)。
其中func是一个字符串,表示被积函数的.m文件名(函数名);a、b分别表示定积分的下限和上限。
信号与系统实验报告
信号与系统实验报告
实验名称:信号与系统实验
一、实验目的:
1.了解信号与系统的基本概念
2.掌握信号的时域和频域表示方法
3.熟悉常见信号的特性及其对系统的影响
二、实验内容:
1.利用函数发生器产生不同频率的正弦信号,并通过示波器观察其时域和频域表示。
2.通过软件工具绘制不同信号的时域和频域图像。
3.利用滤波器对正弦信号进行滤波操作,并通过示波器观察滤波前后信号的变化。
三、实验结果分析:
1.通过实验仪器观察正弦信号的时域表示,可以看出信号的振幅、频率和相位信息。
2.通过实验仪器观察正弦信号的频域表示,可以看出信号的频率成分和幅度。
3.利用软件工具绘制信号的时域和频域图像,可以更直观地分析信号的特性。
4.经过滤波器处理的信号,可以通过示波器观察到滤波前后的信号波形和频谱的差异。
四、实验总结:
通过本次实验,我对信号与系统的概念有了更深入的理解,掌
握了信号的时域和频域表示方法。
通过观察实验仪器和绘制图像,我能够分析信号的特性及其对系统的影响。
此外,通过滤波器的处理,我也了解了滤波对信号的影响。
通过实验,我对信号与系统的理论知识有了更加直观的了解和应用。
信号与系统实验分析及总结
信号与系统实验分析及总结信号与系统实验是信号与系统课程中的重要环节,通过实际操控信号和系统的实验现象,深化学生对信号和系统的理论知识的理解,并培养学生的实际动手能力和解决问题的能力。
本文将对信号与系统实验进行分析和总结,探讨实验的重要性和实验中遇到的问题。
首先,信号与系统实验对于学生理解信号与系统的概念和原理起到了重要作用。
在实验中,学生可以通过操控信号源、滤波器等设备,观察信号的特征和系统的响应。
这样,学生可以将书本中的知识与实际现象相结合,更加直观地感受信号与系统的特性。
例如,在实验中,学生可以通过调节频率、幅度等参数,来观察信号的频谱特征,进而理解频域分析的概念和原理。
其次,信号与系统实验对于培养学生的实际动手能力和解决问题的能力具有重要意义。
在实验过程中,学生需要独立操作仪器设备、进行数据采集、处理和分析。
这样的实践训练,可以提高学生的实际操作技能,培养学生的实验观察能力和数据处理能力。
同时,由于实验中可能会遇到各种问题,如设备故障、数据异常等,学生需要运用所学知识和解决问题的方法来解决这些困难,培养学生的问题解决能力和创新思维能力。
然而,信号与系统实验也存在一些问题和挑战。
首先,实验设备的质量和状态可能会对实验结果产生影响。
如果设备的性能较差或者存在故障,可能会导致实验结果的不准确性,从而影响实验的有效性和可靠性。
解决这个问题的关键在于加强实验设备的维护和管理,定期检查设备状态和性能,及时更新和维修设备。
其次,实验中的数据采集和处理可能存在误差和偏差。
由于实验中操作的局限性和人为因素的影响,采集到的数据可能存在误差,这会对实验结果的分析和结论产生影响。
解决这个问题的关键在于规范实验操作流程,减少人为因素的影响,并运用合理的数据处理方法来减小误差和偏差。
综上所述,信号与系统实验是信号与系统课程中的重要环节,通过实际操控信号和系统的实验现象,深化学生对信号和系统的理论知识的理解,并培养学生的实际动手能力和解决问题的能力。
信号与系统实验报告
信号与系统实验报告目录1. 内容概要 (2)1.1 研究背景 (3)1.2 研究目的 (4)1.3 研究意义 (4)2. 实验原理 (5)2.1 信号与系统基本概念 (7)2.2 信号的分类与表示 (8)2.3 系统的分类与表示 (9)2.4 信号与系统的运算法则 (11)3. 实验内容及步骤 (12)3.1 实验一 (13)3.1.1 实验目的 (14)3.1.2 实验仪器和设备 (15)3.1.4 实验数据记录与分析 (16)3.2 实验二 (16)3.2.1 实验目的 (17)3.2.2 实验仪器和设备 (18)3.2.3 实验步骤 (19)3.2.4 实验数据记录与分析 (19)3.3 实验三 (20)3.3.1 实验目的 (21)3.3.2 实验仪器和设备 (22)3.3.3 实验步骤 (23)3.3.4 实验数据记录与分析 (24)3.4 实验四 (26)3.4.1 实验目的 (27)3.4.2 实验仪器和设备 (27)3.4.4 实验数据记录与分析 (29)4. 结果与讨论 (29)4.1 实验结果汇总 (31)4.2 结果分析与讨论 (32)4.3 结果与理论知识的对比与验证 (33)1. 内容概要本实验报告旨在总结和回顾在信号与系统课程中所进行的实验内容,通过实践操作加深对理论知识的理解和应用能力。
实验涵盖了信号分析、信号处理方法以及系统响应等多个方面。
实验一:信号的基本特性与运算。
学生掌握了信号的表示方法,包括连续时间信号和离散时间信号,以及信号的基本运算规则,如加法、减法、乘法和除法。
实验二:信号的时间域分析。
在本实验中,学生学习了信号的波形变换、信号的卷积以及信号的频谱分析等基本概念和方法,利用MATLAB工具进行了实际的信号处理。
实验三:系统的时域分析。
学生了解了线性时不变系统的动态响应特性,包括零状态响应、阶跃响应以及脉冲响应,并学会了利用MATLAB进行系统响应的计算和分析。
信号与系统实验报告
信号与系统实验报告一、信号的时域基本运算1.连续时间信号的时域基本运算两实验之一实验分析:输出信号值就等于两输入信号相加(乘)。
由于b=2,故平移量为2时,实际是右移1,符合平移性质。
两实验之二心得体会:时域中的基本运算具有连续性,当输入信号为连续时,输出信号也为连续。
平移,伸缩变化都会导致输出结果相对应的平移伸缩。
2.离散时间信号的时域基本运算两实验之一实验分析:输出信号的值是对应输入信号在每个n值所对应的运算值,当进行拉伸变化后,n值数量不会变,但范围会拉伸所输入的拉伸系数。
两实验之二心得体会:离散时间信号可以看做对连续时间信号的采样,而得到的输出信号值,也可以看成是连续信号所得之后的采样值。
二、连续信号卷积与系统的时域分析1.连续信号卷积积分两实验之一实验分析:当两相互卷积函数为冲激函数时,所卷积得到的也是一个冲激函数,且该函数的冲激t值为函数x,函数y冲激t值之和。
两实验之二心得体会:连续卷积函数每个t值所对应的卷积和可以看成其中一个在k值取得的函数与另外一个函数相乘得到的一个分量函数,并一直移动k值直至最后,最后累和出来的最终函数便是所得到的卷积函数。
3.RC电路时域积分两实验之一实验分析:全响应结果正好等于零状态响应与零输入响应之和。
两实验之二心得体会:具体学习了零状态,零输入,全响应过程的状态及变化,与之前所学的电路知识联系在一起了。
三、离散信号卷积与系统的时域分析1.离散信号卷积求和两实验之一实验分析:输出结果的n值是输入结果的k号与另一个n-k的累和两实验之二心得体会:直观地观察到卷积和的产生,可以看成连续卷积的采样形式,从这个方面去想,更能深入地理解卷积以及采样的知识。
2.离散差分方程求解两实验之一实验分析:其零状态响应序列为0 0 4 5 7.5,零输入响应序列为2 4 5 5.5 5.75,全状态响应序列为2 4 9 10.5 13.25,即全状态=零输入+零状态。
两实验之二心得体会:求差分方程时,可以根据全状态响应是由零输入输入以及零状态相加所得,分开来求,同时也加深了自己对差分方程的求解问题的理解。
《信号与系统》实验报告
信号与系统实验报告班级:姓名:信息与通信工程学院实验一 系统的卷积响应实验性质:提高性 实验级别:必做 开课单位:信息与通信工程学院 学 时:2一、实验目的:深刻理解卷积运算,利用离散卷积实现连续卷积运算;深刻理解信号与系统的关系,学习MATLAB 语言实现信号通过系统的仿真方法。
二、实验设备: 计算机,MATLAB 软件 三、实验原理: 1、 离散卷积和: 调用函数:conv ()∑∞-∞=-==i i k f i f f f conv S )()(1)2,1(为离散卷积和,其中,f1(k), f2 (k) 为离散序列,K=…-2, -1, 0 , 1, 2, …。
但是,conv 函数只给出纵轴的序列值的大小,而不能给出卷积的X 轴序号。
为得到该值,进行以下分析:对任意输入:设)(1k f 非零区间n1~n2,长度L1=n2-n1+1;)(2k f 非零区间m1~m2,长度L2=m2-m1+1。
则:)(*)()(21k f k f k s =非零区间从n1+m1开始,长度为L=L1+L2-1,所以S (K )的非零区间为:n1+m1~ n1+m1+L-1。
2、 连续卷积和离散卷积的关系:计算机本身不能直接处理连续信号,只能由离散信号进行近似: 设一系统(LTI )输入为)(t P ∆,输出为)(t h ∆,如图所示。
)t)()(t h t P ∆∆→)()(lim )(lim )(0t h t h t P t =→=∆→∆∆→∆δ若输入为f(t):∆∆-∆=≈∑∞-∞=∆∆)()()()(k t P k f t f t f k得输出:∆∆-∆=∑∞-∞=∆∆)()()(k t hk f t y k当0→∆时:⎰∑∞∞-∞-∞=∆→∆∆→∆-=∆∆-∆==ττδτd t f k t P k f t f t f k )()()()(lim)(lim )(0⎰∑∞∞-∞-∞=∆→∆∆→∆-=∆∆-∆==τττd t h f k t hk f t y t y k )()()()(lim)(lim )(0所以:∆∆-∆=-==∑⎰→∆)()(lim)()()(*)()(212121k t f k fd t f f t f t f t s τττ如果只求离散点上的f 值)(n f ∆])[()()()()(2121∑∑∞-∞=∞-∞=∆-∆∆=∆∆-∆∆=∆k k k n f k f k n f k fn f所以,可以用离散卷积和CONV ()求连续卷积,只需∆足够小以及在卷积和的基础上乘以∆。
信号与系统实验报告 连续信号的时域描述与运算
信号与系统实验报告课程名称:信号与系统实验实验项目名称:连续信号的时域描述与运算专业班级:姓名:学号:完成时间:年月日一、实验目的1.通过绘制典型信号的波形,了解这些信号的基本特征。
2.通过绘制信号运算结果的波形,了解这些信号运算对信号所起的作用。
二、实验原理1.基于MATLAB的信号描述方法如果一个信号在连续时间范围内(除有限个间断点外)有定义,则称该信号为连续时间信号,简称为连续信号。
从严格意义上讲, MATLAB数值计算的方法并不能处理连续信号,但是可利用连续信号在等时间间隔点的采样值来近似表示连续信号,即当采样间隔足够小时,这些离散采样值能够被MATLAB处理,并且能较好地近似表示连续信号。
(1)向量表示法对于连续时间信号f(t),可以定义两个行向量f和t来表示,其中向量t是形如t=t1:Δt:t2的MATLAB命令定义的时间范围向量,t1为信号起始时间,t2为终止时间,Δt为时间间隔;向量f为连续时间信号f(t)在向量t所定义的时间点上的采样值。
(2)符号运算表示法如果信号可以用一个符号表达式来表示,则可用ezplot命令绘制出信号的波形。
2.连续信号的基本运算(1)信号的相加与相乘信号的已知信号f1(t)、f2(t),信号相加和相乘记为f(t)=f1(t)+f2(t)f(t)=f1(t)·f2(t)(2)微分与积分对于连续时间信号,其微分运算是用diff函数来完成的。
其语句格式为:diff(function,’variable’,n);其中function表示需要进行求导运算的信号,或者被赋值的符号表达式;variable为求导运算的独立变量;n为求导的阶数,默认值为求一阶导数。
连续信号的积分运算用int函数来完成。
其语句格式为:int(function,’variable’,a,b);其中function表示被积信号,或者被赋值的符号表达式;variable为积分变量;a,b为积分上、下限,a和b省略时求不定积分。
华工电信学院信号与系统实验一报告参考模板
华工电信学院信号与系统实验信号与系统实验报告(一)实验项目名称:MATLAB 编程基础及典型实例 上机实验题目:信号的时域运算及MA TLAB 实现 一、实验目的学习并掌握使用MATLAB 产生基本信号、绘制信号波形、实现信号的可视化表示,为信号分析和系统设计奠定基础。
二、实验内容1. 利用Matlab 产生下列连续信号并作图。
(1) 51),1(2)(<<---=t t u t x(2) 2000,)8.0cos()1.0cos()(<<=t t t t x ππ 2. 利用Matlab 产生下列离散序列并作图。
(1) ⎩⎨⎧≤≤-=其他,055,1][k k x , 设1515-≤<k 。
(2) )]25.0cos()25.0[sin()9.0(][k k k x k ππ+=,设2020-≤<k 。
3. 已知序列]3,2,1,0,1,2;2,3,1,0,2,1[][--=-=k k x , ]21,0,1,1,1[][=-=k k h 。
(1) 计算离散序列的卷积和][][][k h k x k y *=,并绘出其波形。
(2) 计算离散序列的相关函数][][][n k y k x k R k xy +=∑∞-∞=,并绘出其波形。
(3) 序列相关与序列卷积有何关系?三、实验细节1. 利用Matlab 产生下列连续信号并作图。
(1) 51),1(2)(<<---=t t u t xt=-1:0.01:5;x=-2.*((t-1)>=0); plot(t,x);axis([-1,5,-2.2,0.2])-112345-2-1.5-1-0.5(2) 2000,)8.0cos()1.0cos()(<<=t t t t x ππ t=0:2:200;x=cos(0.1*pi*t).*cos(0.8*pi*t); plot(t,x);20406080100120140160180200-1-0.8-0.6-0.4-0.200.20.40.60.812. 利用Matlab 产生下列离散序列并作图。
信号与系统实验报告实验一 信号与系统的时域分析
实验一信号与系统的时域分析一、实验目的1、熟悉与掌握常用的用于信号与系统时域仿真分析的MA TLAB函数;2、掌握连续时间与离散时间信号的MA TLAB产生,掌握用周期延拓的方法将一个非周期信号进行周期信号延拓形成一个周期信号的MA TLAB编程;3、牢固掌握系统的单位冲激响应的概念,掌握LTI系统的卷积表达式及其物理意义,掌握卷积的计算方法、卷积的基本性质;4、掌握利用MA TLAB计算卷积的编程方法,并利用所编写的MA TLAB程序验证卷积的常用基本性质;掌握MA TLAB描述LTI系统的常用方法及有关函数,并学会利用MATLAB求解LTI系统响应,绘制相应曲线。
基本要求:掌握用MA TLAB描述连续时间信号与离散时间信号的方法,能够编写MATLAB程序,实现各种信号的时域变换与运算,并且以图形的方式再现各种信号的波形。
掌握线性时不变连续系统的时域数学模型用MA TLAB描述的方法,掌握卷积运算、线性常系数微分方程的求解编程。
二、实验原理信号(Signal)一般都就是随某一个或某几个独立变量的变化而变化的,例如,温度、压力、声音,还有股票市场的日收盘指数等,这些信号都就是随时间的变化而变化的,还有一些信号,例如在研究地球结构时,地下某处的密度就就是随着海拔高度的变化而变化的。
一幅图片中的每一个象素点的位置取决于两个坐标轴,即横轴与纵轴,因此,图像信号具有两个或两个以上的独立变量。
在《信号与系统》课程中,我们只关注这种只有一个独立变量(Independent variable)的信号,并且把这个独立变量统称为时间变量(Time variable),不管这个独立变量就是否就是时间变量。
在自然界中,大多数信号的时间变量都就是连续变化的,因此这种信号被称为连续时间信号(Continuous-Time Signals)或模拟信号(Analog Signals),例如前面提到的温度、压力与声音信号就就是连续时间信号的例子。
信号与系统 实验报告
信号与系统实验报告信号与系统实验报告一、引言信号与系统是电子信息工程领域中的重要基础课程,通过实验可以加深对于信号与系统理论的理解和掌握。
本次实验旨在通过实际操作,验证信号与系统的基本原理和性质,并对实验结果进行分析和解释。
二、实验目的本次实验的主要目的是:1. 了解信号与系统的基本概念和性质;2. 掌握信号与系统的采样、重建、滤波等基本操作;3. 验证信号与系统的时域和频域特性。
三、实验仪器与原理1. 实验仪器本次实验所需的主要仪器有:信号发生器、示波器、计算机等。
其中,信号发生器用于产生不同类型的信号,示波器用于观测信号波形,计算机用于数据处理和分析。
2. 实验原理信号与系统的基本原理包括采样定理、重建定理、线性时不变系统等。
采样定理指出,对于带限信号,为了能够完全恢复原始信号,采样频率必须大于信号最高频率的两倍。
重建定理则是指出,通过理想低通滤波器可以将采样得到的离散信号重建为连续信号。
四、实验步骤与结果1. 采样与重建实验首先,将信号发生器输出的正弦信号连接到示波器上,观察信号的波形。
然后,将示波器的输出信号连接到计算机上,进行采样,并通过计算机对采样信号进行重建。
最后,将重建得到的信号与原始信号进行对比,分析重建误差。
实验结果显示,当采样频率满足采样定理时,重建误差较小,重建信号与原始信号基本一致。
而当采样频率不满足采样定理时,重建信号存在失真和混叠现象。
2. 系统特性实验接下来,通过调节示波器和信号发生器的参数,观察不同系统对信号的影响。
例如,将示波器设置为高通滤波器,通过改变截止频率,观察信号的低频衰减情况。
同样地,将示波器设置为低通滤波器,观察信号的高频衰减情况。
实验结果表明,不同系统对信号的频率特性有着明显的影响。
高通滤波器会使低频信号衰减,而低通滤波器则会使高频信号衰减。
通过调节滤波器的参数,可以实现对信号频率的选择性衰减。
五、实验分析与讨论通过本次实验,我们对信号与系统的基本原理和性质有了更深入的理解。
《信号与系统》实验报告
《信号与系统》实验报告目录一、实验概述 (2)1. 实验目的 (2)2. 实验原理 (3)3. 实验设备与工具 (4)二、实验内容与步骤 (5)1. 实验一 (6)1.1 实验目的 (7)1.2 实验原理 (7)1.3 实验内容与步骤 (8)1.4 实验结果与分析 (9)2. 实验二 (10)2.1 实验目的 (12)2.2 实验原理 (12)2.3 实验内容与步骤 (13)2.4 实验结果与分析 (14)3. 实验三 (15)3.1 实验目的 (16)3.2 实验原理 (16)3.3 实验内容与步骤 (17)3.4 实验结果与分析 (19)4. 实验四 (20)4.1 实验目的 (20)4.2 实验原理 (21)4.3 实验内容与步骤 (22)4.4 实验结果与分析 (22)三、实验总结与体会 (24)1. 实验成果总结 (25)2. 实验中的问题与解决方法 (26)3. 对信号与系统课程的理解与认识 (27)4. 对未来学习与研究的展望 (28)一、实验概述本实验主要围绕信号与系统的相关知识展开,旨在帮助学生更好地理解信号与系统的基本概念、性质和应用。
通过本实验,学生将能够掌握信号与系统的基本操作,如傅里叶变换、拉普拉斯变换等,并能够运用这些方法分析和处理实际问题。
本实验还将培养学生的动手能力和团队协作能力,使学生能够在实际工程中灵活运用所学知识。
本实验共分为五个子实验,分别是:信号的基本属性测量、信号的频谱分析、信号的时域分析、信号的频域分析以及信号的采样与重构。
每个子实验都有明确的目标和要求,学生需要根据实验要求完成相应的实验内容,并撰写实验报告。
在实验过程中,学生将通过理论学习和实际操作相结合的方式,逐步深入了解信号与系统的知识体系,提高自己的综合素质。
1. 实验目的本次实验旨在通过实践操作,使学生深入理解信号与系统的基本原理和概念。
通过具体的实验操作和数据分析,掌握信号与系统分析的基本方法,提高解决实际问题的能力。
信号与系统MATLAB实验
2016-2017学年第一学期信号与系统实验报告班级:姓名:学号:成绩:指导教师:实验一 常见信号的MATLAB 表示及运算一.实验目的1.熟悉常见信号的意义、特性及波形2.学会使用MATLAB 表示信号的方法并绘制信号波形 3. 掌握使用MATLAB 进行信号基本运算的指令 4. 熟悉用MATLAB 实现卷积积分的方法二.实验原理信号一般是随时间而变化的某些物理量;按照自变量的取值是否连续,信号分为连续时间信号和离散时间信号,一般用()f t 和()f k 来表示;若对信号进行时域分析,就需要绘制其波形,如果信号比较复杂,则手工绘制波形就变得很困难,且难以精确;MATLAB 强大的图形处理功能及符号运算功能,为实现信号的可视化及其时域分析提供了强有力的工具;根据MATLAB 的数值计算功能和符号运算功能,在MATLAB 中,信号有两种表示方法,一种是用向量来表示,另一种则是用符号运算的方法;在采用适当的MATLAB 语句表示出信号后,就可以利用MATLAB 中的绘图命令绘制出直观的信号波形了;下面分别介绍连续时间信号和离散时间信号的MATLAB 表示及其波形绘制方法;1.连续时间信号所谓连续时间信号,是指其自变量的取值是连续的,并且除了若干不连续的点外,对于一切自变量的取值,信号都有确定的值与之对应;从严格意义上讲,MATLAB 并不能处理连续信号;在MATLAB 中,是用连续信号在等时间间隔点上的样值来近似表示的,当取样时间间隔足够小时,这些离散的样值就能较好地近似出连续信号;在MATLAB 中连续信号可用向量或符号运算功能来表示;⑴ 向量表示法对于连续时间信号()f t ,可以用两个行向量f 和t 来表示,其中向量t 是用形如12::t t p t 的命令定义的时间范围向量,其中,1t 为信号起始时间,2t 为终止时间,p 为时间间隔;向量f 为连续信号()f t 在向量t 所定义的时间点上的样值; 说明: plot 是常用的绘制连续信号波形的函数;严格说来,MATLAB 不能表示连续信号,所以,在用plot 命令绘制波形时,要对自变量t 进行取值,MATLAB 会分别计算对应点上的函数值,然后将各个数据点通过折线连接起来绘制图形,从而形成连续的曲线;因此,绘制的只是近似波形,而且,其精度取决于t 的取样间隔;t 的取样间隔越小,即点与点之间的距离越小,则近似程度越好,曲线越光滑;例如:图1-1是在取样间隔为p=时绘制的波形,而图1-2是在取样间隔p=时绘制的波形,两相对照,可以看出图1-2要比图1-1光滑得多;在上面的f=sint . /t 语句中,必须用点除符号,以表示是两个函数对应点上的值相除;⑵ 符号运算表示法如果一个信号或函数可以用符号表达式来表示,那么我们就可以用前面介绍的符号函数专用绘图命令ezplot 等函数来绘出信号的波形; ⑶ 常见信号的MATLAB 表示对于普通的信号,应用以上介绍的两种方法即可完成计算函数值或绘制波形,但是对于一些比较特殊的信号,比如单位阶跃信号t 、符号函数sgnt 等,在MATLAB 中这些信号都有专门的表示方法; 单位阶跃信号单位阶跃信号的定义为:10()0t t t ε>⎧=⎨<⎩,单位阶跃信号是信号分析的基本信号之一,在信号与系统分析中有着非常重要的作用,通常,我们用它来表示信号的定义域,简化信号的时域表示形式;例如:可以用两个不同延时的单位阶跃信号来表示一个矩形门信号,即:2()(1)(1)G t t t εε=+--在MATLAB 中,可通过多种方法得到单位阶跃信号,下面分别介绍; 方法一: 调用Heavisidet 函数在MATLAB 的Symbolic Math Toolbox 中,有专门用于表示单位阶跃信号的函数,即Heavisidet 函数,用它即可方便地表示出单位阶跃信号以及延时的单位阶跃信号,并且可以方便地参加有关的各种运算过程;首先定义函数Heavisidet 的m 函数文件,该文件名应与函数名同名即;%定义函数文件,函数名为Heaviside,输入变量为x,输出变量为y function y= Heavisidety=t>0; %定义函数体,即函数所执行指令%此处定义t>0时y=1,t<=0时y=0,注意与实际阶跃信号定义的区别; 方法二:数值计算法在MATLAB 中,有一个专门用于表示单位阶跃信号的函数,即stepfun 函数,它是用数值计算法表示的单位阶跃函数()t ε;其调用格式为: stepfunt,t0其中,t 是以向量形式表示的变量,t0表示信号发生突变的时刻,在t0以前,函数值小于零,t0以后函数值大于零;有趣的是它同时还可以表示单位阶跃序列()k ε,这只要将自变量以及取样间隔设定为整数即可;有关单位阶跃序列()k ε的表示方法,我们后面有专门论述,下面通过一个例子来说明如何调用stepfun 函数来表示单位阶跃函数; 符号函数符号函数的定义为:10sgn()1t t t >⎧=⎨-<⎩在MATLAB 中有专门用于表示符号函数的函数sign ,由于单位阶跃信号 t和符号函数两者之间存在以下关系:1122()sgn()t t ε=+,因此,利用这个函数就可以很容易地生成单位阶跃信号;下面举个例子来说明如何利用sign 函数生成单位阶跃信号,并同时绘制其波形; 2.离散时间信号离散时间信号又叫离散时间序列,一般用()f k 表示,其中变量k 为整数,代表离散的采样时间点采样次数;三.实验内容1.分别用MATLAB 的向量表示法和符号运算功能,表示并绘出下列连续时间信号的波形:2[]()cos()()(4)2tf t t t πεε=--syms t;f=sym'costpit/2heavisidet-heavisidet-4'; ezplotf,-2,8; 423()(2)f t t t ε=+syms t;f=sym'2/3theavisidet+2'; ezplotf,-4,8;2.分别用MATLAB 表示并绘出下列离散时间信号的波形: 2[]()()(8)f t k k k εε=-- t=0:8; t1=-10:15;f=zeros1,10,t,zeros1,7; stemt1,faxis-10,15,0,10 4 ()(2)f k k ε=-+t=-20:10;f=ones1,23,zeros1,8; stemt,f stemt,f3.已知信号f t 的波形如下图所示,试用MATLAB 绘出满足下列要求的信号波形; 2(2)f t -t=-1::4; t0=0; t1=1; t2=2;ut=2stepfunt,t0-2stepfunt,t1+stepfunt,t1-stepfunt,t2;plott-2,ut axis-3,2,0,3 4(0.51)f t +t=-1::4;t0=0;t1=1; t2=2;ut=2stepfunt,t0-2stepfunt,t1+stepfunt,t1-stepfunt,t2; plott+1,ut axis-1,2,0,44.已知两信号1()(1)()f t t t εε=+-,2()()(1)f t t t εε=--,求卷积积分12()()()g t f t f t =*,并与例题比较;程序清单:t1=-1::0; t2=0::1; t3=3::5;f1=onessizet1; f2=onessizet2; g=convf1,f2; plott3,g 信号波形:5.已知两信号1()()f t t t ε=,20()()0t tt te t f t t eε-≥⎧=⎨<⎩,求卷积积分12()()()g t f t f t =*;程序代码:t1=0::5;t2=-5::5; t3=-5::10; f1=t1;f2=expt2.t2<0+t2.exp-t2.t2>=0; g=convf1,f2; plott3,g;运行结果截图:6.已知{}{}12()1,1,1,2,()1,2,3,4,5f k f k ==,求两序列的卷积和 ; 程序清单:f1=1,1,1,2,0; f2=1,2,3,4,5; f=convf1,f2; x=0:8;stemx,f,'filled' 信号波形:实验二 LTI 系统的响应一、实验目的1.熟悉连续时间系统的单位冲激响应、阶跃响应的意义及求解方法2.熟悉连续离散时间系统在任意信号激励下响应的求解方法3.熟悉应用MATLAB 实现求解系统响应的方法二、实验原理1.连续时间系统对于连续的LTI 系统,当系统输入为ft ,输出为yt ,则输入与输出之间满足如下的线性常系数微分方程:()()0()()nmi j i j i j a y t b f t ===∑∑,当系统输入为单位冲激信号δt 时产生的零状态响应称为系统的单位冲激响应,用ht 表示;若输入为单位阶跃信号εt 时,系统产生的零状态响应则称为系统的单位阶跃响应,记为gt ,如下图所示;系统的单位冲激响应ht 包含了系统的固有特性,它是由系统本身的结构及参数所决定的,与系统的输入无关;我们只要知道了系统的冲激响应,即可求得系统在不同激励下产生的响应;因此,求解系统的冲激响应ht 对我们进行连续系统的分析具有非常重要的意义;在MATLAB 中有专门用于求解连续系统冲激响应和阶跃响应, 并绘制其时域波形的函数impulse 和step ;如果系统输入为ft ,冲激响应为ht ,系统的零状态响应为yt ,则有:()()()y t h t f t =*;若已知系统的输入信号及初始状态,我们便可以用微分方程的经典时域求解方法,求出系统的响应;但是对于高阶系统,手工计算这一问题的过程非常困难和繁琐;在MATLAB 中,应用lsim 函数很容易就能对上述微分方程所描述的系统的响应进行仿真,求出系统在任意激励信号作用下的响应;lsim 函数不仅能够求出连续系统在指定的任意时间范围内系统响应的数值解,而且还能同时绘制出系统响应的时域波形图;2.离散时间系统LTI 离散系统中,其输入和输出的关系由差分方程描述:00()()n mi ji j a y k i bf k j ==+=+∑∑前向差分方程()()nmi ji j a y k i bf k n j ==-=-+∑∑ 后向差分方程当系统的输入为单位序列δk 时产生的零状态响应称为系统的单位函数响应,用hk 表示;当输入为 εk 时产生的零状态响应称为系统的单位阶跃应,记为:gk ,如下图所示;如果系统输入为ek ,冲激响应为hk ,系统的零状态响应为y k ,则有:()()()y k h k f k =*;与连续系统的单位冲激响应ht 相类似,离散系统的单位函数响应hk 也包含了系统的固有特性,与输入序列无关;我们只要知道了系统的单位函数响应,即可求得系统在不同激励信号作用下产生的响应;因此,求解系统的单位函数响应hk 对我们进行离散系统的分析也同样具有非常重要的意义;MATLAB 中为用户提供了专门用于求解离散系统单位函数响应, 并绘制其时域波形的函数impz ;同样也提供了求离散系统响应的专用函数filter ,该函数能求出由差分方程所描述的离散系统在指定时间范围内的输入序列作用时,产生的响应序列的数值解;当系统初值不为零时,可以使用dlsim 函数求出离散系统的全响应,其调用方法与前面连续系统的lsim 函数相似;另外,求解离散系统阶跃响应可以通过如下两种方法实现:一种是直接调用专用函数dstep ,其调用方法与求解连续系统阶跃响应的专用函数step 的调用方法相似;另一种方法是利用求解离散系统零状态响应的专用函数filter ,只要将其中的激励信号看成是单位阶跃信号εk 即可;三、实验内容1. 已知描述系统的微分方程和激励信号et 分别如下,试用解析方法求系统的单位冲激响应ht 和零状态响应rt ,并用MATLAB 绘出系统单位冲激响应和系统零状态响应的波形,验证结果是否相同;①''()4'()4()'()3()y t y t y t f t f t ++=+;()()t f t e t ε-=程序清单:a=1 4 4;b=1 3; impulseb,a,10 p=;t=0:p:10; x=exp-1t;y=filterb,a,xsubplot2,1,1,impulseb,a,10 subplot2,1,2,lsimb,a,x,t②''()2'()26()'()y t y t y t f t ++=;()()f t t ε= 单位冲激响应程序代码: a=1 2 26;b=1;subplot2,1,1, impulseb,a,4 subplot2,1,2, stepb,a,4 运行结果截图:零状态响应程序代码: a=1 2 26;b=1; p1=;t1=0:p1:10; x1=t1;lsimb,a,x1,t1 运行结果截图:③''()4'()3()()y t y t y t f t ++=;2()()t f t e t ε-=a=1 4 3;b=1; p=;t=0:p:10; x=exp-2ty=filterb,a,xsubplot2,1,1,impulseb,a,10 subplot2,1,2,lsimb,a,x,t④如下图所示的电路中,已知1234()R R R ===Ω,121()L L H ==,且两电感上初始电流分别为12(0)2(),(0)0()i A i A ==,如果以电阻3R 上电压()y t 作为系统输出,请求出系统在激励()12()f t t ε=v 作用下的全响应;程序清单:A=-8 4;4 -8;B=1;0;C=-4 4;D=0; x0=2;0; t=0::10;E=12.onessizet;r,x=lsimA,B,C,D,E,t,x0; plott,r信号波形:阶跃响应程序代码: a=1,-5/6,1/6;b=1,0,-1; k=0:20;x=heavisidek; y=filterb,a,xsubplot2,1,1,stemk,x title'输入序列'subplot2,1,2,stemk,y title'输出序列' 运行结果截图:④一带通滤波器可由下列差分方程描述:()0.81(2)()(2)y k y k f k f k +-=--,其中()f k 为系统输入, ()y k 为系统输出;请求出当激励[]()1010cos(/2)10cos()()f k kn kn k ε=++选取适当的n 值时滤波器的稳态输出;a=1,0,81/100; b=1,0,-1; k=0:20;x=10+10.cos1/2.k+10.cosk; y=filterb,a,xsubplot3,1,1,impzb,a,0:20, subplot3,1,2,dstepb,a,0:20, subplot3,1,3,stemk,y实验三连续时间信号的频域分析一、实验目的1.熟悉傅里叶变换的性质 2.熟悉常见信号的傅里叶变换3.了解傅里叶变换的MATLAB 实现方法二、实验原理傅里叶变换是信号分析 的最重要的内容之一;从已知信号()f t 求出相应的频谱函数()F j ω的数学表示为:()F j ω()j t f t e dt ω∞--∞=⎰()f t 的傅里叶变换存在的充分条件是()f t 在无限区间内绝对可积,即()f t 满足下式: ()f t dt ∞-∞<∞⎰但上式并非傅里叶变换存在的必要条件;在引入广义函数概念之后,使一些不满足绝对可积条件的函数也能进行傅里叶变换;傅里叶反变换的定义为:1()()2j tf t F j e d ωωωπ∞-∞=⎰; 在这一部分的学习中,大家都体会到了这种数学运算的麻烦;在MATLAB 语言中有专门对信号进行正反傅里叶变换的语句,使得傅里叶变换很容易在MATLAB 中实现;在MATLAB 中实现傅里叶变换的方法有两种,一种是利用MATLAB 中的Symbolic Math Toolbox 提供的专用函数直接求解函数的傅里叶变换和傅里叶反变换,另一种是傅里叶变换的数值计算实现法;1.直接调用专用函数法①在MATLAB 中实现傅里叶变换的函数为: F=fourier f对ft 进行傅里叶变换,其结果为FwF =fourierf,v对ft 进行傅里叶变换,其结果为FvF=fourier f,u,v对fu 进行傅里叶变换,其结果为Fv②傅里叶反变换f=ifourier F 对Fw 进行傅里叶反变换,其结果为fx f=ifourierF,U对Fw 进行傅里叶反变换,其结果为fu f=ifourier F,v,u 对Fv 进行傅里叶反变换,其结果为fu由于MATLAB 中函数类型非常丰富,要想了解函数的意义和用法,可以用mhelp 命令;如在命令窗口键入:mhelp fourier 回车,则会得到fourier 的意义和用法; 注意:1在调用函数fourier 及ifourier 之前,要用syms 命令对所有需要用到的变量如t,u,v,w 等进行说明,即要将这些变量说明成符号变量;对fourier 中的f 及ifourier 中的F 也要用符号定义符sym 将其说明为符号表达式;2采用fourier 及fourier 得到的返回函数,仍然为符号表达式;在对其作图时要用ezplot 函数,而不能用plot 函数;3fourier 及fourier 函数的应用有很多局限性,如果在返回函数中含有δω等函数,则ezplot 函数也无法作出图来;另外,在用fourier 函数对某些信号进行变换时,其返回函数如果包含一些不能直接表达的式子,则此时当然也就无法作图了;这是fourier 函数的一个局限;另一个局限是在很多场合,尽管原时间信号ft 是连续的,但却不能表示成符号表达式,此时只能应用下面介绍的数值计算法来进行傅氏变换了,当然,大多数情况下,用数值计算法所求的频谱函数只是一种近似值;三、 实验内容1.编程实现求下列信号的幅度频谱1求出1()(21)(21)f t t t εε=+--的频谱函数F1jω,请将它与上面门宽为2的门函数()(1)(1)f t t t εε=+--的频谱进行比较,观察两者的特点,说明两者的关系;f1t 函数程序代码: syms t w;Gt=sym'Heaviside2t+1-Heaviside2t-1'; Fw=fourierGt,t,w;FFw=maple'convert',Fw,'piecewise'; FFP=absFFw;ezplotFFP,-10pi 10pi;grid; axis-10pi 10pi 0 ; f1t 函数图像: 门函数程序代码: syms t w;Gt=sym'Heavisidet+1-Heavisidet-1'; Fw=fourierGt,t,w;FFw=maple'convert',Fw,'piecewise'; FFP=absFFw;ezplotFFP,-10pi 10pi;grid; axis-10pi 10pi 0 ; 门函数图像:2 三角脉冲21||||1()0||1t t f t t -≤⎧=⎨>⎩程序清单:f2t 函数程序代码:syms t wGt=sym'1+tHeavisidet+1-Heavisidet+1-tHeavisidet-Heavisidet-1'; Fw=fourierGt,t,w;FFw=maple'convert',Fw,'piecewise';FFP=absFFw;ezplotFFP,-10pi 10pi;grid;axis-10pi 10pi 0 ;f2t 函数图像:3 单边指数信号 3()()t f t e t ε-=程序清单:syms t wGt=sym'exp-1theavisidet';Fw=fourierGt,t,w;FFw=maple'convert',Fw,'piecewise';FFP=absFFw; ezplotFFP,-7pi 7pi;grid;axis-7pi 7pi 0信号波形:4 高斯信号23()t f t e -=程序清单:syms t w;Gt=exp-t.^2;Fw=fourierGt,t,wFFP=absFwezplotFw,-30 30;grid;axis-30 30 0 2信号波形:2.利用ifourier 函数求下列频谱函数的傅氏反变换 122()16F j j ωωω=-+ 程序清单:syms t wFw=sym '-j2w/16+w^2';ft=ifourier Fwft=ifourierFw,w,t;运行结果:ft=-jexp-4absxsignx1i222()58()()65j jF jj jωωωωω+-=++syms t wFw=sym'jw^2+5jw-8/jw^2+6jw+5';ft=ifourierFwft=ifourierFw,w,t;运行结果:ft =2pidiracx + piexp-x1i/jsignimag1/j3i/j -piexp-x5i/jsignimag1/j2i/j - piexp-x1i/jsignx3i/j +piexp-x5i/jsignx2i/j/2pi实验四离散信号与系统的时域分析一、实验目的1.学会用MATLAB表示常用离散信号的方法;2.学会用MATLAB实现离散信号卷积的方法;3.学会用MATLAB求解离散系统的单位响应;4.学会用MATLAB求解离散系统的零状态响应;二、实验原理1.离散信号的MATLAB表示表示离散时间信号fk需要两个行向量,一个是表示序号k= ,一个是表示相应函数值f= ,画图命令是stem;2.离散信号的卷积和两个有限长序列f1,f2卷积可调用MATLAB函数conv,调用格式是f=convf1,f2, f是卷积结果,但不显示时间序号,可自编一个函数dconv给出f 和k,并画图;3.离散系统的单位响应MATLAB提供画系统单位响应函数impz,调用格式是impzb,a 式中b和a是表示离散系统的行向量;impzb,a,n 式中b和a是表示离散系统的行向量,时间范围是0~n;impzb,a,n1,n2 时间范围是n1~n2 ;y=impzb,a,n1,n2 由y给出数值序列;4.离散系统的零状态响应MATLAB提供求离散系统零状态响应数值解函数filter,调用格式为filterb,a,x,式中b和a是表示离散系统的向量,x是输入序列非零样值点行向量,输出向量序号同x一样;三、上机实验内容1.验证实验原理中程序离散信号的MATLAB表示例2-1正弦序列信号 正弦序列信号可直接调用MATLAB 函数cos,例)cos(ϕω+k ,当ωπ/2是整数或分数时,才是周期信号;画)8/cos(ϕπ+k ,)2cos(k 波形程序是:k=0:40;subplot2,1,1stemk,coskpi/8,'filled'title'coskpi/8'subplot2,1,2stemk,cos2k,'filled'title'cos2k'2.已知)2(2)1(3)()2()1(2)(2-+-+=-+--k f k f k f k y k y k y ,画单位响应波形;a=2,-2,1;b=1,3,2;impzb,aimpzb,a,60impzb,a,-10:403.已知)()2(25.0)1()(k f k y k y k y =-+-+,输入)()(k t f ε=,画输出波形,范围0~15;a=1 1 ;b=1 ;t=0:15;x=t;y=filterb,a,xsubplot2,1,1stemt,xtitle'输入序列'subplot2,1,2stemt,ytitle'响应序列'实验五 连续信号与系统的S 域分析一、实验目的1. 熟悉拉普拉斯变换的原理及性质2. 熟悉常见信号的拉氏变换3. 了解正/反拉氏变换的MATLAB 实现方法和利用MATLAB 绘制三维曲面图的方法4. 了解信号的零极点分布对信号拉氏变换曲面图的影响及续信号的拉氏变换与傅氏变换的关系二、实验原理拉普拉斯变换是分析连续时间信号的重要手段;对于当t ∞时信号的幅值不衰减的时间信号,即在ft 不满足绝对可积的条件时,其傅里叶变换可能不存在,但此时可以用拉氏变换法来分析它们;连续时间信号ft 的单边拉普拉斯变换Fs 的定义为: 拉氏反变换的定义为:1()()2j st j f t F s e ds j σωσωπ+-=⎰显然,上式中Fs 是复变量s 的复变函数,为了便于理解和分析Fs 随s 的变化规律,我们将Fs 写成模及相位的形式:()()()j s F s F s e ϕ=;其中,|Fs|为复信号Fs 的模,而()s ϕ为Fs 的相位;由于复变量s=σ+jω,如果以σ为横坐标实轴,jω为纵坐标虚轴,这样,复变量s 就成为一个复平面,我们称之为s 平面;从三维几何空间的角度来看,|()|F s 和()s ϕ分别对应着复平面上的两个曲面,如果绘出它们的三维曲面图,就可以直观地分析连续信号的拉氏变换Fs 随复变量s 的变化情况,在MATLAB 语言中有专门对信号进行正反拉氏变换的函数,并且利用 MATLAB 的三维绘图功能很容易画出漂亮的三维曲面图;①在MATLAB 中实现拉氏变换的函数为:F=laplace f 对ft 进行拉氏变换,其结果为FsF=laplace f,v对ft 进行拉氏变换,其结果为FvF=laplace f,u,v对fu 进行拉氏变换,其结果为Fv②拉氏反变换f=ilaplace F对Fs 进行拉氏反变换,其结果为ftf=ilaplaceF,u 对Fw 进行拉氏反变换,其结果为fuf=ilaplaceF,v,u 对Fv 进行拉氏反变换,其结果为fu注意: 在调用函数laplace 及ilaplace 之前,要用syms 命令对所有需要用到的变量如t,u,v,w 等进行说明,即要将这些变量说明成符号变量;对laplace 中的f 及ilaplace 中的F 也要用符号定义符sym 将其说明为符号表达式;三、实验内容1.求出下列函数的拉氏变换式,并用MATLAB 绘制拉氏变换在s 平面的三维曲面图① 3()2()5()t t f t e t e t εε--=+函数程序代码:syms t sft=sym'2exp-tHeavisidet+5exp-3tHeavisidet';Fs=laplaceft运算结果:绘制三维曲面图的程序代码:syms x y ss=x+iy;FFs=2/s+1+5/s+3;FFss=absFFs;ezmeshFFss;ezsurfFFss;colormaphsv;运算结果截图:② ()()(2)f t t t εε=--函数程序代码:syms t sft=sym'Heavisidet-Heavisidet-2';Fs=laplaceft运算结果:绘制三维曲面图的程序代码:syms x y ss=x+iy;FFs=1/s-exp-2s/s;FFss=absFFs;ezmeshFFss;ezsurfFFss;colormaphsv;运算结果截图:③ 3()sin()()t f t e t t ε-=函数程序代码:syms t sft=sym'exp-3tsintHeavisidet';Fs=laplaceft运算结果:绘制三维曲面图的程序代码:syms x y ss=x+iy;FFs=1/s+3^2+1;FFss=absFFs;ezmeshFFss;ezsurfFFss;colormaphsv;运算结果截图:④ []()sin()()(2)f t t t t πεε=--函数程序代码:syms t sft=sym'sinpitHeavisidet-Heavisidet-2';Fs=laplaceft运算结果:绘制三维曲面图的程序代码:syms x y ss=x+iy;FFs= pi/s^2+pi^21/s-exp-2s/s;FFss=absFFs;ezmeshFFss;ezsurfFFss;colormaphsv;运算结果截图:2. 已知信号的拉氏变换如下,请用MATLAB 画出其三维曲面图,观察其图形特点,说出函数零极点位置与其对应曲面图的关系,并且求出它们所对应的原时间函数f t ①22(3)(3)()(5)(16)s s F s s s -+=-+ 函数程序代码:syms x y ss=x+iy;FFs=2s-3s+3/s-5s^2+16;FFss=absFFs;ezmeshFFss;ezsurfFFss;colormaphsv;运行结果截图:求原函数的程序代码:syms t sFs =sym'2s-3s+3/s-5s^2+16'ft=ilaplaceFs原函数:ft = 50cos4t/41 + 32exp5t/41 + 125sin4t/82 ②(1)(3)()(2)(5)s s F s s s s ++=++ 函数程序代码:syms x y ss=x+iy;FFs=s+1s+3/ss+2s+5;FFss=absFFs;ezmeshFFss;ezsurfFFss;colormaphsv;运行结果截图:求原函数的程序代码:syms t sFs =sym's+1s+3/ss+2s+5'ft=ilaplaceFs原函数:ft = exp-2t/6 + 8exp-5t/15 + 3/103. 已知连续时间信号[]()s(2)()(4)f t co t t t πεε=--,请分别求出该信号的拉氏变换()F s 及其傅里叶变换()F j ω,并用MATLAB 绘出()F s 的曲面图及振幅频谱()F j ω的波形,观察()F s 的曲面图在虚轴上的剖面图,并将它与信号的振幅频谱曲线进行比较,分析两者的对应关系;1拉氏变换:程序代码:syms t sft=sym'cos2pitHeavisidet-Heavisidet-4';Fs=laplaceft运算结果:Fs=laplaceHeavisidet, t, s-pi2i/2 + laplaceHeavisidet, t, s+pi2i/2 - laplaceHeavisidet - 4, t, s - pi2i/2 - laplaceHeavisidet - 4, t, s + pi2i/22傅里叶变换:程序代码:syms t wGt=sym'cos2pitHeavisidet-Heavisidet-4';Fw=fourierGt运算结果:Fw = fouriercos2pitHeavisidet, t, w - fourierHeavisidet - 4cos2pit, t, w四、总结报告由于平时都是在上课过程中学习理论知识,而这次实验是在理论知识的基础上来进行实验操作,但并是不全是上课时学习的理论知识,也存在许多的新知识;所以对于这次把上课时的理论知识并结合新知识一起应用于实践操作来说是有点困难的; 信号与系统的实验不同于大物实验,一开始说可以多人合作完成的实验,到最后是一个人单独完成;在为数不多的四次实验中,我深深感受到了团队合作在实验中的重要性;在自己对自己写出的代码,运行出现错误的时候,两个人或者多个人对实验的共同理解是实验高效、误差小完成的基础;参考文献:MATLAB应用大全MATLAB无师自通。
信号与系统的实验报告(2)
信号与系统实验报告——连续时间系统的复频域分析班级:05911101学号:**********姓名:***实验五连续时间系统的复频域分析——1120111487 信息工程(实验班)蒋志科一、实验目的①掌握拉普拉斯变换及其反变换的定义,并掌握MA TLAB 实现方法 ②学习和掌握连续时间系统系统函数的定义及其复频域分析方法③掌握系统零极点的定义,加深理解系统零极点分布与系统特性的关系。
二、实验原理与方法 1、拉普拉斯变换连续时间信号x(t)的拉普拉斯变换定义为:X s =x (t )e −st dt +∞−∞拉普拉斯反变换为:x t =12πj X (s )e st ds σ+j ∞σ−j ∞在MA TLAB 中可以采用符号数学工具箱中的laplace 函数和ilaplace 函数进行拉氏变换和拉氏反变换。
L=laplace(F)符号表达式F 的拉氏变换,F 中时间变量为t ,返回变量为s 的结果表达式。
L=laplace(F,t)用t 替换结果中的变量s 。
F=ilaplace(L)以s 为变量的符号表达式L 的拉氏反变换,返回时间变量t 的结果表达式。
F=ilaplace(L,x)用x 替换结果中的变量t 。
2、连续时间系统的系统函数连续时间系统的系统函数是系统单位冲激响应的拉氏变换H s =ℎ(t )e −st dt +∞−∞此外,连续时间系统的系统函数还可以由系统输入和输出信号的拉氏变换之比得到H s =Y(s)/X(s) 单位冲激响应h(t)反映了系统的固有性质,而H(s)从复频域反映了系统的固有性质。
对于H(s)描述的连续时间系统,其系统函数s 的有理函数H s =b M s M +b M−1s M−1+⋯+b 0a n s n +a n −1s M−1+⋯+a 03、连续时间系统的零极点分析系统的零点指使式H s 的分子多项式为零的点,极点指使分母多项式为零的点,零点使系统的值为零,极点使系统函数的值无穷大。
信号与系统实验报告-1(常用信号的分类与观察)
信号与系统实验报告-1(常用信号的分类与观察)实验一:信号的时域分析一、实验目的1.观察常用信号的波形特点及产生方法2.学会使用示波器对常用波形参数的测量二、实验仪器1.信号与系统试验箱一台(型号ZH5004)2.40MHz双踪示波器一台3.DDS信号源一台三、实验原理对于一个系统特性的研究,其中重要的一个方面是研究它的输入输出关系,即在一特定的输入信号下,系统对应的输出响应信号。
因而对信号的研究是对系统研究的出发点,是对系统特性观察的基本手段与方法。
在本实验中,将对常用信号和特性进行分析、研究。
信号可以表示为一个或多个变量的函数,在这里仅对一维信号进行研究,自变量为时间。
常用信号有:指数信号、正弦信号、指数衰减正弦信号、复指数信号、Sa(t)信号、钟形信号、脉冲信号等。
1、信号:指数信号可表示为f(t)=Ke at。
对于不同的a取值,其波形表现为不同的形式,如下图所示:图1―1 指数信号2、信号:其表达式为f(t)=Ksin(ωt+θ),其信号的参数:振幅K、角频率ω、与初始相位θ。
其波形如下图所示:图1-2 正弦信号3、指数衰减正弦信号:其表达式为其波形如下图:图1-3 指数衰减正弦信号4、Sa(t)信号:其表达式为:。
Sa(t)是一个偶函数,t= ±π,±2π,…,±nπ时,函数值为零。
该函数在很多应用场合具有独特的运用。
其信号如下图所示:图1-4 Sa(t)信号5、钟形信号(高斯函数):其表达式为:其信号如下图所示:图1-5 钟形信号6、脉冲信号:其表达式为f(t)=u(t)-u(t-T),其中u(t)为单位阶跃函数。
其信号如下图所示:f(t)…………0 t图1-6 脉冲信号7、方波信号:信号为周期为T,前T/2期间信号为正电平信号,后T/2期间信号为负电平信号,其信号如下图所示U(t)…………0 t图1-7 方波信号四、实验内容及主要步骤下列实验中信号产生器的工作模式为111、指数信号观察通过信号选择键1,设置A组输出为指数信号(此时信号输出指示灯为000000)。