扭摆法测定物体转动惯量 计算程序

本代码适用于大学物理实验“扭摆法测定物体转动惯量”

因实验环境不同，使用前需调整参数

#include<stdio.h> 
int main()
{
	double m1=0.7297; //塑料圆柱体(高)-质量（kg） 
	double m2=0.7142; //金属圆筒-质量（kg）
	double m3=0.5367;//塑料圆柱体(矮)-质量（kg） 

	double D1=10.00867; //塑料圆柱体(高)-平均直径（cm） 
	double D21=9.99600; //金属圆筒-平均外径（cm）
	double D22=9.33333; //金属圆筒-平均内径（cm）
	double D3=10.01733;//塑料圆柱体(矮)-平均直径（cm）
	
	double T0=0.77500; //载物盘-平均周期（s）
	double T1=1.26003; //塑料圆柱体(高)-平均周期（s） 
	double T2=1.5723; //金属圆筒-周期平均（s）
	double T3=1.1517; //塑料圆柱体(矮)-平均周期（s）
	
	double Pai=3.14;
	
	double I10=(m1*D1*D1)/8;//塑料圆柱体(高)-理论值 
	double I20=(m2*((D21*D21)+(D22*D22)))/8;//金属圆筒-理论值
	double I30=(m3*D3*D3)/8;//塑料圆柱体(矮)-理论值
	double K=(4*Pai*Pai*I10)/((T1*T1)-(T0*T0));//利用 塑料圆柱体(高) 所求公式系数 
	double I0=(I10*T0*T0)/((T1*T1)-(T0*T0));//载物盘-理论值 
	double I1=((K*T1*T1)/(4*Pai*Pai))-I0;//塑料圆柱体(高)-实际值 
	double I2=((K*T2*T2)/(4*Pai*Pai))-I0;//金属圆筒-实际值
	double I3=((K*T3*T3)/(4*Pai*Pai))-I0;//塑料圆柱体(矮)-实际值
	
	printf("I10=%lf\n\nI20=%lf\n\nI30=%lf\n\nK=%lf\n\nI0=%lf\n\nI1=%lf\n\nI2=%lf\n\nI3=%lf",I10,I20,I30,K,I0,I1,I2,I3);
 }