《FORTRAN数值方法及其在物理学中应用省公共课一等奖全国赛课获奖课件.pptx》由会员分享，可在线阅读，更多相关《FORTRAN数值方法及其在物理学中应用省公共课一等奖全国赛课获奖课件.pptx（76页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、 第二章 物理图形、图象与计算机模拟物理图形、图象与计算机模拟 2.1简谐振动及其合成曲线模拟简谐振动及其合成曲线模拟2.2阻尼运动和阻尼振动模拟阻尼运动和阻尼振动模拟2.3驻波模拟驻波模拟2.5波干涉和衍射图形模拟波干涉和衍射图形模拟2.4点电荷与点电荷系等势线和电场模拟点电荷与点电荷系等势线和电场模拟1第1页2.1简谐振动及其合成曲线模拟简谐振动及其合成曲线模拟一、简谐振动一、简谐振动 曲线和曲线和 曲线曲线例例1：画出：画出曲线及对应曲线及对应曲线，其中曲线，其中2第2页Implicitreal*8(a-h,o-z)open(1,file=x-t.dat)open(2,file=v-t.
2、dat)write(*,*)inputA,w,phi,Nread(*,*)A,w,phi,Npi=3.1415926do10I=1,Nt=2.*pi/wt=t*float(I)/Nx=A*cos(w*t+phi)v=A*w*cos(w*t+phi+pi/2.)write(1,*)t,x10 write(2,*)t,vend计算程序计算程序3第3页图形模拟图形模拟Origin介绍介绍4第4页二、简谐振动合成二、简谐振动合成1.同方向简谐振动合成同方向简谐振动合成 同频率情况同频率情况 其中其中 5第5页例例2：试给出两个同方向同频率简谐振动合成程序：试给出两个同方向同频率简谐振动合成程序open
3、(1,file=x1.dat)open(2,file=x2.dat)open(3,file=x.dat)write(*,*)inputA1,A2,w,phi1,phi2=?read(*,*)A1,A2,w,phi1,phi2pi=3.1415926do10I=1,1000t=2.*pi/wt=t*float(I)/1000 x1=A1*cos(w*t+phi1)x2=A2*cos(w*t+phi2)x=x1+x2write(1,*)t,x1write(2,*)t,x210write(3,*)t,xend 6第6页图形模拟图形模拟两个同方向同频率简谐振动合成两个同方向同频率简谐振动合成 7第7页
4、EX2-1:编程完成例编程完成例2。两个同方向同频率简谐振动合成两个同方向同频率简谐振动合成 8第8页 不一样频率情况不一样频率情况 若若 ，会出现，会出现拍现象拍现象。合振动不再是简谐振动，利用合振动不再是简谐振动，利用旋转矢量法旋转矢量法能够求得合振动振幅为能够求得合振动振幅为v 中间经历时间中间经历时间 称为周期，显然称为周期，显然 ，频率：频率：v 振幅在振幅在 和和 间周期性地改变，属振动调制。间周期性地改变，属振动调制。合振动振幅从一次极大到相邻另一次极大。合振动振幅从一次极大到相邻另一次极大。9第9页假设两个分振动振幅都为假设两个分振动振幅都为，圆频率，圆频率相差较小，取它们初相
5、位相差较小，取它们初相位此时合成运动位移可写成：此时合成运动位移可写成：改变主要取决于改变主要取决于，振幅按改变，振幅按改变。都是零，则能够分别表示为：都是零，则能够分别表示为：其中其中。因为圆频率因为圆频率 远大于圆频率远大于圆频率 ，10第10页图形模拟图形模拟两个同方向频率近似简谐振动合成两个同方向频率近似简谐振动合成 11第11页2.两个相互垂直方向简谐振动合成两个相互垂直方向简谐振动合成若若 ，则有合振动方程：，则有合振动方程：（一、三象限直线方程一、三象限直线方程）（二、四象限直线方程二、四象限直线方程）（椭圆方程椭圆方程）12第12页质点轨迹曲线质点轨迹曲线n下列图所表示为两个频
6、率相同、振幅相等、相互垂直而相下列图所表示为两个频率相同、振幅相等、相互垂直而相位差位差两个频率相同、振幅相等、相互垂直简谐振动合成两个频率相同、振幅相等、相互垂直简谐振动合成 为为 下质点轨迹曲线。下质点轨迹曲线。13第13页若若 ，但但满满足一定整数倍数比关系足一定整数倍数比关系时时，则则会会 下列图所表示为两个频率不一样（满足下列图所表示为两个频率不一样（满足 ）、振幅）、振幅相等、相等、利萨如图形示意图利萨如图形示意图 出现出现利萨如图形利萨如图形：相互垂直而相位差为相互垂直而相位差为0，下质点轨迹曲线。下质点轨迹曲线。14第14页2.2阻尼运动和阻尼振动模拟阻尼运动和阻尼振动模拟一、
7、一、阻尼情况下物体运动阻尼情况下物体运动 曲线曲线例例3：质量为：质量为摩托快艇以速度摩托快艇以速度行驶，它受到摩擦阻力行驶，它受到摩擦阻力与与速速度度成成正正比比，设设百百分分比比系系数数为为，则则，试试求求关关闭发动机后，对闭发动机后，对改变规律。改变规律。（取（取解解：物理分析与数学模型：物理分析与数学模型可用函数作图法可用函数作图法 方法方法1:15第15页n方法方法2：用用 （差商法）（差商法）函数近似值作图法函数近似值作图法 16第16页open(1,file=vt.dat)write(*,*)inputa,v0,t=?read(*,*)a,v0,tv1=v0t0=0.v10=v0
8、write(1,*)t0,v0,v10dt=t/1000.do10j=1,1000tt=t*float(j)/1000.v=v0*exp(-1.)*a*tt)v1=v1-a*v1*dt10write(1,*)tt,v,v1end模拟程序模拟程序17第17页二、二、阻尼振动阻尼振动n问题问题：弹簧振子阻尼振动方程为：弹簧振子阻尼振动方程为n 阻尼因子阻尼因子n 弹簧振子角频率弹簧振子角频率试用函数近似法作出位移与时间函数改变曲线。试用函数近似法作出位移与时间函数改变曲线。18第18页问题分析问题分析:n解：将解：将二阶二阶微分方程微分方程化为一阶化为一阶微分方程微分方程即即而而 19第19页实例
9、说明：实例说明：n例例4 4：画出当：画出当 ，时，时，秒下秒下 曲线。曲线。计算程序：计算程序：open(1,file=v-t.dat)open(2,file=x-t.dat)write(*,*)inputB,w0,x0,v0,t=?read(*,*)B,w0,v0,x0,tdt=t/1000.v=v0 x=x0tt0=0.0write(1,*)tt0,v0write(2,*)tt0,x0do10j=1,1000tt=float(j)*dtf=-2.*B*v-w0*2*xv=v+f*dtx=x+v*dtwrite(1,*)tt,v10write(2,*)tt,x end20第20页阻尼振动曲
10、线示意图阻尼振动曲线示意图 图形模拟图形模拟21第21页不一样阻尼情况振动曲线示意图不一样阻尼情况振动曲线示意图 图形模拟图形模拟22第22页EX2-2:编程完成例编程完成例3。EX2-3:一石子从空中静止下落，已知一石子从空中静止下落，已知式中式中为常数，试绘制石子下落为常数，试绘制石子下落曲线。其中曲线。其中。作业作业23第23页2.3驻波模拟驻波模拟 定义：定义：两列振幅、振两列振幅、振动动方向和方向和频频率都相同而率都相同而传输传输方向相方向相 反两列同反两列同类类波相干叠加形成波相干叠加形成驻驻波波。设有两列设有两列振动方向相同、振幅相同、频率相同振动方向相同、振幅相同、频率相同平面
11、余弦波，平面余弦波，按叠加原理，合成驻波波函数为按叠加原理，合成驻波波函数为:轴正、负方向传输。轴正、负方向传输。分别沿分别沿24第24页n在在值满足下式各点，值满足下式各点，振幅为零振幅为零驻波波节处驻波波节处相邻两波节距离为相邻两波节距离为半波长半波长，即：，即：讨论讨论:p 因子因子是时间是时间 余弦函数，说明形成驻波后，余弦函数，说明形成驻波后，各各质点都在作同频率谐振动质点都在作同频率谐振动。p 另一因子另一因子是坐标余弦函数，说明是坐标余弦函数，说明各质点各质点振幅按余弦函数规律分布振幅按余弦函数规律分布。25第25页在在值满足下式各点，值满足下式各点，振幅最大振幅最大驻波驻波波腹
12、波腹处处相邻两波腹间距离也为相邻两波腹间距离也为半波长半波长，即：，即：l l波节处波节处质点振动质点振动振幅为零振幅为零，一直处于静止；，一直处于静止；l l波腹处波腹处质点振动质点振动振幅最大振幅最大，等于，等于。l l其它各处质点振动振幅则在零与最大之间。其它各处质点振动振幅则在零与最大之间。l l两两相邻波节相邻波节或两或两相邻波腹相邻波腹之间之间相距半波长相距半波长。l l波腹和相邻波节波腹和相邻波节间间距离为距离为 ，波腹和波节交替作等距离排列。，波腹和波节交替作等距离排列。特征特征:讨论讨论:26第26页驻波模拟程序流图驻波模拟程序流图27第27页implicitreal*8(a
13、-h,o-z)open(1,file=zhubo.dat)open(2,file=zhubo1.dat)open(3,file=zhubo2.dat)write(*,*)inputA,x,wavelengthread(*,*)A,x,wapi=3.1415926T=2.*pitime=0.5*Tdx=x/1000do10I=1,1000 x=dx*float(i)y1=A*cos(2*pi)*(time/T-x/wa)y2=A*cos(2*pi)*(time/T+x/wa)y=2*A*cos(2*pi*x/wa)*cos(2*pi*time/T)write(1,*)x,ywrite(2,*)x
14、,y110write(3,*)x,y2end模拟程序模拟程序28第28页图形模拟图形模拟，29第29页30第30页每每一一时时刻刻，驻驻波波都都有有一一定定波波形形，此此波波形形既既不不向向右右移移，也也不不向向左左移移，各各点点以以各各自自确确定定振振幅幅在在各各自自平平衡衡位位置置附附近近振振动动，没有振动状态或相位传输，因而称为没有振动状态或相位传输，因而称为驻波驻波。又沿相反方向同时经过平衡位置。又沿相反方向同时经过平衡位置。在驻波进行过程中，在驻波进行过程中，没有振动状态（相位）和波形定向没有振动状态（相位）和波形定向传输，能够证实也无能量定向传输传输，能够证实也无能量定向传输，这也
15、是，这也是行波和驻波行波和驻波主要区分主要区分所在。所在。将相邻两波节之间各点称为一段，将相邻两波节之间各点称为一段，每一段各点每一段各点 具具有相同符号有相同符号，而，而相邻两端符号总是相反相邻两端符号总是相反，这说明在驻，这说明在驻波中波中同一段上各质点振动相位相同同一段上各质点振动相位相同，而，而相邻两段中各点相邻两段中各点振动相位相反振动相位相反。l 同一段内各点沿相同方向同时到达各自振动位移最大值，同一段内各点沿相同方向同时到达各自振动位移最大值，又沿相同方向同时经过平衡位置；又沿相同方向同时经过平衡位置；l 波节两侧各点同时沿相反方向到达振动位移正、负最大值，波节两侧各点同时沿相反
16、方向到达振动位移正、负最大值，31第31页 当当波在自由端反射时，则无相位突变波在自由端反射时，则无相位突变，形成驻波时，在，形成驻波时，在 半波损半波损失失 在弦线上进行驻波试验，在弦线上进行驻波试验，反射点反射点处弦线是固定不动，处弦线是固定不动，这一点这一点只能是波节只能是波节。这说明。这说明反射波和入射波相位在反射反射波和入射波相位在反射点恰好相反点恰好相反，即入射波在反射点反射时相位有，即入射波在反射点反射时相位有 突变。突变。依据相位差依据相位差 与波程差关系（与波程差关系（），相位差为），相位差为 说，说，反射波和入射波之间存在着半个波长波程差反射波和入射波之间存在着半个波长波程
17、差，这种，这种相位突变相位突变 称为称为半波损失半波损失。就相当于半个波长波程差，这说明对固定端反射点来就相当于半个波长波程差，这说明对固定端反射点来自由端出现波腹自由端出现波腹。32第32页EX2-4:一沿一沿方向传输入射波波函数为方向传输入射波波函数为在在处发生反射，反射点为一节点。处发生反射，反射点为一节点。求（求（1）反射波波函数及在）反射波波函数及在时图形。时图形。（2）合成波（驻波）波函数及在）合成波（驻波）波函数及在 时图形。时图形。作业作业33第33页2.4点电荷与点电荷系等势线和电场模拟点电荷与点电荷系等势线和电场模拟一、一、等势线方程等势线方程要求要求：作出满足等势线方程：
18、作出满足等势线方程 等势线。等势线。例例5：一个点电荷：一个点电荷在在处产生电势处产生电势等势线方程为：等势线方程为：即即圆心在圆心在，半径为，半径为圆。圆。由上式可得：由上式可得：画图时可采取参数方程：画图时可采取参数方程：34第34页例例6：两个点电荷电势分布。：两个点电荷电势分布。等势线方程为：等势线方程为：35第35页二、隐函数曲线绘制（二、隐函数曲线绘制（）（这里引入参数（这里引入参数）n指定指定 (如如 )，计算计算出出 和和 ，由上式，由上式可得可得 和和 。n若点若点 处电势为处电势为 ，则，则 处电势也为处电势也为 。36第36页编程步骤编程步骤:把上述点把上述点 连起来，就
19、是连起来，就是 等势线。等势线。输入输入给定给定（或（或）Step1:求求根根令令如取如取Step2:37第37页三、三、等势线作图步骤等势线作图步骤n1 1关于关于 轴对称。轴对称。2.对于给定对于给定其等势线出发点选在其等势线出发点选在轴上，轴上，即由即由给出给出 由由点出发。点出发。38第38页作图步骤作图步骤:n3.而而一样有：一样有：所以有：所以有：取：39第39页解解：令：令则有：则有：例例7：(异号异号)，画出画出时等势线。时等势线。(忽略忽略)注意若求得初若求得初值值有多个，需有多个，需对对每一个初每一个初值进值进行迭代行迭代计算才会确保曲线完整性。计算才会确保曲线完整性。40
20、第40页open(1,file=a1.dat)write(*,*)inputq1,q2,x1,x2=?read(*,*)q1,q2,x1,x2x=x*y=0.0write(1,*)x,ydt=0.005do10j=1,100r1=(x-x1)*2+y*2r1=sqrt(r1)r2=(x-x2)*2+y*2r2=sqrt(r2)dx=y*(q1*r2*3+q2*r1*3)*dtx=x+dxdy=(-1.)*(q1*(x-x1)*r2*3+q2*(x-x2)*r1*3)*dty=y+dy10write(1,*)x,yend模拟程序模拟程序41第41页图形模拟图形模拟 ,（异号异号）两个点电荷电势等
21、势线示意图两个点电荷电势等势线示意图（同号同号）42第42页q1=1,q2=1,x1=-2,x2=2 q1=-1,q2=-1,x1=-2,x2=2同号情况同号情况43第43页q1=-1,q2=1,x1=-2,x2=2 q1=1,q2=-1,x1=-2,x2=2异号情况异号情况44第44页四、点电荷系四、点电荷系电场线电场线图像模拟图像模拟单个点电荷电场：单个点电荷电场：点电荷系电场：点电荷系电场：例例8：如图所表示，：如图所表示，A、B、C三个点电荷组成一个点电荷系，其中三个点电荷组成一个点电荷系，其中A(0,0)点带电为点带电为，B(0,1)点带电亦为点带电亦为，最终，最终C(0,-1)点带
22、电点带电，试模拟，试模拟平面内电场线分布。平面内电场线分布。45第45页点电荷系在某点产生电场强度等于各点电荷单独点电荷系在某点产生电场强度等于各点电荷单独在该点产生电场强度矢量和，这称为在该点产生电场强度矢量和，这称为电场强度叠电场强度叠加原理加原理。所以，在所以，在XOY平面内电场强度为：平面内电场强度为：电场线上每一点切线方向反应了该点场强方向，有：电场线上每一点切线方向反应了该点场强方向，有：解上述方程可得：解上述方程可得：46第46页此超越方程（等值线）即为这个点电荷系电场线方程，此超越方程（等值线）即为这个点电荷系电场线方程，点电荷系电场线模拟点电荷系电场线模拟常数常数C选取不一样
23、值对应不一样电场线。选取不一样值对应不一样电场线。(C=0.13,间隔间隔0.1）47第47页例例9：电偶极振子电场模拟：电偶极振子电场模拟设电偶极子电偶极矩为：设电偶极子电偶极矩为：在球坐标系下，任意时刻在球坐标系下，任意时刻t空间任意处空间任意处r辐射电场为：辐射电场为：辐射电场线满足方程辐射电场线满足方程:48第48页解上述方程即可得电场线满足方程：解上述方程即可得电场线满足方程：当当C 取不一样值时可得到不一样辐射电场线。取不一样值时可得到不一样辐射电场线。电偶极子辐射电场线模拟电偶极子辐射电场线模拟 49第49页五、带电粒子在电磁场中运动五、带电粒子在电磁场中运动 在均匀磁场中带电粒
24、子受洛伦茨力作用在均匀磁场中带电粒子受洛伦茨力作用 时：时：50第50页粒子回转周期与频率分别为：粒子回转周期与频率分别为：当粒子速度与磁场夹角为当粒子速度与磁场夹角为 时：时：带电粒子作带电粒子作螺旋运动螺旋运动：51第51页其轨迹为一圆柱螺旋线，参数方程为：其轨迹为一圆柱螺旋线，参数方程为：其中其中 为角速度，为角速度，h为螺距螺距。52第52页均匀磁场中带电粒子运动轨迹模拟均匀磁场中带电粒子运动轨迹模拟 53第53页设原子核电量为设原子核电量为，位置坐标为，位置坐标为，其质量远大于，其质量远大于 其中：其中：在直角坐标系下：在直角坐标系下：粒子质量粒子质量m。粒子电量为粒子电量为2e，其
25、位置用，其位置用表示，表示，则则 粒子受斥力为：粒子受斥力为：六、六、粒子散射试验粒子散射试验54第54页对应加速度对应加速度a：粒子速度为：粒子速度为：所以所以 粒子坐标为：粒子坐标为：55第55页散射粒子运动轨迹模拟散射粒子运动轨迹模拟 粒子由粒子由位置入射，位置入射，选选取不一取不一样样 对应不一样粒子轨迹对应不一样粒子轨迹56第56页2.5波干涉和衍射图形模拟波干涉和衍射图形模拟一、一、波干涉波干涉图图形模形模拟拟 物理分析物理分析（在均匀介质中传输）：（在均匀介质中传输）：n波源波源：n 点点：这里这里 57第57页点合振动：点合振动：其中其中 (隐函数等值线方程隐函数等值线方程)若
26、若则有则有设有两列相干波，假如初相相等，则干涉条件可简化为用设有两列相干波，假如初相相等，则干涉条件可简化为用波程差波程差 来表示，即当满足：来表示，即当满足：处空间各点，合振幅最大，这时合振幅处空间各点，合振幅最大，这时合振幅 ，强度加强；强度加强；58第58页当满足：当满足：处空间各点，合振幅最小，这时合振幅处空间各点，合振幅最小，这时合振幅当波程差界于上述二者之间，则合振幅也处于上述二者当波程差界于上述二者之间，则合振幅也处于上述二者波干涉图象波干涉图象(两波源相距为两波源相距为4m,=1.5m，)(隐函数等值线方程隐函数等值线方程)大小之间。大小之间。强度减弱；强度减弱；59第59页二
27、、等厚干涉（二、等厚干涉（牛顿环牛顿环）以下列图所表示：以下列图所表示：d表示入射点处膜厚度表示入射点处膜厚度,当垂直入射单当垂直入射单色平行光透过平凸透镜色平行光透过平凸透镜B后后,在空气层上、下表面发生反在空气层上、下表面发生反射形成两束向上相干光射形成两束向上相干光。这两束相干光在平凸透镜下表。这两束相干光在平凸透镜下表面相遇而发生干涉。面相遇而发生干涉。等厚干涉试验及牛顿环干涉图样等厚干涉试验及牛顿环干涉图样 60第60页两束相干光光程差为：两束相干光光程差为：其中：其中：两束相干光相位差为：两束相干光相位差为：牛顿环干涉光强分布为：牛顿环干涉光强分布为：直角坐标系中直角坐标系中 半波
28、损失半波损失61第61页模拟牛顿环干涉图样光强分布模拟牛顿环干涉图样光强分布 62第62页明纹位置：明纹位置：暗纹位置：暗纹位置：讨论：讨论：(1)测透镜球面半径测透镜球面半径R：已知：已知，测，测m、rk+m、rk，可得，可得R(2)测波长测波长：已知：已知R，测出，测出m、rk+m、rk，可得，可得(3)检测透镜曲率半径误差及其表面平整度。检测透镜曲率半径误差及其表面平整度。(4)透射图样与反射图样互补。透射图样与反射图样互补。63第63页三、三、波波衍射衍射图图形模形模拟拟 均匀光源夫琅禾费多缝衍射，其衍射场强度分布为：均匀光源夫琅禾费多缝衍射，其衍射场强度分布为：式中式中 ，其中，其中
29、 是狭缝宽度，是狭缝宽度，是光栅常数，是光栅常数，是是参数选择是任意，但必须服从参数选择是任意，但必须服从物理规律物理规律：缝宽缝宽 必须大于波长必须大于波长 ；光栅常数光栅常数 与与 同数量级；同数量级；衍射角，衍射角，是入射波长，是入射波长，是狭缝有效数目。是狭缝有效数目。64第64页单缝衍射单缝衍射 显然当衍射角显然当衍射角 时，时，在衍射屏上某一点在衍射屏上某一点 处强度为处强度为单缝衍射强度分布曲线单缝衍射强度分布曲线 其中其中 ，为为 点合振动振幅。点合振动振幅。此时此时 ，光强最大，光强最大，为中央明纹中心处强度。为中央明纹中心处强度。65第65页单缝衍射图样模拟单缝衍射图样模拟
30、66第66页讨论讨论:当当 时，时，n由此得由此得 ，解得：，解得：次极大光强比中央明纹中心强度小得多，且随次极大次极大光强比中央明纹中心强度小得多，且随次极大，此时为暗纹条件。在相邻暗纹之间，有一次极大，此时为暗纹条件。在相邻暗纹之间，有一次极大，出现次极大条件为：出现次极大条件为：增加很快降低。增加很快降低。n 由上式可知，由上式可知，67第67页孔径为孔径为R衍射屏衍射屏相对光强曲线相对光强曲线f试验衍射屏放置在试验衍射屏放置在XOY平面上，在屏上平面上，在屏上O点有二点有二分之一径为分之一径为R圆孔，一束平行于圆孔，一束平行于Z轴光线垂直照射圆孔，轴光线垂直照射圆孔，在透镜后方焦平面上
31、可得衍射图样。在透镜后方焦平面上可得衍射图样。四、圆孔夫琅禾费衍射四、圆孔夫琅禾费衍射68第68页圆孔衍射光强公式为：圆孔衍射光强公式为：其中其中 是一阶贝塞尔函数是一阶贝塞尔函数 最终利用矩阵存放空间各点光强大小，再将该矩阵映最终利用矩阵存放空间各点光强大小，再将该矩阵映射成灰度图：射成灰度图：69第69页圆孔夫琅禾费衍射模拟圆孔夫琅禾费衍射模拟 70第70页五、矩形孔夫琅禾费衍射五、矩形孔夫琅禾费衍射 试验衍射屏放置在试验衍射屏放置在XOY平面上，在屏上平面上，在屏上O点有一边长为点有一边长为a、b矩形孔，一束平行于矩形孔，一束平行于Z轴光线垂直照射，在透镜后方轴光线垂直照射，在透镜后方焦
32、平面焦平面L上可得衍射图样。上可得衍射图样。其中：其中：a、b为矩形孔长和宽，为矩形孔长和宽，、为二维衍射角，为二维衍射角，矩形孔衍射光强公式为矩形孔衍射光强公式为:对于直角坐标系有对于直角坐标系有 71第71页矩形孔夫琅禾费衍射光强分布矩形孔夫琅禾费衍射光强分布 72第72页矩形孔夫琅禾费衍射矩形孔夫琅禾费衍射73第73页作业作业EX2-5：编程完成例编程完成例7图形模拟。图形模拟。EX2-6:编程完成波干涉图形模拟。编程完成波干涉图形模拟。EX2-7:有一块有一块光栅，狭缝宽度光栅，狭缝宽度mm，光栅常数光栅常数mm，入射波，入射波，试给出，试给出衍射角衍射角在在范围内光强范围内光强随随分
33、布图形。分布图形。j74第74页EX2-7：图形模拟：图形模拟:讨论讨论：在某一确定衍射角下，光强随入射频率改变。：在某一确定衍射角下，光强随入射频率改变。75第75页本章小结本章小结简谐振动及其合成曲线模拟简谐振动及其合成曲线模拟模拟程序模拟程序（同方向、垂直方向合成）（同方向、垂直方向合成）图形模拟图形模拟阻尼运动和阻尼振动模拟阻尼运动和阻尼振动模拟 差商法差商法（微分方程数值解：高阶降为一阶）（微分方程数值解：高阶降为一阶）程序程序和图形模拟和图形模拟驻波模拟驻波模拟模拟模拟、特征分析、半波损失、特征分析、半波损失隐函数曲线模拟隐函数曲线模拟隐函数曲线编程隐函数曲线编程步骤步骤点电荷系等势线、波干涉图形等编程及曲线模拟点电荷系等势线、波干涉图形等编程及曲线模拟76第76页