Matlab 最小二乗解の求め方 ②

過剰条件の連立一次方程式

3X1 + 4X2 = 1000
1X1 + 7X2 = 1000
2X1 + 8X2 = 1000

（ X1 と X2  は製品A, Bの個数 ）

における最小二乗解を Matlab で解いたメモ

>> A = [3,4;1,7;2,8]

A =

     3     4
     1     7
     2     8

>> y = [1000,1200,1500]'

y =

        1000
        1200
        1500

>> x = (A'*A)^(-1)*A'*y

x =

  128.7435
  154.2169

>> x = inv(A'*A)*A'*y

x =

  128.7435
  154.2169

>> x = A\y

x =

  128.7435
  154.2169


実際に X1  と X2 を代入して確認すると


3 * 128 + 4 * 154 = 1000
1 * 128 + 7 * 154 = 1206
2 * 128 + 8 * 154 = 1488

となり、１２０６は１２００をオーバーしているのでだめ。

X1 = 128
X2 = 153

が経営者が取るべき現実的な解となる。

ーーーーー

ここで扱った問題は以下の論文から引用させて頂きました。

・ 過剰条件の連立一次方程式における最小二乗解, 近藤寿紀 
    「2003年度卒業論文要旨集」(2004.3.7 南山大学数理情報学部 数理科学科)
    http://www.seto.nanzan-u.ac.jp/ise/gr-thesis/ms/2003/torii/00mm043.pdf

昼ごはん

Phillbrook のダイニングホールにて。
UNHの学食はビュッフェ形式です。

今日のメニュー

・ サラダ（きゅうり、たまねぎ、トマト、塩・レモンで味付け）
・ じゃがいもの炒めたやつ
・ チキンナゲット
・ ゆで卵
・ ホットドッグ
・ オレンジジュース
・ フルーツ（オレンジ、パイナップル、メロン、レーズン）＆ヨーグルト

ゆで卵が一番うまい。次はじゃがいも。
 

これで ５．０３８ドル なのだから安いよね。
（１５０回で７５５．７０ドル）

自分の健康はこの学食によって保たれている。

Matlab 複数行（or 列）の削除

A の第i 列を削除した行列X を作る場合
X = A;
X(:,i) = []

同様に，A の第i 行を削除した行列X を作る場合
X = A;
X(i,:) = []

また，A の第i 列から第j 列までを削除した行列X を作る場合
X = A;
X(:,i:j) = []

－－－－－－－
以下の HP から引用させて頂きました

島根大学 総合理工学部 電子制御システム工学科
制御工学研究室　吉田和信　Matlab 教材
http://www.ecs.shimane-u.ac.jp/~kyoshida/matlab%282001%29.pdf

Matlab 関数の掛け算

超基礎的なことだが、今日はじめて知る。

y = exp(-2*t).*sin(t)

*（かける） の前に . （ドット）をつけること

家政婦のミタ

の５話から最後までを一気に見た。

ボロ泣きしてしまった。

どういう表現が好ましいか？

１．事実であること
２．シンプルであること
３．直感的であること
４．美しいこと

Googleのプライバシーポリシーに関するメールを参照。

殺傷事件発生

土曜日の今日、CCOMに来てメールをチェックするとこんなメールが。

－－－－

Sat 1/28/2012 2:17 AM

STABBING OCCURRED AT ADAMS TOWER WEST.   POLICE ARE ON SCENE AND INVESTIGATING.

Sent by Erin Storace to All users (E-mail, Pagers, Cell phones) through UNH Alert

... powered by the *** Alert Network

ーーーー

STABBING　という単語が分からないので辞書で調べると

「突き刺すこと、殺傷」

 とある。 アダムズタワーというとこで殺傷事件が起きたらしい。

まじで！怖え！

メールが送られてきた時間から深夜。

  その後も、続報がいっぱい来ている。

ーーーー

Sat 1/28/2012 2:52 AM

UNH POLICE ARE LOOKING FOR 2 WHITE COLLEGE AGED MALES BOTH 6 FEET TALL. ONE WEARING A RED FLANNEL JACKET WITH BROWN HAIR, THE OTHER WITH A BLACK SHIRT AND BLONDE HAIR. DURHAM AND UNH POLICE DETECTIVES ARE NOW ON SCENE. ANY INFORMATION ON THESE SUBJECTS PLEASE CALL UNH POLICE AT 603-862-1427.

ーーーー

容疑者は２人組で、年をとった(AGED)、白人(WHITE)で、１８０ｃｍ(6 feet) らしい。

どういうジャケットを着ているとか、詳しい情報が書かれている。

結局、その後、無事に解決されたらしい。

デーブに聞くと、殺傷事件は９年間で初めてとのこと。

UNH ではメールアドレスを登録するとこのような危険情報を知らせてくれる。

吹雪のような自然災害も。なかなか便利。

日本の大学ではこういうサービスあるのかな？