FORTRAN数值方法与其在物理学中应用2.pptVIP

第二章 物理图形、图象与计算机模拟 ;§2.1 简谐振动及其合成曲线模拟; Implicit real*8(a-h,o-z) open(1,file=’x-t.dat’) open(2,file=’v-t.dat’) write(*,*)’input A,w,phi,N’ read(*,*)A,w,phi,N pi=3.1415926 do 10 I=1,N t=2.*pi/w t=t*float(I)/N x=A*cos(w*t+phi) v=A*w*cos(w*t+phi+pi/2.) write(1,*)t,x 10 write(2,*)t,v end;图形模拟;二、简谐振动的合成;例2：试给出两个同方向同频率简谐振动的合成程序 ;图形模拟;EX2-1: 编程完成例2。 ; 不同频率情况 ;假设两个分振动振幅都为 ，圆频率 相差较小，取它们的初相位;图形模拟;2. 两个相互垂直方向简谐振动的合成;质点轨迹曲线;若 ，但满足一定整数倍数比关系时，则会 ;§2.2 阻尼运动和阻尼振动的模拟 ;方法2：用 （差商法） ; open(1,file=’vt.dat’) write(*,*)’input a,v0,t=?’ read(*,*) a,v0,t v1=v0 t0=0. v10=v0 write(1,*)t0,v0,v10 dt=t/1000. do 10 j=1,1000 tt=t*float(j)/1000. v=v0*exp((-1.)*a*tt) v1=v1-a*v1*dt 10???write(1,*)tt,v,v1 end;二、阻尼振动;问题分析:;实例说明：;阻尼振动曲线示意图 ;不同阻尼情况振动曲线示意图 ;EX2-2: 编程完成例3。;§2.3 驻波的模拟 ;在 值满足下式的各点，振幅为零－驻波波节处;在 值满足下式的各点，振幅最大－驻波波腹处;驻波模拟程序流图; implicit real*8(a-h,o-z) open(1,file=zhubo.dat) open(2,file=zhubo1.dat) open(3,file=zhubo2.dat) write(*,*)input A,x,wavelength read(*,*)A, x, wa pi=3.1415926 T=2.*pi time=0.5*T dx=x/1000 ;图形模拟;趁脸层涤埠许函朝襟恃降糊穆对海怪氏构久摸藻君碗感棺攻欢务刺蒜咆唬FORTRAN数值方法与其在物理学中应用2FORTRAN数值方法与其在物理学中应用2;每一时刻，驻波都有一定的波形，此波形既不向右移，也不向左移，各点以各自确定的振幅在各自的平衡位置附近振动，没有振动状态或相位的传播，因而称为驻波。; 当波在自由端反射时，则无相位突变，形成驻波时，在 ;EX2-4:一沿;§2.4 点电荷与点电荷系的等势线和电场模拟 ;例6：两个点电荷的电势分布。;二、隐函数曲线的绘制（ ）;编程步骤:;三、 等势线作图步骤;作图步骤:;解：令; open(1,file=’a1.dat’) write(*,*)input q1,q2,x1,x2=? read(*,*)q1,q2,x1,x2 x=x* y=0.0 write(1,*)x