Skip to content

Instantly share code, notes, and snippets.

View naoyat's full-sized avatar

naoya_t naoyat

116 followers · 104 following
View GitHub Profile
@naoyat
naoyat / gist:4065596
Created November 13, 2012 12:48
http://d.hatena.ne.jp/n_shuyo/20091018/gaussian のコードのOctave訳
% PRML's synthetic data set
curve_fitting.x = [0.000000 0.111111 0.222222 0.333333 0.444444 0.555556 0.666667 0.777778 0.888889 1.000000];
curve_fitting.t = [0.349486 0.830839 1.007332 0.971507 0.133066 0.166823 -0.848307 -0.445686 -0.563567 0.261502];
function [product] = outer(a, b, fn)
R = length(a);
C = length(b);
product = zeros(R, C);
for r = 1:R
for c = 1:C
@naoyat
naoyat / gist:4065608
Created November 13, 2012 12:50
PRML図6.8の再現
1;
clear
global theta = [0.33 25 1 0.03];
beta = 50
x = linspace(0, 1, 100);
data = load("curvefitting.txt");
data = data(1:7, :); %% 一番右から3つを捨てる
@naoyat
naoyat / gist:4065614
Created November 13, 2012 12:51
PRML図6.8のパラメータをsimulated annealingで求めたい
1; %% これがなくて関数定義から始まるとOctaveではファイル名と関数名が違うぞと怒られるので、スクリプトは1;から始める
function [score] = score68(theta)
global data_x;
global data_t;
global N;
global mu;
M = 6;
x = linspace(0, 1, M);
@naoyat
naoyat / gist:4065620
Created November 13, 2012 12:53
PRML図6.4図6.5の再現
1;
clear
fig_id = 6;
thetas_6_5 = [ 1 4 0 0;
9 4 0 0;
1 64 0 0;
1 0.25 0 0;
1 4 10 0;
@naoyat
naoyat / gist:4065623
Created November 13, 2012 12:53
PRML図6.9の再現
1;
clear
%global eta = [1 1]; %% 左図
global eta = [1 0.01]; %% 右図
global theta0 = 0.5;
N = 30; %% N=30は割と時間がかかる
xl = linspace(-1, 1, N);
@naoyat
naoyat / Makefile.diff
Created November 17, 2012 14:20
Mac OS X LionでGNU Octave最新版を使いたい ref: http://qiita.com/items/9bfa6c4e3d5119bcd5ee
--- octave-3.6.3/Makefile 2012-11-17 22:31:59.000000000 +0900
+++ octave-3.6.3-naoyat/Makefile 2012-11-17 21:54:30.000000000 +0900
@@ -242,7 +242,8 @@
ETAGS = etags
CTAGS = ctags
DIST_SUBDIRS = $(SUBDIRS)
-DISTFILES = $(DIST_COMMON) $(DIST_SOURCES) $(TEXINFOS) $(EXTRA_DIST)
+# DISTFILES = $(DIST_COMMON) $(DIST_SOURCES) $(TEXINFOS) $(EXTRA_DIST)
+DISTFILES = $(DIST_COMMON) $(DIST_SOURCES) $(EXTRA_DIST)
distdir = $(PACKAGE)-$(VERSION)
@naoyat
naoyat / evernote_markdown.scpt
Created November 22, 2012 16:57
Evernoteでmadeverでmarkdown、但しmdソースは温存 ref: http://qiita.com/items/1dfb0e0ff7ac6b1bb2bd
property src_notebook : "markdown.src"
property dest_notebook : "markdown.html"
global madever
global CSS
set madever_dir to (path to scripts folder from user domain as text) & "madever"
set madever to load script file (madever_dir & ":madever.scpt")
set loadtheme to load script file (madever_dir & ":loadtheme.scpt")
set CSS to loadtheme's main()
@naoyat
naoyat / kobitonote.py
Created November 23, 2012 09:38
kobitonote.py - Kobitoで編集したアイテムをEvernoteに同期保存するpythonスクリプト ref: http://qiita.com/items/4c465380b4020185a224
#!/usr/bin/env python
# -*- coding: utf-8 -*-
##
## kobitonote - KobitoのアイテムをEvernoteにも保存する
##
## (c)2012 by @naoya_t
##
import os
import re
import sqlite3
@naoyat
naoyat / fig37.matlab
Created November 28, 2012 03:49
PRML図3.7の再現コード ref: http://qiita.com/items/c31d850e516646f712d3
1;
function g = gaussian(x, mu,sigma)
g = 1/sqrt(2*pi*sigma^2) * exp(-1/(2*sigma^2)*(x-mu)^2);
endfunction
function g = gaussian2(x, mu,Sigma)
g = 1/(2*pi*sqrt(det(Sigma))) * exp(-1/2*(x-mu)'*inv(Sigma)*(x-mu));
endfunction
@naoyat
naoyat / file0.txt
Created November 29, 2012 10:21
TeX記法ちょっと練習 ref: http://qiita.com/items/29e1fddf566a80a78525
1;
%
% コーシー分布 Cauthy distribution のpdf
%
function p = CauthyDistributionPDF(y)
p = arrayfun(@(x) 1.0/(pi*(1+x^2)), y);
endfunction
%
% コーシー分布に従う乱数