验证可控硅通过15KHz信号实验报告（5篇）

第一篇：验证可控硅通过15KHz信号实验报告

验证可控硅可通过15KHz信号实验报告

实验目的：验证普通可控硅器件是否可控制15KHz频率信号通过。

实验器材：可控硅JCT2515、MTC55-12各一只，直流电源一台，直流开关电源GSM4012改成的信号源，电阻150Ω/0.25w和1KΩ各一只，导线若干、100W白炽灯一只。测试设备：福禄克电力分析仪。

实验过程:1.将直流开关电源GSM4012中滤波电容焊下，经福禄克观察，通电后虽然开关电源成保护状态，但可间歇输出一个19.98KHz正向的锯齿波信号，由于间歇出现，福禄克保存不下来。但可明确产生19.98KHz幅值在4~10V左右的锯齿波信号。用它来代替15KHz信号源。

2.由于负载为100W白炽灯，负载电流很小，先用jct2515小功率单向硅做实验。首先用jct2515的阳极加与门极串入1k电阻的方式触发可控硅，阳极端加19.98KHz信号，可控硅阴极输出端未见信号输出，后门极触发电阻改为150Ω/0.25w电阻，可见输出19.98KHz信号。分析也可能为通过门极串过的信号，就改用直流电源进行有源触发方式。在门极串入150Ω电阻，硅的门极为电源正极入，硅的阴极为电源付端，接入直流电压为5v。通电后可控硅受控，可通断控制19.98KHz信号，通过的信号幅值应该有些衰减，衰减值大约1v左右。

3.小功率硅通过测试后，改用MTC55-12硅经测试也可受控。

操作者：朱东岳2014.6.6

第二篇：实验报告：验证机械能守恒定律

实验报告：验证机械能守恒定律

高一()班姓名实验时间

一、实验原理

1．机械能守恒定律

(1)当只有重力做功时，物体的但机械能的持不变．

(2)做自由落体运动的物体，只受重力作用，其机械能是守恒的．

2．实验原理

(1)如右图所示，借助打点计时器打出的纸带，测出物体自由下落的高度h和该时刻的速度

v，以纸带上的第n个点为例，如下图中的纸带，打第n个计数点时的瞬时速度等于以该时刻为中间时刻的某一段时间内的平均速度．即vn＝

(2)物体下落的高度为h时速度为v，则物体的重力势能

11减小量为mgh，动能增加量为mv2，如果mghv2，2

21即＝v2，就验证了机械能守恒定律． 2

二、实验器材

铁架台,，重物(带纸带夹子)，纸带，复写纸，导线，．

三、注意事项

1．安装打点计时器时，必须使两个限位孔的中线严格，以减小摩擦阻力．

2．应选用质量和密度较可使空气阻力减小．

3．实验时，应先接通，让打点计时器工作正常后再松开纸带让重物下落．

4．本实验中的两种验证方法中，均不需要测

5．速度不能用v＝或v＝计算，而应用vn＝进行测量并计算． ．．．

四、探究步骤,数据记录及处理

1．安装置：将打点计时器固定在铁架台上；用导线将打点计时器与

2．接电源，打纸带：把纸带的一端在重物上用夹子固定好，另一端穿过打点计时器的限位孔，用手提着纸带使重物停靠在打点计时器附近，接通电源，待打点稳定后松开纸带，让重物自由下落．重复几次，打下3～5条纸带．

3．选纸带：选取0,1,2,3….hn＋1－hn－14．数据处理：测出0到点

1、点

2、点3„的距离，即为对应的下落高度h1、h2、h3„；利用公式vn＝，2T

六、实验作业

1．在做“验证机械能守恒定律”的实验时，发现重锤减少的重力势能总是大于重锤增加的动能，造成这种现象的原因是()

A．选用重锤的质量过大B．空气对重锤的阻力和打点计时器对纸带的阻力

C．选用重锤的质量过小D．实验时操作不仔细，实验数据测量不准确

2.(2011年广州高一检测)某同学做验证机械能守恒定律实验时，不慎将一条挑选出的纸带一部分损坏，损坏的是前端部分．剩下的一段纸带上各相邻点间的距离已测出标在图7－9－4中，单位是cm.已知打点计时器工作频率为50 Hz，重力加速度g取9.8 m/s2.(1)重物在2点的速度v2＝________，在5点的速度v5＝________，此过程中动能增加量ΔEk＝________，重力势能减少量ΔEp＝________.(2)比较得ΔEk________ΔEp(填“大于”“等于”“小于”)，原因是__________________．由

以上可得出实验结论____________________．

3．某研究性学习小组在做“验证机械能守恒定律”的实验中，已知打点计时器所用电源的频率为50Hz．查得当地的重力加速度g＝9.80m/s2．测得所用重物的质量为1.00kg．实验中得到一条点迹清晰的纸带，把第一个点记作O，每两个计数点之间有四点未画出，另选连续的3个计数点A、B、C作为测量的点，如图所示．经测量知道A、B、C各点到O点的距离分别为50.50cm、86.00cm、130.50cm．根据以上数据，计算出打B点时重物的瞬时速度vB＝_____m/s；重物由O点运动到B点，重力势能减少了_______J，动能增加了_______J（保留3位有效数字）． 根据所测量的数据，还可以求出物体实际下落的加速度为_______ m/s2，则物体在下落的过程中所受到的阻力为_______N．

4.(2010年高考课标全国卷)如左图为验证机械能守恒定律的实验装置示意图．现有的器材为：带铁夹的铁架台、电磁打点计时器、纸带、带铁夹的重锤、天平．回答

下列问题：

(1)为完成此实验，除了所给的器材，还需要的器材有________．(填入正确选项前的字母)

A．米尺B．秒表C．0～12 V的直流电源D．0～12

V的交流电源

(2)实验中误差产生的原因有(写出两个原因)

第三篇：验证相对论关系实验报告

验证快速电子的动量与动能的相对论关系实验报告

摘要：实验利用β磁谱仪和NaI(Tl)单晶γ闪烁谱仪，通过对快速电子的动量值及动能的同时测定来验证动量和动能之间的相对论关系。同时介绍了β磁谱仪测量原理、NaI(Tl)单晶γ闪烁谱仪的使用方法及一些实验数据处理的思想方法。

关键词：电子的动量 电子的动能 相对论效应 β磁谱仪 闪烁记数器。

引言：

经典力学总结了低速的宏观的物理运动规律，它反映了牛顿的绝对时空观，却在高速微观的物理现象分析上遇见了极大的困难。随着20世纪初经典物理理论在电磁学和光学等领域的运用受阻，基于实验事实，爱因斯坦提出了狭义相对论，给出了科学而系统的时空观和物质观。为了验证相对论下的动量和动能的关系，必须选取一个适度接近光束的研究对象。的速度几近光速，可以为我们研究高速世界所利用。本实验我们利用源90Sr—90Y射出的具有连续能量分布的粒子和真空、非真空半圆聚焦磁谱仪测量快速电子的动量和能量，并验证快速电子的动量和能量之间的相对论关系。

实验方案：

一、实验内容测量快速电子的动量。测量快速电子的动能。验证快速电子的动量与动能之间的关系符合相对论效应。

二、实验原理

经典力学总结了低速物理的运动规律，它反映了牛顿的绝对时空观：认为时间和空间是两个独立的观念，彼此之间没有联系；同一物体在不同惯性参照系中观察到的运动学量(如坐标、速度)可通过伽利略变换而互相联系。这就是力学相对性原理：一切力学规律在伽利略变换下是不变的。

19世纪末至20世纪初，人们试图将伽利略变换和力学相对性原理推广到电磁学和光学时遇到了困难；实验证明对高速运动的物体伽利略变换是不正确的，实验还证明在所有惯性参照系中光在真空中的传播速度为同一常数。在此基础上，爱因斯坦于1905年提出了狭义相对论；并据此导出从一个惯性系到另一惯性系的变换方程即“洛伦兹变换”。

洛伦兹变换下，静止质量为m0，速度为v的物体，狭义相对论定义的动量p为：

p

m02vmv(5—1)

1mm0/,v/c。相对论的能量E为： 式中

Emc2(5—2)

这就是著名的质能关系。mc2是运动物体的总能量，当物体静止时v=0，物体的能量为E0=m0c2称为静止能量；两者之差为物体的动能Ek，即

Ekmc2m0c2m0c2（当β« 1时，式（5—3）可展开为

1

1)

(5—3)

2p1v112

Ekm0c2(1m0v22)m0c2c22m0(5—4)

即得经典力学中的动量—能量关系。

由式(5—1)和(5—2)可得：

E2c2p2E02(5—5)

这就是狭义相对论的动量与能量关系。而动能与动量的关系为：

EkEE0c2p2m02c4m0c2

(4─6)

这就是我们要验证的狭义相对论的动量与动能的关系。对高速电子其关系如图所示，图中pc用MeV作单位，电子的m0c2=0.511MeV。式(5—4)可化为：

p2c21p2c2

Ek

2m0c220511.以利于计算。

三、验仪器的介绍及方法：

1、实验装置主要由以下部分组成：①真空、非真空半圆聚焦β磁谱仪②β-放射源

Sr-90Y（强度≈1毫居里），定标用γ放射源137Cs和60Co（强度≈2微居里）③200μm厚Al窗NaI(Tl)闪烁探测器④数据处理软件⑤高压电源、放大器、多道脉冲幅度分析器。图2是实验装置图。

图1 实验装置图图2β磁谱仪的结构简图

β源：β源是放射高速运动β电子的源，高能β粒子的速度可接近光速。如其PC为1MeV时，v=0.89c，PC为2MeV时，v=0.97c。实验所使用的90Sr-90Yβ粒子源强度约为1.5 毫居里，在0~2.27MeV的范围内形成一连续的β谱。γ放射源：γ射线是一种波长很短的电磁波。γ射线与物质的相互作用要比带电粒子弱得多，因而它具有较强的穿透本领。我们实验中采用的137Cs与60Co两个γ放射源是作为定标源的，它们的强度约2微居里。

β磁谱仪：图2是β磁谱仪的结构简图，图中间长方形区域是一均匀磁场区，它是由垂直纸面的上、下两层产生均匀磁场的材料组成。而中间的空间可放一与其空间相吻合的真空盒，真空盒与真空泵、真空表相联结。均匀磁场方向是垂直纸面穿过真空盒的,真空盒的放入可使高速电子运动的区域为真空区。在磁场外左侧有一固定架可放置β源。β源放射的电子在保持磁场区B均匀不变的情况下，各个不同动量的电子将以不同半径R的半圆周运动被分离，这也称为磁分离技术。而闪烁探头与多道分析器是进行能量探测与能量幅度甄别的，与计算机相联后，探测到的粒子能量与粒子数将即时地在计算机上显示并图示。这里，闪烁探头是由碘化钠晶体和光电倍增管组成的，碘化钠晶体可把入射的高速B粒子动能转化成可见光脉冲；然后光电倍增管把这些光脉冲转化为电脉冲。磁谱仪长方形区域的右侧小区是放置γ放射源，进行定标与其他实验应用的。

微机与多道分析器：由光电倍增管产生的电脉冲经线性放大器放大后，由微机与多道分析器对它们进行幅度分析，按电脉冲幅度大小微机与多道分析器将其可分成512道或1024道（相当于一阶梯），即不同幅度的电脉冲计入不同的道（阶梯），电脉冲幅度越高，则所处的道数应越大。电脉冲幅度与阶梯道数关系的线性度与斜率可通过调节光电倍增管的高压与增益改变。所以每次实验测量前，需对微机与多道分析器的512道（或1024道）进行定标。定标实验可采用两个γ放射源的已知能谱图进行。

2、实验方法：β源射出的高速β粒子经准直后垂直射入一均匀磁场中(VB)，粒子因受到与运动方向垂直的洛伦兹力的作用而作圆周运动。如果不考虑其在空气中的能量损失(一般情况下为小量)，则粒子具有恒定的动量数值而仅仅是方向不断变化。粒子作圆周运动的方程为：

dp

evBdt(5—7)

e为电子电荷，v为粒子速度，B为磁场强度。由式(5—1)可知p=mv，对某一确定的动量数

值P，其运动速率为一常数，所以质量m是不变的，故

dvv2dpdv

m,dtR dtdt且

所以

peBR(5—8)

式中R为β粒子轨道的半径，为源与探测器间距的一半。

在磁场外距β源X处放置一个β能量探测器来接收从该处出射的β粒子，则这些粒子的能量(即动能)即可由探测器直接测出，而粒子的动量值即为：peBReBX/2。由于β

90源38

Sr9039Y(0～2.27MeV)射出的β粒子具有连续的能量分布(0～2.27MeV)，因此探测器在不同位置(不同就可测得一系列不同的能量与对应的动量值。这样就可以用实验方法确定测量范围内动能与动量的对应关系，进而验证相对论给出的这一关系的理论公式的正确性。

四、实验步骤:

1、检查仪器线路连接是否正确，然后开启高压电源，开始工作； 打开

动闪烁探测器使其狭缝对准

源的出射孔并开始记数测量；

定标源的盖子，移

2、调整加到闪烁探测器上的高压和放大数值，使测得的理的位置

3、选择好高压和放大数值后，稳定 10～20 分钟；

4、正式开始对 NaI(Tl)闪烁探测器进行能量定标，首先测量的1.33MeV 峰位道数在一个比较合的γ能谱，等1.33MeV 光

电峰的峰顶记数达到1000 以上后（尽量减少统计涨落带来的误差），对能谱进行数据分析，记录下1.17 和1.33MeV 两个光电峰在多道能谱分析器上对应的道数CH3、CH4；

5、移开探测器，关上使其狭缝对准

定标源的盖子，然后打开

定标源的盖子并移动闪烁探测器

源的出射孔并开始记数测量，等0.661MeV 光电峰的峰顶记数达到1000

后对能谱进行数据分析，记录下0.184MeV 反散射峰和0.661 MeV 光电峰在多道能谱分析

器上对应的道数CH1、CH2；

6、关上

定标源，打开机械泵抽真空（机械泵正常运转2～3 分钟即可停止工作）；

7、盖上有机玻璃罩，打开源的盖子开始测量快速电子的动量和动能，探测器与源的距

离ΔX最近要大于9cm、最远要小于24cm，保证获得动能范围0.4～1.8MeV 的电子；

8、选定探测器位置后开始逐个测量单能电子能峰，记下峰位道数 CH 和相应的位置坐标X；

9、全部数据测量完毕后关闭

10、实验完毕后，需要洗手

源及仪器电源，进行数据处理和计算。

五、数据记录与处理

1、定标数据：高压电源为667kv；放大倍数为0.3倍；放射源位置41.8

表格一

2、使用β源进行探测，β源位置为10.0cm处

表格二（坐标和道数的数据已经取平均值）

备注：选择四个孔分别为第2、4、6、8个

将表中的数据填入到数据处理软件进行数据处理，得到拟合曲线如附图所示以及得到的信息如下表格：

五、实验注意事项

1.闪烁探测器上的高压电源、前置电源、信号线绝对不可以接错； 2.装置的有机玻璃防护罩打开之前应先关闭β源； 3.应防止β源强烈震动，以免损坏它的密封薄膜； 4.移动真空盒时应格外小心，以防损坏密封薄膜；

六、实验结论

通过实验，我不仅巩固了放射源、闪烁探测器的正确使用。同时了解了β磁谱仪测量原理，并在实验中互相组合，从而验证快速电子的动量与动能的相对论关系。在实验过程中，林老师不仅向我们解释了一些实验原理中的难点，更向我们讲了物理学习乃至以后的物理教学过程中应注意的问题，以及解决方法，使我们懂得很多实验以外的知识。对实验现象和实验原理的理解不能只敷衍与表面，要深入探究。对于不理解的知识要提前查资料弄明白。林老师更是严格的指出了我预习报告中诸多的不足之处。本次实验我受益匪浅，相信必定在接下来乃至将来的工作上起到重要的意义。

七、参考文献：

[1]林根金等.近代物理实验讲义[M]．浙江师范大学数理信息学院近代物理实验室，2010 [2]林木欣．近代物理实验教程[M].北京：科学出版社，1999 [3]张天喆、董有尔.近代物理实验[M].北京：科学出版社，2004

第四篇：高频小信号放大器实验报告

南京信息工程大学滨江学院

高频电子线路实验报告

作者 徐飞 学号 20092334925 系部 电子工程系 专业班级 通信三班

实验一 高频小信号放大器实验

一、实验原理

高频小信号放大器的作用就是放大无线电设备中的高频小信号，以便作进一步变换或处

理。所谓“小信号”，主要是强调放大器应工作在线性范围。高频与低频小信号放大器的基 本构成相同，都包括有源器件（晶体管、集成放大器等）和负载电路，但有源器件的性能及负载电路的形式有很大差异。高频小信号放大器的基本类型是以各种选频网络作负载的频带 放大器，在某些场合，也采用无选频作用的负载电路，构成宽带放大器。

频带放大器最典型的单元电路如图所示，由单调谐回路做法在构成晶体管调谐放大器。

图电路中，晶体管直流偏置电路与低频放大器电路相同，由于工作频率高，旁路电

容Cb.、Ce可远小于低频放大器中旁路电容值。调谐回路的作用主要有两个：

晶体管单调谐回路调谐放大器

第一、选频作用，选择放大ff0的信号频率，抑制其它频率信号。

第二、提供晶体管集电极所需的负载电阻，同时进行阻抗匹配变换。

高频小信号频带放大器的主要性能指标有：

（1）中心频率 f0：指放大器的工作频率。它是设计放大电路时，选择有源器件、计算

谐振回路元件参数的依据。

（2）增益：指放大器对有用信号的放大能力。通常表示为在中心频率上的电压增益和

功率增益。

电压增益 AVOVO/Vi

功率增益 APOPO/Pi

式中 VO、Vi分别为放大器中心频率上的输出、输入电压幅度，PO、Pi分别为放大器中心频率上的输出、输入功率。增益通常用分贝表示。

（3）通频带：指放大电路增益由最大值下降 3db 时对应的频带宽度。它相当于输入不

变时，输出电压由最大值下降到 0.707 倍或功率下降到一半时对应的频带宽度。（4）选择性：指放大器对通频带之外干扰信号的衰减能力。通常有两种表征方法： 其一，用矩形系数说明邻近波道选择性的好坏。

其二，用抑制比来说明对带外某一特定干扰频率 fn信号抑制能力的大小，其定义为中心频率上功率增益 APf0与特定干扰频率fn上的功率增益 APfn之比：

df0

ApfnAp还有其它一些性能指标参数，如工作稳定性，噪声系数等。

高频小信号谐振放大电路如图所示：

高频小信号谐振放大器

晶体管基极为正偏，工作在甲类，负载为 LC 并联谐振回路，调谐在输入信号的频率

465khz 上。该放大电路能够对输入的高频小信号进行反向放大。

在 Multisim 7 电路窗口中，创建如图所示的高频小信号放大电路图，其中晶体管

Q1 选用虚拟晶体管。单击“防真”按钮，就可以从示波器中观察到输入与输出的信号波形。

二、实验内容

（一）频带放大器的测量

1.观察高频小信号放大器输入输出信号的波形，注意幅度变化和相位关系。

高频小信号放大器输入输出信号

2.高频小信号的选频作用

观察输入输出波形，分析产生此种现象的原因

3.高频小信号放大电路的通频带和矩形系数

利用 Multisim 7 仿真软件中所提供的波特图仪观察上述高频小信号放大电路的通频

带，将波特图仪接入高频小信号谐振放大电路，观察幅频特性。

4.观察双调谐回路高频小信号放大器输入与输出波形，分析幅频特性。

（二）宽带放大器的测量

观察输入输出信号的波形，分析幅频特性。

第五篇：信号与系统实验报告,

实验三

常见信号得ＭATLAB 表示及运算 一、实验目得 １。熟悉常见信号得意义、特性及波形 2.学会使用 MATLAB 表示信号得方法并绘制信号波形 ３、掌握使用ＭATLAB 进行信号基本运算得指令 4、熟悉用ＭATＬAB 实现卷积积分得方法 二、实验原理 根据ＭAＴLAB 得数值计算功能与符号运算功能，在 MATＬAＢ中，信号有两种表示方法,一种就是用向量来表示,另一种则就是用符号运算得方法。在采用适当得 MＡＴLAB 语句表示出信号后，就可以利用 MATＬAＢ中得绘图命令绘制出直观得信号波形了。

１、连续时间信号

从严格意义上讲,ＭATLＡＢ并不能处理连续信号。在ＭAＴLＡB 中，就是用连续信号在等时间间隔点上得样值来近似表示得,当取样时间间隔足够小时,这些离散得样值就能较好地近似出连续信号。在 MATＬAB 中连续信号可用向量或符号运算功能来表示。

⑴

向量表示法 对于连续时间信号，可以用两个行向量 f 与 t 来表示,其中向量 t 就是用形如得命令定义得时间范围向量,其中，为信号起始时间，为终止时间,p 为时间间隔。向量 f 为连续信号在向量 ｔ所定义得时间点上得样值． ⑵

符号运算表示法 如果一个信号或函数可以用符号表达式来表示,那么我们就可以用前面介绍得符号函数专用绘图命令 ezplot（）等函数来绘出信号得波形。

⑶

得 常见信号得 M ＡTLA Ｂ表示

单位阶跃信号 单位阶跃信号得定义为:

方法一：

调用 H ｅａviside(t)函数 首先定义函数 Heａｖisｉdｅ（ｔ)得ｍ函数文件,该文件名应与函数名同名即Ｈeaviside、m.％定义函数文件，函数名为 Hｅavｉsiｄｅ,输入变量为 x,输出变量为ｙ function y= Hｅaviside（t）

y＝(t＞0）;

%定义函数体,即函数所执行指令 %此处定义ｔ>0 时 y=1,t<＝0 时ｙ=0,注意与实际得阶跃信号定义得区别.方法二：数值计算法 在ＭATLAB 中,有一个专门用于表示单位阶跃信号得函数，即 s ｔe ｐｆun()函数，它就是用数值计算法表示得单位阶跃函数．其调用格式为: st ｅpfun(t,t0)

其中,t 就是以向量形式表示得变量，t0 表示信号发生突变得时刻,在ｔ０以前,函数值小于零,ｔ０以后函数值大于零。有趣得就是它同时还可以表示单位阶跃序列,这只要将自变量以及

取样间隔设定为整数即可。

符号函数 符号函数得定义为:

在 MATLAB 中有专门用于表示符号函数得函数 s ｉgn(),由于单位阶跃信号(t）与符号函数两者之间存在以下关系：,因此，利用这个函数就可以很容易地生成单位阶跃信号.2、离散时间信号 离散时间信号又叫离散时间序列,一般用 表示，其中变量 k 为整数，代表离散得采样时间点（采样次数)。

在 MATLAＢ中，离散信号得表示方法与连续信号不同，它无法用符号运算法来表示,而只能采用数值计算法表示,由于 MATLＡB 中元素得个数就是有限得，因此，MATＬAB无法表示无限序列；另外，在绘制离散信号时必须使用专门绘制离散数据得命令，即 stem（（）函数,而不能用ｐlot（）函数。

单位序列

单位序列）得定义为

单位阶跃序列 单位阶跃序列得定义为 ３、卷积积分 两个信号得卷积定义为:

MATＬAB 中就是利用 conv 函数来实现卷积得.功能：实现两个函数与得卷积.格式：ｇ=cｏnv(f1,ｆ2)

说明:ｆ１=f 1（t),f2=f 2(t）

表示两个函数,g＝ｇ(t)表示两个函数得卷积结果。

三、实验内容 1、分别用 MATLAＢ得向量表示法与符号运算功能，表示并绘出下列连续时间信号得波形:

⑴

⑵

（1）

t=－１:0、01:１0;t1=-１：０、01:-0、０1;ｔ2=０：0、01:10； f1=［zeros(1,lｅngth（t１)），ｏneｓ（1,lｅngｔh(t2))］;f=(2—exp（－２*ｔ））、＊f1； pｌoｔ（t，f)axiｓ（[-１,10,0,2、１]）

syｍs t;ｆ=ｓym（’(2-ｅxp(—２＊t))*hｅavisiｄe（ｔ)“)； ｅzｐｌoｔ(f,[-1,10]）；

(2)t=—2:0、01:8； f=0、*(t<０）+ｃos(pi*t／2)、*（t>０&ｔ〈４)+0、*（t〉4);ｐlot(ｔ,ｆ)

syms ｔ;f=sym(”ｃos（ｐi*t/２）*［hｅavｉside(t)—heavisｉｄｅ（t—4）] “）;ezplｏt(f，[-2，8］);

２、分别用 MATLAB 表示并绘出下列离散时间信号得波形:

⑵

⑶

（2）

t=0：8； ｔ１=—10:1５； ｆ=［zeros(1，10),ｔ,zeros(1,7)］;stem（ｔ1,f)axｉs(［—10，１5,0，1０]);

(3)t=0:5０;t１＝—10:５0； f=[ｚeｒos(1，10),sin(t*pi/４)];sｔem(t1,f）

ａxis（［—10，50,—2，２])

3、已知两信号，求卷积积分，并与例题比较。

t1=—1:0、01:0； t２=0:0、01:1;ｔ3＝—1:0、01：１； f1=ｏｎｅs(ｓiｚe(t1));f2＝ｏnes(size（t2）);g＝coｎｖ(ｆ1,f2)； sｕｂploｔ(3,1,１),pｌoｔ(t1,f1)； subplot(３，1,２),plot(t2，f2);subｐｌot(３,1，3),ｐlot(t3,ｇ);

与例题相比较,ｇ(t)得定义域不同,最大值对应得横坐标也不同。

4、已知，求两序列得卷积与 .N＝４;Ｍ=５； L=N+M—1； ｆ１＝[１，１，1，2］； f２=[1,2，3,4，５];g=ｃｏｎv(ｆ1,f2)； kf1=0:N-1； kf２=0：Ｍ－1;kg=0：L—１;ｓｕbｐlot（1,３,1）,stem(kf１，f1,’*k’);xｌabel（”k“)； yｌaｂｅl（’f1（k)”);grｉd on sｕbpｌｏt(1，3，2)，stem（ｋf2，f2，’*k“)；xlabeｌ（＇k’);ｙlabｅl（”f2(k)’);grｉd ｏｎ ｓubplot（1,3,3);steｍ（kｇ，g，＇＊k’)；xlabeｌ(＇k“)； ylabel（”g（ｋ)＇）;grid ｏn

实验心得:第一次接触 Mutlab 这个绘图软件,觉得挺新奇得,同时 ,由于之前不太学信号与系统遇到一些不懂得问题，结合这些图对信号与系统有更好得了解。

实验四

连续时间信号得频域分析 一、实验目得 １。熟悉傅里叶变换得性质 ２.熟悉常见信号得傅里叶变换 ３。了解傅里叶变换得ＭＡTＬAB 实现方法 二、实验原理 从已知信号求出相应得频谱函数得数学表示为：

傅里叶反变换得定义为：

在 MAＴLAＢ中实现傅里叶变换得方法有两种,一种就是利用 MATLAB 中得 Ｓy ｍbo ｌｉｃ Math Too ｌｂｏx 提供得专用函数直接求解函数得傅里叶变换与傅里叶反变换，另一种就是傅里叶变换得数值计算实现法.1、直接调用专用函数法 ①在 MATLAB 中实现傅里叶变换得函数为:

F=fｏurｉer(f)

对ｆ(t)进行傅里叶变换,其结果为 F(w)

Ｆ=fourier(f,v)

对 f（t)进行傅里叶变换,其结果为Ｆ（v)

F=fouriｅr(f,u，v）

对ｆ(u)进行傅里叶变换,其结果为 F(ｖ)②傅里叶反变换

f=ifourｉer（F)

对 F(w)进行傅里叶反变换，其结果为 f(ｘ)

f=ｉｆoｕriｅr(F,U)

对Ｆ（w)进行傅里叶反变换，其结果为ｆ（u)

f＝ifourｉeｒ(F,v，u）

对Ｆ（ｖ)进行傅里叶反变换,其结果为 f(ｕ）

注意：

（1)在调用函数 fｏuｒier（）及 ifouriｅｒ（）之前，要用 syms 命令对所有需要用到得变量(如 t,u,v，w)等进行说明,即要将这些变量说明成符号变量。对ｆourｉer（）中得 f 及ifourier()中得 F 也要用符号定义符 sｙm 将其说明为符号表达式。

(2)采用 fｏuｒier()及 fourｉeｒ（）得到得返回函数,仍然为符号表达式。在对其作图时要用 ezplｏｔ（）函数,而不能用ｐｌｏt（）函数.(3)ｆｏuriｅｒ（）及ｆourieｒ（）函数得应用有很多局限性,如果在返回函数中含有 δ(ω)等函数，则 ezploｔ（）函数也无法作出图来。另外,在用 fourier()函数对某些信号进行变换时,其返回函数如果包含一些不能直接表达得式子，则此时当然也就无法作图了。这就是ｆouriｅr（）函数得一个局限。另一个局限就是在很多场合，尽管原时间信号 f(t)就是连续得,但却不能表示成符号表达式，此时只能应用下面介绍得数值计算法来进行傅氏变换了,当然，大多数情况下，用数值计算法所求得频谱函数只就是一种近似值。

2、傅里叶变换得数值计算实现法 严格说来，如果不使用 symboｌic 工具箱,就是不能分析连续时间信号得。采用数值计算方法实现连续时间信号得傅里叶变换，实质上只就是借助于ＭATLAB 得强大数值计算功能,特别就是其强大得矩阵运算能力而进行得一种近似计算。傅里叶变换得数值计算实现法得原理如下: 对于连续时间信号 f(t),其傅里叶变换为：

其中 τ 为取样间隔,如果 f(t)就是时限信号,或者当|t｜大于某个给定值时，f(t)得值已经衰减得很厉害,可以近似地瞧成就是时限信号,则上式中得ｎ取值就就是有限得,假定为 N，有:

若对频率变量 ω 进行取样,得:

通常取：,其中就是要取得频率范围，或信号得频带宽度。采用 MAＴLAB 实现上式时，其要点就是要生成 f(t)得Ｎ个样本值得向量，以及向量,两向量得内积(即两矩阵得乘积）,结果即完成上式得傅里叶变换得数值计算。

注意：时间取样间隔 τ 得确定,其依据就是 τ 必须小于奈奎斯特(Nyquist)取样间隔。如果 f（t）不就是严格得带限信号，则可以根据实际计算得精度要求来确定一个适当得频率为信号得带宽。

三、实验内容 1、编程实现求下列信号得幅度频谱（1)

求出得频谱函数 F 1（jω)，请将它与上面门宽为 2 得门函数得频谱进行比较,观察两者得特点，说明两者得关系。

(2)三角脉冲

(3)单边指数信号

(4）

高斯信号

(1）

syｍs t w

Ｇt＝sym（“Heaｖiside（２*t+1）—Heaviｓiｄe(2*ｔ-１）’);

Fｗ=fｏurｉｅr(Gｔ,t,ｗ）；

FFｗ=maｐｌe(’conveｒt’，Fw,’pieｃewise”)；

ＦFP＝abs(FFｗ);

eｚplot（FＦＰ,［—10*pi 10*pｉ］);grｉd;

ａxｉs(［-１０*ｐｉ １0＊pi 0 2、2])

与得频谱比较,得频谱函数 F 1(jω)最大值就是其得１／2．（２)sｙmｓ ｔ w;Ｇt=sym（“（1＋t）＊（Hｅaviｓｉｄe(ｔ+1)—Heａviｓｉdｅ（t)）+(1-t）＊（Ｈeaviside（t)—Heaｖiside（t—1）)”);Fｗ=fourｉer(Gt，t，w）；

ＦFw=ｍaple(“coｎｖeｒt＇,Fw，’piｅcewｉｓe”）;

FFP＝abｓ（ＦＦw)；

ｅzpｌot（FFP，[—1０*pi 10＊pｉ］）；grid；

axis([—10*pi １０*pｉ 0 ２、２］)

(3)ｓymｓ t w

Ｇt＝sｙm(’exｐ(-ｔ)＊Ｈeavisiｄe(t）’);

Fｗ＝fouｒieｒ(Ｇt,t,w)；

ＦＦｗ=maple(“ｃoｎvｅrt”,Fw,’ｐiecｅｗｉse’);

ＦFＰ＝aｂs(ＦFw);

eｚplｏt(FＦP,［—10*pi １0*pi])；grid；

axis([—1０＊pi 10*ｐi —1 ２］）

(4)symｓ t w

Gt=ｓym(’eｘp(-t^2)“）;

Fｗ=foｕrｉｅｒ（Gt,t,w);

FFw＝mapｌｅ（＇convert’，Fw，’pｉecｅｗise’)；

eｚpｌｏｔ(ＦＦw，[－30 30］）;grid；

axis（[—30 30 —1 2]）

２、利用 iｆouriｅr（）函数求下列频谱函数得傅氏反变换（１)

(2）

(1)syms t w

Fw=sｙｍ(’-i*2＊w/(16＋w^２)’）;

ｆt=ifourier(Ｆｗ,w，t）;

ｆt 运行结果: ｆt = —ｅxp(４*ｔ)*ｈeａｖｉsｉde(—t)+exp（—４＊t）＊heaｖｉｓｉdｅ（t)(2）

syms t w

Fw＝sym(”(（i*w）^2+5*i＊w-8）／((i＊w)^2+6*i*w+5）’)；

ft=ｉｆoｕriｅｒ（Ｆw,w，t）;

ft 运行结果: ft = diｒac(t)+(-3＊exp(-t)＋2＊exp(-5*t)）*ｈｅａvisiｄe（t)实验 心得 maｔlab 不但具有数值计算能力,还能建模仿真,能帮助我们理解不同时间信号得频域分析。

实验五 连续时间系统得频域分析 一、实验目得 1.学习由系统函数确定系统频率特性得方法.2.学习与掌握连续时间系统得频率特性及其幅度特性、相位特性得物理意义.3.通过本实验了解低通、高通、带通、全通滤波器得性能及特点。

二、实验原理及方法 频域分析法与时域分析法得不同之处主要在于信号分解得单元函数不同。在频域分析法中，信号分解成一系列不同幅度、不同频率得等幅正弦函数,通过求取对每一单元激励产生得响应,并将响应叠加，再转换到时域以得到系统得总响应。所以说,频域分析法就是一种变域分析法.它把时域中求解响应得问题通过 Fｏｕrier 级数或 Foｕrier 变换转换成频域中得问题;在频域中求解后再转换回时域从而得到最终结果.在实际应用中,多使用另一种变域分析法：复频域分析法,即 Laplａce 变换分析法。

所谓频率特性，也称频率响应特性,就是指系统在正弦信号激励下稳态响应随频率变化得情况,包括幅度随频率得响应与相位随频率得响应两个方面.利用系统函数也可以确定系统频率特性，公式如下:

幅度响应用表示,相位响应用表示。

本实验所研究得系统函数 H（ｓ)就是有理函数形式,也就就是说,分子、分母分别就是 m、n 阶多项式。

要计算频率特性，可以写出

为了计算出、得值，可以利用复数三角形式得一个重要特性：

而，则 利用这些公式可以化简高次幂,因此分子与分母得复数多项式就可以转化为分别对实部与虚部得实数运算,算出分子、分母得实部、虚部值后，最后就可以计算出幅度、相位得值了。

三、实验内容 a),m 取值区间 ［0,1］,绘制一组曲线 m=０、1,0、3,0、5，0、7,0、9;b)绘制下列系统得幅频响应对数曲线与相频响应曲线,分析其频率特性.(1)

(2)

(3)

a)% deｓign2、m

figuｒe

ａlpha＝［0、1，0、３,0、5,0、7,0、9］;

ｃolorｎ=[＇ｒ’ ’g’ ’b“ ’y” “ｋ＇];

%

r ｇ b y m c k(红,绿，蓝,黄，品红,青,黑)

fｏｒ ｎ=１:5

b=[0 ａlpha(n)］;

% 分子系数向量

a=［alpha(n）-alpｈa（n）＾２ 1]；

% 分母系数向量

printsys(ｂ,a,”s“)

[Hz,ｗ］=freｑs(ｂ，a)；

ｗ=w、/pi;

mａgｈ=ａbｓ(Hz）；

ｚeroｓIｎdx＝fｉｎd(ｍagh==０);

mａgh(zｅｒｏsIndx）=1;

ｍagh＝20＊loｇ10(ｍaｇｈ）；

magh（zeroｓIndx）=-ｉnｆ；

angh＝angle（Hz）;

ａngh=ｕnｗｒap（anｇh)＊１80/ｐi;

sｕbpｌot（１,2,1)

ｐlot(w，ｍａｇh,coｌorn（n））；

ｈolｄ ｏn

subｐlｏt(1,２,2）

pｌoｔ（w,angh，coｌorn（n)）;

ｈold on

end

ｓubpｌoｔ（1，２,１）

ｈoｌd ofｆ

xlａｂｅl(”特征角频率（timｅspi rａd/sample)“)

title(＇幅频特性曲线 |H(ｗ）|(dB)”);

ｓｕbpｌot（1,2,2)

ｈold ｏｆf

xlａbel(’特征角频率(tiｍｅsｐｉ rad/sａｍple）’）

ｔitle（“相频特性曲线 theta(w)（degｒees）’）;

ｂ)(1)% deｓｉgn1、ｍ ｂ=［1,０]；

% 分子系数向量 a=[1，1］;

% 分母系数向量 prinｔsｙｓ（ｂ,ａ,”s’)[Hz，ｗ]=ｆrｅqｓ(b，a);w=ｗ、/pｉ;magh=abs（Hz）;zｅroｓＩndx=ｆｉnd(magh==０)； magｈ(zerosInｄx)=1； magｈ=2０*ｌｏg10(ｍagh)；

% 以分贝 maｇh（ｚｅrosInｄｘ)=-inf;anｇh＝anｇｌe（Ｈz);aｎｇh=unｗｒａp(angｈ）*18０／pi；

% 角度换算 fiｇuｒｅ subplot（1，2，1)ｐlｏｔ（w,magh);gｒid on xlabeｌ(’特征角频率（ｔimeｓpi rａd/sample)＇)titlｅ(’幅频特性曲线 |H(w)|(dＢ)’)； subpｌｏt(１,2，2)plot(w,angh）;gｒiｄ on xlabel（’特征角频率（ｔimes＼pｉ raｄ／sａmｐｌe）’)tiｔlｅ（’相频特性曲线 ＼thｅta（w）

(deｇreeｓ)’)；

（2）

% desｉgn1、ｍ b=[０,１,0］;

% 分子系数向量 a=[１，3，２]；

% 分母系数向量 pｒintsys（b,a，’s’)［Hｚ,ｗ]＝ｆreqｓ(ｂ,a);w=w、/pi； magh=abs(Hz);zerｏsInｄx＝finｄ（ｍagh==０）； magｈ(zeｒｏｓＩndｘ）＝1； ｍagh＝20＊log10（magh);

% 以分贝 magh（zerosIndx）＝-iｎｆ； angh=ａngle(Ｈz);angｈ=ｕnｗrａp(angh)*180/pi;

％ 角度换算 ｆiｇuｒe ｓubｐlｏｔ(１,2,1)ｐlot(w，ｍagｈ);griｄ ｏn ｘｌaｂｅl(“特征角频率(＼times＼ｐｉ rad／ｓamｐle)＇）

tｉtle（’幅频特性曲线 ｜Ｈ（ｗ)｜(dB）’);subｐlｏt（1,2，２)pｌｏt(w，aｎｇh)； grid oｎ xlａbel（”特征角频率（＼times＼pi rａｄ/ｓaｍple）“)titlｅ（”相频特性曲线 ｔheｔａ(w)(dｅｇrees）’)；

(3）

% ｄesiｇn1、m ｂ=［１,-1］;

％ 分子系数向量 ａ=[１,1］;

% 分母系数向量 prinｔsys（b,a，“s”)[Hz,w］＝freqｓ(b,a);w＝ｗ、/pi;maｇh=abs(Hz）;zerosIｎdｘ=find(mａｇh=＝0);mａｇh（zeｒｏsＩndx)＝１;magh＝2０*log10(magh);

% 以分贝 ｍａgｈ(ｚerosIndx)=-inｆ;angh=angle(Hz);ａngｈ=unwｒap（ａｎgh)*18０／pｉ；

％ 角度换算 ｆiｇure ｓubplｏt(1,2,１）

plｏt(w，ｍagh)； gｒiｄ on xlabel(’特征角频率(ｔimｅsｐi ｒad／ｓaｍｐle）“）

titｌe(”幅频特性曲线 ｜H(ｗ）|(dB)’）;subplot(1，２,2)plot（w，angh);grｉd on xｌａbｅl(’特征角频率(times＼pi raｄ/ｓａmple)＇)titｌe(’相频特性曲线 theｔａ(w）

(degｒeｅs)“）；

实验心得: :虽然之前用公式转换到频域上分析,但就是有时会觉得挺抽象得,不太好理解。根据这些图像结合起来更进一步对信号得了解。同时，这个在编程序时，虽然遇到一些问题，但就是总算解决了。

实验六

离散时间系统得 Z 域分析 一、实验目得 1.学习与掌握离散系统得频率特性及其幅度特性、相位特性得物理意义。

2.深入理解离散系统频率特性与对称性与周期性。

3.认识离散系统频率特性与系统参数之间得系统 4.通过阅读、修改并调试本实验所给源程序，加强计算机编程能力。

二、

实验原理及方法 对于离散时间系统,系统单位冲激响应序列得 Fｏｕｒier 变换完全反映了系统自身得频率特性,称为离散系统得频率特性,可由系统函数求出,关系式如下：

（6 – １)由于就是频率得周期函数,所以系统得频率特性也就是频率得周期函数,且周期为,因此研究系统频率特性只要在范围内就可以了.       n n nj jn n h j n n h e n h e H)sin()()cos()()()(  

（6 – 2)容易证明，其实部就是得偶函数，虚部就是得奇函数,其模得得偶函数,相位就是得奇函数。因此研究系统幅度特性、相位特性,只要在范围内讨论即可。

综上所述,系统频率特性具有周期性与对称性，深入理解这一点就是十分重要得。

当离散系统得系统结构一定，它得频率特性将随参数选择得不同而不同,这表明了系统结构、参数、特性三者之间得关系,即同一结构,参数不同其特性也不同。

例如，下图所示离散系统，其数学模型由线性常系数差分方程描述：

系统函数: 系统函数频率特性：

幅频特性: 相频特性：

容易分析出,当时系统呈低通特性,当时系统呈高通特性;当时系统呈全通特性.同时说明,在系统结构如图所示一定时，其频率特性随参数 a 得变化而变化.三、实验内容 ａ)。

b)c)ａ)% dｅsigｎ1、m b=[1，０,-1];

％ 分子系数向量 ａ=[1,0,—0、81］；

% 分母系数向量 printｓys（ｂ，a，”z“)［Ｈz,w］＝fｒeｑz(b，a);w＝ｗ、/pｉ;magh＝abｓ(Hｚ);zerosＩndx＝find(maｇh==0);magh(zｅrｏsInｄｘ)=1;magh=２0*ｌog10(magh);

% 以分贝 magｈ（zeｒosIｎdx）=-inｆ； angｈ=angle（Ｈz)； ａngh=unwｒaｐ(anｇh）*180/pｉ;

% 角度换算 ｆigurｅ ｓubplot（1,2,１）

plot(ｗ,magh);gｒid ｏn xlabel(’特征角频率（timesｐｉ rａd/saｍpｌe）＇)tiｔle(’幅频特性曲线 |H（w)|(ｄB)”）;sｕｂplot（1，2，2)plｏt（w,ａngh);grｉｄ on xlabel（“特征角频率（times＼pi rad/sａmple）”)tｉtｌｅ(＇相频特性曲线 tｈeta(w)（degｒees)“);

带通

ｂ)% ｄesｉgｎ１、m b=［0、１,—0、3,0、3，-0、1]；

％ 分子系数向量 ａ=［1,0、6,０、4，0、1］；

% 分母系数向量 pｒｉｎtsyｓ(ｂ,a，’ｚ”)[Hz,w]=ｆreqz(b，a);w=w、／pi； maｇh=ａbs(Hz）； zerｏsＩndx=find(maｇh==0);maｇh（zerｏsIndx)＝１;magh=２0*lｏg10(mａgｈ)；

％ 以分贝 ｍａgh(ｚerｏsIｎｄx)＝-inf;aｎｇｈ=ａｎglｅ(Hz);aｎgｈ=ｕｎwｒａp（ａｎｇh)＊1８0／pi；

% 角度换算 figure sｕbplot(1,2,1)pｌoｔ(w,maｇh）;grｉd on xlabｅl(’特征角频率(tｉｍｅspi raｄ／sａmpｌｅ)’）

titｌe(“幅频特性曲线 |H(w)|(ｄB)”）;ｓubｐlot（1，2，2)ploｔ（w,angh）;ｇｒid on

xlabel(“特征角频率(＼timｅspi raｄ/saｍｐle）’)title(”相频特性曲线 theta（w）

(dｅｇrees）’）；

高通

ｃ)% ｄesign1、m b=[１，—1,0];

% 分子系数向量 a=［1，０,０、8１];

% 分母系数向量 prinｔｓyｓ（b，ａ,“z’)［Hｚ,w]=freqｚ(b,a）;w=w、/pｉ； magｈ=aｂs(Hz）； zerosIｎdx=fiｎd(mａgｈ＝＝0）;magh（zerosIndx）＝1;mａgh=20*log1０（magh）;

％ 以分贝 ｍagh(ｚｅrosIｎdx）=—ｉnｆ;anｇh=ａngｌe（Hｚ)； angh＝unwrap（aｎgh)*１80/pｉ；

％ 角度换算 fｉgure suｂpｌoｔ（1,２,1)pｌｏt（w,magｈ);ｇrid ｏｎ xｌabｅｌ(”特征角频率(＼timeｓ＼pi rａｄ/saｍple）＇)title（“幅频特性曲线 |H（w)｜(dB)”）;ｓubplｏt(１,2，2）

pｌot(w,aｎgh);

griｄ on xｌaｂｅl（’特征角频率(＼tiｍespi rad/sａmple)")tiｔle（’相频特性曲线 theｔa（ｗ）

(dｅｇreｅs）’)；

带通

实验心得: :本来理论知识不就是很强得,虽然已经编出程序得到相关图形,但就是不会辨别相关通带,这让我深刻地反省。