p082105 - 経済統計練習帳

package iron;
 
public class pro {
 
       int s;
  double[] p=new double[101];
double[] y=new double[101];
    double[] pi=new double[101];
    double maxy,miny;
    double h;
    int[] k=new int[101];
    double[] ys=new double[101];
     double[][] x=new double[1000][10];
     double[] sol=new double[100];
     double[][] v=new double[1000][161];
    int[][] gotok=new int[1000][151];
 
    int sx,number;
int[] year=new int[500];
int[] month=new int[500];
 
int[] y1=new int[500];
int[] m1=new int[500];
 
double[] delta=new double[500];
 
 
int[] vv=new int[500];
int[] pp=new int[500];
 
 
 
 
 
int[] op=new int[161];
 
public static void main(String[] args) {
        pro t=new pro();
    }
    pro(){
 
 
        int s;
 
        catchdata sub=new catchdata();
        sub.makedata();
        y=sub.y;
        p=sub.p;
        year=sub.year;
        month=sub.month;
 
 
 
 
 
      pi[1]=0;
      for(s=2;s<85;s++){
      pi[s]=p[s]-p[s-1];
      }
 
 
     maxy=0;
      miny=999;
 
      for(s=1;s<85;s++){
      if(y[s]>maxy)maxy=y[s];    
      if(y[s]<miny)miny=y[s];
      }  
 
      System.out.println(miny+","+maxy);
 
 
 
 
 
      for(s=1;s<85;s++){
      k[s]=100;
      }
 
      int tr;
 
      for(tr=1;tr<11;tr++){
 
 
 
      for(s=1;s<85;s++){
      x[s][1]=1;
      x[s][2]=y[s]-90-0.1*k[s];
      }
 
      ols sub3=new ols();
      sub3.datanumber=84;
      sub3.number=2;
      sub3.x=x;
      sub3.y=pi;
      sub3.makedata();
      sol=sub3.sol;
 
      System.out.println(sol[1]+","+sol[2]);
 
 
 for(s=1;s<150;s++){
 double p1=sol[1]+sol[2]*(y[1]-90-0.1*s);
 v[1][s]=(p1-pi[1])*(p1-pi[1]);
   }
 
   int t,sx,sxs;
  double u1,v1,minv;
 
 for(t=2;t<85;t++){
 for(s=1;s<150;s++){  
 double p1=sol[1]+sol[2]*(y[t]-90-0.1*s);
u1=(p1-pi[t])*(p1-pi[t]);
minv=9999;    
sxs=0;
for(sx=1;sx<150;sx++){   
    v1=u1+v[t-1][sx];
    double e=(s-sx)*(s-sx);
    if(e>50)v1=9999;
    if(v1<minv)sxs=sx;
    if(v1<minv)minv=v1;
    }
v[t][s]=minv;
gotok[t][s]=sxs;
}
}   
 
 
 
     minv=9999;
        sx=0;
         for(s=1;s<150;s++){
         if(v[84][s]<minv)sx=s;
         if(v[84][s]<minv)minv=v[84][s];
         }
 
        op[84]=sx;
         for(s=1;s<84;s++){
         sx=op[84-s+1];
         op[84-s]=gotok[84-s+1][sx];
         }
 
 
          for(s=1;s<85;s++){
          k[s]=op[s];
          }
 
 
 
 
 
    }      
 
 
 
      int n;  
 
      for(n=1;n<6;n++){    
 
         double[] sig=new double[5];
 
      for(s=1;s<85;s++){
      x[s][1]=1;
      x[s][2]=y[s]-90-0.1*k[s];
      }
 
      double[] delta=new double[100];
 
      for(s=1;s<84;s++){
      delta[s]=0.1*(k[s+1]-k[s]);
      }
      delta[84]=0;
 
      ols sub4=new ols();
      sub4.datanumber=84;
      sub4.number=2;
      sub4.x=x;
      sub4.y=delta;
      sub4.makedata();
      sig=sub4.sol;
 
           System.out.println(sig[1]+","+sig[2]);
 
 
 
      for(tr=1;tr<11;tr++){
 
 
 
      for(s=1;s<85;s++){
      x[s][1]=1;
      x[s][2]=y[s]-90-0.1*k[s];
      }
 
      ols sub31=new ols();
      sub31.datanumber=84;
      sub31.number=2;
      sub31.x=x;
      sub31.y=pi;
      sub31.makedata();
      sol=sub31.sol;
 
      System.out.println(sol[1]+","+sol[2]);
 
 
 for(s=1;s<150;s++){
 double p1=sol[1]+sol[2]*(y[1]-90-0.1*s);
 v[1][s]=(p1-pi[1])*(p1-pi[1]);
   }
 
   int t,sx,sxs;
  double u1,v1,minv;
 
 for(t=2;t<85;t++){
 for(s=1;s<150;s++){  
double p1=sol[1]+sol[2]*(y[t]-90-0.1*s);
u1=(p1-pi[t])*(p1-pi[t]);
double u2=0.1*(k[t]-k[t-1])-sig[1]-sig[2]*(y[t]-90-0.1*k[t-1]);
u1=u1+u2*u2;
minv=9999;    
sxs=0;
for(sx=1;sx<150;sx++){   
    v1=u1+v[t-1][sx];
    double e=(s-sx)*(s-sx);
    if(e>250)v1=9999;
    if(v1<minv)sxs=sx;
    if(v1<minv)minv=v1;
    }
v[t][s]=minv;
gotok[t][s]=sxs;
}
}   
 
 
 
     minv=9999;
        sx=0;
         for(s=1;s<150;s++){
         if(v[84][s]<minv)sx=s;
         if(v[84][s]<minv)minv=v[84][s];
         }
 
        op[84]=sx;
         for(s=1;s<84;s++){
         sx=op[84-s+1];
         op[84-s]=gotok[84-s+1][sx];
         }
 
 
          for(s=1;s<85;s++){
          k[s]=op[s];
          }
 
 
 
     }      
 
 
      }
 
 
 
         for(s=1;s<85;s++){
         System.out.println(year[s]+"年"+month[s]+"月,"+(90+0.1*k[s])+","+y[s]);
          }
 
 
 
 
 
 
 
 
 
 
 
    }
 
 
     int numx(String str){
 
    int nx;
 
    nx=0;
 
 try{
 
    nx=Integer.parseInt(str);
 
 }catch (Exception e) {}
 
 return nx;
 
 
 
}
 
 
 
 
 
double num(String str){
 
    double nx;
 
    nx=0;
 
 try{
 
    nx=Double.parseDouble(str);
 
 }catch (Exception e) {}
 
 return nx;
 
 
 
}
 
 
 
}