密立根油滴实验-平衡法-数据处理-C程序源代码.doc

《密立根油滴实验-平衡法-数据处理-C程序源代码.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、/*此程序用于CCD显微镜密立根油滴实验实验**********只针对静态平衡法******************参数设置只需修改预处理命令中的参数值******************最后的相对误差计算未取绝对值***************该程序已经进行过相当的测验，尚未发现问题。即便如此，也请先进行测试，再使用。2011-10-18*/#include"stdio.h"#include"math.h"#defineb0.00823//修正常数#definep101325//大气压强#defineq979//油滴密度#defineg9.794//重力加速度#definePi3.1415

2、926//圆周率#defined0.005//极板距离#definej0.0000183//空气粘滞系数#definel0.0016//下落距离#definee1.602*0.0000000001*0.000000001//单位电荷量10E-19longdoubleex[6],sume=0;//用于存储各个油滴的计算电荷量和所有油滴电荷量总和的存储doubleQ=0,a=0,U=0,t=0,tp=0,sumt=0,tx[6];//Q——油滴计算带电量、a——中间量……longdoublee0=0;inti,N=0,k;intflag=0;charz;voidmain(){e0=e;print

3、f("是否需要计算平均下落时间(yorn)?");z=getchar();if(z=='y'){for(i=1;i<=5;i++){sumt=0;printf("****第%d滴油滴****",i);for(k=1;k<=5;k++){printf("请输入第%d次下落时间t=",k);scanf("%lf",&tp);sumt=sumt+tp;}tx[i]=sumt/5;printf("平均下落时间：t=%lf",tx[i]);}flag=1;}if((z=='n')

4、

5、flag==1){for(i=1;i<=5;i++){printf("****第%d滴油滴****"

6、,i);printf("请输入电压平均值U：U=");scanf("%lf",&U);if(flag==0){printf("请输入下落时间平均值t：t=");scanf("%lf",&t);}if(flag==1)t=tx[1];a=sqrt(9*j*l/2/q/g/t);printf("中间量a=%E",a);Q=18*Pi/sqrt(2*q*g)*d/U;Q=Q*pow((j*l/t/(1+b/p/a)),1.5);N=(int)((longdouble)Q/(longdouble)e0);if(((longdouble)Q/(longdouble)e0)-N>=0.5)N

7、++;//四舍五入ex[i]=(longdouble)Q/N;sume=sume+ex[i];printf("电荷量q=%E电荷个数N=%d计算单位电荷量ei=%e",Q,N,ex[i]);}printf("平均电子电荷e=%e相对误差E=%e",sume/5,(sume/5-e0)/e0);}}