読者です 読者をやめる 読者になる 読者になる

室内音響指標C値の計算 (R Advent Calendar 2013)

あけましておめでとうございます。R Advent Calendar 2013の32日目の記事です。僕の専門は音響と心理との学際領域で、ふだんは音に関する主観評価実験の準備・分析にMatlabやRを使用しています。今回は室内音響分野で使われる計算をRでやってみます(その分…

シェッフェの一対比較法(中屋の変法)

先日の浦の変法に続き、シェッフェの一対比較法(中屋の変法)を実装してみました。佐藤信『統計的官能検査法』(日科技連出版社)に書かれている方法をそのままコードにして、今度はp値の計算結果も表示するようにしただけです。詳しいことは同書を参照して…

シェッフェの一対比較法(浦の変法)

自分の仕事でよく使うシェッフェの一対比較法(浦の変法)を、Rの勉強がてら実装してみました。佐藤信『統計的官能検査法』(日科技連出版社)に書かれている方法をそのままコードにしただけなので、詳しいことは同書を参照して下さい。使用例は以下の通り。…

数字を読み上げる音声ファイルを作る

主観評価実験をするときに、刺激を逐次呈示して、それに対する評価を記述してもらう方法があります。10個程度の刺激ならあまり問題にはならないのですが、刺激が多くなると、いま何番目の刺激なのか、いま使っている回答欄は正しいのか、など不安になること…

Excel 2011で95%信頼区間を計算する

普段はMatlabやRやSPSSで統計処理をしているのですが、Excelしか使えない場合もあるので、Excelでも簡単な統計処理はできるようにしています。そんな中で、何年か前に、信頼区間の計算でつまづいたことがあったので備忘メモです。Excel 2008にはCONFIDENCE関…

東京事変『林檎の唄』のクラビを聞く

東京事変のファーストアルバム『教育』の第1曲目にクラビネットが入っている『林檎の唄』があります。右方向から音が出ていますが、ヒイズミマサユ機がどんな演奏をしているのか詳しく聞きたかったので、クラビネットを強調するプログラムを書きました。 [x,…

Octaveによる主成分分析の計算

Matlabにはprincompという主成分分析を計算する関数があるけれど、Octaveには見あたらなかったので、英語版Wikipediaに書かれていたやり方どおりに作ってみました。 [pc, score, latent] = princomp(A); みたいに使えます。pcが主成分の係数、scoreが主成分…

情動の種類

メモ的に書いておきます。Shaver (1987)による階層構造になった情動の分類。基本5情動というのはLove、Joy、Anger、Sadness、Fearですね。 快 Love Fondness Infatuation Joy Bliss Contentment Pride 不快 Anger Annoyance Hostility Contempt Jealousy Sad…

AHP

昨日まで参加していた音楽音響研究会に、ファジーAHPというものを用いた研究が発表されていました。AHPとはAnalytic Hierarchy Processの略で、日本語だと階層分析法と言うらしいです。複数の対象から一つを選ぶという意志決定を行うときに、対象から一つを…

Matlabで2次元プロット

論文に載せる2次元プロットを作成するとき、僕はMatlabを使って以下のような感じで書いています。多次元尺度法の分析結果を図にするときに相性がいいですよ。背景になる部品から順番に重ねていく必要があるので、まずは縦軸と横軸を書きます。grid onでグリ…

vDSP関連サイト

Matlabでやっていた信号処理をC/C++/Objective-Cあたりでガリガリやってみようと思い立ちました。そこでAppleのAcceleration Frameworkを使ってみようと思ったのですが、Appleのサイトにはリファレンスとサンプルプログラムくらいしか頼りになるものがなくて…

sound qualityとsound characterの邦訳

Francis Rumsey ``Spatial Quality Evaluation for Reproduced Sound: Terminology, Meaning, and a Scene-Based Paradigm'' Journal of Audio Engineering Society, v.50, n.9, pp.651--666, Sep 2002を再読中。日本語で説明しようとして、sound qualityとs…

Mac OS X用の数値計算ライブラリ

Mac OS Xで簡単に使えそうなBLASやLAPACKみたいな数値計算ライブラリを探していたら、Accelerate.frameworkなんてのが出てきました。なんと、Mac OS XにデフォルトでBLASもLAPACKも入っているではないですか。しかもvImageやvDSPという画像処理・信号処理の…

プログラム修正

過去に書いたMatlabプログラムのうち、どうもうまく動いていないようだった部分をまとめて修正しました。最大のミスは、音ファイルの周波数特性を描き、それを1/3オクターブ幅で平均してくれるプログラムのはずが、1/3オクターブ刻みのそれぞれにおいて、バ…

MUSHRAをMax/MSPで実行する

MUSHRA (MUltiple Stimuli with Hidden Reference and Anchors)という実験手法があります。これは知覚符号化音源など、違いの小さな複数の音刺激の音質評価を行うために開発された手法です、複数(10個くらい)の刺激をボタンで切り替えながら試聴し、それぞれ…

音の周波数特性を描くMatlabプログラム

音の周波数特性を描くMatlabプログラムを改善したので、メモ的に置いておきます。今回の目玉は「1/3オクターブごとの平均値を描くようにした」です。500-2000Hz帯の平均値が0dBになるようにスケーリングしています。

R言語で分散分析(と多重比較)

例えば、A〜Eの5種類の肥料をそれぞれ別の畑で使ったときのジャガイモの収穫量を表す以下のようなデータがあったとき、肥料による収穫量の明確な違いがあるかどうかを調べようと思います。一元配置の分散分析(対応なし)を行います。(ジャガイモ収穫量の表…

MIT全教員の論文がオープンアクセス

MIT全教員の学術論文が無料でアクセスできるようになりました。自分が過去に書いた論文なんかは全部公開してしまいたいんですが、頒布権が学会に属しているものもあるので、なかなかうまくいきません。僕の分野の学会で言えば、Acoustical Society of Americ…

シェッフェの一対比較法(中屋変法)のMatlabコード

研究で必要になったので、佐藤信『統計的官能検査法』(日科技連 1985年)を参考に、シェッフェの一対比較法(中屋変法)による分散分析を行うMatlabスクリプトを作りました。実のところ「参考にした」どころではなく、書いてあった手計算のアルゴリズムをそ…

日本音響学会2009年春季研究発表会

今日の目当てはHRTFの特別セッション。千葉工業大学の飯田さんの発表が面白かったです。僕が会津大学にいた時に、修士生が取り組んでいた問題に似たものですが、基礎的な部分で解決しているので、とても応用範囲が広そうです。関連文献は以下(ElsevierのApp…

でたらめ選択がどのくらい有り得るものなのか(後編)

昨日は、「2つ聞こえた音は同じものですか?」という問いに対して「はい」か「いいえ」で答えさせるタイプの実験を行ったとして、2項分布を使った分析方法について紹介しました。これにより、例えばWAVとMP3を聞き比べて、それらが同じように聞こえるかどう…

でたらめ選択がどのくらい有り得るものなのか(前編)

音響心理実験で「2つ聞こえた音は同じものですか?」という問いに対して「はい」か「いいえ」で答えさせるタイプの実験があります。たとえば、CDからリッピングしただけのWAVファイルと、それをMP3化したものの音質が聞き分けられるかを調べたりする実験なん…

AESポッドキャスト

いつのまにかAESがポッドキャストをはじめていてびっくり。

ピッチとチューニングについての論文

授業でピッチやチューニングについて取り上げるときのために、忘れないうちにメモ。 大串健吾「音の高さに関する生理・心理学と楽器のチューニング」日本音響学会誌59巻3号(2003), pp.153-158 森太郎「ピアノ調律技術」日本音響学会誌60巻5号(2004), pp.…

SICPに相当するもの

id:higeponさんの「計算機プログラムの構造と解釈(SICP)」を読み終えてを読んで、僕の『計算機プログラムの構造と解釈』は本棚の肥やしになっているなぁ、と思い出したのでした。日本語版は翻訳が悪いだの何だの言われていたので、意気込んで英語版を買った…

ロジットとプロビット

ロジットしか使ったことなかったのですが、プロビットも同じ目的で使えるのですね。どちらも独立変数が連続値で従属変数がブール値のときに、シグモイド関数を重回帰分析ではめる分析手法です。SPSSのページによると、違いは「ロジスティック関数を使うか正…

アストロノーカ

山本弘『神は沈黙せず』という小説を読み、遺伝的アルゴリズムに興味を持ちました。遺伝的アルゴリズムなどの人工知能的な物には以前から興味があったのですが、特に勉強するでもなく過ごしてきました。僕が博士課程にいるときに、同期が研究していた遺伝的…

ワインの品質を計算する方程式

経済学者Orley Ashenfelterによる「ワインの品質を計算する方程式」というのがあって、それは ワインの質=12.465+0.00117×冬の降雨+0.0614×育成期平均気温−0.00386×収穫期降雨 なのだそうです。その年の気候データを入れて計算すると、ワインの品質がどの…

ダミーヘッドがあるものだから、つい調子に乗って

今日から学生は登校禁止になっていたのですね。すっかり忘れて研究や面談の用事を入れてしまいまっていました。借りていたダミーヘッドを使って集中的にスタジオやヘッドホンのインパルス応答もできた、充実した一日でした。でも、何時間も小部屋にこもって…

インパルス応答録音

朝から夕方まで集中的にコンサートホールのインパルス応答を測定し、合計で1400個以上の録音ファイルが出来ました。天井高を変化させたり、さまざまなマイクロホン位置で録音したりした、とても面白そうなデータです。これからみっちり分析しないといけませ…

実験見せてもらった

ヤマハに勤めている友人が来学してくれて、いま取りくんでいる研究を見せてもらいました。けっこうスタンダードなことをやってるんだなーと思っていたら、いろいろと細かいところで興味深い処理をやっています。プロトタイプはすべてMax/MSPで行っているそう…

有楽町マルイのアンケート

有楽町マルイで1週間にわたってBGMアンケートを行いました。ご協力いただいたお客様、ありがとうございます。6月からはBGMががらりと変わります。7月、8月と、お買い物のしやすさを向上してゆくために、アンケートの結果を分析していきます。目標は「対流時…

インパルス応答測定

授業の中で、吸音パネルがどのくらい吸音するのか測定してみました。とは言っても、吸音素材の厳密な測定はかなりの難問なので、「音がどのくらい透過するか」ということを中心に計りました。測定に使った信号はPink TSPです。TSP信号の時間引き延ばしを線形…

音響計測関連の文献

TSPなんかを使って音響計測をしていたんですが、もっと良い方法を探して、文献あさりをしていました。そんな中で見つけたのがAudio Precisionという計測機器メーカーが無料公開している、計測マニュアル。計測手法そのものの説明ではなく、より一般的なこと…