clear all;
clc;
y=[];
y3=[];
data=[
 2.0408
    2.8807
    4.1667
    5.7082
    9.6491
   10.6195
   12.3596
   15.4734
   19.9041
   25.6281
];
x1=data(:,1);
y8=5./(100+x1).*1000;
y8
t=900;
t1=20;
a=0.00013;
x=x1+x1*a*(t-t1);
y=[y
    x];
c=y(:,1);
y1=5./(100+c).*1000;
y1
y6=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9];
y2=y6';
y3=[y3 y2];
y3=[y3 y1];
x2=y3(:,1);
y5=y3(:,2);
A(:,1)=ones(size(x2));A(:,2)=x2;
c1=A\y5;
xx=0:0.02:1;
g=c1(1)*ones(size(xx))+c1(2)*xx;
data1=[
      0   49.8
      0.1 49.0
      0.2 48.6
      0.3 48.0
      0.4 47.3
      0.5 45.6
      0.6 45.2
      0.7 44.5
      0.8 43.3
      0.9 41.7];
x7=data1(:,1);
y7=data1(:,2);
figure(1),subplot(1,2,1),plot(x7,y8,'o',x7,y8,'k');axis([0 0.9 39.0 49.0]);title('電流變化曲線');
legend('離散電流值','電流變化曲線')
xlabel('電阻燒蝕長度（單位：mm）')
ylabel('電流值（單位：mA)')
subplot(1,2,2),plot(x2,y5,'o',x2,y5,'k');axis([0 0.9 39.0 49.0]);
title('修正后電流變化曲線');
legend('離散電流值','電流變化曲線')
xlabel('電阻燒蝕長度（單位：mm）')
ylabel('電流值（單位：mA)')
figure(2),plot(x2,y5,'sk',xx,g,':k');axis([0 0.9 39 49]);
hold on
plot(x2,y5,'o',x2,y5,'k');axis([0 0.9 39 49]);
title('函數組合曲線擬合');
legend('電流變化曲線','擬合后函數曲線')
xlabel('電阻燒蝕長度（單位：mm）')
ylabel('電流值（單位：mA)')
% figure(3),plot(x7,y8,'o',x7,y8,'k');axis([0 0.9 44.0 49.0]);title('電流變化曲線');
% legend('離散電流值','電流變化曲線')
% xlabel('電阻燒蝕長度（單位：mm）')
% ylabel('電流值（單位：mA)')