clear
close all
yyy=load("Thermocouple_alldata.txt");
indsgood=find(yyy(:,1)<100);
yyy=yyy(indsgood,:);
all_tempdata=yyy(:,1);
all_thermnumber=yyy(:,2);
all_tablenumber=yyy(:,3);
all_tiptemp=yyy(:,5);
all_roomtemp=yyy(:,6);
all_teamnumber=yyy(:,7);
teamlist=unique(all_teamnumber);
all_basetemp=65;
theta=all_tempdata-all_roomtemp;
theta_b=all_basetemp-all_roomtemp;
ThetaRatio=theta./theta_b;
for tnow=1:9
indsnow1=find(tnow==all_thermnumber);
XX=ThetaRatio(indsnow1);
[MUHAT,SIGMAHAT,MUCI,~] = normfit(XX,0.05);
mean_allthermnumber(tnow)=MUHAT;
MUCI1_allthermnumber(tnow)=MUCI(1);
MUCI2_allthermnumber(tnow)=MUCI(2);
all_thermnumber_list(tnow)=tnow;
end
h=39;
L=.45;
k=167;
rad=0.01;
P=2*pi*rad;
Ac=pi*rad^2;
m=sqrt(h*P/(k*Ac));
x=[0 0.0055 .01 .015 .02 .025 .03 .035 .04]*10;
figure(1)
plot(all_thermnumber(:),all_tempdata(:),'r.')
figure(2)
hold on
plot(all_thermnumber(:),ThetaRatio(:),'k.')
plot(all_thermnumber_list,mean_allthermnumber,'k.', 'LineWidth', 2)
plot(all_thermnumber_list,MUCI1_allthermnumber,'k.', 'LineWidth',2)
plot(all_thermnumber_list,MUCI2_allthermnumber,'k.', 'LineWidth', 2)
ylim([0 1.5])
figure(3)
ThetaL=all_tiptemp - all_roomtemp;
for mm=1:4
subplot(1,4,mm)
hold on
if mm==1; Theta_ThetaB=(cosh(m*(L-x))+(h/(m*k))*sinh(m*(L-x)))/(cosh(m*L)+(h/
(m*k))*sinh(m*L)); title('Convective'); end
if mm==2; Theta_ThetaB=(cosh(m*(L-x)))/cosh(m*L); title('Adiabatic'); end
