先日のTokyoRで「(続)Rで力学系」というタイトルで発表させていただきました。 発表では残念ながら時間切れでStandard mapのホモクリニック軌道を重み付け(lyapunov weighted sampling)によって探索する手法の紹介までには至りませんでした。やや不完全な部…

多次元データの回帰、判別などの多変量解析においてベクトルとして表現されるデータを 処理しやすい空間(もとのデータより次元の高い場合が多い)に写像する方法としてカーネル法が知られています。 カーネル法は直感的にはベクトルの内積を拡張したものと見…

レコード大賞には間にあいませんが、紅白のまえにごり押しました。 こちらの画像を使いました。 http://akb48matomemo.com/wp-content/uploads/a14738fa6ac358c1f3e90347f2e32e91.jpg 手順 まず単純に48等分します。 split.cpp #include <iostream> #include <string> #include <stdio.h></stdio.h></string></iostream>…

アンサイクロペディアのPythonの項に感銘を受けたので書いてみました。 ramanujan1 = head [y|(y,i,j)<-pre, (x,k,l)<-pre, y==x,i/=k ,j/=l,i/=l] where pre=[(i^3+j^3,i,j)|i<-ll, j<-ll] where ll=[1..] しかし無限リストにすると数分待っても計算が終わ…

C++の行列演算ライブラリのEigenを使うと比較的簡単にアルゴリズムを実装することができます。 MDS(多次元尺度法)を実装してみました。 #include <eigen3/Eigen/Dense> #ifdef _MDS_WITH_REDSVD #include <redsvd/redsvd.hpp> #else #include <eigen3/Eigen/SVD> #endif using namespace Eigen; //d: distance matrix o</eigen3/eigen/svd></redsvd/redsvd.hpp></eigen3/eigen/dense>…