clear all
a=0:10:100;
b=[52.6 52.2 51.9 51.5 51.1 50.8 50.5 50.1 49.8 49.5 49.1];
plot(a,b)
hold on
st1=75.1
st2=50.8 %单位是摄氏度
c=0.389*10^3 %单位是J/(Kg*℃)
hp=[8.50 8.44 8.52 8.54 8.50];hp=hp*10^-3;
dp=[130.10 130.16 130.18 130.24 130.24];dp=dp*10^-3;
hb=[9.56 9.62 9.66 9.62 9.60];hb=hb*10^-3;
db=[129.92 129.92 129.92 129.88 129.92];db=db*10^-3;
mp=[1117.73 1117.73 1117.73 1117.76 1117.75];mp=mp*10^-3;
hpb=sum(hp)/5
dpb=sum(dp)/5
hbb=sum(hb)/5
dbb=sum(db)/5 %单位是米
mpb=sum(mp)/5 %单位是千克
x1=linspace(0,100);
a1=40:10:60;
b1=[51.1 50.8 50.5];%在50.8摄氏度找切线斜率
xl=polyfit(a1,b1,1);
y1=xl(1)*x1+xl(2);
plot(x1,y1)
xl(1)
k=mpb*c*xl(1)*(dpb+4*hpb)/(dpb+2*hpb)*hbb/(st1-st2)*2/pi/dbb^2
ubmp=0.01*10^-3/3^0.5
ub=0.02*10^-3/3^0.5
uahp=0;
uadp=0;
uahb=0;
uadb=0;
for i=1:5
uahp=uahp+(hp(i)-hpb)^2;
end
uahp=(uahp/20)^0.5
for i=1:5
uadp=uadp+(hp(i)-dpb)^2;
end
uadp=(uadp/20)^0.5
for i=1:5
uahb=uahb+(hp(i)-hbb)^2;
end
uahb=(uahb/20)^0.5
for i=1:5
uadb=uadb+(db(i)-dbb)^2;
end
uadb=(uadb/20)^0.5
uhp=(uahp^2+ub^2)^0.5
udp=(uadp^2+ub^2)^0.5
uhb=(uahb^2+ub^2)^0.5
udb=(uadb^2+ub^2)^0.5
ukk=((ubmp/mpb)^2+((1/(dpb+4*hpb)-1/(dpb+2*hpb))*udp)^2+((4/(dpb+4*hpb)-2*(dpb+2*hpb))*uhp)^2+(uhb/hbb)^2+(2*udb/dbb)^2)^0.5
uk=k*ukk