薄膜干涉实验的计算机模拟分析

薄膜干涉实验的计算机模拟分析
薄膜干涉实验的计算机模拟分析

目录

1 绪论 (1)

2 薄膜干涉的相关介绍 (2)

2.1薄膜干涉的理论基础 (2)

2.1.1 薄膜干涉公式推导 (2)

2.1.2 薄膜干涉的附加光程差和条纹特点 (4)

2.1.3 薄膜干涉中的半波损失与薄膜厚度 (6)

2.2 薄膜干涉实验种类 (8)

2.2.1 劈尖干涉 (8)

2.2.2 牛顿环干涉 (9)

2.2.3 迈克尔逊干涉 (9)

3 薄膜干涉实验的计算机模拟与分析 (10)

3.1 劈尖干涉实验的模拟与分析 (11)

3.2 牛顿环实验的模拟与分析 (14)

3.3 迈克尔逊干涉实验的模拟与分析 (16)

3.4 薄膜干涉实验仿真结果对比分析 (18)

结论 (20)

参考文献 (21)

致谢 (22)

附录Ⅰ、Ⅱ、Ⅲ (23)

1 绪论

近几十年来,现代光学的最重要的进展之一就是光学信息处理与数字光计算的飞速发展。光学信息处理是以傅里叶分析方法为核心,研究光学成像和光学变换的理论与技术。它以光子传递信息,利用光学或者光电子器件进行操作运算,用光的折射、干涉和衍射等特性来实现对输入信息的各种变换与处理。光学信息处理的历史可以追溯到1873年阿贝提出的二次成像理论。傅里叶光学的核心是基于标量衍射理论的菲涅尔衍射和夫琅和费衍射、透镜成像性质以及利用傅里叶变换分析光学成像系统等[1]。这些处理过程通常都需要在光学暗室中进行,而且很容易受到外界环境的影响。

在计算机高速发展的今天,计算机仿真作为虚拟实验手段已经得到了长足的发展。其具有良好的可控性(参数可以根据需要调整)、无破坏性(不会导致器件的损坏)、可重复性(受一些随机因素影响较小)、易观察性(不至于稍纵即逝)和经济性(不需要贵重仪器)。所以利用仿真结果指导实验,可以减少或者避免仪器不必要的损伤。国外的光学信息处理仿真是在模拟设计和优化光学系统的基础上发展起来的。最有代表性的是美国劳伦斯利弗莫尔实验室的Prop92和法国的光学传输软件Miro。俄罗斯也有一套较成熟的商业光传输软件菲涅尔。我国在这方面起步较晚,SG99刚在神光一Ⅲ上运行的可行性论证[2]。

事实上,利用MATLAB在科学计算上的优点借助计算机对真实实验的模仿,使用者可以通过计算机提示的信息直接参与模拟操作。可以与真实仪器相比较,而且比真实实验直观准确,减少仪器特别是贵重仪器的损伤率。在进行仪器设计之前,还能够对其性能作出评估以及可行性论证,节约时间和经费。长期以来,光学课程的抽象性较强,理论教学对实验的依赖性较大,给学生学习该门课程带来了诸多困难[3]。将MA TLAB这一工具引入到光学实验中,利用计算机对一些光学实验进行模拟,对计算结果可视化,化抽象思维为形象思维,从而更好地洞察含义、理解概念、发现规律。还可以在分析中改变参数值,便于理论研究,也给初学者带来方便。与传统的实验方式相辅相成,实现更好的教学效果[4]。

2 薄膜干涉的相关介绍

光学是一门古老的科学,到了l9世纪波动光学理论经过杨、菲涅尔和麦克斯韦等人的工作变得日臻完善。在大学物理的课程教学中,要求学生对波动光学的掌握是非常重要的一项教学目标,干涉作为波动光学中的一部分显得尤为重要。

2.1薄膜干涉的理论基础 当两列相干光相遇时,其合振动的强度为()??++=cos 22121I I I I I ,式中1I 、2I 是两列光波的强度,??是这两列光相遇时它们的相位差。而??取决于两列波的光程差δ。

2.1.1 薄膜干涉公式推导

设薄膜的折射率为n ,厚度为h ,薄膜上方媒质折射率为1n ,薄膜下方媒质折射率为2n ,并设n >1n ,n >2n 。如图2-1所示。

图2-1 折射图

单色入射光a 由薄膜上表面A 点和下表面B 点反射,形成相互平行的两束光线1a 和2a 。根据波在传播过程中各点之间的相位关系,光线1a 上C '点相位比A 点相位落后

πλπ

-'C A 12

其中π为A 点反射时的半波损失现象所引起的附加相位落后,π前亦可取正号;同理,光线2a 上C 点相位比A 点落后

()BC AB +膜λπ

2

B 点反射无半波损失现象。两光线1a 和2a 的相位差(在

C '点及C 点以后不会再引起相位差)为

()πλπλπ

?+-+=?'221

AC BC AB 膜 (2-1) ()?

??+=干涉减弱,干涉加强,ππ122k k 设真空中光的波长为λ,则有

n ,11λλλλ==

膜n

(2-1)式变为 ()

????????+-+=?212'1λλπ?AC n BC AB n ()???+=干涉减弱

,干涉加强,ππ122k k (2-2) 经简单整理得

()21λ

+'-+C A n BC AB n ()???+=干涉减弱

,干涉加强,ππ122k k (2-3) (2-1)式中出现了折射率×几何路程,即干涉结果可用折射率×几何程之差来表示。因此把折射率×几何路程定义为光程,式中()BC AB n +光线2a 从A →B →C '的光程;C A n '1为光线a 1从A →C '的光程。2

λ为薄膜上表面A 点反射时因有半波损失现象所引起的附加光程,用δ表示这两条光线的光程差,于是

δ=()=+'-+21λC A n BC AB n ()???+干涉减弱

,干涉加强,ππ122k k (2-4) 再利用图中的几何关系及简单的数学推导,可得:

=+

-=2sin 22212λδi n n h ()???+干涉减弱,干涉加强,ππ122k k (2-5) (2-5)式中2

λ这项可有可无,依具体情况而定。如果在薄膜上表面和下表面反射其

中只有一处有半波损失现象,则此项一定存在。

2.1.2 薄膜干涉的附加光程差和条纹特点

薄膜等倾干涉是光学的重要内容之一,在薄膜干涉中,往往存在附加光成差,教材一般只对干涉条纹的特点给出定性的结论。实际上可根据菲涅尔公式,从半波损失的角度,较简单地判断出是否存在附加光成差。另外,对干涉条纹也可以通过简单的数学推导给出结果,本文就讨论这两个问题。

(1)附加光程差

设薄膜折射率为2n ,薄膜上方介质折射率为1n ,其下方介质折射率为3n ,即由上自下的折射率分别为1n 、2n 、3n ,如图2-2所示。

图2-2 光程差① 当1n >2n >3n 或1n <2n <3n 时。前者由于入射光线由光密介质射入光疏介质并在界面上产生反射。根据菲涅尔对半波损失的解释,上界面的反射无半波损失;在下界面上,由于2n >3n ,光线也是由光密介质射入光疏介质并在下界面上产生反射的,故也无半波损失;总的过程无半波损失。所以,当1n >2n >3n 时无半波损失,即由上表面出射的两束相干光的附加光程差δ=0。对于后者,上、下界面反射均存在半波损失

2

λ,两者刚好相互抵消,所以,在1n <2n <3n 的条件下,附加光程差也为零。 ②1n >2n >3n 或1n <2n <3n ,前者光线在上界面反射存在半波损失,而在下界面反射时无半波损失,所以,由上界面出射的两柬相干光的附加光程差

δ=2λ;同理后者也有附加光程差δ=2

λ。 (2)条纹特点

① 干涉条纹的形状是明暗相间的同心圆环,如图2-3所示。先考虑扩展光源

上一点1S 的情形,从1S 发出

2

sin n -d 21221220λδ-=i n 可知,不同的倾角对应不同的级次,凡同一倾角的光线,如S1发出的光线l 和2,经薄膜上下表面反射后,形成一个薄层圆锥面,如图中1' 、2'、2'' 所示;再经透镜会聚后光线仍关于透镜主轴对称,必然在透镜焦平面上形成一个圆。扩展光源上不同点形成的同一级次的干涉条纹总是重合的。因此,于涉条纹的形状是明暗相间的同心圆环。

② 干涉条纹的分布足不均匀的,内疏外密。根据光程差公式,第j 级条纹为

22122122221sin sin d j i n i n n λ??? ??+==- (2-6) j 越大,2cos i 越大,即2i 越小,也就是说,越靠近环心,干涉级次越高;越向外,干涉级次越低,这和牛顿环是不同的。第j+l 级条纹为

022*********sin sin d j i n i n n λ??? ?

?+='=- (2-7) 由(2)减(1)式,得22cos cos i i -'022d n λ=,由微分学知22

cos cos i i -'22sin i i ?-≈,估()222sin 2i n i λ

=?。由此可见,干涉条纹越向外,由于2i 增大。2sin i 增大,2

i ?减小,因此,越向外条纹越密,越不易辨认[5]。

图2-3 等倾干涉特点

2.1.3 薄膜干涉中的半波损失与薄膜厚度

在解释薄膜厚度问题时,无法回避半波损失的概念。所谓光波的半波损失是指:当光波被光密介质反射时,反射波的周相要发生π的突变或接近π的突变。

(1) 增透膜中的半波损失

图2-4中为增透膜示意图,其中0n 、1n 、2n 分别表示空气、膜层和玻璃的折射率。当入射光线SA 入射到薄膜时,将分裂成功率递减的1、2、3…各条反射光线和Ⅰ、Ⅱ、Ⅲ…各条透射光线。若薄膜前后两面相互平行,则反射光线和透射光线将组成两平行光线组。

图2-4 增透膜示意图

① 反射光组中的半波损失

对增透膜而言,0n <1n <2n (如空气折射率0n =1、MgF 2折射率1n =1.38、冕牌玻璃2n =1.52),光线在空气进入薄膜(光线由光疏到光密媒质)和由薄膜进入玻璃时,上下两面反射均出现半波损失,即产生π的突变。因此,反射光组1、2、3…中由于反射引起的周相差为零。如果相邻两反射光相互抵消,则有光程差:

01)2

1(cos 2λγ+==?k dn (2-8) (其中0λ为光在空气中的波长,d 为薄膜的厚度,γ为光在薄膜中的折射角,k=0、1、2、…),当SA 正射薄膜时,γ=0,则

01212λ??? ??+==?k dn ,2

2101λ??? ??+=k d n (2-9) 又101λλ=n ,所以薄膜厚度14

12λ+=k d ,当k=0时;41λ=d ,即增透膜薄膜厚度

为入射光在薄膜中波长的1/4。

② 透射光组中的半波损失

由于反射光组1、2、3…的互相抵消,就大大减少了反射光组的能量。同样因为0n <1n <2n ,所以光在薄膜中连续两次反射时,只在薄膜下表面出现半波损失,产生π的突变。直接由薄膜下表面透射的光与在薄膜中连续两次反射后再透射的光,由于反射引起的周相差为π,同理可知在透射光组I 、Ⅱ、Ⅲ…中相邻两透射光线由于反射引起的周相差均为π。当反射光组1、2、3…互相抵消时,透射光组I 、Ⅱ、Ⅲ…中相邻的两透射光线的光程差也为:

01)21(cos 2λγ+==?k dn 加上周相差π,所以出射时,相邻的两透射光线将是互相加强,从而使透射光的

能量也得以加强。

在增透膜实现增透的过程中不难发现:反射光组1、2、3…的互相抵消与透射光组I 、Ⅱ、Ⅲ…的互相加强形成了互补效应,体现出于涉效果只是将能量重新分配到不同地方,而总能量仍然守恒[6]。

(2) 皂液薄膜干涉中的半波损失

① 理解上的偏差

增透膜是薄膜干涉知识的应用,我们在学完增透膜的内容后,对薄膜干涉现象的认识将会出现反复。在增透膜的增透过程中,由教科书可知:“当薄膜的厚度是人射光在薄膜中波长的1/4时,在薄膜的两个面上反射的光,路程差恰好等于半个波长,因互相抵消”。在双缝干涉的学习中,在屏上距两个狭缝的路程差是半个波长的奇数倍处出现暗条纹。图2-5为皂液薄膜的干涉模型,我们在判断出现暗纹处的薄膜厚度d 时,依据上述两情形的结论极易作出如下推理:当皂液薄膜前后两表面反射回来的两反射光的路程差:λ??? ?

?+=?21k 时(其中λ为光在皂液中的波长,k=0、1、2…),两反射光也互相抵消。从而得出“当薄膜厚度λ4

12+=k d 处(如图2-5示),两反射光互相抵消,形成暗条纹”这一错误结论。

图2-5 皂液薄膜的干涉模型 比较皂液薄膜与增透膜相干的过程,可以看出:由于MgF 2薄膜前后两面反射均出现半波损失,使反射引起的周相差为零。而皂液薄膜前后两面反射时只在前一面出现半波损失,使反射引起的周相差为π。虽然薄膜厚度均为1/4个波长的奇数倍,但现象却是一个互相抵消,一个互相加强,结果截然相反。教材在前后处理双缝干涉和增透膜的机理中,引入了路程差的概念,对相干问题加以量化处理,接连两次涉及出现暗纹的条件,强化了“路程差为半个波长的奇数倍”这一条件,而在皂液薄膜干涉及空气薄膜干涉中却回避了路程差的问题[7]。究其原因,教科书在编写上难辞其咎,既然在皂液薄膜干涉中由于半波损失而回避了路程差及薄膜厚度问题,那么在增透膜中也应回避此类问题,这样才能保证教材前后的连续性和完整性,才能有利于学生系统地把握知识,更有利于学生思维的正向迁移[6]。

2.2 薄膜干涉实验种类

在教学中,常见的几种薄膜干涉有劈尖干涉、牛顿环干涉以及迈克尔逊干涉仪所产生的干涉现象。

2.2.1 劈尖干涉

如图2-6所示,将一劈尖放置于空气中,劈尖材料折射率为n ,劈尖上下两表面左端点距离为H ,夹角为θ。一束波长为λ的平行光自上而下垂直地射向劈尖,上下两表面的反射光相遇后会形成稳定的干涉条纹。干涉条纹的定域是一个非常复杂的问题,教学中考虑的是定域在劈尖的上表面。设所考察的定域点距离上表面左端的距离为x ,则该处劈尖的厚度为H x sin +θ ,因此两反射光的光程差为:

()2sin 2λθδ+

+=H x n

图2-6 劈尖干涉装置

2.2.2 牛顿环干涉

牛顿环实验装置如图2-7所示,

图2-7 牛顿环装置

上部为一曲率半径为R 的平凸透镜,下部为一平板玻璃,中间是折射率为n 的介质。设透镜底部与玻璃板的距离为H ,单色平行入射光波长为λ。图2-7中所示半径为r 处两反射光的光程差为

()2222λ

δ++--=H r R R n

2.2.3 迈克尔逊干涉

迈克尔逊干涉仪光路如图2-8所示,1G 、2G 分别是分光板和补偿板,1M 、2M 是两块平面镜用于将光反射回去。其中1M 是固定的,2M 可左右移动并旋转,

'2M 是2M 关于1G 所成的虚像。

为了便于讨论,我们将光路图简化为图2-9。A 点代表观察者的位置,O 点是A 点在1M 平面上的垂足,z =AO 。P 代表两反射光的相遇点,α表示AP 的连线与Z 轴的夹角,实际上这个角等于在P 点相遇的两光的入射角[7]。设图中Z 轴与两平面交点之间的距离为d ,则P 点处,沿Z 轴方向上两平面的距离为

θtan x d +,

因此,P 点处两反射光的光程差为()αθδcos tan 2x d +=。三角形AOP 是一直角三角形,AP 是斜边,则222cos z y x z AP AO ++==

α,最终相遇于P 点

的两反射光的光程差可表示为: ()222tan 2z y x z

x d +++=θδ

图2-9 迈克尔逊光路简化图

3 薄膜干涉实验的计算机模拟与分析

常见的干涉实验有劈尖干涉、牛顿环、迈克尔逊干涉。作为薄膜干涉的典型代表,对它们进行计算机的模拟对教学实验具有一定的指导意义。

3.1 劈尖干涉实验的模拟与分析

根据光的干涉理论,应用MATLAB 软件编程,模拟了劈尖干涉实验光强分布图形,通过改变劈尖的厚度、入射波长和介质折射率等参数,观察到劈尖干涉条纹的相应变化,反过来,通过干涉条纹的相应变化,也可以测量介质薄膜的厚度、折射率和入射波长等。这种计算机仿真方法扩展了劈尖干涉问题的研究途径与方法,具有一定的应用价值,尤其对于教学具有直观作用,有助于更加深刻地理解劈尖干涉的特征和规律。

设劈尖的长度L=5cm 、厚度d=5m μ,入射波长λ=500nm ,介质折射率n=1,逐一改变上述参数,可以得到等厚干涉条纹强度的各种变化情况。

图3-1 是劈尖干涉实验计算机模拟的程序流程图(程序见附录Ⅰ):

图3-1 劈尖干涉实验计算机模拟的程序流程图

(1) 劈尖厚度对干涉条纹的影响

改变劈尖的厚度d ,可以得到不同干涉图样,如图3-2(a)-(c)所示,其劈尖厚

度分别为5m μ、10m μ和15m μ,从下列图中可以看出,随着劈尖厚度的增加,即劈尖的倾角增大,条纹变密,其原因由式3-1可知,

()、、、暗纹,明纹,3,2,1,0k 21k 2k 2nd 2=???+=+=ππλ? (3-1)

随着d 的增大,在相同的区域,能显示的条纹级数增多。通过这种方法,可以测定纸张的厚度、头发丝的直径等。

劈尖干涉强度分布

1000

2000

3000

4000

5000

6000

7000强度曲线分布

劈尖干涉强度分布02000400060008000

10000

强度曲线分布劈尖干涉强度分布

0500100015002000250030003500

4000强度曲线分布

(a) (b) (c)

图3-2 劈尖厚度不同的干涉图样

(2) 入射波长对干涉条纹的影响

改变入射波长λ,可以得到不同干涉图样,如图3-3(a)~(c)所示,其入射波长分别为380nm 、500nm 和780nm ,从下列图中可以看出,随着入射波长的增加,条纹逐渐变疏,其原因由式3-1可知,随着λ的增大,在相同的区域,能显示的条纹级数减少。通过这种方法,可以测定入射波的波长。

劈尖干涉强度分布

-0.4-0.200.20.4

500

1000

1500

2000

2500

3000

024

强度曲线分布劈尖干涉强度分布

050010001500200025003000

强度曲线分布劈尖干涉强度分布

-0.4-0.200.20.4050010001500200025003000

024

强度曲线分布

(a) (b) (c)

图3-3 入射波长不同的干涉图样

(3) 介质折射率对干涉条纹的影响

改变介质的折射率n ,可以得到不同干涉图样,如图3-4(a)~(c)所示,其折射率分别为1(空气)、1.133(水)和1.165(液溴),从下列图中可以看出,随着介质折射率的增加,条纹逐渐的变密,其原因由式3-1可知,随着n 的增大,在相同的区域,能显示的条纹级数增多。通过这种方法,可以测定介质的折射率。

劈尖干涉强度分布

-0.4-0.200.20.4010002000

3000

4000

50006000

024强度曲线分布劈尖干涉强度分布

050010001500200025003000

强度曲线分布劈尖干涉强度分布

-0.4-0.200.20.40100020003000400050006000

024

强度曲线分布

(a) (b) (c)

图3-4 介质折射率不同的干涉图样

3.2 牛顿环实验的模拟与分析

牛顿环干涉现象是光干涉的一个重要内容,对于入射光波长、透镜半径及透镜与平面玻璃板间距改变时牛顿环的移动方向的解释比较抽象,很难理解和接受,而具体实现这些变化的实验也比较困难。用计算机对这些现象进行模拟,能直观看到牛顿环移动方向,很容易的掌握干涉条纹形成规律,并加深了对光干涉现象的理解。

设牛顿环装置透镜半径R=2m ,入射光波长λ=6328?,建立一个大小为8mm ?8mm 具有400?400个像素点的视场,视场各点(x ,y)到视场中心的距离则可表示为22y x r +=。逐一改变程序中的透镜半径和入射光波长,可以得到干涉条纹强度的各种变化情况。

图3-5 是牛顿环实验计算机模拟的程序流程图(程序见附录Ⅱ):

图3-5 牛顿环实验计算机模拟程序流程图

根据附录Ⅱ中的程序,取k 为不同的整数,运行后观察仿真结果,得到的结论有:

(1)透镜半径对牛顿环的影响

我们清晰看到,视场中干涉条纹在不停的向外扩张,并随透镜半径的增大而展得越宽,条纹间距也随之越大。图3-6为透镜半径分别为R =2m ,R =3m ,R =4m 时取得的图样,从图形看出,视场范围内可以观察到的干涉条暗纹级数分别为a k =12,b k =8,c k =6,充分反映了波长不变情况下,距接触中心处为r 的某点的干涉暗纹的级次与透镜半径成反比。

(a) R=2m (b) R=3m (c) R=4m

图3-6 视场中牛顿环随透镜半径变化图 (2)入射波长对牛顿环的影响

保持透镜半径R =4m 不变,改变入射光波长λ,获取波长从400nm 分变化到700nm 的每个画面帧。可以看到,视场中牛顿环也在不停的向外扩张,条纹变得越来越宽,条纹间距也随之越大。图3-7为波长分别为400 nm ,500 nm ,600 nm 取得的图样,视场中4 mm ?4 mm 范围内能够看到的干涉暗纹级数分别为a k =20,b k =16,c k =13,同样的反映了视场范围大小一定的情况下,能够看到的最大的干涉暗纹级次与波长成反比。同时可以利用公式λkR r =暗计算某一

级次暗纹的半径。

(a)λ=400 nm (b)λ=500 nm (c)λ=600 nm

图3-7 视场中牛顿环随波长变化图

(3)透镜与玻璃平面板间距变化对干涉圆环的影响

假设透镜半径R =2 m ,入射光波长λ=400nm ,透镜与平面玻璃板的距离逐渐拉开,透镜中心与平面玻璃板的距离由0逐渐增加到1200nm ,获取过程中的每个画面帧。可以看到,视场中干涉圆环不停的收缩,当距离每拉开半个波长,牛顿环中心处“陷入”一个圆环,当计算机运行完一个循环时,先后看到有六条牛顿环被“陷入”。当透镜中心与平面玻璃板距离L=100mm 时,干涉条纹中心出现亮纹,说明中心处相遇的的两列光相位差为πk 2(k 取自然数),出现了相长干涉;当距离L=200mm 时,干涉条纹中心又出现暗纹,说明中心处相遇的的两列光相位差为()π12+k (k 取非负整数),出现了相消干涉。这与光的干涉理论完全相符。

3.3 迈克尔逊干涉实验的模拟与分析

迈克耳逊干涉仪是利用光的分振幅干涉原理制成的一种精密光学仪器,在科学技术上有着非常广泛的重要应用,同时也是大学物理及实验教学的重要内容之

一。根据薄膜干涉原理分析了迈克耳孙干涉仪实验中各级等倾干涉环的产生条件,得到了接收屏上的光强分布与入射角的关系。在此基础上采用计算机模拟的方法,用Matlab 编写程序,对迈克耳孙干涉仪等倾干涉条纹的图样进行了模拟和仿真。

用钠光灯作为单色光源,其波长λ=5893?,设张角θ=0.15rad ,透镜到接受

屏的距离D=0.2m ,1M 和'2M 之间的等效空气膜厚度4k d λ=,改变参数k 的取

值,可以得到等厚干涉条纹强度的各种变化情况。

图3-8 是迈克尔逊干涉实验计算机模拟的程序流程图(程序见附录Ⅲ):

图3-8 迈克尔逊干涉实验计算机模拟程序流程图

根据附录Ⅲ中的程序,取k为不同的整数,运行后观察仿真结果,得到的结论有:

(1) 干涉图样是一组明暗相间的圆环,

(2) 当k取偶数时,中心为一亮斑;k取奇数时,中心为一暗斑。图4-1(a)中干涉图样对应k=3000,图4-2(b)中干涉图样对应k=3001。

(a) (b)

图4-1 干涉环中心的明暗和空气厚度的关系 (3) k 越大,即薄膜越厚时,条纹越窄、越密。图4-2(a)中干涉图样对应k=2000;图4-2(b)中干涉图样对应k=4000;图4-2(c)中干涉图样对应k=6000。

(a) (b) (c)

图4-2 干涉环的疏密和空气膜厚度的关系 光线“2”也可看做是从2M 在半透明铬层中的虚像'2M 反射回来的,在研究

干涉时,'2M 与2M 是等效的;当1M 和'2M 互相平行时(即2M 和1M 恰好互相垂

直)将观察到圆形条纹(等倾干涉);当1M 和'2M 交很小角度时,将观察到直

线形干涉仪条纹(等厚干涉)。在实际应用中,主要利用圆形和直线形干涉条纹。

3.4 薄膜干涉实验仿真结果对比分析

利用MATLAB 提供的强大的图形图像处理功能,可以轻松实现各种二维图形的绘制,对计算机模拟得出的仿真图形进行分析,比较所得结果与理论实验是否相符合;如有不符,讨论引起差别的原因。

(1) 劈尖干涉与牛顿环干涉比较分析

迈克尔逊干涉仪产生的是等倾干涉条纹,条纹的明暗变化,和入射角度有关,相同入射角的位置干涉条纹明暗情况一致,条纹间距,条纹粗细都不等,影响条

纹干涉变化的主要原因是光源入射角度的问题。牛顿环是等厚干涉条纹,条纹的明暗变化,和上下表面中间加的空气厚度有关,相同厚度的位置,干涉条纹明暗一致(迈克耳逊干涉仪等倾干涉形成的空气劈尖是上下表面平行的),条纹间距和粗细也都不等,但是两者条纹粗细,条纹间隔计算公式完全不一样。另外牛顿环中心是零级干涉,迈克耳逊干涉仪中心是最大级干涉。

迈克尔逊干涉仪所观察到的是单色点光源入射时形成的等倾干涉,而牛顿环干涉是等厚干涉,区别在于空气薄膜的厚度,等倾干涉中薄膜厚度是均匀的,而等厚干涉中薄膜厚度是有变化的。如果在迈克尔逊干涉中调节M1与M2不是相互垂直的,空气的厚度就会不均匀,从而产生等厚干涉。迈克耳逊干涉仪干涉条纹宽度不一,干涉级次中心最大,边缘最小,牛顿环干涉条纹宽度几乎一致,干涉级次中心最小,边缘最大,迈克耳逊干涉仪通过调节镜子,中心亮暗间隔交替出现,牛顿环中心始终是暗纹。

(2) 迈克尔逊干涉与牛顿环干涉比较分析

牛顿环干涉条纹为不等间距的内粗外细的一组同心圆。光学劈尖干涉条纹为一组等间距的平行条纹。

相同点:都是薄膜干涉,且这个薄膜一般都是空气。

不同点:劈尖干涉,形成的是等厚的直条纹,且劈尖夹角θ越大,条纹越密;顿环,形成的同心干涉环,中心必是暗点,且曲率半径R 越大,环越稀疏若观察到光的稳定的干涉条纹,都必须采用分光的方法得到两列相干光源。劈尖是利用契状空气薄膜来分光的,因同一条条纹对应的空气膜的厚度相等(又称等厚线)。

(3) 劈尖干涉、牛顿环干涉和迈克尔逊干涉的比较分析

在迈克尔逊干涉仪实验中,使用He -Ne 激光器做光源,当1M 与'2M 严格平

行时,出现等倾干涉,干涉图样为同心圆环,1M 于'2M 拉开时,圆环涌出,变

小;反之,陷入,变大。当1M 与'2M 有微小夹角时,出现劈尖干涉,1M 与'2M 由

远及近时,干涉花样由圆弧趋向直线,继续移动1M 时,干涉花样由直线变成圆

弧,当1M 与2M 的交线在视场中间时,近似等厚干涉,干涉图样为直线。

计算机仿真与建模实验报告

中南大学 计算机仿真与建模 实验报告 题目:理发店的服务过程仿真 姓名:XXXX 班级:计科XXXX班 学号:0909XXXX 日期:2013XXXX

理发店的服务过程仿真 1 实验案例 (2) 1.1 案例:理发店系统研究 (2) 1.1.1 问题分析 (3) 1.1.2 模型假设 (3) 1.1.3 变量说明 (3) 1.1.4 模型建立 (3) 1.1.5 系统模拟 (4) 1.1.6 计算机模拟算法设计 (5) 1.1.7 计算机模拟程序 (6) 1实验案例 1.1 案例:理发店模拟 一个理发店有两位服务员A和B顾客随机地到达该理发店,每分钟有一个顾客到达和没有顾客到达的概率均是1/2 , 其中60%的顾客理发仅用5分钟,另外40%的顾客用8分钟. 试对前10分钟的情况进行仿真。 (“排队论”,“系统模拟”,“离散系统模拟”,“事件调度法”)

1.1.1 问题分析 理发店系统包含诸多随机因素,为了对其进行评判就是要研究其运行效率, 从理发店自身利益来说,要看服务员工作负荷是否合理,是否需要增加员工等考 虑。从顾客角度讲,还要看顾客的等待时间,顾客的等待队长,如等待时间过长 或者等待的人过多,则顾客会离开。理发店系统是一个典型的排队系统,可以用 排队论有关知识来研究。 1.1.2 模型假设 1. 60%的顾客只需剪发,40%的顾客既要剪发,又要洗发; 2. 每个服务员剪发需要的时间均为5分钟,既剪发又洗发则花8分钟; 3. 顾客的到达间隔时间服从指数分布; 4. 服务中服务员不休息。 1.1.3 变量说明 u :剪发时间(单位:分钟),u=5m ; v: 既剪发又理发花的时间(单位:分钟),v=8m ; T : 顾客到达的间隔时间,是随机变量,服从参数为λ的指数分布,(单位: 分钟) T 0:顾客到达的平均间隔时间(单位:秒),T 0=λ 1; 1.1.4 模型建立 由于该系统包含诸多随机因素,很难给出解析的结果,因此可以借助计算机 模拟对该系统进行模拟。 考虑一般理发店的工作模式,一般是上午9:00开始营业,晚上10:00左 右结束,且一般是连续工作的,因此一般营业时间为13小时左右。 这里以每天运行12小时为例,进行模拟。 这里假定顾客到达的平均间隔时间T 0服从均值3分钟的指数分布, 则有 3小时到达人数约为603 603=?人, 6小时到达人数约为1203 606=?人, 10小时到达人数约为2003 6010=?人, 这里模拟顾客到达数为60人的情况。 (如何选择模拟的总人数或模拟总时间)

关于计算机实验报告的参考范文

关于计算机实验报告的参考范文 篇一 一、实验题目 文件和文件夹的管理 二、实验目的 1.熟悉Windows XP的文件系统。 2.掌握资源管理器的使用方法。 3.熟练掌握在Windows XP资源管理器下,对文件(夹)的选择、新建、移动、复制、删除、重命名的操作方法。 三、实验内容 1.启动资源管理器并利用资源管理器浏览文件。 2.在D盘创建文件夹 3.在所创建文件夹中创建Word文件。 4.对所创建文件或文件夹执行复制、移动、重命名、删除、恢复、创建快捷方式及设置共享等操作。 四、实验步骤 (一)文件与文件夹管理 1.展开与折叠文件夹。右击开始,打开资源管理器,在左窗格中点击“+”展开,点击“—”折叠 2.改变文件显示方式。打开资源管理器/查看,选择缩略、列表,排列图标等

3.建立树状目录。在D盘空白处右击,选择新建/文件夹,输入经济贸易学院,依次在新建文件夹中建立经济类1103 4..创建Word并保存。打开开始/程序/word,输入内容。选择文件/另存为,查找D盘/经济贸易学院/1103班/王帅,单击保存 5.复制、移动文件夹 6.重命名、删除、恢复。右击文件夹,选择重命名,输入新名字;选择删除,删除文件 7.创建文件的快捷方式。右击王帅文件夹,选择发送到/桌面快捷方式 8.设置共享文件。右击王帅,选择属性/共享/在网络上共享这个文件/确定 9.显示扩展名。打开资源管理器/工具/文件夹选项/查看/高级设置,撤销隐藏已知文件的扩展名 (二)控制面板的设置。 1.设置显示属性。右击打开显示属性/桌面、屏幕保护程序 2.设置鼠标。打开控制面板/鼠标/按钮(调整滑块,感受速度)、指针 3.设置键盘。打开控制面板/键盘/速度(调整滑块,感受速度)、硬件 4.设置日期和时间打开控制面板/日期和时间

计算机实验报告范文.doc

计算机实验报告范文 课程:大学计算机基础班级:***** 学号:** 姓名:*** 组别:同组者姓名:仪器编号:实验日期: 实验windows 2000 操作 [实验目的] 1. 掌握windows 2000的启动和退出。熟悉windows 2000的桌面环境,掌握“回收站”、任务栏、帮助菜单的使用。 2. 掌握windows 2000的窗口、菜单以及对话框的操作。掌握“资源管理器”和“我的电脑”的使用。 3. 掌握文件与文件夹的创建、移动、复制等基本操作。 4. 掌握应用程序的安装与删除、移动与退出,快捷方式的创建与删除。 5. 掌握windows 2000系统的设置,了解windows2000 附件的使用。 [实验环境] 硬件:pentium 以上的计算机。 软件:windows2000 操作系统。 [实验内容]

见附件 [实验结果] 1.建立了如下图所示目录树:D:\ USER NEW1.COD A USER2 B BBB NEW2.DOC BBB 2.完成了“计算机”快捷方式的创建。 3.完成了控制面板中显示、区域选项等属性的设置。 实验指导教师(签名)实验成绩实验报告批改日期: 实验内容: 1.启动资源管理器选中D盘在右窗格空白处单击鼠标右键选择“新建/文件夹”命名为UESR 双击打开该文件夹(在当前目录下重复上面的操作,分别创建USER、A、B、C、USER2文件夹)双击打开USER1文件夹右击空白处选择“新建/WORD 文档”命名为“OLD1.DOC”,重复该操作在相应的文件夹内分别创建“OLD2.DOC”,“TOM.DOC”。 2.(1)打开B文件夹选中BBB 单击右键后选择“复制”命令打开USER文件夹在空白处单击右键后选择“粘贴”命令。

通信对抗原理大作业题目

通信对抗原理仿真大作业题目 基本要求:仿真大作业分组完成,每个组3~5人,至少选择4个题目,并且在每一类中至少选择一个题目。利用MATLAB完成计算机仿真,并且撰写仿真实验报告。大作业完成情况将作为评价平时成绩的依据。 第一类:测频方法仿真 1.FFT法数字测频技术仿真。仿真模拟通信信号或者数字通信信号三种以上, 基于FFT法进行载波频率测量。画出信号的时域、频域波形,给出FFT法测量的结果。进一步在0-20dB信噪比条件给出不同信噪比下的测量曲线,分析信噪比与测量误差的关系。 2.互相关法数字测频技术仿真。仿真模拟通信信号或者数字通信信号三种以上, 基于互相关法进行载波频率测量。画出信号的时域、频域波形,给出测量结果。进一步在0-20dB信噪比条件给出不同信噪比下的测量曲线,分析信噪比与测量误差的关系。 3.相位差分法数字测频技术仿真。仿真模拟通信信号或者数字通信信号三种以 上,基于相位差分法法进行载波频率测量。画出信号的时域、频域波形,给出测量结果。进一步在0-20dB信噪比条件给出不同信噪比下的测量曲线,分析信噪比与测量误差的关系。 第二类:测向方法仿真 4.相位干涉仪测向方法仿真。仿真模拟通信信号或者数字通信信号两种以上, 基于相位干涉仪测向方法,对不同方向到达的通信信号进行测向。画出信号的时域、频域波形,给出到达方向测量结果。进一步在0-20dB信噪比条件给出不同信噪比下的测量曲线,分析信噪比、到达角与测量误差的关系。 5.到达时差测向方法仿真。仿真模拟通信信号或者数字通信信号两种以上,基 于到达时差测向方法,对不同方向到达的通信信号进行测向。画出信号的时域、频域波形,观察相关函数,给出测量结果。进一步在0-20dB信噪比条件给出不同信噪比下的测量曲线,分析信噪比、到达角与测量误差的关系。6.多普勒测向方法仿真。仿真模拟通信信号或者数字通信信号两种以上,基于 多普勒测向方法,对不同方向到达的通信信号进行测向。画出信号的时域、频域波形,给出测量结果。进一步在0-20dB信噪比条件给出不同信噪比下的测量曲线,分析信噪比、到达角与测量误差的关系。 7.沃森-瓦特测向方法仿真。仿真模拟通信信号或者数字通信信号两种以上,基 于沃森-瓦特测向方法,对不同方向到达的通信信号进行测向。画出信号的时域、频域波形,给出测量结果。进一步在0-20dB信噪比条件给出不同信噪比下的测量曲线,分析信噪比、到达角与测量误差的关系。 第三类:信号处理技术仿真 8.信号带宽和幅度测量方法仿真。仿真模拟通信信号或者数字通信信号三种以 上,基于FFT法进行信号带宽、信号相对幅度测量。画出信号的时域、频域

计算机仿真实验-基于Simulink的简单电力系统仿真

实验七 基于Simulink 的简单电力系统仿真实验 一. 实验目的 1) 熟悉Simulink 的工作环境及SimPowerSystems 功能模块库; 2) 掌握Simulink 的的powergui 模块的应用; 3) 掌握发电机的工作原理及稳态电力系统的计算方法; 4)掌握开关电源的工作原理及其工作特点; 5)掌握PID 控制对系统输出特性的影响。 二.实验内容与要求 单机无穷大电力系统如图7-1所示。平衡节点电压0 44030 V V =∠? 。负荷功率10L P kW =。线路参数:电阻1l R =Ω;电感0.01l L H =。发电机额定参数:额定功率100n P kW =;额定电压440 3 n V V =;额定励磁电流 70 fn i A =;额定频率50n f Hz =。发电机定子侧参数:0.26s R =Ω, 1 1.14 L mH =,13.7 md L mH =,11 mq L mH =。发电机转子侧参数:0.13f R =Ω,1 2.1 fd L mH =。发电机阻尼绕组参数:0.0224kd R =Ω, 1 1.4 kd L mH =,10.02kq R =Ω,11 1 kq L mH =。发电机转动惯量和极对数分别 为224.9 J kgm =和2p =。发电机输出功率050 e P kW =时,系统运行达到稳态状态。在发电机输出电磁功率分别为170 e P kW =和2100 e P kW =时,分析发电机、平衡节点电源和负载的电流、电磁功率变化曲线,以及发电机转速和功率角的变化曲线。

G 发电机节点 V 负 荷 l R l L L P 图 7.1 单机无穷大系统结构图 输电线路 三.实验步骤 1. 建立系统仿真模型 同步电机模块有2个输入端子、1个输出端子和3个电气连接端子。模块的第1个输入端子(Pm)为电机的机械功率。当机械功率为正时,表示同步电机运行方式为发电机模式;当机械功率为负时,表示同步电机运行方式为电动机模式。在发电机模式下,输入可以是一个正的常数,也可以是一个函数或者是原动机模块的输出;在电动机模式下,输入通常是一个负的常数或者是函数。模块的第2个输入端子(Vf)是励磁电压,在发电机模式下可以由励磁模块提供,在电动机模式下为一个常数。 在Simulink仿真环境中打开Simulink库,找出相应的单元部件模型,构造仿真模型,三相电压源幅值为4403,频率为50Hz。按图连接好线路,设置参数,建立其仿真模型,仿真时间为5s,仿真方法为ode23tb,并对各个单元部件模型的参数进行修改,如图所示。

计算机网络实验报告 答案讲解

计算机网络实验报告 专业计算机科学与技术 班级计102 学号109074057 姓名王徽军 组号一组D 指导教师毛绪纹 安徽工业大学计算机学院 二○一二年十二月

目录 实验总体说明 (3) 实验一以太网帧的构成 (3) 实验三路由信息协议RIP (8) 实验四传输控制协议TCP (10) 实验五邮件协议SMTP、POP3、IMAP (12) 实验六超文本传输协议HTTP (14)

实验总体说明 1.实验总体目标 配合计算机网络课程的教学,加强学生对计算机网络知识(TCP/IP协议)的深刻理解,培养学生的实际操作能力。 2.实验环境 计算机网络协议仿真实验室: 实验环境:网络协议仿真教学系统(通用版)一套 硬件设备:服务器,中心控制设备,组控设备,PC机若干台 操作系统:Windows 2003服务器版 3.实验总体要求 ●按照各项实验内容做实验,记录各种数据包信息,包括操作、观察、记录、分析, 通过操作和观察获得直观印象,从获得的数据中分析网络协议的工作原理; ●每项实验均提交实验报告,实验报告的内容可参照实验的具体要求,但总体上应包 括以下内容:实验准备情况,实验记录,实验结果分析,算法描述,程序段,实验过程中遇到的问题以及对思考问题的解答等,实验目的、实验原理、实验步骤不需要写入实验报告中。 实验一以太网帧的构成 实验时间:_____________ 成绩:________________ 实验角色:_____________ 同组者姓名:______________________________

练习一:领略真实的MAC帧 q....U 00000010: 85 48 D2 78 62 13 47 24 58 25 00 00 00 00 00 00 .H襵b.G$X%...... 00000020: 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 ................ 00000030: 00 00 00 00 00 00 00 00 00 00 00 00 ............ 练习二:理解MAC地址的作用 ●记录实验结果 表1-3实验结果 本机MAC地址源MAC地址目的MAC地址是否收到,为什么 主机B 8C89A5-7570BB 8C89A5-757113 8C89A5-7570C1 是,主机A与主机B接在同一共享模块 主机D 8C89A5-771A47 8C89A5-757113 8C89A5-7570C1 是,主机C与主机D接在同一共享模块 主机E 8C89A5-757110 无无否,与主机A、C都不在同一共享模块 主机 F 8C89A5-7715F8 无无否,与主机A、C都不在同一共享模块 练习三:编辑并发送MAC广播帧 ●结合练习三的实验结果,简述FFFFFF-FFFFFF作为目的MAC地址的作用。 答:该地址为广播地址,作用是完成一对多的通信方式,即一个数据帧可发送给同一网段内的所有节点。 练习四:编辑并发送LLC帧 ●实验结果 帧类型发送序号N(S)接受序号N(R) LLC 001F 0 ●简述“类型和长度”字段的两种含义 答:一是如果字段的值小于1518,它就是长度字段,用于定义下面数据字段的长度;二是如果字段的值大于1536,用于定义一个封装在帧中的PDU分组的类型。 思考问题: 1.为什么IEEE802标准将数据链路层分割为MAC子层和LLC子层? 答:出于厂商们在商业上的激烈竞争,IEEE的802委员会未能形成一个统一的、最佳的局域网标准,而是被迫制定了几个不同标准,如802.4令牌总线网、802.5令牌环网等。为了使数据链路层能更好地适应多种局域网标准,802委员会就将局域网的数据链路层拆成两个子层,即逻辑链路控制

江南大学数媒0902基于虚拟现实技术大作业报告

课程:虚拟现实题目:沸腾的水壶 班级:数媒0902 学号:0305090206 姓名:沈玉婷 日期:2012.12

1、绪论 1.1 虚拟现实动画简介 虚拟现实动画就是用虚拟现实的技术以动画的形式表现出来(这是建立在虚拟现实及动画技术的基础上出现的)。我们以了解什么是虚拟现实及动画的意思后就能全面理解虚拟现实动画的概念。 1.2 关于虚拟现实技术 虚拟现实(Virtual Reality,简称VR;又译作灵境、幻真)是近年来出现的高新技术,也称灵境技术或人工环境。虚拟现实是利用电脑模拟产生一个三维空间的虚拟世界,提供使用者关于视觉、听觉、触觉等感官的模拟,让使用者如同身历其境一般,可以及时、没有限制地观察三度空间内的事物。 VR是一项综合集成技术,涉及计算机图形学、人机交互技术、传感技术、人工智能等领域,它用计算机生成逼真的三维视、听、嗅觉等感觉,使人作为参与者通过适当装置,自然地对虚拟世界进行体验和交互作用。使用者进行位置移动时,电脑可以立即进行复杂的运算,将精确的3D世界影像传回产生临场感。该技术集成了计算机图形(CG)技术、计算机仿真技术、人工智能、传感技术、显示技术、网络并行处理等技术的最新发展成果,是一种由计算机技术辅助生成的高技术模拟系统。概括地说,虚拟现实是人们通过计算机对复杂数据进行可视化操作与交互的一种全新方式,与传统的人机界面以及流行的视窗操作相比,虚拟现实在技术思想上有了质的飞跃。 2、需求分析 随着CAD技术的发展,人们就开始研究立体声与三维立体显示相结合的计算机系统。目的在于建立一种新的用户界面,使用户可以置身于计算机所表示的三维空间资料库环境中,并可以通过眼、手、耳或特殊的空间三维装置在这个环境中"环游",创造出一种"亲临其境"的感觉。 虚拟现实是人们通过计算机对复杂数据进行可视化、操作以及实时交互的环境。与传统的计算机人――机界面(如键盘、鼠标器、图形用户界面以及流行的Windows等)相比,虚拟现实无论在技术上还是思想上都有质的飞跃。传统的人――机界面将用户和计算机视为两个独立的实体,而将界面视为信息交换的媒介,由用户把要求或指令输入计算机,计算机对信息或受控对象作出动作反馈。虚拟现实则将用户和计算机视为一个整体,通过各种直观的工具将信息进行可视化,形成一个逼真的环境,用户直接置身于这种三维信息空间中自由地使用各种信息,并由此控制计算机。目前,虚拟现实技术已经遍布我们生活中的每一个行业,城市规划中的应用、旅游景观的应用、医学中应用、娱艺教中的应用、军事与航天中的应用、室内设计中的应用、房产开发中的应用、工业仿真中的应用、应急推演中的应用。由此可知,虚拟

计算机仿真实验

计算机仿真实验报告 专业:电气工程及其自动化班级:09电牵一班学号:22 姓名:饶坚指导老师:叶满园实验日期:2012年4月30日 一、实验名称 三相桥式SPWM逆变电路仿真 二、目的及要求 1.了解并掌握三相逆变电路的工作原理; 2.进一步熟悉MA TLAB中对Simulink的使用及构建模块; 3.掌握SPWM原理及构建调制电路模块; 4.复习在Figure中显示图形的程序编写和对图形的修改。 三、实验原理与步骤、电路图 1、实验原理图

2、电路原理(采用双极性控制方式) U、V和W三相的PWM控制通常公用一个三角波载波Uc,三相的调制信号Uru、Urv和Urw依次相差120°。 电路工作过程(U相为例):当Uru>Uc时,上桥臂V1导通,下桥臂V4关断,则U相相对于直流电源假想中点N’的输出电压Uun’=Ud/2。当Uru

对电路模型进行封装如下图示:

其中Subsystem1为主电路,Subsystem2为负载,Subsystem3为检测电路,Subsystem4为输入信号,Subsystem5为调制电路,Scope 为示波器,Repeating Sequence为三角载波。 各子系统电路分别如下所示: Subsystem1 Subsystem2 Subsystem3

计算机网络实验报告(以太网帧格式分析)

计算机网络实验报告 学院计算机与通信工程学院专业网络工程班级1401班 学号20姓名实验时间:2016.5.13 一、实验名称: FTP协议分析实验 二、实验目的: 分析FTP 报文格式和FTP 协议的工作过程,同时学习 Serv-U FTP Server服务软件的基本配置和FTP 客户端命令的使用。 三、实验环境: 实验室局域网中任意两台主机PC1,PC2。 四、实验步骤及结果: 步骤1:查看实验室PC1和PC2的IP地址,并记录,假设PC1的IP 地址为10.64.44.34,PC2的IP地址为10.64.44.35。 步骤2:在PC1上安装Serv-U FTP Server,启动后出现图1-20所示界面。 点击新建域,打开添加新建域向导,完成如下操作。 添加域名:https://www.360docs.net/doc/af13595887.html,;设置域端口号:21(默认);添加域IP地址:10.28.23.141;设置密码加密模式:无加密,完成后界面如图1-21所示。 完成上述操作后,还需要创建用于实验的用户帐号。点击图1.20中

浮动窗口中的“是”按钮,打开添加新建用户向导:添加用户名:test1;添加密码:123;设置用户根目录(登陆文件夹);设置是否将用户锁定于根目录:是(默认);访问权限:只读访问,完成后界面如图1-22所示。 新建的用户只有文件读取和目录列表权限,为完成实验内容,还需要为新建的用户设置目录访问权限,方法为点击导航——〉目录——〉目录访问界面,然后点击添加按钮, 按照图1-23所示进行配置。 步骤3:在PC1 和PC2 上运行Wireshark,开始捕获报文。 步骤4:在PC2 命令行窗口中登录FTP 服务器,根据步骤2中的配置信息输入用户名和口令,参考命令如下: C:\ >ftp ftp> open To 10.28.23.141 //登录ftp 服务器 Connected to 10.28.23.141 220 Serv-U FTP Server v6.2 for WinSock ready... User(none): test1 //输入用户名 331 User name okay, need password. Password:123 //输入用户密码 230 User logged in, proceed. //通过认证,登录成功

武汉理工大学-计算机仿真实验作业答案

五、(10分)已知系统的传递函数为6 168682)(232+++++=s s s s s s G 。 语言建立系统传递函数模型,并求: ⑴ 该系统的单位阶跃响应;(2分) ⑵ 输入函数为u(t)时的响应;(3分) (u(t)正弦信号,周期2秒,仿真时间8秒,采样周期0.1); (3) 输入函数为u(t)时的响应;(3分) (u(t)方波输入信号,周期10秒,仿真时间20秒,采样周期0.05) (4) 绘出系统的波德图(Bode )。(2分) 解答: num=[2 8 6]; den=[1 8 16 6]; sys=tf(num,den); t=0:0.1:8; y1=step(sys,t); u=sin(t*pi); y2=lsim(sys,u,t); subplot(2,2,1);plot(t,y1); grid; title('阶跃响应曲线'); xlabel('响应时间'); ylabel('响应值'); hold on; subplot(2,2,2);plot(t,y2); grid on ; title('对sin(t)的响应曲线'); xlabel('响应时间'); ylabel('响应值'); t=0:0.05:20 u=square(pi/5*t) y3=lsim(sys,u,t); subplot(2,2,3);plot(t,y3) grid on ; title('对方波信号的响应曲线'); xlabel('响应时间'); ylabel('响应值'); subplot(2,2,4);bode(sys);

grid ;title('bode 图'); 运行结果: 六、(10分)设二阶动力学系统的传递函数如下,假设将无阻尼固有频 率固定为ωn =1 rad/s ,将阻尼比的值分别设置成ζ=0,0.1,0.2,0.3,…, MATLAB 语言编程,分析在这些阻尼比ζ的取值下该系统的阶跃响应。 2222)(n n n s s s G ω?ωω++= 解答:wn=1; kesi=[0:0.1:1,2,3,4,5]; figure('color',[1 1 1]); hold on for i=kesi num=wn.^2 den=[1,2*i*wn,wn.^2]; step(num,den);

计算机仿真实验报告7

山东工商学院计算机仿真及应用实验报告 实验七 MATLAB的基本应用(二)及Simulink仿真 (验证性实验) 学院: 专业班级: 实验时间: 学号: 姓名:

一、实验目的 1、掌握连续信号的仿真和傅里叶分析方法 2、掌握连续系统的分析方法(时域分析法,拉氏变换法和傅里叶分析法); 3、掌握离散信号的仿真和分析运算方法 4、掌握离散系统的分析方法(时域分析法); 5、掌握符号运算方法; 6、掌握Simulink仿真工具; 二、实验原理 1、连续信号的仿真和分析法,参考教材第6.1节,重点: 单位冲激信号的仿真方法;单位阶跃信号的仿真方法;复指数信号的仿真方法 2、连续系统的分析方法,参考教材第6.1节,重点: 例6.2,LTI系统的零输入响应的求解方法; 例6.3,LTI系统的冲激响应的求解方法 例6.5,LTI系统的零状态响应的求解方法 例6.6,系统中有重极点时的计算 3、系统的频域分析方法,参考教材第6.2节,重点: 例6.7,方波分解为多次正弦波之和 例6.8:全波整流电压的频谱 例6.10:调幅信号通过带通滤波器 例6.12:用傅里叶变换计算滤波器的响应和输出 4、离散信号的仿真和分析法,参考教材第6.3节,7.1节,重点: 单位脉冲序列impseq,单位阶跃序列stepseq 例7.1:序列的相加和相乘 例7.2:序列的合成与截取 例7.3:序列的移位和周期延拓运算 三、实验内容(包括内容,程序,结果) 以自我编程练习实验为主,熟悉各种方法和设计,结合课堂讲授,实验练习程序代码。 1、根据教材第6.1节的内容,练习连续信号和系统的时域分析和拉氏变换方法。 q602 clear,clc a=input('输入分母系数向量a=[a1,a2,...]= '); n=length(a)-1; Y0=input('输入初始条件向量Y0=[y0,Dy0,D2y0,...]= '); p=roots(a);V=rot90(vander(p));c=V\Y0'; dt=input('dt= ');tf=input('tf= '); t=0:dt:tf;y=zeros(1,length(t));

计算机仿真期末大作业Mersenne Twister随机数发生器及随机性测试

Mersenne Twister随机数发生器及随机性测试 一、实验目的 用MATLAB实现Mersenne Twister随机数发生器,并对其随机性进行测试。二、实验原理 伪随机数的产生,首先是选取种子,然后是在此种子基础上根据具体的生成算法计算得到一个伪随机数,然后利用此伪随机数再根据生成算法递归计算出下二个伪随机数,直到将所有不重复出现的伪随机数全部计算出来。这个伪随机数序列就是以后要用到的伪随机数序列。上面的计算过程可以一次性计算完毕,也可以使用一次递归计算一次,每次生成的伪随机数就是这个伪随机数序列中的一个,不过不管怎么样,只要确定了种子,确定了生成算法,这个序列就是确定的了。所谓种子,就是一个对伪随机数计算的初始值。 Mersenne Twister算法是一种随机数产生方法,它是移位寄存器法的变种。该算法的原理:Mersenne Twister算法是利用线性反馈移位寄存器(LFSR)产生随机数的,LFSR的反馈函数是寄存器中某些位的简单异或,这些位也称之为抽头序列。一个n位的LFSR能够在重复之前产生2^n-1位长的伪随机序列。只有具有一定抽头序列的LFSR才能通过所有2^n-1个内部状态,产生2^n - 1位长的伪随机序列,这个输出的序列就称之为m序列。为了使LFSR成为最大周期的LFSR,由抽头序列加上常数1形成的多项式必须是本原多项式。一个n阶本原多项式是不可约多项式,它能整除x^(2*n-1)+1而不能整除x^d+1,其中d能整除2^n-1。例如(32,7,5,3,2,1,0)是指本原多项式x^32+x^7+x^5+x^3+x^2+x+1,把它转化为最大周期LFSR就是在LFSR小邓第32,7,5,2,1位抽头。利用上述两种方法产生周期为m的伪随机序列后,只需要将产生的伪随机序列除以序列的周期,就可以得到(0,1)上均匀分布的伪随机序列了。 伪代码如下: // 建立624位随机序列数组 int[0..623] MT int index = 0 //初始化随机序列数组 function initializeGenerator(int seed) { MT[0] := seed for i from 1 to 623 { MT[i] := last 32 bits of(1812433253 * (MT[i-1] xor(right shift by 30 bits(MT[i-1]))) + i) // 0x6c078965 } }

电路计算机仿真实验报告

电路计算机仿真分析 实验报告

实验一直流电路工作点分析和直流扫描分析 一、实验目的 1、学习使用Pspice软件,熟悉它的工作流程,即绘制电路图、元件类别的选择及其参数的赋值、分析类型的建立及其参数的设置、Probe窗口的设置和分析的运行过程等。 2、学习使用Pspice进行直流工作点分析和直流扫描分析的操作步骤。 二、原理与说明 对于电阻电路,可以用直观法(支路电流法、节点电压法、回路电流法)列写电路方程,求解电路中各个电压和电流。PSPICE软件是采用节点电压法对电路进行分析的。 使用PSPICE软件进行电路的计算机辅助分析时,首先在capture环境下编辑电路,用PSPICE 的元件符号库绘制电路图并进行编辑、存盘。然后调用分析模块、选择分析类型,就可以“自 动”进行电路分析了。需要强调的是,PSPICE软件是采用节点电压法“自动”列写节点电 压方程的,因此,在绘制电路图时,一定要有参考节点(即接地点)。此外,一个元件为一 条“支路”(branch),要注意支路(也就是元件)的参考方向。对于二端元件的参考方向定 义为正端子指向负端子。 三、示例实验 应用PSPICE求解图1-1所示电路个节点电压和各支路电流。 图1-1 直流电路分析电路图

4.000V R2 1 2.000A 0V Idc2 4Adc 4.000A 6.000V R1 1 4.000A Idc1 2Adc 2.000A R3 3 2.000A 图1-2 仿真结果 四、选做实验 1、实验电路图 (1)直流工作点分析,即求各节点电压和各元件电压和电流。 (2)直流扫描分析,即当电压源Us1的电压在0-12V之间变化时,求负载电阻R L中电流I RL随电压源Us1的变化曲线。 R4 3 Is3 2Adc 0Vs2 10Vdc RL 1 Is1 1Adc Is2 1Adc R1 4 I Is5 3Adc R2 2 12Vdc IPRINT Vs3 5Vdc Vs4 7Vdc 图1-3 选做实验电路图 2、仿真结果

郑大计算机实验报告答案

习题及实验(一) 第一部分习题 一、简答题 1计算机的发展阶段: 四个发展阶段: 第一个发展阶段:1946-1956年电子管计算机的时代。1946年第一台电子计算机问世美国宾西法尼亚大 学,它由冯·诺依曼设计的。占地170平方,150KW。运算速度慢还没有人快。是计算机发展历史上的一个里程碑。(ENIAC)(electronic numerical integator and calculator)全称叫“电子数值积分和计算机”。 第二个发展阶段:1956-1964年晶体管的计算机时代:操作系统。 第三个发展阶段:1964-1970年集成电路与大规模集成电路的计算机时代 (1964-1965)(1965-1970) 第四个发展阶段:1970-现在:超大规模集成电路的计算机时代。 第一代计算机 1946 1957 电子管运算速度较低,耗电量大存储容量小。 第二代计算机 1958 1964 晶体管体积小,耗电量较少,运算速度高,价格下降。第三代计算机 1965 1971 中小规模集成电路体积功能进一步减少,可靠性及速度进一步提高。 第四代计算机 1972年至今大规模及超大规模集成电路性能到规模提高,价格大幅度降低,广泛应用于社会生活的各个领域,走进办公室和家庭 2.主要应用:计算机的应用极其广泛,早期的计算机主要体现在科学计算机,数据处理,计算机控制等几个方面.随着微型计算机的发慌和迅速普及,计算机的应用

已渗透到国民经济各个总门及社会生活的各个方面现代计算机除了传统的应用外,还应用于以下几个大方面. 1.办化自动化 2.计算机辅助系统 3.虚拟现实 4.人工智能 5.电子商务 3. 1.管理系统中的各种资源,包括硬件资源和软件资源。 1)监视资源 2)决定分配资源策略 3)分配资源 4)回收资源 2.为用户提供友好的界面。 1)命令行界面 2)图形化界面 4.操作系统大致可分为6种类型。 简单操作系统。分时系统。实时操作系统。网络操作系统。分布操作系统。智能操作系。目前微机上常见的操作系统有DOS、OS/2、UNIX、XENIX、LINUX、Windows、Netware等。 5. 系统软件,应用软件。 系统软件:用以实现计算机系统的管理、控制、运行、维护,并完成应用程序的装入、编译等任务的程序。系统软件是开发和运行应用软件的平台,系统软件的核心是操作系统。

计算机仿真实训实验报告实验1-4

实验一 熟悉MATLAB 工作环境 16电气5班 周树楠 20160500529 一、实验目的 1.熟悉启动和退出MATLAB 软件的方法。 2.熟悉MATLAB 软件的运行环境。 3.熟悉MATLAB 的基本操作。 二、实验设备及条件 计算机一台(带有MATLAB6.0以上的软件境)。 三、实验内容 1.练习下面指令: cd,clear,dir,path,help,who,whos,save,load 。 2.建立自己的工作目录MYBIN 和MYDATA ,并将它们分别加到搜索路径的前面或者后面。 3.求23)]47(*212[÷-+的算术运算结果。 4.M 文件的建立,建立M 文件,求出下列表达式的值: ?? ????-+=++=+= 545.0212),1ln(21 185sin 2222 1i x x x z e z o 其中

5.利用MATLAB的帮助功能分别查询inv、plot、max、round函数的功能和用法。 四、运行环境介绍及注意事项 1.运行环境介绍 打开Matlab软件运行环境有图1-1所示的界面

图1-1 MATLAB的用户界面 操作界面主要的介绍如下: 指令窗( Command Window ),在该窗可键入各种送给 MATLAB 运作的指令、函数、表达式,并显示除图形外的所以运算结果。 历史指令窗( Command History ),该窗记录已经运行过的指令、函数、表达式;允许用户对它们进行选择复制、重运行,以及产生 M 文件。 工作空间浏览器( Workspace Browser ),该窗口罗列出 MATLAB 工作空间中所有的变量名、大小、字节数;并且在该窗中,可对变量进行观察、编辑、提取和保存。 其它还有当前目录浏览器( Current Directory Browser )、 M 文件编辑 / 调试器(Editor/Debugger )以及帮助导航/ 浏览器(Help Navigator/Browser )等,但通常不随操作界面的出现而启动。 利用 File 菜单可方便对文件或窗口进行管理。其中 File | New 的各子菜单, M-file ( M 文件)、 Figure (图形窗口)、或 Model ( Simulink 编辑界面)分别可创建对应文件或模块。 Edit 菜单允许用户和 Windows 的剪切板交互信息。 2.在指令窗操作时应特别注意以下几点 1)所有输入的指令、公式或数值必须按下回车键以后才能执行。例如: >>(10*19+2/4-34)/2*3 (回车) ans= 234.7500 2)所有的指令、变量名称都要区分字母的大小写。 3)%作为MATLAB注释的开始标志,以后的文字不影响计算的过程。 4)应该指定输出变量名称,否则MATLAB会将运算结果直接存入默认的输出变量名ans。 5)MATLAB可以将计算结果以不同的精确度的数字格式显示,可以直接在指令视窗键入不同的数字显示格式指令。例如:>>format short (这是默认的) 6)MATLAB利用了↑↓二个游标键可以将所输过的指令叫回来重复使用。按下↑则前一次输入的指令重新出现,之后再按Enter键,即再执行前一次的指令。

控制系统数字仿真大作业.

《控制系统数字仿真》课程 大作业 姓名: 学号: 班级: 日期: 同组人员:

目录 一、引言 (2) 二、设计方法 (2) 1、系统数学模型 (2) 2、系统性能指标 (4) 2.1 绘制系统阶跃响应曲线、根轨迹图、频率特性 (4) 2.2 稳定性分析 (6) 2.3 性能指标分析 (6) 3、控制器设计 (6) 三、深入探讨 (9) 1、比例-微分控制器(PD) (9) 2、比例-积分控制(PI) (12) 3、比例-微分-积分控制器(PID) (14) 四、设计总结 (17) 五、心得体会 (18) 六、参考文献 (18)

一、引言 MATLAB语言是当今国际控制界最为流行的控制系统计算机辅助设计语言,它的出现为控制系统的计算机辅助分析和设计带来了全新的手段。其中图形交互式的模型输入计算机仿真环境SIMULINK,为MATLAB应用的进一步推广起到了积极的推动作用。现在,MATLAB语言已经风靡全世界,成为控制系统CAD领域最普及、也是最受欢迎的软件环境。 随着计算机技术的发展和应用,自动控制理论和技术在宇航、机器人控制、导弹制导及核动力等高新技术领域中的应用也愈来愈深入广泛。不仅如此,自动控制技术的应用范围现在已发展到生物、医学、环境、经济管理和其它许多社会领域中,成为现代社会生活中不可或缺的一部分。随着时代进步和人们生活水平的提高,在人类探知未来,认识和改造自然,建设高度文明和发达社会的活动中,控制理论和技术必将进一步发挥更加重要的作用。作为一个自动化专业的学生,了解和掌握自动控制的有关知识是十分必要的。 利用MATLAB软件及其SIMULINK仿真工具来实现对自动控制系统建模、分析与设计、仿真,能够直观、快速地分析系统的动态性能和稳态性能,并且能够灵活的改变系统的结构和参数,通过快速、直观的仿真达到系统的优化设计,以满足特定的设计指标。 二、设计方法 1、系统数学模型 美国卡耐尔基-梅隆大学机器人研究所开发研制了一套用于星际探索的系统,其目标机器人是一个六足步行机器人,如图(a)所示。该机器人单足控制系统结构图如图(b)所示。 要求: (1)建立系统数学模型; (2)绘制系统阶跃响应曲线、根轨迹图、频率特性; (3)分析系统的稳定性,及性能指标; (4)设计控制器Gc(s),使系统指标满足:ts<10s,ess=0,,超调量小于5%。

控制系统计算机仿真实验报告

计算机仿真试验报告 自动化1201 ** 3120502007 [实验目的] (1).掌握采样控制系统数字仿真的特点。 (2).了解数字控制器对系统动态性能的影响。 (3).学会编制双重循环法的仿真程序。 (1). 复习采样控制系统的仿真原理及特点。 (2).根据理论分析,初步估计系统在给定条件下可能出现的动态过程。

(1).按实验目的、要求和已知条件,建立系统的Simulink模型,并且编制双重循环法的仿真程序。 1) Simulink模型建立: 根据题目给出的条件,数字控制系统的结构图如下图所示: 其中的其中数字控制器为: 根据上面结构图,所建立 2) 编制双重循环法的仿真程序 根据数字控制系统的结构图与条件(1)式,我们可以得到得到被控对象的状态空间模型:

[]112212()()0010()() ()110()()01()x t x t u t x t x t x t y t x t ?????????=+?????????-? ????????? ?? ?=??? ??? (2) 按连续系统离散相似算法将(2)式离散化。为了保证精度,其离散化时的步长h (虚 拟采样周期)应比数字控制器的实际采样周期T (=1s)小得多。为简化起见,取h=T/N=T/100=0.01T=0.01s 。 利用MATLAB 控制系统工具箱提供的将连续系统转换成离散系统的函数c2d ,把连续状态空间模型(2)变换为离散状态空间模型。 离散程序如下(程序1): clear; h=0.01; A=[0 0;1 -1]; B=[10;0]; [G,H]=c2d(A,B,h) 运行后的结果为: 即: 1 0(())0.010.99T ??==???? G Φ (3) 0.1(())0.0005T ?? ==? ? ?? H Γ (4) 故连续系统被控对象(2)的等价离散化状态方程为: []112212(1)()100.1()(1)()0.010.990.0005()()01()x k x k u k x k x k x k y k x k ?+???????? =+????????? +???? ?????? ?? ?=??? ??? (5)

计算机实验报告

实验指导 实验一 Visual C++开发环境使用 大气科学专业实验日期 4月 18日姓名:学号 1.实验目的 (1)熟悉Visual C++集成开发环境。 (2)掌握C语言程序的书写格式和C语言程序的结构。 (3)掌握C语言上机步骤,了解C程序的运行方法。 (4)能够熟练地掌握C语言程序的调试方法和步骤 2. 实验内容 输入如下程序,实现两个数的乘积。 #include ; int main() { x=10,y=20 p=prodct(x,t) printf("The product is : ",p) int prodct(int a ,int b ) int c c=a*b return c } (1)在编辑状态下照原样键入上述程序。 (2)编译并运行上述程序,记下所给出的出错信息。 (3)再编译执行纠错后的程序。如还有错误,再编辑改正,直到不出现语法错误为止。 3.分析与讨论 (1)记下在调试过程中所发现的错误、系统给出的出错信息和对策。分析讨论成功或失败的原因。(2)总结C程序的结构和书写规则。 实验心得: 通过本次实验,我了解到C语言的特点,初步认识程序设计方法和程序设计一般步骤,掌握C语言程序编译、链接和运行过程,为我进一步学好C语言打下了基础。

实验二数据类型、运算符和表达式 大气科学专业实验日期4 月 25日姓名:刘园园学号327 1.实验目的 (1)理解常用运行符的功能、优先级和结合性。 (2)熟练掌握算术表达式的求值规则。 (3)熟练使用赋值表达式。 (4)理解自加、自减运算符和逗号运算符 (5)掌握关系表达式和逻辑表达式的求值 2.实验内容 (1)整数相除 #include<> int main() { int a=5,b=7,c=100,d,e,f; d=a/b*c; e=a*c/b; f=c/b*a; printf("d=%d , e=%d ,f=%d\n",d,e,f); return 0; } (2)自加、自减运算 #include<> int main() { int a=5,b=8; printf("a++=%d\n",a++);

相关文档
最新文档