2013-10-01から1ヶ月間の記事一覧
SUMOのデータ構造。 ネットワーク（道路と交差点の集合） http://sumo-sim.org/userdoc/Networks/SUMO_Road_Networks.html 道路はネットワークである。 ノードがジャンクション（交差点） エッジが道路。 道路は、レーン（車線）が定義されている。 レーンに…
SUMOの交通シミレーションを参考に道路などｦ作った後に｡｡ http://d.hatena.ne.jp/arupaka-_-arupaka/20131008 http://sumo-sim.org/userdoc/Demand/Introduction_to_demand_modelling_in_SUMO.html demandモデリングを参照して｡｡車の数を増やす。 p=-0.001 p…
RとCの連携の基本は、 Windows上でのRとCの連携 http://d.hatena.ne.jp/arupaka-_-arupaka/20131020/1382252302 を参考に。 具体的な手順 (1)まず、準備。mingwを利用して、cygwin上でcygwinに依存しないgslを作る。 (Rはcygwinのライブラリを参照できないの…
(1)まず、gslのソースをもってくる http://ftp.jaist.ac.jp/pub/GNU/gsl/ (2)次にコンパイル ./configure --host=x86_64-w64-mingw32 --prefix=/cygdrive/c/gsl CC=x86_64-w64-mingw32-gcc.exe CXX=i686-pc-mingw32-g++.exe make make install configure の…
まず、Rには、C言語との連携に .C関数と.Call 関数が用意されている。 .C関数は簡便用で、.Call関数は、C言語からRのオブジェクトを操作できる。ここでは、.C関数を利用する。以下のサイトを参考にしました。http://www.okada.jp.org/RWiki/?R%A4%AB%A4%E9%C…
http://www.okada.jp.org/RWiki/?R%A4%AB%A4%E9%C2%BE%B8%C0%B8%EC%CD%F8%CD%D1を参考に (1)Mingw を導入 (2) PATHを通す. システム→システムの詳細設定→環境変数 PATH ... (3)Rtoolsを導入 http://essrc.hyogo-u.ac.jp/cran/ をダウンロード (4)ダウンロー…
gslでBスプラインでデータの平滑化. Bスプイラインは自然スプラインより計算が安定。 Bスプライン（拡張スプライン）同じところにノットを重ねることで 不連続を表現できる、３次だと、１回重ねると２回微分、２回重ねると１回微分、 ３回重ねると関数そのも…
汎用ポインタ *void の関数を構造体と組み合わせて使う。構造体をポインタとして操る。 struct points *a; ポインタ状態のときは->演算子で操作。 a->x=3; aにはアドレスが入っている。アドレスを渡す。 void my_func(void *a2){キャスト struct points* a=(…
例：GSLで最適化 http://www.gnu.org/software/gsl/manual/html_node/Multimin-Examples.html#Multimin-Examples と http://www.gnu.org/software/gsl/manual/html_node/Providing-a-function-to-minimize.html#Providing-a-function-to-minimize を参考に。…
C++で自然スプライン補完 #include<stdio.h> #include<gsl/gsl_linalg.h> void f(double* x,int size,double* x_knots,int knots_size,double* a,double* b,double* c,double* d,double* output){ //x should be sorted; int j=0; double r3; double r2; double r1; double r0; int k=0;</gsl/gsl_linalg.h></stdio.h>…
gslで線形方程式、連立一次方程式をとく。 #include<stdio.h> #include<gsl/gsl_linalg.h> void solve(double** m,double* v,int len){ using namespace std; int s; double *mb; int k=0; mb=(double*)malloc(sizeof(double)*len*len); for(int i=0;i</gsl/gsl_linalg.h></stdio.h>
http://inet-lab.naist.jp/sumo_define_net_and_route/ を全面に参考 道路生成 netgenerate --grid --grid.length 100 --grid.number 10 --default.lanenumber 2 --default.speed 11.1111 --tls.guess --output-file manhattan.net.xml 出発地と目的地の生成…
DoE.baseで直交表の自動生成(実験計画法) library("DoE.base") name<-c("Taro","Jiro","Hanako","Honda") gender<-c("male","female") year<-c("1964","1852") amada<-c("1","2","3") oaTable<-oa.design(factor.names=list(name=name,genger=gender,year=ye…
重みつきスコア法によるポアソンパラメータのロバスト推定 参考（べき分布の場合） http://d.hatena.ne.jp/arupaka-_-arupaka/20130616/1371343477 #Saiyuuhou z<-rpois(10000,lambda=5) f<-function(lambda){ -sum(log(dpois(z,lambda))) } optimize(f,uppe…
c1 #c1 #c1 #qcsv("data/sikousakugo/siokouji_s.csv")[,3] n z x #x x1 yy #yy #y1 #plot(x1,yy) knot f12 #y# plot(x,y) h len v m for(i in 1:(len-1)){ if(i==1){ m[i,] } if(i>1 & i*1/6/diff(x_knots) d c f sum for(i in 1:(length(x_knots)-1)){ sum…
c1 #c1 #c1 #qcsv("data/sikousakugo/siokouji_s.csv")[,3] n z x #x x1 yy #yy #y1 #plot(x1,yy) knot f12 #y# plot(x,y) h len v m for(i in 1:(len-1)){ if(i==1){ m[i,] } if(i>1 & i*1/6/diff(x_knots) d c f sum for(i in 1:(length(x_knots)-1)){ sum…
