cdbhgfx2.asv

第一次上载程序

clear

%%%%%%%%%%%%%%%           DGBstu2

%%%   每月回归

%%%%% 岳城水库1962年--1990年水位和降雨量数据。


% jssjly2=[
% 1	12.5	1.5	9.7	21.3	60.3	189.5	108.8	136.7	34	56.5	0.6
% 0	0.3	22.1	33.1	117.5	39.2	189	325	95.3	2.6	16	9.5
% 14.9	12.1	17.9	66.1	73.3	30.7	238.2	86.2	118	72.3	7.4	2.1
% 2.3	6.2	11.4	64.8	20.3	34	87.5	54	27.8	27.6	11.4	0.7
% 1.1	6.9	19.4	31.8	14.1	126.8	234.7	133.4	31.8	30.7	8.1	1.8
% 10.9	13.9	22.7	30.1	25	88.2	109.2	209.1	140.7	5.6	46.6	0
% 3.5	0.2	1.7	28.7	21	30.9	160.6	115	79.8	111.4	28.4	7.9
% 5.3	7.5	12.4	69.6	21.8	30.2	137.2	77	153.3	15.8	24.4	0
% 0.6	10	9.6	36.9	75.9	65.3	155.8	86.5	38.7	19.8	1.9	1.1
% 10.4	7.2	12.1	34.3	28.7	183.2	151.8	207.2	77.3	17.4	41.4	12
% 16.1	11.1	8.6	9.7	17.1	34.5	164.9	130.4	40.7	23.4	15	3.6
% 7.6	3.7	5.7	19.8	33.3	149.2	165.3	185.9	76.4	98.9	1.6	0
% 7.7	6.2	17.3	8.3	22.1	27.4	141.7	88.6	69.4	40.6	40.5	19.2
% 1.5	3.2	7.8	51.6	11.3	57.6	215.7	99	133.2	57.1	6.6	14.7
% 0	30.6	8.4	55.8	4.9	35.7	236.2	188.1	62	27.9	16	4.9
% 3.1	0.4	3.6	28.9	53.1	67.3	194.2	134.3	24.9	30	13.8	9.1
% 0.1	15.4	9.2	2.7	34.5	30.2	213.2	85	70.4	41.5	8.9	4.4
% 3.1	12.3	37.6	30.6	11	79.9	138.2	69.6	52.7	4.4	2.2	8.7
% 1.9	0.9	34.8	16.6	56.6	101.8	90.1	91.3	50.7	31.3	0.9	0
% 4.9	1.8	13.9	14.6	7.2	50.5	148.6	142.3	38.2	11.5	8.9	0.3
% 0.6	19.3	20.5	9.5	24.9	49.6	174.8	242	59.1	5.4	17.9	0
% 0.6	3.2	22.8	25.8	104.1	29.1	89.3	74.8	99.9	65.9	2.9	0
% 0.5	2.1	10.3	10.4	69.1	94.3	102.1	81.7	91.8	20.9	16.3	17.4
% 3.5	10.3	6.8	4	108.1	55.2	125.3	104.5	173.2	49.9	0.1	1.5
% 0.1	2.3	14	14.6	55.6	47.8	70.5	50	35	43.2	1.8	14.1
% 3.1	7.9	22.9	18.6	48.2	132.8	63.7	125.7	27.9	60	17	0
% 0.6	6.2	15.5	4.7	70.2	42.8	278.7	184.7	38.6	7.6	0	2
% 18.5	7.7	22.1	11.6	42.8	84.3	150.1	123.6	53	20.6	14.8	13.7
% 14.2	29.8	58.9	51.3	66.7	105.3	144.6	98.7	66.1	9.9	16.8	1.6
%     
% ];
% 
% sksw1=[
% 130.71	131.76	132.16	128.33	114.32	114.83	112.36	113.67	111.2	112.23	115.51	119.31
% 126.67	130.24	131.21	131.67	129.69	127.95	116.71	138.57	137.82	142.74	140.61	135.53
% 136.13	139.21	140.17	133.93	119	113.91	112.66	120.95	133.64	133.68	132.95	139
% 139.17	139.13	139.73	138.39	137.22	123.57	117.22	115.75	115.45	115.23	118.03	117.23
% 117.43	121.24	120.81	116.62	115.28	115.24	119.73	135.47	143.74	142.77	142.57	140.82
% 141.53	142.61	141.59	139.82	136.36	131.44	126.83	132.27	144.56	145.33	143.73	139.65
% 141.53	143.27	144.77	143.16	135.78	123.49	120.54	128.87	132.59	138.13	144.51	145.58
% 146.07	147.05	147.91	146.55	144.78	136.25	124.75	129.97	133.79	139	141.55	142.77
% 144.11	144.97	142.37	135.55	131.75	128.4	125.92	133.88	137.15	135.8	135.67	134.3
% 132.76	132.41	130.76	128.25	123.49	122.09	133.87	139.84	147.89	145.71	146.74	148.05
% 148.86	148.78	148.82	145.16	136.99	130.11	128.47	130.2	134.78	135.79	136.72	137.73
% 138.74	138.7	134.76	128.99	127.07	127.33	135.62	140.92	147.99	145.79	145.57	146.2
% 146.41	147.85	146.99	142.12	137.54	132.42	120.03	122.17	126.85	128.05	129.89	131.49
% 132.65	133.01	124.96	116.62	116.69	116.51	119.01	134.87	135.42	139.36	143.75	142.29
% 142.71	143.32	143.82	139.84	132.48	118.34	122.03	135.96	147.95	148.1	148.36	147.32
% 147.58	148.28	143.88	136.77	133.05	128.28	129.48	138.53	142.41	140.17	141.24	141.55
% 141.69	141.81	139.94	132.08	123.9	122.32	122.63	129.16	129.68	124.55	125.56	125.52
% 125.63	126.11	126.16	126.34	125.39	121.07	123.09	123.14	122.36	122.82	123.06	122.37
% 119.99	118	118	118	118	118	118	119.49	121.72	122.34	122.75	122.44
% 122.7	123.06	119.42	118	118	118	119.27	123.35	127.32	126.93	127.03	126.97
% 126.26	125.42	124.14	121.81	121.81	120.17	119.27	144.55	148.17	148.17	146.94	144.33
% 144.21	144.76	144.23	140.31	140.31	138.09	132	127.66	128.36	128.76	130.35	130.97
% 131.24	131.56	131.63	131.62	131.62	131.21	130.82	129.65	128.88	129.32	130.62	131.19
% 131.76	132.17	132.39	132.45	132.45	132.12	129.76	129.77	133.22	140.37	142.29	143
% 143.4	143.61	142.56	141.21	141.21	136.83	133.76	132.49	132.06	132	131.84	131.17
% 130.57	129.95	129.06	128.07	128.07	125.86	125.35	123.97	126.05	126.14	126.2	126.12
% 125.66	125	123.97	123.03	123.03	123	125.5	138.41	146.43	145.9	145.35	143.47
% 142.96	143.65	143.73	142.16	142.16	136.89	131.67	132.45	136.17	135.04	134.77	134.37
% 134.65	135.79	136.95	138.11	138.11	141.09	133.28	132.8	135.54	136.25	136.59	136.97
%     
% ];
% 
% 
% skjs3=[
% 0	4.9	1.9	7	0.2	12.9	110.2	108	143	30.9	56.1	0.2
% 0	1.2	25.8	18.1	131	48.2	147.8	635.6	61.3	0.5	8.6	18.2
% 14.2	10.8	8.9	125.3	70.8	9.7	288.8	109.3	48	73.6	4	2.1
% 2	8.9	3.6	47.4	12.7	5	65.2	30.9	31.1	23.4	11.8	0
% 0	6.5	13.1	23.6	6.8	32.5	222.9	50.6	15.6	17.3	17.4	0.1
% 8.2	14	16.6	24	32.9	73.7	76.6	127.3	86.9	0	40.8	0
% 0.8	0.1	0	31.6	15.5	49.9	132.2	86.2	57.3	88.6	27	19.2
% 1.9	4.8	5.3	66.8	23.7	31.9	129.7	183.4	148.1	10.7	12.4	0
% 0.3	1.4	2	21.6	72.7	42.3	111.3	62.6	19.5	16.7	3.5	0
% 0	7.5	10.8	32	21.8	161.4	109.4	143.7	11.2	17.3	36.8	14.7
% 13.1	6.1	2.4	7.6	48.6	0.4	251.6	146.4	49.5	31	35.6	0.2
% 8.4	4.5	0	27.6	11.5	112.1	239.4	98.9	62.9	79.9	6	0
% 5.8	5.6	13.9	13.9	43.7	112	63.5	121.3	64.7	43.5	44.7	21.5
% 0	0.3	0.6	32.6	0	41.9	161.7	160.7	53.3	44.9	4.6	6.9
% 0	35.9	12.1	29.2	1	36.7	343.6	99.8	84.1	12.7	30.9	1.5
% 0.9	0	4.7	23.6	19.5	76	184.1	26.7	2.7	32.3	20.9	0.5
% 0	8.7	8.9	0	47.5	36.4	178.7	16.7	10.6	34.9	7	2.5
% 6.1	12.1	19	32.1	8.7	76.3	141.2	8.7	47.5	0.4	3	6.9
% 0.7	0.2	17.1	17	59.4	74.4	97.7	38.9	73.9	41.7	0.6	0
% 4.4	0	17.9	1.3	1.9	34.5	99.4	63.9	38.1	7.9	5.4	0
% 0	16.7	9.1	6.8	54.7	12.2	149	517.9	19	14.8	11.7	0
% 0	1.4	29.1	29.8	74.3	16.8	73.5	46.3	112.8	70.2	0.5	0.2
% 0	0	6.8	11.9	29.7	103.7	80.9	248.9	49.5	8.9	8.8	13.2
% 2.1	6.6	3	3	61.5	34.9	71.2	120.4	141.9	67	0	0.8
% 0	0	7.7	3.6	35.7	81.6	42	28.5	17.1	33	1.1	3.1
% 4.1	4.2	25	13.3	30.2	100	58.3	113.9	29	75.4	13.8	0
% 0	0.7	11.8	0	78.3	12.5	182.7	124.2	3.7	9.1	0	0.6
% 19.5	3.3	22.3	11.1	47.5	69.8	168.7	47.1	11	3.8	14.9	10.9
% 18.8	4.2	44.3	62	85.8	57.2	82.3	116.9	35.4	2.4	31.2	0
%     
% ];
% 
% 
% tqll4=[
% 16.4	22.6	21.3	17.8	18.3	16.7	61.6	80.9	44.2	68.4	42.3	35.5
% 28.1	22.7	26.5	29.2	30.6	84.4	109	156	130	84.8	51	39.8
% 32.8	32.1	38.3	35.9	34.3	52	91.4	75.6	74.6	57.7	54.7	47.2
% 30.9	28.7	20	26.7	40	22.7	21.9	10.7	6.74	6.76	14.8	7.04
% 7.35	11.5	5.9	3.36	2.54	5.09	87.2	114	45.2	20.9	16.4	20.7
% 13.5	14	14.7	21.4	29.4	13.2	18.4	105	111	66.2	28.5	26.8
% 17.7	16.2	16.3	8.05	11.8	16.3	33	34.5	34.5	71	34.4	18.2
% 10.7	10.8	17.2	20.3	28.7	5.65	27.1	33.3	20.1	19.8	9.7	14.6
% 12.6	9.42	7.71	4.29	19.4	15.2	72.1	64.1	9.66	5.79	5.32	4.68
% 4.92	3.25	8.62	6.71	5.63	22.8	111	135	118	19.7	30.1	27.9
% 22.1	19.1	25.6	23.5	6.38	2.02	32.7	21	27.6	13.7	8.5	16.5
% 6.43	4.78	2.66	6.3	4.51	30.6	50.5	46.8	75.9	112	33.8	27.1
% 13.4	16.2	15.7	16.1	6.57	4.05	10.2	20.4	4.08	4.17	2.92	2.3
% 2.54	3.71	2.86	2.12	4.5	3	59.2	67.8	25.3	32.2	32.3	14.5
% 9.88	5.66	7.81	20.2	7.78	3.84	71.2	164	87.1	23.9	20.9	26.6
% 13.4	12.4	8.02	5.58	6.81	8.64	37	126	30.2	6.37	7.08	9.19
% 8.62	4.45	3.59	2.89	1.31	1.17	23	12.2	44.3	8.8	6.04	3.9
% 5.19	3.92	3.75	3.73	1.29	3.95	16.1	5.84	2.99	2.72	1.21	1.27
% 1.44	2.56	1.72	1.95	1.2	5.77	5.22	10.8	2.36	1.99	1.89	2.36
% 2.07	1.85	1.55	1.61	1.62	2.48	8.88	23.5	5.83	3.12	1.35	1.4
% 1.69	1.86	1.59	1.21	1.01	2.55	6.08	65.4	28.4	9.83	3.42	2.54
% 1.79	1.74	1.67	7.4	4.49	10.7	2.47	4.14	5.77	3.78	2.23	1.2
% 1.19	1.61	1.82	1.47	1.45	6.35	5.5	2.95	3.59	2.59	2.25	1.91
% 2.4	2.39	2.19	1.24	1.83	1.38	4.49	5.02	39.4	25.2	7.47	3.31
% 2.76	2.67	2.37	12	2.12	1.88	3.82	2.51	2.06	2.32	2.75	1.8
% 2.99	1.8	1.91	2.84	2.23	2.86	3.12	7.42	2.31	1.91	1.86	1.85
% 2.38	1.81	1.65	2.08	1.29	1.66	22.5	111	29.3	11	3.24	1.91
% 2.25	2.62	2.73	2.59	7.08	7.69	33.6	36.2	12.7	5.15	2.62	3.11
% 3.71	4.79	7.44	11.2	19.5	20.5	30.1	31.9	8	6.7	5.61	3.28
%     
% ];
% 
% 
% 
% 
% y1ys=[];
% x2ys=[];
% x3ys=[];
% z4ys=[];
% 
% for i=1:1:length(sksw1)
%     y1ys=[y1ys sksw1(i,:)];
%     x2ys=[x2ys jssjly2(i,:)];
%     x3ys=[x3ys skjs3(i,:)];
%     z4ys=[z4ys tqll4(i,:)];
% end
% 
% 
% i=1:1:length(y1ys);y1(i)=(y1ys(i)-min(y1ys)+mean(y1ys))/(max(y1ys)-min(y1ys));
% i=1:1:length(x2ys);x2(i)=(x2ys(i)-min(x2ys)+mean(x2ys))/(max(x2ys)-min(x2ys));
% i=1:1:length(x3ys);x3(i)=(x3ys(i)-min(x3ys)+mean(x3ys))/(max(x3ys)-min(x3ys));
% i=1:1:length(z4ys);z4(i)=(z4ys(i)-min(z4ys)+mean(z4ys))/(max(z4ys)-min(z4ys));
% 
% 
% for i=1:1:length(sksw1)
%     j=1:1:12;y1jz(i,j)=y1((i-1)*12+j);
%     j=1:1:12;x2jz(i,j)=x2((i-1)*12+j);
%     j=1:1:12;x3jz(i,j)=x3((i-1)*12+j);
%     j=1:1:12;z4jz(i,j)=z4((i-1)*12+j);
% end

%%%%%%%%%%%%%%%%%%%%%%%%%%%
corrcoef(y1,z4)   ;  %  xiang guan xi shu


%%%%%%%%%%%%%%%%%%%%%%%%%%  线性回归分析
eee= ones(1,length(z4jz(:,1)))  ;

for ij=1:1:12
  YSJ = y1jz(:,ij)';
  

% XX =[eee'  x2jz(:,ij)  x3jz(:,ij)];  

XX =[eee'  x2jz(:,ij)  x3jz(:,ij)  z4jz(:,ij)];

%%%%%

HGXS1 =(YSJ*XX)/(XX'*XX);
y1ycjz(:,ij)= XX*HGXS1';

%%%%%%%%%%%%%%%%%%%%%%%%%  线性回归分析完毕

% i=1:1:length(y1);plot(i,y1,i,y1yc,'r'),grid on


%%%%%%%%%%%%%%%%%%%%%%%%%%  二阶非线性回归分析
XX22=[];

m=length(XX(1,:));
M= m +m*(m-1)/2;  % 二阶多项式项数

XX22(:,1)=eee';
for i=1:1:length(XX(:,1))
      for j=1:1:m-1
          XX22(i,j+1)=XX(i,j+1);
      end    
      for j=1:1:m-1
          XX22(i,j+m)=XX(i,j+1)^2;
      end    
end
for i=1:1:length(XX(:,1))
      l=2*m;
            for j=1:1:m-1
              for k=(j+1):1:m-1
                XX22(i,l)=XX(i,j+1)*XX(i,k+1);
                l=l+1;
              end
           end
 end            
            
HGXS2 =(YSJ*XX22)/(XX22'*XX22);

y1yc22jz(:,ij)= XX22*HGXS2';

%%%%%%%%%%%%%%%%%%%%%%%%%%  二阶非线性回归分析完毕

end
%%%%%%%%%
y1yc=[];
for i=1:1:length(y1jz(:,1))
    y1yc= [y1yc y1ycjz(i,:)];
end

for i=1:1:length(y1)
    wch1(i)= 100*(y1(i)-y1yc(i))/y1(i);
end
%%%%%%%%%

y1yc22=[];
for i=1:1:length(y1jz(:,1))
    y1yc22= [y1yc22 y1yc22jz(i,:)];
end

for i=1:1:length(y1)
    wch2(i)= 100*(y1(i)-y1yc22(i))/y1(i);
end

%%%%%%%%%


i=1:1:length(y1);plot(i,y1,i,y1yc22,'r'),grid on            
    
i=1:1:length(y1);plot(i,wch1,'g',i,wch2,'r'),grid on


i=1:1:length(y1);plot(i,y1,i,y1yc,'g',i,y1yc22,'r'),grid on   


[mean(abs(wch1))  mean(abs(wch2))]