スポンサーサイト

上記の広告は1ヶ月以上更新のないブログに表示されています。
新しい記事を書く事で広告が消せます。

2度目のZipf-Mandelbrot分布

今日は休んで、午前は反応拡散パターン探索、午後からは、あるデータの解析を続けていた。

$$

修正ZM分布版(完成。稀なデータでズレ)
setwd("F:/Seimei")
data <- read.table("seimei-08Feb2014.csv", header=T, sep=",")
data[,2] <- as.numeric(gsub("位","",data[,2]))
x <- c(1:length(data[,1]))
y <- data[,3]/sum(data[,3])
resid <- function(par){
xx <- c(1:1000000)
C <- sum((par[3]^x)/(xx+par[2])^par[1])
yhat <- (par[3]^x)/((x+par[2])^par[1])/C
sum((y-yhat)^2) # 残差平方和を返す
}
result <- optim(c(1,1,1),resid,method="Nelder-Mead")
C <- sum((result$par[3]^x)/(x+result$par[2])^result$par[1])
plot(x, (result$par[3]^x)/((x+result$par[2])^result$par[1])/C)
points(x, y, col="red")
ey <- (result$par[3]^x)/((x+result$par[2])^result$par[1])/C
plot(log(x), log(ey))
points(log(x), log(y), col="red")

> pdf("seimei-distribution_26Feb2014.pdf", height=7, width=10, pointsize=14)
> plot(log(x), log(ey))
> points(log(x), log(y), col="red")
> dev.off()

seimei-distribution_26Feb2014.png


岸辺露伴は動かない (ジャンプコミックス)岸辺露伴は動かない (ジャンプコミックス)
(2013/11/19)
荒木 飛呂彦

スポンサーサイト

PCLによるMax Planck点群データの表面復元

Poisson Surface ReconstructionアルゴリズムがPCLに入っている。これを使いMax Planck点群データからメッシュを構築する。

元の点群データ:

Max_Planck.png

点群データは下記サイトからダウンロードし使用している。
http://www.cs.unc.edu/~ravishm/robustPointVisibility_smi_10.html

PoissonSurfaceReconstruction_22Feb2014_2.cppとして保存した。

Max_Planck_mesh.png


ミミイカの動画

YouTubeで探した。下は、もの凄くよく撮れている動画。
最後の砂に潜る姿がかわいい。




イカ・タコガイドブックイカ・タコガイドブック
(2002/04)
土屋 光太郎、阿部 秀樹 他

for文によるインターネットからのファイルの繰り返しダウンロード

import urllib.request
for i in range(1,481):
  u = '%s%03d%s%03d%s' % ('http://www.abcdefg/profile/', i, '.jpg/$File/', i, '.jpg')
  v = '%s%03d%s' % ('C:\\Users\\hijklmn\\Pictures\\download\\', i, '.jpg')
  urllib.request.urlretrieve(u,v)

パブ 夢精

この記事を閲覧するにはパスワードが必要です
パスワード入力

アンドロイド

この記事を閲覧するにはパスワードが必要です
パスワード入力

Visual StudioでQ_OBJECTマクロを宣言する方法

この記事を閲覧するにはパスワードが必要です
パスワード入力
Profile

zoa

Author:zoa
Metropolis (1927)

Calendar
01 | 2014/02 | 03
- - - - - - 1
2 3 4 5 6 7 8
9 10 11 12 13 14 15
16 17 18 19 20 21 22
23 24 25 26 27 28 -
Labels
Admin
Previous Posts
Recent Comments
Recent Trackback
Archive
このページのトップへ
上記広告は1ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。