FORTRAN數(shù)值方法及其在物理學(xué)中應(yīng)用.pptVIP

1,第二章 物理圖形、圖象與計算機(jī)模擬,2,2.1 簡諧振動及其合成曲線模擬,一、簡諧振動的 曲線和 曲線,例1：畫出 曲線及對應(yīng)的 曲線，其中,3,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,計算程序,4,圖形模擬,Origin簡介,5,二、簡諧振動的合成,1.同方向簡諧振動的合成,同頻率情況,其中,6,例2：試給出兩個同方向同頻率簡諧振動的合成程序,open(1,file=x1.dat) open(2,file=x2.dat) open(3,file=x.dat) write(*,*)input A1,A2,w,phi1,phi2=? read(*,*)A1,A2,w,phi1,phi2 pi=3.1415926 do 10 I=1,1000 t=2.*pi/w t=t*float(I)/1000 x1=A1*cos(w*t+phi1) x2=A2*cos(w*t+phi2) x=x1+x2 write(1,*)t,x1 write(2,*)t,x2 10 write(3,*)t,x end,7,圖形模擬,兩個同方向同頻率簡諧振動的合成,8,EX2-1: 編程完成例2。,兩個同方向同頻率簡諧振動的合成,9,不同頻率情況,若 ，會出現(xiàn)拍的現(xiàn)象。,合振動不再是簡諧振動，利用旋轉(zhuǎn)矢量法可以求得合振動的振幅為,中間經(jīng)歷的時間 稱為周期，顯然 ， 頻率：,振幅在 和 間周期性地變化，屬振動調(diào)制。 合振動振幅從一次極大到相鄰的另一次極大。,10,假設(shè)兩個分振動振幅都為 ，圓頻率 相差較小，取它們的初相位,此時合成運動的位移可寫成：,變化主要取決于 ，振幅按變化 。,都是零，則可以分別表示為：,由于圓頻率 遠(yuǎn)大于圓頻率 ，,11,圖形模擬,兩個同方向頻率近似的簡諧振動的合成,12,2. 兩個相互垂直方向簡諧振動的合成,若 ，則有合振動方程：, （一、三象限直線方程）,（二、四象限直線方程）, （橢圓方程）,13,質(zhì)點軌跡曲線,下圖所示為兩個頻率相同、振幅相等、相互垂直而相位差,兩個頻率相同、振幅相等、相互垂直簡諧振動的合成,為 下的質(zhì)點軌跡曲線。,14,若 ，但滿足一定整數(shù)倍數(shù)比關(guān)系時，則會,下圖所示為兩個頻率不同（滿足 ）、振幅相等、,利薩如圖形示意圖,出現(xiàn)利薩如圖形：,相互垂直而相位差為0， 下的質(zhì)點軌跡曲線。,15,2.2 阻尼運動和阻尼振動的模擬,一、阻尼情況下物體運動的 曲線,例3：質(zhì)量為 的摩托快艇以速度 行駛，它受到的摩擦阻力 與速度成正比，設(shè)比例系數(shù)為 ，則 ，試求關(guān) 閉發(fā)動機(jī)后，對 的變化規(guī)律。 （取,解：物理分析與數(shù)學(xué)模型,方法1:,16,方法2：用 （差商法）,17,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,模擬程序,18,二、阻尼振動,問題：彈簧振子阻尼振動的方程為 阻尼因子 彈簧振子的角頻率,試用函數(shù)近似法作出位移與時間的函數(shù)變化曲線。,19,問題分析:,解：將二階微分方程化為一階微分方程,即,而,20,實例說明：,例4：畫出當(dāng) ， ， 時， ， ， ， 秒下的 曲線。,計算程序： open(1,file=v-t.dat) open(2,file=x-t.dat) write(*,*)input B,w0,x0,v0,t=? read(*,*)B,w0,v0,x0,t dt=t/1000. v=v0 x=x0 tt0=0.0 write(1,*)tt0,v0 write(2,*)tt0,x0,do 10 j=1,1000 tt=float(j)*dt f=-2.*B*v-w0*2*x v=v+f*dt x=x+v*dt write(1,*)tt,v 10 write(2,*)tt,x end,21,阻尼振動曲線示意圖,圖形模擬,22,不同阻尼情況振動曲線示意圖,圖形模擬,23,EX2-2: 編程完成例3。,EX2-3:一石子從空中靜止下落，已知,曲線。其中,。,作業(yè),24,2.3 駐波的模擬,定義：兩列振幅、振動方向和頻率都相同而傳播方向相 反的兩列同類波相干疊加形成駐波。,設(shè)有兩列振動方向相同、振幅相同、頻率相同的平面余弦波，,按疊加原理，合成的駐波的波函數(shù)為:,25,在 值滿足下式的各點，振幅為零駐波波節(jié)處,相鄰兩波節(jié)的距離為半波長，即：,討論:,各質(zhì)點都在作同頻率的諧振動。,的振幅按余弦函數(shù)規(guī)律分布。,26,在 值滿足下式的各點，振幅最大駐波波腹處,相鄰兩波腹間的距離也為半波長，即：,l 波節(jié)處的質(zhì)點振動的振幅為零，始終處于靜止； l 波腹處的質(zhì)點振動的振幅最大，等于 。 l 其他各處質(zhì)點振動的振幅則在零與最大之間。 l 兩相鄰波節(jié)或兩相鄰波腹之間相距半波長。 l 波腹和相鄰波節(jié)間的距離為 ，波腹和波節(jié)交替作等距離排列。,討論:,27,駐波模擬程序流圖,28,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,do 10 I=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,y write(2,*)x,y1 10 write(3,*)x,y2 end,模擬程序,29,圖形模擬,，,30,31,每一時刻，駐波都有一定的波形，此波形既不向右移，也不向左移，各點以各自確定的振幅在各自的平衡位置附近振動，沒有振動狀態(tài)或相位的傳播，因而稱為駐波。,又沿相反方向同時通過平衡位置。,有相同的符號，而相鄰的兩端符號總是相反的，這說明在駐,波中同一段上各質(zhì)點的振動相位相同，而相鄰兩段中的各點,振動相位相反。,同一段內(nèi)各點沿相同方向同時達(dá)到各自振動位移的最大值，,又沿相同方向同時通過平衡位置；,波節(jié)兩側(cè)各點同時沿相反方向達(dá)到振動位移的正、負(fù)最大值，,32,當(dāng)波在自由端反射時，則無相位突變，形成駐波時，在,半波損失,在弦線上進(jìn)行的駐波實驗，反射點處弦線是固定不動的，,這一點只能是波節(jié)。這說明反射波和入射波的相位在反射,點正好相反，即入射波在反射點反射時相位有 的突變。,根據(jù)相位差 與波程差的關(guān)系（ ），相位差為,說，反射波和入射波之間存在著半個波長的波程差，這種,相位突變 稱為半波損失。,就相當(dāng)于半個波長的波程差，這說明對固定端的反射點來,自由端出現(xiàn)波腹。,33,作業(yè),34,2.4 點電荷與點電荷系的等勢線和電場模擬,一、等勢線方程,要求：作出滿足等勢線方程 的等勢線。,例5：一個點電荷 在 處產(chǎn)生的電勢,等勢線方程為：,即,圓心在 ，半徑為 的圓。,由上式可得：,畫圖時可采用參數(shù)方程：,35,例6：兩個點電荷的電勢分布。,等勢線方程為：,36,二、隱函數(shù)曲線的繪制（ ）,（這里引入?yún)?shù) ）,指定 (如 )，計算出 和 ，由上式可得 和 。,若點 處電勢為 ，則 處的電勢也為 。,37,編程步驟:,把上述點 連起來，就是 的等勢線。,Step2:,38,三、 等勢線作圖步驟,39,作圖步驟:,3.,而,同樣有：,因此有：,取：,40,計算才會保證曲線的完整性。,41,open(1,file=a1.dat) write(*,*)input q1,q2,x1,x2=? read(*,*)q1,q2,x1,x2 x=x* y=0.0 write(1,*)x,y dt=0.005,do 10 j=1,100 r1=(x-x1)*2+y*2 r1=sqrt(r1) r2=(x-x2)*2+y*2 r2=sqrt(r2) dx=y*(q1*r2*3+q2*r1*3)*dt x=x+dx dy=(-1.)*(q1*(x-x1) *r2*3+ q2*(x-x2)*r1*3)*dt y=y+dy 10 write(1,*)x,y end,模擬程序,42,圖形模擬, （異號）,兩個點電荷電勢等勢線示意圖,（同號）,43,同號情況,44,q1=-1, q2=1, x1=-2, x2=2,q1=1, q2=-1, x1=-2, x2=2,異號情況,45,四、點電荷系電場線的圖像模擬,單個點電荷的電場：,點電荷系的電場：,例8：如圖所示，A、B、C三個點電荷組成一個點電荷系，其中A(0,0)點帶電為 ， B(0,1)點帶電亦為 ，最后C(0,-1)點帶電 ，試模擬 平面內(nèi)電場線分布。,46,點電荷系在某點產(chǎn)生的電場強(qiáng)度等于各點電荷單獨在該點產(chǎn)生的電場強(qiáng)度的矢量和，這稱為電場強(qiáng)度疊加原理。,所以，在XOY平面內(nèi)電場強(qiáng)度為：,電場線上每一點的切線方向反映了該點的場強(qiáng)方向，有：,解上述方程可得：,47,此超越方程（等值線）即為這個點電荷系的電場線方程，,點電荷系的電場線模擬,常數(shù)C選取不同的值對應(yīng)不同的電場線。(C=0.13,間隔0.1）,48,例9：電偶極振子電場的模擬,設(shè)電偶極子的電偶極矩為：,在球坐標(biāo)系下，任意時刻t空間任意處r的輻射電場為：,輻射的電場線滿足方程:,49,解上述方程即可得電場線滿足的方程：,當(dāng)C 取不同值時可得到不同的輻射電場線。,電偶極子輻射電場線模擬,50,五、帶電粒子在電磁場中的運動,在均勻磁場中帶電粒子受洛倫茨力的作用,時：,51,粒子回轉(zhuǎn)周期與頻率分別為：,當(dāng)粒子速度與磁場夾角為 時：,帶電粒子作螺旋運動：,52,其軌跡為一圓柱螺旋線，參數(shù)方程為：,53,均勻磁場中帶電粒子運動軌跡模擬,54,其中：,在直角坐標(biāo)系下：,六、 粒子散射實驗,55,對應(yīng)的加速度a：,粒子速度為：,所以 粒子坐標(biāo)為：,56,散射粒子運動軌跡模擬,57,2.5 波的干涉和衍射圖形模擬,一、波的干涉圖形模擬,物理分析（在均勻介質(zhì)中傳播）：,波源：,點：,這里,58,點的合振動：,其中,(隱函數(shù)等值線方程),設(shè)有兩列相干波，如果初相相等，則干涉條件可簡化為用,波程差 來表示，即當(dāng)滿足：,處的空間各點，合振幅最大，這時合振幅 ，,強(qiáng)度加強(qiáng)；,59,處的空間各點，合振幅最小，這時合振幅,當(dāng)波程差界于上述兩者之間，則合振幅也處于上述兩者,波的干涉圖象 (兩波源相距為4m, =1.5m， ),(隱函數(shù)等值線方程),大小之間。,強(qiáng)度減弱；,60,二、等厚干涉（牛頓環(huán)）,如下圖所示：d表示入射點處膜的厚度, 當(dāng)垂直入射的單色平行光透過平凸透鏡B后, 在空氣層的上、下表面發(fā)生反射形成兩束向上的相干光。這兩束相干光在平凸透鏡下表面相遇而發(fā)生干涉。,等厚干涉實驗及牛頓環(huán)干涉圖樣,61,兩束相干光的光程差為：,其中：,兩束相干光的相位差為：,牛頓環(huán)干涉光強(qiáng)分布為：,直角坐標(biāo)系中,62,模擬牛頓環(huán)干涉圖樣光強(qiáng)分布,63,明紋位置：,暗紋位置：,討論：,(1) 測透鏡球面的半徑R：已知，測m、rk+m、rk，可得R,(2) 測波長：已知R，測出m、rk+m、rk，可得,(3) 檢測透鏡的曲率半徑誤差及其表面平整度。,(4) 透射圖樣與反射圖樣互補(bǔ)。,64,三、波的衍射圖形模擬,均勻光源的夫瑯禾費多縫衍射，其衍射場強(qiáng)度分布為：,式中 ， ，其中 是狹縫寬度， 是光柵常數(shù)， 是,衍射角， 是入射波長， 是狹縫的有效數(shù)目。,65,單縫衍射,顯然當(dāng)衍射角 時， ，,在衍射屏上某一點 處的強(qiáng)度為,單縫衍射強(qiáng)度分布曲線,其中 ， 為 點合振動的振幅。,此時 ，光強(qiáng)最大， 為中央明紋中心處強(qiáng)度。,66,單縫衍射圖樣模擬,67,討論:,當(dāng) 時， ，,由此得 ，解得：,，此時為暗紋條件。在相鄰暗紋之間，有一次極大，,出現(xiàn)次極大條件為：,的增加很快減少。,由上式可知，,68,實驗衍射屏放置在XOY平面上，在屏上O點有一半徑為R的圓孔，一束平行于Z軸的光線垂直照射圓孔，在透鏡后方的焦平面上可得衍射圖樣。,四、圓孔的夫瑯禾費衍射,69,圓孔衍射光強(qiáng)公式為：,其中 是一階