matlab: なんちゃって構造体

空の構造体は

>> a = struct([])
a =
0x0 struct array with no fields.

とすれば簡単に作れる．では，「フィールドを持たないけど空ではない構造体」はどうすれば作れるか？

すぐに思いつくのは，rmfield を使って，フィールドを消して行くやり方だろう．

>> b.foo = []; a = rmfield(b, 'foo')
a =
1x1 struct array with no fields.

これはかっこ悪いので却下．だが実現は出来ることがわかる．
では，rmfield はどうやってこれを実現しているか？
rmfield は実は組み込み関数でなく .m 関数なので，type rmfield して中身をのぞくと，どうやら cell2struct を呼び出しているらしい．
ならばはじめから cell2struct を使えばいいのである．しかも cell2struct は組み込み関数だから，わりと根源的なやり方といえるだろう．

>> a = cell2struct(cell(0,1),{})
a =
1x1 struct array with no fields.

おお，できた．これでいいや．

しかし，恒例の whos でみると

>> whos a
  Name      Size                    Bytes  Class
  a         1x1                         0  struct array
Grand total is 0 elements using 0 bytes

のように，メモリが割り当てられていないではないか．しかも a の repmat で b を作ってみても，b も Bytes は 0 のままである．

実際にメモリを割り当てるにはフィールドを定義しないといけないらしい．AMD64 な Linux 上で AMD64 な MatlabR2006a を使った場合，

>> b = repmat(a,[100,1]);
>> whos b
  Name      Size                    Bytes  Class
  b       100x1                         0  struct array
Grand total is 0 elements using 0 bytes
>> b(1).foo = '';
>> whos b
  Name      Size                    Bytes  Class
  b       100x1                       952  struct array
Grand total is 100 elements using 952 bytes

となるが，同じことを IA32 な Windows 上で IA32 な Matlab6.5.1 を使って行った場合，

>> b=repmat(a,[100,1]);
>> whos b
  Name      Size                   Bytes  Class
  b       100x1                        0  struct array
Grand total is 0 elements using 0 bytes
>> b(1).foo = '';
>> whos b
  Name      Size                   Bytes  Class
  b       100x1                      520  struct array
Grand total is 100 elements using 520 bytes

となる．なんで違うんでしょう．

oblivion: なんとなく screenshot: 1

Sophia ※がある人を Cloud Ruler Temple に連れて行くと，集合した Blades の前で演説が行われた．

※Sophia = プレイヤー

matlab: なんちゃってcomplex

matlab (今回は 6.5.1 でしか試してません)は複素数について直感に反する挙動を示すことがある．例を見てみよう．

>> A = repmat(complex(0),[10,1]);
>> whos A
  Name      Size                   Bytes  Class
  A        10x1                      160  double array (complex)
Grand total is 10 elements using 160 bytes
>> A(end) = complex(1,0);
>> whos A
  Name      Size                   Bytes  Class
  A        10x1                       80  double array
Grand total is 10 elements using 80 bytes

complex じゃなくなってる．およ？一応確認すると，A に代入した値は実数でなく複素数と認識されているはずだ．

>> clear a; a=complex(1,0); whos a
  Name      Size                   Bytes  Class
  a         1x1                       16  double array (complex)
Grand total is 1 element using 16 bytes

となるからだ．ただし，入力引数に対する inputname のように，出力引数についての情報取得機構が関数 complex 内で働いているとすると，よくわからない話になるかもしれない．

一つ考えられる仮説は，「虚部が全て 0 の配列は実数型にされてしまう」というものだ．実際，その方がメモリ効率やら計算効率やらよくなる気はする．

念のため，虚部が 0 でない複素数を代入すれば，配列は実数型に変換されることは無い．

>> A = repmat(complex(0),[10,1]);
>> A(end)=complex(0,1);
>> whos A
  Name      Size                   Bytes  Class
  A        10x1                      160  double array (complex)
Grand total is 10 elements using 160 bytes

そして，要素を虚部が 0 の複素数に順に置き換えていくと，

>> A=complex([1,2],[3,4]);
>> A(1)=complex(real(A(1))); whos A
  Name      Size                   Bytes  Class
  A         1x2                       32  double array (complex)
Grand total is 2 elements using 32 bytes
>> A(2)=complex(real(A(2))); whos A
  Name      Size                   Bytes  Class
  A         1x2                       16  double array
Grand total is 2 elements using 16 bytes

となった．ここまでは前述の仮説で当たっていそうである．

しかし，もちろん話はこれで終わらない．次の例を見てみよう．

>> A=zeros(1,2);
>> A(:)=complex(0);
>> whos A
  Name      Size                   Bytes  Class
  A         1x2                       32  double array (complex)
Grand total is 2 elements using 32 bytes

なんと配列は複素数のままであり，これは前述の仮説と矛盾する．ところがさらに，代入の方法を変えると，

>> A=zeros(1,2);
>> A(1:end)=complex(0);
>> whos A
  Name      Size                   Bytes  Class
  A         1x2                       16  double array
Grand total is 2 elements using 16 bytes

となって前述の仮説が復活する．どうなってるのという話である．

ここで新たに，｢全代入(whole substition とか？)｣という概念を考えてみる．
全代入と考えるのは次の４つ．

1. A = B;
2. A(:) = B;
3. A = B(:);
4. A(:) = B(:);
   

変数間の代入で虚部 0 でも complex を維持するには，どうやら厳密にこの４形式しか許されないようだ．
たとえば２次元配列で A(:,:) としてもダメだったし，列ベクトルで A(:,1) としてもダメだった．
この概念を利用して立てた仮説がこれだ：
「全代入では愚直な代入が行われ，それ以外では値の取り出し・代入いずれかで可能なら虚部が 0 の複素数は実数に落ちる」
いまのところこの仮説で説明できるような気がする．
この他に，左辺が全参照（全代入と同じ気分）で右辺が「適当な」式なら，やはり実数に落ちないようだ．以下の例を参照のこと．

>> A=complex(rand(3));
>> A(:)=complex(0); whos A
  Name      Size                   Bytes  Class
  A         3x3                      144  double array (complex)
Grand total is 9 elements using 144 bytes
>> A(:)=complex(zeros(3)); whos A
  Name      Size                   Bytes  Class
  A         3x3                      144  double array (complex)
Grand total is 9 elements using 144 bytes

いや～．重箱の隅つつくといろいろ出てくるなぁ．
ちなみに，大方の予想通り

>> A=complex(rand(3));
>> A(:,:)=complex(0); whos A

  Name      Size                   Bytes  Class
  A         3x3                       72  double array
Grand total is 9 elements using 72 bytes
>> A=complex(rand(3));
>> A(:)=complex(zeros(3)) + 1; whos A
  Name      Size                   Bytes  Class
  A         3x3                       72  double array
Grand total is 9 elements using 72 bytes

となる．もういいって？
今回の現象は，subsref, subsasgn, subsindex と型キャストの実装のコンボで生じるんだろうけど，意図された仕様なんですかねぇ．

MM: 携帯動画の変換（mp4へ）

以前の記事にも書いた携帯で録画した動画を何とかするのの続き．
YAMB は MP4Box のフロントエンドであるらしく，Options->General->MP4Box->View or Edit MP4Box Command Line. にチェックを入れると，どうやって MP4Box を呼び出しているかがわかる．
これを調べると，次のようなバッチファイルでいけそうだ．

"C:\Program Files\YAMB\MP4Box.exe" -avi 1 "%~f1"
"C:\Program Files\YAMB\MP4Box.exe" -raw 2 "%~f1"
"C:\Program Files\YAMB\MP4Box.exe" -add "%~d1%~p1%~n1_track1.avi" -add "%~d1%~p1%~n1_track2.amr" "%~f1.mp4"
del "%~d1%~p1%~n1_track1.avi"
del "%~d1%~p1%~n1_track2.amr"

この内容でバッチファイルを作り，そのショートカットを Send To フォルダに入れておく．
すると，hogehoge.3g2 を右クリックし，「送る」で hogehoge.3g2.mp4 に変換できるようになる．
このバッチファイルは，w43h で録画すると出来る .3g2 を前提にしており，track2 が amr 形式の音声トラックだと仮定している．

バッチパラメータ（ %~f1 等）についてはykr さんのページや，Windows 一般に共通なのかわからないけどMicrosoft TechNetのページなどを参照．

ちなみに，track 1 の抽出を -raw 1 とすると .cmp ファイルができる．これでも .mp4 へ mux できたが，なぜか後ろの方が再生できなかったりした．.avi 経由だと問題なかった．なんで？？

さて，あとは mpeg2 に変換して DVD に焼くところだ．

matlab: 呼び出しの調査

matlab 上で関数を呼び出すと疑似コンパイル(p-code)され，メモリに残る．
メモリ上の p-code buffer に存在する m-file の名前を見るには，inmem という関数を使う．
これで，実際に呼ばれた関数の一覧が得られる．

関数のコードを走査して再帰的に依存を調べるには，depfun という関数を使う．
depfun の賢さは限定的なので，もれなくわかると期待してはいけない．depdir も同様．

profile なんかも結構役に立つ．単に echo してみるのも面白いかもしれない．

MM: 携帯で録画した動画

MM はまるちめでぃあのつもり．

携帯（W43HにmicroSD 2GB）で録画した動画（3gpp2 ってやつ？）を，USB で PC にもってきた．
PRIVATE\AU\DF\D_MA に入っている .3g2 な動画を再生してみると，
結構長時間連続録画したはずなのに，mpc でも vlc でも，再生時間が15秒と認識されてしまう．
いわゆる断片化ということに関連した問題らしい．

いじれそうなツールとしてYAMBを使ってみた．
YAMB からは長時間の動画として認識されたので，
動画を avi に，音声を amr に抽出して，mp4 に remux.
無事に全時間分再生できた．

もちっとワンタッチなやり方があるといいんだけど．
MP4Boxを直に使ってバッチ処理とか出来るんだろうか．

ゲーム: Oblivion

Elder Scrolls IV: Oblivion というゲームがあります．

RPGです．洋ゲーです．
4Gamer.net のレビュー，開発元，販売元

いくつかのプラットフォームで出ています．
今度日本語版も出るそうです．

日本語の wiki は Oblivion Wiki JP というのです．
プレイ日記も結構たくさんあるのです．
Simらしのなく頃に+Oblivion日記
あたりが面白いのです．

ですです．です？です‥‥ですっ！

ということで，通販で買ってしまいました．PC 英語版．
PC なのは MOD を入れたいから．英語版なのは趣味で．
このために6年ぶりに PC を刷新したといってもイイ（いずれ書くかも）．
メモリも積んだしCPUも変えた！
キーボードまで diNovo Edge（engadget japanese, スタパブログ, マイコミジャーナルのレビュー）にしちまった（いずれ書くかも）！

そして Oblivion をプレイしてみた．

やばいです．楽しいです．

下手すると今後の記事がひたすらプレイ日記になってしまうかもしれない．
どうしよう．