有哪位哥可以讲一下这个函数嘛(matlab笨笨瑟瑟发抖)

load abc.mat

x = abc.t;

y = abc.T;

for n=1:7500
if(abs(y(n))>1000.0)
y(n)=20.0;
end
end

y1=smoothdata(y,'movmedian',20);
plot(x,y1);

temp1=ones(1,40);
time1=(1:1:40)0.1;
for n=1:40
temp1(n)=y1(n10+2300);
end
plot(time1,temp1)

y2=ones(1,20);
for n=1:20
y2(n)= (temp1(n)+28)/(32+28);
end
t=(1:1:20)*0.1;
plot(t,y2);

y3=ones(1,17);
t1=(1:1:17)*0.1;
for n=1:17
y3(n)= log(1-y2(n));
end
plot(t1,y3,'.')%标点
p=polyfit(t1,y3,1);
X=0.1:1.9;
Y=p(1)*X+p(2);%拟合
hold on
plot(X,Y)
hold off

tau1=1/p(1);

%1.从热水中拿出
temp2=ones(1,40);
time2=(1:1:40)0.1;
for n=1:40
temp2(n)=y1(n10+5400);
end
plot(time2,temp2)

r1=ones(1,20);
for n=1:20
r1(n)= (temp2(n)-13)/(28-13);
end
t=(1:1:20)*0.1;
plot(t,r1);

r2=ones(1,20);
t2=(1:1:20)*0.1;
for n=1:20
r2(n)= log(1-r1(n));
end
plot(t2,r2,'.')%标点
r=polyfit(t2,r2,1);
X=0.1:1.9;
Y=r(1)*X+r(2);%拟合
hold on
plot(X,Y)
hold off
tau2=1/r(1);

%2放入冷水中
temp3=ones(1,30);
time3=(1:1:30)0.1;
for n=1:30
temp3(n)=y1(n10+5850);
end
plot(time3,temp3)

s1=ones(1,20);
for n=1:20
s1(n)= (temp3(n)-26)/(46-26);
end
t=(1:1:20)*0.1;
plot(t,s1);

s2=ones(1,20);
t2=(1:1:20)*0.1;
for n=1:20
s2(n)= log(1-s1(n));
end
plot(t2,s2,'.')%标点
s=polyfit(t2,s2,1);
X=0.1:1.9;
Y=s(1)*X+s(2);%拟合
hold on
plot(X,Y)
hold off
tau3=1/s(1);

%在冷水中
temp4=ones(1,50);
time4=(1:1:50)0.1;
for n=1:50
temp4(n)=y1(n10+6300);
end
plot(time4,temp4)

k1=ones(1,40);
for n=1:40
k1(n)= (temp4(n)-34)/(43-34);
end
t=(1:1:40)*0.1;
plot(t,k1);

k2=ones(1,30);
t2=(1:1:30)*0.1;
for n=1:30
k2(n)= log(1-k1(n));
end
plot(t2,k2,'.')%标点
k=polyfit(t2,k2,1);
X=0.1:2.9;
Y=k(1)*X+k(2);%拟合
hold on
plot(X,Y)
hold off
tau4=1/k(1);