电热法测量焦耳热功当量实验MATLAB源码分享

NUM=1:30;

x=60.*NUM;

c0=4.18.*1000;

m0=0.21611;%水比热与质量 J/kgK kg

c1=389;

m1=0.07975;%铜杯比热与质量 J/kgK kg

c2m2=64.35;%功率电阻和铝搅拌器总

CM=c0.*m0+c1.*m1+c2m2;

K0=263.15;

theta0=26+K0;

A=0.0039083;

B=-0.0000005775;

V=15;

I=0.407;

t=60.*31;

U=[1.1104,1.1119,1.1131,1.1143,1.1155,1.1168,1.1180,1.1192,1.1204,1.1217,1.1228,1.1240,1.1252,1.1263,1.1275,1.1287,1.1297,1.1306,1.1320,1.1332,1.1342,1.1355,1.1365,1.1377,1.1398,1.1409,1.1419,1.1429,1.1441,1.1452,1.1462];

Rt=1000.*U;

R0=1000;%欧姆

theta=(-A+sqrt(A.*A-4.*B.*(1-Rt./R0)))./(2.*B);

theta=theta+K0;

dtheta=theta-theta0;

Q=CM.*dtheta;

for i=NUM;

 theta1(i)=(theta(i)+theta(i+1))./2;

end

%plot(x,theta1) %作散点图

Lxx=sum((x-mean(x)).^2);

Lxy=sum((x-mean(x)).*(theta1-mean(theta1)));

b1=Lxy/Lxx;

b0=mean(theta1)-b1.*mean(x);

theta2=b1.*x+b0;

xlabel('x（时间/秒））','fontsize', 12) %横坐标名

ylabel('y（温度/K）', 'fontsize',12) %纵坐标名

plot(x,theta2);

J1=V.*I.*t./Q(1,31);

J2=V.*I./(b1.*CM);

r=Lxy./sqrt(Lxx)./sqrt(sum((theta1-mean(theta1)).^2));

tao=sqrt(Lxx)./(sqrt(30.*29));

taoba=tao./sqrt(30);

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删