📄 tar.m
字号:
load sunspot.dat;
x=(sunspot(:,2))';
maxx=max(x);
minx=min(x);
N=length(x);
hist(x);
%------------------------
x1=zeros(1,N); x2=zeros(1,N); x3=zeros(1,N); x4=zeros(1,N);
x5=zeros(1,N); x6=zeros(1,N); x7=zeros(1,N); x8=zeros(1,N);
x9=zeros(1,N); x10=zeros(1,N); x11=zeros(1,N); x12=zeros(1,N);
x13=zeros(1,N); x14=zeros(1,N); x15=zeros(1,N); x16=zeros(1,N);
x17=zeros(1,N); x18=zeros(1,N); x19=zeros(1,N); x20=zeros(1,N);
x1d=zeros(1,N); x2d=zeros(1,N); x3d=zeros(1,N); x4d=zeros(1,N);
x5d=zeros(1,N); x6d=zeros(1,N); x7d=zeros(1,N); x8d=zeros(1,N);
x9d=zeros(1,N); x10d=zeros(1,N); x11d=zeros(1,N); x12d=zeros(1,N);
x13d=zeros(1,N); x14d=zeros(1,N); x15d=zeros(1,N); x16d=zeros(1,N);
x17d=zeros(1,N); x18d=zeros(1,N); x19d=zeros(1,N); x20d=zeros(1,N);
xd1=zeros(1,N); xd2=zeros(1,N); xd3=zeros(1,N); xd4=zeros(1,N);
xd5=zeros(1,N); xd6=zeros(1,N); xd7=zeros(1,N); xd8=zeros(1,N);
xd9=zeros(1,N); xd10=zeros(1,N); xd11=zeros(1,N); xd12=zeros(1,N);
xd13=zeros(1,N); xd14=zeros(1,N); xd15=zeros(1,N); xd16=zeros(1,N);
xd17=zeros(1,N); xd18=zeros(1,N); xd19=zeros(1,N); xd20=zeros(1,N);
%----------------------------------------
j=zeros(1,20);j1=zeros(1,20);j2=zeros(1,20);
for i=1:N-1
if x(i)<=minx+10
j(1)=j(1)+1;x1(i)=x(i+1);
elseif x(i)<=minx+20&x(i)>minx+10
j(2)=j(2)+1;x2(i)=x(i+1);
elseif x(i)<=minx+30&x(i)>minx+20
j(3)=j(3)+1;x3(i)=x(i+1);
elseif x(i)<=minx+40&x(i)>minx+30
j(4)=j(4)+1;x4(i)=x(i+1);
elseif x(i)<=minx+50&x(i)>minx+40
j(5)=j(5)+1;x5(i)=x(i+1);
elseif x(i)<=minx+60&x(i)>minx+50
j(6)=j(6)+1;x6(i)=x(i+1);
elseif x(i)<=minx+70&x(i)>minx+60
j(7)=j(7)+1;x7(i)=x(i+1);
elseif x(i)<=minx+80&x(i)>minx+70
j(8)=j(8)+1;x8(i)=x(i+1);
elseif x(i)<=minx+90&x(i)>minx+80
j(9)=j(9)+1;x9(i)=x(i+1);
elseif x(i)<=minx+100&x(i)>minx+90
j(10)=j(10)+1;x10(i)=x(i+1);
elseif x(i)<=minx+110&x(i)>minx+100
j(11)=j(11)+1;x11(i)=x(i+1);
elseif x(i)<=minx+120&x(i)>minx+110
j(12)=j(12)+1;x12(i)=x(i+1);
elseif x(i)<=minx+130&x(i)>minx+120
j(13)=j(13)+1;x13(i)=x(i+1);
elseif x(i)<=minx+140&x(i)>minx+130
j(14)=j(14)+1;x14(i)=x(i+1);
elseif x(i)<=minx+150&x(i)>minx+140
j(15)=j(15)+1;x15(i)=x(i+1);
elseif x(i)<=minx+160&x(i)>minx+150
j(16)=j(16)+1;x16(i)=x(i+1);
elseif x(i)<=minx+170&x(i)>minx+160
j(17)=j(17)+1;x17(i)=x(i+1);
elseif x(i)<=minx+180&x(i)>minx+170
j(18)=j(18)+1;x18(i)=x(i+1);
elseif x(i)<=minx+190&x(i)>minx+180
j(19)=j(19)+1;x19(i)=x(i+1);
else x(i)<=maxx&x(i)>minx+190
j(20)=j(20)+1;x20(i)=x(i+1);
end
end
%--------------------------------
for i=1:N-2
if x(i)<=minx+10
j1(1)=j1(1)+1;x1d(i)=x(i+2);
elseif x(i)<=minx+20&x(i)>minx+10
j1(2)=j1(2)+1;x2d(i)=x(i+2);
elseif x(i)<=minx+30&x(i)>minx+20
j1(3)=j1(3)+1;x3d(i)=x(i+2);
elseif x(i)<=minx+40&x(i)>minx+30
j1(4)=j1(4)+1;x4d(i)=x(i+2);
elseif x(i)<=minx+50&x(i)>minx+40
j1(5)=j1(5)+1;x5d(i)=x(i+2);
elseif x(i)<=minx+60&x(i)>minx+50
j1(6)=j1(6)+1;x6d(i)=x(i+2);
elseif x(i)<=minx+70&x(i)>minx+60
j1(7)=j1(7)+1;x7d(i)=x(i+2);
elseif x(i)<=minx+80&x(i)>minx+70
j1(8)=j1(8)+1;x8d(i)=x(i+2);
elseif x(i)<=minx+90&x(i)>minx+80
j1(9)=j1(9)+1;x9d(i)=x(i+2);
elseif x(i)<=minx+100&x(i)>minx+90
j1(10)=j1(10)+1;x10d(i)=x(i+2);
elseif x(i)<=minx+110&x(i)>minx+100
j1(11)=j1(11)+1;x11d(i)=x(i+2);
elseif x(i)<=minx+120&x(i)>minx+110
j1(12)=j1(12)+1;x12d(i)=x(i+2);
elseif x(i)<=minx+130&x(i)>minx+120
j1(13)=j1(13)+1;x13d(i)=x(i+2);
elseif x(i)<=minx+140&x(i)>minx+130
j1(14)=j1(14)+1;x14d(i)=x(i+2);
elseif x(i)<=minx+150&x(i)>minx+140
j1(15)=j1(15)+1;x15d(i)=x(i+2);
elseif x(i)<=minx+160&x(i)>minx+150
j1(16)=j1(16)+1;x16d(i)=x(i+2);
elseif x(i)<=minx+170&x(i)>minx+160
j1(17)=j1(17)+1;x17d(i)=x(i+2);
elseif x(i)<=minx+180&x(i)>minx+170
j1(18)=j1(18)+1;x18d(i)=x(i+2);
elseif x(i)<=minx+190&x(i)>minx+180
j1(19)=j1(19)+1;x19d(i)=x(i+2);
else x(i)<=maxx&x(i)>minx+190
j1(20)=j1(20)+1;x20d(i)=x(i+2);
end
end
%------------------------------
for i=1:N-3
if x(i)<=minx+10
j2(1)=j(1)+1;xd1(i)=x(i+3);
elseif x(i)<=minx+20&x(i)>minx+10
j2(2)=j(2)+1;xd2(i)=x(i+3);
elseif x(i)<=minx+30&x(i)>minx+20
j2(3)=j(3)+1;xd3(i)=x(i+3);
elseif x(i)<=minx+40&x(i)>minx+30
j2(4)=j(4)+1;xd4(i)=x(i+3);
elseif x(i)<=minx+50&x(i)>minx+40
j2(5)=j(5)+1;xd5(i)=x(i+3);
elseif x(i)<=minx+60&x(i)>minx+50
j2(6)=j(6)+1;xd6(i)=x(i+3);
elseif x(i)<=minx+70&x(i)>minx+60
j2(7)=j(7)+1;xd7(i)=x(i+3);
elseif x(i)<=minx+80&x(i)>minx+70
j2(8)=j(8)+1;xd8(i)=x(i+3);
elseif x(i)<=minx+90&x(i)>minx+80
j2(9)=j(9)+1;xd9(i)=x(i+3);
elseif x(i)<=minx+100&x(i)>minx+90
j2(10)=j2(10)+1;xd10(i)=x(i+3);
elseif x(i)<=minx+110&x(i)>minx+100
j2(11)=j(11)+1;xd11(i)=x(i+3);
elseif x(i)<=minx+120&x(i)>minx+110
j2(12)=j(12)+1;xd12(i)=x(i+3);
elseif x(i)<=minx+130&x(i)>minx+120
j2(13)=j(13)+1;xd13(i)=x(i+3);
elseif x(i)<=minx+140&x(i)>minx+130
j2(14)=j(14)+1;xd14(i)=x(i+3);
elseif x(i)<=minx+150&x(i)>minx+140
j2(15)=j(15)+1;xd15(i)=x(i+3);
elseif x(i)<=minx+160&x(i)>minx+150
j2(16)=j(16)+1;xd16(i)=x(i+3);
elseif x(i)<=minx+170&x(i)>minx+160
j2(17)=j(17)+1;xd17(i)=x(i+3);
elseif x(i)<=minx+180&x(i)>minx+170
j2(18)=j(18)+1;xd18(i)=x(i+3);
elseif x(i)<=minx+190&x(i)>minx+180
j2(19)=j(19)+1;xd19(i)=x(i+3);
else x(i)<=maxx&x(i)>minx+190
j2(20)=j(20)+1;xd20(i)=x(i+3);
end
end
%-------------------
sum=zeros(1,20);sum1=zeros(1,20);sum2=zeros(1,20);
ex=zeros(1,20);ex1=zeros(1,20);ex2=zeros(1,20);
%--------------------------
for i=1:N
if(x1(i)~=0)
sum(1)=sum(1)+x1(i); ex(1)=sum(1)/j(1);
end
if(x2(i)~=0)
sum(2)=sum(2)+x2(i); ex(2)=sum(2)/j(2);
end
if(x3(i)~=0)
sum(3)=sum(3)+x3(i); ex(3)=sum(3)/j(3);
end
if(x4(i)~=0)
sum(4)=sum(4)+x4(i); ex(4)=sum(4)/j(4);
end
if(x5(i)~=0)
sum(5)=sum(5)+x5(i); ex(5)=sum(5)/j(5);
end
if(x6(i)~=0)
sum(6)=sum(6)+x6(i); ex(6)=sum(6)/j(6);
end
if(x7(i)~=0)
sum(7)=sum(7)+x7(i); ex(7)=sum(7)/j(7);
end
if(x8(i)~=0)
sum(8)=sum(8)+x8(i);ex(8)=sum(8)/j(8);
end
if(x9(i)~=0)
sum(9)=sum(9)+x9(i); ex(9)=sum(9)/j(9);
end
if(x10(i)~=0)
sum(10)=sum(10)+x10(i);ex(10)=sum(10)/j(10);
end
if(x11(i)~=0)
sum(11)=sum(11)+x11(i); ex(11)=sum(11)/j(11);
end
if(x12(i)~=0)
sum(12)=sum(12)+x12(i); ex(12)=sum(12)/j(12);
end
if(x13(i)~=0)
sum(13)=sum(13)+x13(i); ex(13)=sum(13)/j(13);
end
if(x14(i)~=0)
sum(14)=sum(14)+x14(i); ex(14)=sum(14)/j(14);
end
if(x15(i)~=0)
sum(15)=sum(15)+x15(i); ex(15)=sum(15)/j(15);
end
if(x16(i)~=0)
sum(16)=sum(16)+x16(i); ex(16)=sum(16)/j(16);
end
if(x17(i)~=0)
sum(17)=sum(17)+x17(i); ex(17)=sum(17)/j(17);
end
if(x18(i)~=0)
sum(18)=sum(18)+x18(i);ex(18)=sum(18)/j(18);
end
if(x19(i)~=0)
sum(19)=sum(19)+x19(i); ex(19)=sum(19)/j(19);
end
if(x20(i)~=0)
sum(20)=sum(20)+x20(i);ex(20)=sum(20)/j(20);
end
end
%-----------------------
for i=1:N
if(x1d(i)~=0)
sum1(1)=sum1(1)+x1d(i); ex1(1)=sum1(1)/j1(1);
end
if(x2d(i)~=0)
sum1(2)=sum1(2)+x2d(i); ex1(2)=sum1(2)/j1(2);
end
if(x3d(i)~=0)
sum1(3)=sum1(3)+x3d(i);ex1(3)=sum1(3)/j1(3);
end
if(x4d(i)~=0)
sum1(4)=sum1(4)+x4d(i);ex1(4)=sum1(4)/j1(4);
end
if(x5d(i)~=0)
sum1(5)=sum1(5)+x5d(i); ex1(5)=sum1(5)/j1(5);
end
if(x6d(i)~=0)
sum1(6)=sum1(6)+x6d(i);ex1(6)=sum1(6)/j1(6);
end
if(x7d(i)~=0)
sum1(7)=sum1(7)+x7d(i);ex1(7)=sum1(7)/j1(7);
end
if(x8d(i)~=0)
sum1(8)=sum1(8)+x8d(i);ex1(8)=sum1(8)/j1(8);
end
if(x9d(i)~=0)
sum1(9)=sum1(9)+x9d(i);ex1(9)=sum1(9)/j1(9);
end
if(x10d(i)~=0)
sum1(10)=sum1(10)+x10d(i); ex1(10)=sum1(10)/j1(10);
end
if(x11d(i)~=0)
sum1(11)=sum1(11)+x11d(i); ex1(11)=sum1(11)/j1(11);
end
if(x12d(i)~=0)
sum1(12)=sum1(12)+x12d(i); ex1(12)=sum1(12)/j1(12);
end
if(x13d(i)~=0)
sum1(13)=sum1(13)+x13d(i);ex1(13)=sum1(13)/j1(13);
end
if(x14d(i)~=0)
sum1(14)=sum1(14)+x14d(i);ex1(14)=sum1(14)/j1(14);
end
if(x15d(i)~=0)
sum1(15)=sum1(15)+x15d(i); ex1(15)=sum1(15)/j1(15);
end
if(x16d(i)~=0)
sum1(16)=sum1(16)+x16d(i);ex1(16)=sum1(16)/j1(16);
end
if(x17d(i)~=0)
sum1(17)=sum1(17)+x17d(i);ex1(17)=sum1(17)/j1(17);
end
if(x18d(i)~=0)
sum1(18)=sum1(18)+x18d(i);ex1(18)=sum1(18)/j1(18);
end
if(x19d(i)~=0)
sum1(19)=sum1(19)+x19d(i);ex1(19)=sum1(19)/j1(19);
end
if(x20d(i)~=0)
sum1(20)=sum1(20)+x20d(i); ex1(20)=sum1(20)/j1(20);
end
end
%-------------------------
for i=1:N
if(xd1(i)~=0)
sum2(1)=sum2(1)+xd1(i);ex2(1)=sum(1)/j2(2);
end
if(xd2(i)~=0)
sum2(2)=sum2(2)+xd2(i);ex2(2)=sum2(2)/j2(2);
end
if(xd3(i)~=0)
sum2(3)=sum2(3)+xd3(i);ex2(3)=sum2(3)/j2(3);
end
if(xd4(i)~=0)
sum2(4)=sum2(4)+xd4(i); ex2(4)=sum2(4)/j2(4);
end
if(xd5(i)~=0)
sum2(5)=sum2(5)+xd5(i); ex2(5)=sum2(5)/j2(5);
end
if(xd6(i)~=0)
sum2(6)=sum2(6)+xd6(i); ex2(6)=sum2(6)/j2(6);
end
if(xd7(i)~=0)
sum2(7)=sum2(7)+xd7(i);ex2(7)=sum2(7)/j2(7);
end
if(xd8(i)~=0)
sum2(8)=sum2(8)+xd8(i);ex2(8)=sum2(8)/j2(8);
end
if(xd9(i)~=0)
sum2(9)=sum2(9)+xd9(i);ex2(9)=sum2(9)/j2(9);
end
if(xd10(i)~=0)
sum2(10)=sum2(10)+xd10(i);ex2(10)=sum2(10)/j2(10);
end
if(xd11(i)~=0)
sum2(11)=sum2(11)+xd11(i);ex2(11)=sum2(11)/j2(11);
end
if(xd12(i)~=0)
sum2(12)=sum2(12)+xd12(i);ex2(12)=sum2(12)/j2(12);
end
if(xd13(i)~=0)
sum2(13)=sum2(13)+xd13(i);ex2(13)=sum2(13)/j2(13);
end
if(xd14(i)~=0)
sum2(14)=sum2(14)+xd14(i); ex2(14)=sum2(14)/j2(14);
end
if(xd15(i)~=0)
sum2(15)=sum2(15)+xd15(i); ex2(15)=sum2(15)/j2(15);
end
if(xd16(i)~=0)
sum2(16)=sum2(16)+xd16(i);ex2(16)=sum2(16)/j2(16);
end
if(xd17(i)~=0)
sum2(17)=sum2(17)+xd17(i);ex2(17)=sum2(17)/j2(17);
end
if(xd18(i)~=0)
sum2(18)=sum2(18)+xd18(i);ex2(18)=sum2(18)/j2(18);
end
if(xd19(i)~=0)
sum2(19)=sum2(19)+xd19(i);ex2(19)=sum2(19)/j2(19);
end
if(xd20(i)~=0)
sum2(20)=sum2(20)+xd20(i);ex2(20)=sum2(20)/j2(20);
end
end
xd0=[5:10:195];
%--------------------------
p1=polyfit(xd0(1:10),ex(1:10),1);
p2=polyfit(xd0(11:16),ex(11:16),1);
x_1=min(xd0(1:10)):0.1:max(xd0(1:10));j2(1)=j(1)+1;
x_2=min(xd0(11:17)):0.1:max(xd0(11:17));
y1=polyval(p1,x_1);
y2=polyval(p2,x_2);
figure(2);plot(xd0,ex,'o',x_1,y1,x_2,y2);xlabel('Xt-1'),ylabel('E^(Xt/Xt-1)');
%------------------------------
p1_1=polyfit(xd0(1:10),ex1(1:10),1);
p2_1=polyfit(xd0(11:16),ex1(11:16),1);
x_1=min(xd0(1:10)):0.1:max(xd0(1:10));
x_2=min(xd0(11:17)):0.1:max(xd0(11:17));
y1_1=polyval(p1_1,x_1);
y2_1=polyval(p2_1,x_2);
figure(3);plot(xd0,ex1,'o',x_1,y1_1,x_2,y2_1);xlabel('Xt-2'),ylabel('E^(Xt/Xt-2)');
%-----------------------------
p1_2=polyfit(xd0(2:9),ex2(2:9),1);
p2_2=polyfit(xd0(11:16),ex2(11:16),1);
x_1=min(xd0(1:10)):0.1:max(xd0(1:10));
x_2=min(xd0(11:16)):0.1:max(xd0(11:16));
y1_2=polyval(p1_2,x_1);
y2_2=polyval(p2_2,x_2);
figure(4);plot(xd0,ex2,'o',x_1,y1_2,x_2,y2_2);xlabel('Xt-3'),ylabel('E^(Xt/Xt-3)');
%----------------------------
%
%
%
%
%⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -