贵州大学 大学物理实验3 气体比热容比的测定

 

#include<stdio.h>
#include<math.h>
float P(float pl, float m, float r) {//计算气压 
    float p = 0.00000;
    p = pl + m * 10 / 3.14 / pow(r, 2);
    return p;
}

float gama(float m, float v, float t, float p, float d) {//计算γ，其中t，d为平均值 
    float gama = 0.00000;
    gama = 64 * m * v / pow(t, 2) / p / pow(d, 4);
    return gama;
}

float um() {//计算um 
    float um = 0.00000;
    um = 5 * pow(10, -5) / sqrt(3);
    return um;
}

float ut(float t1, float t2, float t3, float t4, float t5) {
    float st = 0.00000;
    float avt = (t1 + t2 + t3 + t4 + t5) / 5;
    st = sqrt((pow(t1 - avt, 2) + pow(t2 - avt, 2) + pow(t3 - avt, 2) + pow(t4 - avt, 2) + pow(t5 - avt, 2)) / sqrt(5) / 2);
    float uat = 0.00000;
    uat = 1.14 * st;
    float ubt = 0.01 / sqrt(3);
    float ut = 0.00000;
    ut = sqrt(pow(uat, 2) + pow(ubt, 2));
    return ut;
}

float ud(float d1, float d2, float d3, float d4, float d5) {
    float avd = 0.00000;
    avd = (d1 + d2 + d3 + d4 + d5) / 5;
    float sd = 0.00000;
    sd = sqrt((pow(d1 - avd, 2) + pow(d2 - avd, 2) + pow(d3 - avd, 2) + pow(d4 - avd, 2) + pow(d5 - avd, 2)) / sqrt(5) / 2);
    float ua = 0.00000;
    ua = 1.14 * sd;
    float ub = 0.00000;
    ub = 4 * pow(10, -6) / sqrt(3);
    float ud = 0.00000;
    ud = sqrt(pow(ua, 2) + pow(ub, 2));
    return ud;
}

float ur(float gama, float ud, float um, float ut, float d1, float d2, float d3, float d4, float d5, float t1, float t2, float t3, float t4, float t5, float m) {
    float avd = 0.00000;
    avd = (d1 + d2 + d3 + d4 + d5) / 5;
    float avt = 0.00000;
    avt = (t1 + t2 + t3 + t4 + t5) / 5;
    float er = 0.00000;
    er = sqrt(pow(um / m, 2) + pow(ut * 2 / avt, 2) + pow(4 * ud / avd, 2));
    float ur = 0.00000;
    ur = er * gama;
    return ur;
}

int main() {
    printf("请输入d1 d2 d3 d4 d5 t1 t2 t3 t4 t5 m v p且转化为国际单位制,每输入一个值按一下回车");
    float d1, d2, d3, d4, d5, t1, t2, t3, t4, t5, m, v, p;
    scanf("%f %f %f %f %f %f %f %f %f %f %f %f %f", &d1, &d2, &d3, &d4, &d5, &t1, &t2, &t3, &t4, &t5, &m, &v, &p);
    float tp, tgama, tum, tut, tud, tur;
    tp = P(p, m, (d1 + d2 + d3 + d4 + d5) / 5);
    tgama = gama(m, v, (t1 + t2 + t3 + t4 + t5) / 5, tp, (d1 + d2 + d3 + d4 + d5) / 5);
    tum = um();
    tut = ut(t1, t2, t3, t4, t5);
    tud = ud(d1, d2, d3, d4, d5);
    tur = ur(tgama, tud, tum, tut, d1, d2, d3, d4, d5, t1, t2, t3, t4, t5, m);
    printf("ur=%f gama=%f", tur, tgama);
}