fortran数值方法及其在物理学中应用3.ppt

例6：已知圆环带电量为q ，杆的线密度为? ，长为L。 求：杆对圆环的作用力。 例8：求载流圆线圈轴线上一点 P 的磁感应强度 例9：求绕轴旋转的带电圆盘轴线上的磁场 载流圆线圈轴线上磁感应强度分布 x q r O R P 解： 带电圆盘电荷面密度为： 所取环带上带电量为： 等效电流为： q r O R P 圆环dq在P点产生的磁感应强度dB： 对上式积分： 取σ=1，ω=600，R=1，x=1 其中 由上式可得B=4.536746E-05 a=0，b=1，N=100 矩形法 B=4.5071927E-05 梯形法 B=4.5738361E-05 辛普森法 B=4.5736691E-05 绕轴旋转的带电圆盘轴线上的磁场分布 三、分子物理中数值积分的应用 理想气体在平衡态下分子的速率分布函数: 是分子质量， 是气体温度， 为波耳兹曼常数。 麦克斯韦速率分布律 气体中速率在 到 间的分子数的比率为: 例10：利用气体分子麦克斯韦速率分布律，求0度，127度 下氮分子、氧分子运动的速率分布曲线，并求当速 率在300～500m/s范围内的分子所占的比例，讨论温 度 及分子量 对速率分布曲线的影响。 解：根据麦克斯韦速率分布律可知，任一速率间隔 到 中 对于氮气，分子质量分别为 kg/ ，而氧气分子 的分子数所占的比率可用积分法求出，即 质量 kg/ ，其中 =6.022为阿伏伽德罗常数。 气体分子运动麦克斯韦速率分布曲线 ? 当温度升高和分子量减小时，速率小的分子数减小而速率大的分子 ? 随着温度升高和分子量减小，分布曲线将变得越来越平坦。 氮气分子所占比例为0.3957 (T=273K)，0.3085(T=400K)； 氧气分子所占比例为0.4189 (T=273K)，0.3410(T =400K)。 数增多，分布曲线的极大值随着温度升高和分子量减小而向右移动。 利用辛普森法可得，当速率在300～500m/s范围内： EX3-5: 结合理想气体平衡态下的麦克斯韦速率分布律，利用辛普森法计算0°C下氮气、氧气分子运动的平均速率 和 方均根速率 ，并与解析解 ， 做比较，其中 是气体的摩尔质量， 为摩尔气体常数。 EX3-3 p.67 EX3-4 p.67 作业 dimension y(0:47) coe=2*(4.561*1.0e-5)*9.81/0.0956 data y/0.10135,0.20064,0.27303,0.31095,0.33094,0.33991,0.34474, * 0.35577,0.31508,0.29578,0.27717,0.26131,0.24545,0.23097,0.21718, * 0.20339,0.19167,0.17995,0.16823,0.15789,0.14824,0.13927,0.13288, * 0.12548,0.11859,0.11238,0.10687,0.10204,0.09215,0.09308,0.08894, * 0.08480,0.08067,0.07722,0.07377,0.07032,0.06757,0.06481,0.06205, * 0.05929,0.05654,0.05378,0.05102,0.04826,0.04550,0.04274,0.04067, * 0.03861/ C method 3：辛普森法 y3=0.0 do 30 j=0,47 If(j.eq.0.or.j.eq.47) then y3=y3+1./3.*y(j)*0.0127 else k=j-2*int(j/2) if(k.eq.0) then y3=y3+2./3.*y(j)*0.0127 else y3=y3+4./3.*y(j)*0.0127 end if end if 30 continue write(*,*) sqrt(coe*y3) end 小结： 1．矩形法： 2．梯形法： 其中： 3．抛物线法(辛普森法)： 其中： 问题： 1.二重积分问