はてなブックマーク

（2020年8月15日：Julia 1.5で、REPLではlocalスコープ内からglobal変数を参照できるようになりました。Julia 0.6以前の挙動に戻ったことになります） Juliaの変数スコープはMatlabやC言語と少し違うので混乱した話。（2019年8月25日：global/localスコープの説明が公式ドキュメントと違ったので修正しました） たとえば、1から10までの整数の和を求めたいとき、C言語では以下のようにやるかと思います。 #include <stdio.h> int main(int argc, char *argv[]) { int total = 0; int x; for (x=1; x<=10; x++) { total += x; } printf("Total: %d\n", total); return(0); } totalをfor文の内側から更

Excelでもできる程度の簡単な統計的仮説検定ならSciPyの範囲内でできるもよう。以下の実行環境はMac OS X 10.11.6 El Capitan上でHomebrewのpyenvからインストールしたAnaconda 4.3.1（Python 3.6.0、NumPy 1.11.3、Scipy 0.18.1）です。 下記の例は、いずれも伝統的な（あるいは頻度論に基づく）統計的仮説検定をPythonでやる方法をざっくり書いたものですので、統計については別途教科書などを参照して下さい……。 あらかじめ以下のようにscipy.statsモジュールを読み込んでおきます。 >>> from scipy import stats 2項検定 コイン投げのように成功率1/2（=0.5）のもとで、20回中16回表を出すことに成功したとき。 >>> p = stats.binom_test(16, 20,

さきほど気づいたのですが、ターミナルからSwiftを実行するときにxcrun swiftをしなくてもよくなっていました。僕の環境は「Xcode 6.2 (6C131e)」ですが、/usr/bin/swiftと/usr/bin/swiftcという二つの実行ファイルが用意されています（Yosemiteからなのかもしれません）。これを使えば、これまでより楽にSwiftでシェルスクリプトを書くことができそうですし、ターミナルから実行ファイルも作れそうです。早速試してみました。 シェルスクリプトっぽく使ってみる 次のようなhello.swiftというファイルを作成しました。 #!/usr/bin/env swift let name = "丸井綜研" println("ようこそ" + name + "へ！") これを実行すると $ swift hello.swift ようこそ丸井綜研へ！ となります

PMML 僕が学部生のときに何から何までお世話になった会津大学の西村憲先生が開発されたPMML (Practical Music Macro Language)という言語がありました。PMMLはテキストベースで楽譜を書いていくコンピュータ言語です。PMMLは単に音符をテキストに置き換えただけではなく制御構造やマクロなどのプログラミング言語の要素を多く取り入れられていたため、楽譜上での繰り返し記号であったり、アルペジエータであったりの記述がとても簡単に行えました。しかも短時間での音楽制作が可能になるだけでなく、プログラミング言語としても利用が可能だったのが大きな特徴でした。僕も標準MIDIファイルから（VRMLを使って）3Dアニメーションを自動生成するプログラムをPMMLで作成し、音高やヴェロシティといった音楽要素をアニメ内の任意の要素（色や位置や形状）に変換することについての卒業研究を行い

IncanterはLISP+R on JVMという理想環境 Incanterは「Clojureで動くRっぽい雰囲気の統計環境」とのことで、けっこう期待しているんですが、このところ表立った動きが見えてこないので寂しいです。JVM上でLISPを使って統計ができるのだから最強だと思うんだけどなぁ。（何ができるのか概要を知るためには『統計解析アプリ「Incanter」入門 - あんちべ！』や『Web系女子がLispと出会って統計学に目覚めるまでのお話 - あんちべ！』が良いと思います） で、本家サイトのトップページで「最新ビルドをダウンロード」みたいなところからダウンロードするとバージョン1.2.xのバイナリパッケージが入手できます。ただ、GitHubで行われている開発版を見てみるとバージョンは1.5.xで、すでに1.2の頃から大分と変更がなされているようです。最新版をできるだけ簡単に使うにはHo

Mac OS X 10.7 Lionでは、これまでのマウスホイールとは上下逆方向の挙動をします。これまではスクロールバーを上下に動かすのと同様の方向でしたが、Lionからはそれが逆です。ウェブページを閲覧しているときに、ウィンドウからはみ出て下に隠れている部分を見たければ、スクロールホイールを上向きに転がすことになります。 これは、ウェブページが巻物のように紙として印刷されていると考えると納得のいく操作です。もっと下のほうを読みたければ、紙をつまんで上にスライドさせればいいのです。すでにiPhoneやiPadではおなじみのスクロール方向だったので、徐々にiOSとOS Xを近づけようとするAppleにとっては一つのステップだったのでしょう。 ただ、そうなるとWindowsなどとの操作の不一致が気になります。可能なことならMac以外のマウスホイールも上下逆にしてしまいたくなります。 そこで少し

僕は言語を作るきっかけとか思想とかに興味があるので、日頃から新言語を見つけてはニヤニヤしてるんですが、つい昨日Juliaという新言語のリリース記事を読んで、面白そうだったので、紹介がてら粗く翻訳してみました。 なぜ僕らはJuliaを作ったか（原文：Why We Created Julia） 2012年2月14日(火) | Viral Shah, Jeff Bezanson, Stefan Karpinski, Alan Edelman 端的に言えば、僕らは欲張りだからだ。 僕らはMatlabのパワーユーザーだ。LispハッカーやPython使いやRuby使いもPerlハッカーもいる。髭が生える前からMathematicaを使っていたのもいるし、未だに髭が生えてない仲間もいる。常識的な人にはオススメしないくらい多くのグラフをR言語で描いてきた。そしてC言語は僕らのユートピアだ。 いま挙げた言

若手の音響研究者たちが作成しているOpenAcousticsというドキュメント＋ソースコード群があります。このサイトでは音響数値解析のソースコードだけでなく、充実した日本語の入門書が無料公開されています。だいぶ前に教えてもらっていたのですが、重い腰を上げて、空き時間に読み進めることにしました。 チュートリアルに使われているプログラミング言語はPythonです。Pythonはコードが理解しやすいので好きなのですが、より使い慣れているMatlabを使って自分用に書き直してみました。オリジナル版はGmshを使ってメッシュデータを読み込んでくるなど、より便利に使えるようになっています。下のコードではその部分は省略し、とりあえず動かしています。今度はドキュメントの数式を理解すべく、一行一行読んでいこうと思います。早くBEM、FDTD、CIPと進みたい。 （右のグラフは波の位置・周波数を変えて作ったも

CocoaのFoundationフレームワークで使うNSLog()のメモ。 printfとほとんど同じ使い方。 NSLog(@"hogehoge %d %f %@", intNum, floatNum, anObject); 行頭に日付時刻とプロセス名が、行末に\nが、自動的に付加されます。%〜の意味は以下の表の通り。(Hillegass本第3版より) 記号 表示する型 %@ id %d, %D, %i long %u, %U unsigned long %hi short %hu unsigned short %qi long long %qu unsigned long long %x, %X unsigned long (16進数で表示) %o, %O unsigned long (8進数で表示) %f, %e, %E, %g, %G double %c unsigned char

コンピュータの中にあるWAVファイルどうしの音量比較をしたいときはけっこうあります。音量の指標として「ラウドネス（sone）」があり、以前その計算プログラムも書いたのですが、その音を鳴らしたときの音圧値（dBSPL値）が分からないと正確な値が得られないというものでした。そのため、実際にスピーカーで鳴らして、レベルメーターで計測して、という手間がかかります。複数のWAVファイルの音量をある程度そろえたいときに、いちいち実験室をセッティングするのはめんどうです。どうせWAVファイルはすべてコンピュータ内にあるので、コンピュータ内で完結させられればいいなぁ、なんて思ってました。 そんなときに見つけたのがITU-R BS.1770-1勧告。複数ファイルの相対的な音量を計算するための方法が書かれた国際基準です。単位はLKFSで、0 dBFS、1 kHzのサイン波が約-3.01 LKFSとなります。規

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2023-04-01 14095.9 2.99 156.86 40.80 3.84 2023-07-13 14218.0 2.73 163.00 44.73 3.64 2023-09-25 14380.5 3.60 176.94 45.14 3.92 2023-11-25 14529.5 2.93 162.12 50.85 3.19 2024-02-07 14644.2 2.80 165.00 40.96 4.03 2024-03-04 14775.6 2.93 164.85 44.85 3.68 2024-05-02 14914.0 2.82 165.96 49.08 3.38 2024-06-02 14998.1 2.01 165.17 41.89 3.94 2024-09-2