php dsge 13 data - プログラム倉庫

<?php  
$handle = fopen("php131.txt", "r");
$num = 1;
while (($csvdata = fgetcsv($handle, 100000, ",")) !== FALSE) {  
    for ($rank=0; $rank < 6; $rank++) {
$data[$num] = $csvdata[$rank];
$num = $num +1;      
}
}
fclose($handle);
for ($i=1; $i < 11001; $i++) {
$num=6*($i-1);
$s=$data[$num+1];
$m=$data[$num+2];
$n=$data[$num+3];
$lx[$s][$m][$n]=$data[$num+4];
$cx[$s][$m][$n]=$data[$num+5];
$px[$s][$m][$n]=$data[$num+6];
}
$handle = fopen("php132.txt", "r");
$num = 1;
while (($csvdata = fgetcsv($handle, 100000, ",")) !== FALSE) {  
    for ($rank=0; $rank < 5; $rank++) {
$data[$num] = $csvdata[$rank];
$num = $num +1;      
}
}
fclose($handle);
for ($i=1; $i < 110001; $i++) {
$num=5*($i-1);
$s=$data[$num+1];
$m=$data[$num+2];
$n=$data[$num+3];
$m1=$data[$num+4];
$gotos[$s][$m][$n][$m1]=$data[$num+5];
}
for ( $m=1; $m<=10; $m++ ){
$th[$m] = 0.95 + 0.01 * $m;
}
$phi=0.9;
$beta = 0.95;
$a = 0.33;
$ls = ((1 - $a) * $phi) / ((1 - $a) * $phi + 1);
$ks = $ls * pow((1 / $beta - 1) / ($a * $phi) ,1 / ($a - 1));  
$h = 2 * $ks / 100;
for ($n = 1; $n<=100; $n++){
$k[$n] = $n * $h;
}
$t=1;
$kt[$t]=$k[40];
$mp=5;
for ( $t=1; $t<=100; $t++ ){
$m = rand(1, 10);
$tht[$t] = $th[$m];
$s1=$gotos[$s2][$mp][$n2][$m];
$s2=floor($s1);
$n1 = $kt[$t] / $h;
$n2 = floor($n1);
$n3 = $n2 + 1 ;
$s2=floor($s1);
$s3=$s2+1;
$dn=($n1 - $n2) * ($cx[$s2][$m][$n3] - $cx[$s2][$m][$n2]);
$ds=($s1 - $s2) * ($cx[$s3][$m][$n2] - $cx[$s2][$m][$n2]);
$ct[$t] = $cx[$s2][$m][$n2] +$dn+$ds;
$dn= ($n1 - $n2) * ($lx[$s2][$m][$n3] - $lx[$s2][$m][$n2]);
$ds=($s1 - $s2) * ($lx[$s3][$m][$n2] - $lx[$s2][$m][$n2]);
$lt[$t] = $lx[$s2][$m][$n2] +$dn+$ds;
$dn= ($n1 - $n2) * ($px[$s2][$m][$n3] - $px[$s2][$m][$n2]);
$ds=($s1 - $s2) * ($px[$s3][$m][$n2] - $px[$s2][$m][$n2]);
$pt[$t]=$px[$s2][$m][$n2] +$dn+$ds;
$kt[$t+1] = $kt[$t] + $tht[$t]*pow($kt[$t] , $a) * pow($lt[$t] ,1 - $a) - $ct[$t];
$mp=$m;
}
for ( $t=1; $t<=98; $t++ ){
print($pt[$t]);
print(",");
}
print($pt[99]);
?>