Juliaで任意精度演算

Juliaで階乗を計算する関数は、その名もfactorial()です。たとえば20!を計算するには以下のようにすれば一発です。 julia> factorial(20) 2432902008176640000 ただ、デフォルトでは64ビット符合付き整数での演算なので、その範囲を超える場合にはエラーにな…

ピンクノイズを作るPython関数

Voss法でピンクノイズを作るプログラムをMatlabからPythonに移植しました。ベクトルではなくループで書いているところなんて、Python的・numpy的にも汚い書き方だと思うので、適当に修正して使って下さい。ピンクノイズについては「DSP Generation of Pink N…

Kindleの進化

僕が持っているKindleは、アメリカ国外にも初めて販売開始したKindle 2 Internationalというもの。100カ国ぐらいで3G回線を使って電子書籍のダウンロードができるモデル。当時はiPadが出る前だったので「おおタブレットですか、かっこいい」などと言われたも…

昔のウェブブラウザで

昔のウェブブラウザでURL欄に「about:version」と入力すると、そのブラウザのバージョン情報が表示されるとかいうのがあった気がするのだけど、Safariでそれをやっても検索結果が出てしまう。「about:blank」は動いたが、「about:plugins」とかもあった気が…

Swiftを試してみた(簡単な分布関数)

昨日の続きで(チラシの裏のように)書いてみています。今日は二項分布と正規分布。あとは、累積分布関数の逆関数(norminvとbinoinv)に加えてt分布とF分布があると日常的な統計処理には使えそう。(たぶん仕事はSwiftでなくRでやると思うけど) 関数名はMa…

Swiftを試してみた(簡単な記述統計関数)

新しく使う言語を試すときには音響プログラムを書いてみたり、統計プログラムを書いてみたりします。Pythonが簡単で面白かったので、Swiftはどんなだろうと試してみます。今回はSwiftで簡単な統計関数を作ってみました。と言っても平均や分散程度ですし、エ…

Twitterのような近況報告

1時間だけの数学の講義資料を作るのに10時間くらいかかった。効率が低いままだから、いつまでたっても週末も仕事漬けなのだ。とは言っても昨日は学会会場までバイク乗ったし、今日は10kmほどランニングしたし、仕事だけではないが。これからOblivionを30分ほ…

本日の給油 & 南東北の旅

タイトルに「本日の給油」としてあるけど、ゴールデンウィーク中に行った南東北のカブ旅の給油記録。宿泊込みの単独ツーリングはお伊勢参りのときに続いて4年ぶり2度目。 給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 201…

碁盤の大きさ(19路・13路・9路)

碁盤の大きさとして主なものに何種類かあります。プロの対局にも使われている一般的な大きさは19路盤(じゅうきゅうろばん)ですが、入門用として13路盤と9路盤もよく使われています。他にも6路盤、7路盤、15路盤などもあると聞いたことがありますが、特によ…

Python(とSciPy)で簡単な統計的仮説検定

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再開してみています

15年くらい前には中途半端にPythonを使っていたのですが、最近のデータサイエンス業界でのPythonの盛り上がりが気になったので、Pythonに再入門してみました。(15年前は主としてCとJavaをやっていたのでCPythonではなくJythonを使っていましたが、いろいろ…

本日の給油

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2017-03-18 8559.6 2.76 125.72 45.94 2.74 今日は埼玉県南埼玉郡まで。前回給油が去年10月29日に茨城県つくば市に行ったときだったので、それからは近所への買い物くらいし…

本日の給油

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2016-02-13 6714.7 2.96 108.78 49.26 2.21 2016-04-03 6847.3 2.57 108.95 51.60 2.11 2016-04-10 7025.7 3.42 111.99 52.16 2.15 2016-05-18 7113.8 1.68 108.93 52.44 2.0…

48ヶ月点検

スーパーカブ110プロ(JA07)の4年目の定期点検。春頃から、以前と音が変わったかなーと思ったまま乗り続けていたところ、それがチェーンにけっこうなダメージを与えていたらしいとのこと。今年は調整で何とかしてもらい、来年の点検時にはスプロケット込み…

本日の給油

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2016-02-13 6714.7 2.96 108.78 49.26 2.21 2016-04-03 6847.3 2.57 108.95 51.60 2.11 2016-04-10 7025.7 3.42 111.99 52.16 2.15 2016-05-18 7113.8 1.68 108.93 52.44 2.0…

MaxのGenExprを使ってみた

先日Pure Dataのexternを書きました。 Pure Dataのexternalを作ってみた(前編) - 丸井綜研 Pure Dataのexternalを作ってみた(後編) - 丸井綜研 これを、Maxでも動かせるようにしようと思って、Max SDKのドキュメントを読み始めました。ただMaxにはgen~も…

今後10年の言語(翻訳)

もともとはClojureのことが読みたくてBruce A. Tateの「Seven Languages in Seven Weeks(訳書)」を買ったのですが、読んでいくとErlangが気になりはじめました(ほぼ同時期に発掘して読んだMatz氏のブログ記事もErlangに言及していたので)。その後、2013…

Pure Dataのexternalを作ってみた(後編)

はじめに 先日つくった正弦波を生成するsingen~を拡張して、もう少し複雑なexternにしていきます。正弦波の生成ができたので、それらを複数重ねて音作りをする加算合成シンセにしてみましょう。 とは言ってもいきなり難しいものは無理なので、下図のように基…

Pure Dataのexternalを作ってみた(前編)

はじめに MacでPure Dataのオーディオexternal objectを作ってみたメモ。日本語で書かれた情報が少なかったので、「無いよりはマシ」程度で書き残しておきます。いろいろと間違いがあると思います。指摘していただければ幸いです。 まず参考にしたのはPd-ext…

Pure DataのexternalをMacでビルドするシェルスクリプト

MacでPure Dataの拡張オブジェクト(external object)をビルドする方法はいくつかあり、下記のMakefileを使うものなどが便利だったりします。(永野さん提案の、Xcodeを使う方法もあります) 参考:Pure Data の拡張オブジェクトを作る gcc on Mac OS X 版 …

本日の給油

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2016-02-13 6714.7 2.96 108.78 49.26 2.21 2016-04-03 6847.3 2.57 108.95 51.60 2.11 2016-04-10 7025.7 3.42 111.99 52.16 2.15 2016-05-18 7113.8 1.68 108.93 52.44 2.0…

MFiゲームコントローラはほぼ全てがアナログ入力

iPhoneに接続できるMFiゲームコントローラを探していて、Horipad Ultimateが良さそうだったので購入。Macに接続して調べてみたら、メニューボタン以外の全ボタンがアナログ入力に対応していた。十字キーまでとは驚き。 MFi要件がそうなっているのだろうか。…

ANSI Common Lispの規格書(に近いもの)のPDFを入手する【改】

『ANSI Common Lispの規格書(に近いもの)のPDFを入手する(http://g000001.cddddr.org/3676298293)』を自分の環境で試してみました(元記事に感謝!)。OS X 10.11.5、MacTeX 2016むけに、いくつか変更。ファイルのダウンロードにwgetではなくcurlを使用。P…

本日の給油

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2016-02-13 6714.7 2.96 108.78 49.26 2.21 2016-04-03 6847.3 2.57 108.95 51.60 2.11 2016-04-10 7025.7 3.42 111.99 52.16 2.15 2016-05-18 7113.8 1.68 108.93 52.44 2.0…

Arai CT-Zを購入

いま使っているバイク用ヘルメットはShoeiのMultitecです。いわゆるシステムヘルメットというもので、チンガード全体がシールドのように跳ね上げられる構造になっています。信号待ちで水分補給できるなど便利です。しかし、これは2007年にバイクの新車購入と…

本日の給油

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2016-02-13 6714.7 2.96 108.78 49.26 2.21 2016-04-03 6847.3 2.57 108.95 51.60 2.11 2016-04-10 7025.7 3.42 111.99 52.16 2.15 2016-05-18 7113.8 1.68 108.93 52.44 2.0…

本日の給油

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2016-02-13 6714.7 2.96 108.78 49.26 2.21 2016-04-03 6847.3 2.57 108.95 51.60 2.11 2016-04-10 7025.7 3.42 111.99 52.16 2.15 2016-05-18 7113.8 1.68 108.93 52.44 2.08

本日の給油

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2016-02-13 6714.7 2.96 108.78 49.26 2.21 2016-04-03 6847.3 2.57 108.95 51.60 2.11 2016-04-10 7025.7 3.42 111.99 52.16 2.15

本日の給油

給油日 オドメーター (km) 給油量 (L) 単価 (円/L) 燃費 (km/L) 距離単価 (円/km) 2016-02-13 6714.7 2.96 108.78 49.26 2.21 2016-04-03 6847.3 2.57 108.95 51.60 2.11

親指シフト始めました(NICOLA D型配列?)

昨年5月にDVORAK & JISカナ入力を始めたのですが、ようやく最近になって会議のメモが普通に取れるくらいの速度になってきました。ただし、小指を酷使するのと、指が動く範囲が広いことで、なかなか疲れるものがあります。そこで最近、さらに疲れない入力方法…