第三章物理学中定积分数值计算方法.docVIP

第三章 物理学中定积分的数值计算方法 一、填空题 1、库仑常数k等于 9×109mV/C ，真空中的介电常数ε0等于 8.85×10-12F/m 。 2、对于电量为Q的点电荷，在距离r处产生的电场强度为。 3、已知定积分，被积分函数为，积分区间为。将该区间N等分，步长，用曲线下的虚矩形面积和近似替代积分值，该方法称为矩形法。积分近似计算公式为。 4、毕奥—萨伐尔定律所描述的公式为。 5、玻尔兹曼常数是 k=1.38×1023 J/K 。 6、麦克斯韦速率分布律公式。 7、在计算物理中求解定积分的方法有 辛普森法 、 龙贝格法 、 高斯求积法等。 二、简答 1、写出库仑常数、真空中的介电常数和玻尔兹曼常数的值。 答：库仑常数k= 9×109mV/C ，真空中的介电常数ε0= 8.85×10-12F/m，玻尔兹曼常数是 k=1.38×1023 J/K。 2、什么是矩形法？ 答：已知定积分，被积分函数为，积分区间为。将该区间N等分，步长，用曲线下的虚矩形面积和近似替代积分值，该方法称为矩形法。积分近似计算公式为。 3、毕奥—萨伐尔定律和麦克斯韦速率分布律公式。 答：毕奥—萨伐尔定律所描述的公式为。麦克斯韦速率分布律公式。 三、计算、编程 1、将区间（0，π/2）二等分，分别利用矩形法、梯形法、抛物线法计算积分。 解：将区间（0，π/2）二等分。 矩形法计算结果: 梯形法计算结果： 抛物线法计算结果： 若将区间10等分，有G1=1.0764，G2=0.9974，G3=1.000003。 计算程序： open(1,file=’int.dat’) write(*,*)’input a,b,N=?’ read(*,*)a,b,N ! method 1：矩形法 y1=0.0 do 10 j=0,N-1 x1=a x1=x1+float(j)*(b-a)/float(N) 10 y1=y1+f(x1)*(b-a)/float(N) write(1,*)N,y1 write(*,*)N,y1 ! method 2：梯形法 y2=0.0 do 20 j=0,N x2=a x2=x2+float(j)*(b-a)/float(N) if(j.eq.0.or.j.eq.N) then y2=y2+0.5*f(x2)*(b-a)/float(N) else y2=y2+f(x2)*(b-a)/float(N) end if 20 continue write(1,*)N,y2 write(*,*)N,y2 ！ method 3：抛物线法（N: 偶数） y3=0.0 do 30 j=0,N x3=a x3=x3+float(j)*(b-a)/float(N) if(j.eq.0.or.j.eq.N) then y3=y3+1./3.*f(x3)*(b-a)/float(N) else k=j-2*int(j/2) if(k.eq.0) then y3=y3+2./3.*f(x3)*(b-a)/float(N) else y3=y3+4./3.*f(x3)*(b-a)/float(N) end if end if 30 continue write(1,*)N,y3 write(*,*)N,y3 end function f(x) f=cos(x) end 2、设有一长直导线均匀带电，线电荷密度为λ，长度为2l。求空间任意一点P的电势, 试利用变步长辛普森求积分法求解。(取,a=l,b=-l,x0=0.125,y0=0,λ=5×10-9,l=0.05，ε0=8.85×10-12F/m)解：在导线上取一小段dx，视为点电荷，其电量为λdx，它在P点产生的电势为： 因此有： 利用解析法可得以上定积分结果为： 采用数值法进行计算，则有 计算程序： subroutine simp(a,b,ep,s2,f)