【代码】【matlab】准稳态法测比热和导热系数

%准稳态法测比热和导热系数

%% glass data

vt(1,:)=[0.075 0.091 0.100 0.103 0.106 0.107 0.109 0.109 0.110 0.110 0.110 0.111 0.112 0.111 0.112 0.111 0.112 0.111];

v(1,:)=[ 0.025 0.036 0.054 0.076 0.096 0.119 0.139 0.161 0.182 0.205 0.226 0.248 0.268 0.290 0.310 0.331 0.351 0.371];

for i=1:17

vc(1,i)=v(1,i+1)-v(1,i);

end

%% rubber data

vt(2,:)=[0.097 0.135 0.158 0.170 0.177 0.182 0.184 0.186 0.188 0.188 0.190 0.192 0.194 0.195 0.195 0.195 0.195 0.195];

v(2,:)=[ 0.819 0.137 0.209 0.226 0.247 0.268 0.290 0.312 0.335 0.356 0.381 0.405 0.429 0.453 0.476 0.498 0.521 0.543];

for i=1:17

vc(2,i)=v(2,i+1)-v(2,i);

end

%% other data

ro(1)=1374; %kg*m^-3

ro(2)=1196; %kg*m^-3

R=0.01; %m

a=0.85;

f=a*0.09*0.09;

r=110; %omen

%% formula

for k=1:2

t(k)=vt(k,17)/0.04; %k

q=18^2/(2*r*f);

veo_rise(k)=vc(k,17)/(60*0.04);

lamda(k)=q*R/(2*t(k));

c(k)=q/(ro(k)*R*veo_rise(k));

end