% 一橋大学「金融経済論I」（2010年度、担当：塩路)
%
% 2010年10月19日
% ファイル名＝kiso2.m
%
% 目的：Matlabに親しむ

% (練習2) 行列の足し算、逆行列など。

disp('行列演算の練習１') 
x=[1 2; 3 4]   %2かける2行列の入力
y=x'           %xの転置行列
z=inv(x)       %xの逆行列
w=x*z          %行列の掛け算
v=x+y          %行列の足し算
u=[x x']       %2つの行列を横に結合
uu=[x ; x']    %2つの行列を縦に結合
size(x)        %行列の行と列の数
det(x)         %行列式
eig(x)         %固有値

disp('行列演算の練習２') 
x=[1 2 3; 4 5 6; 7 8 9] 
xx1=x(3,2)     %行列の要素の取り出し(第3行第2列を取り出せ)
xx2=x(1:2,2:3) %部分行列の取り出し(第1行から第2行、第2列から第3列を取り出せ)
xx3=x(1:2,:)   %部分行列の取り出し(第1行から第2行、すべての列を取り出せ)

disp('行列演算の練習３：固有値分解') 
[P,L]=eig(x)   %Lの各対角要素がxの固有値、Pの各列がそれに対応する固有ベクトル
x2=P*L*inv(P)  %本当に固有値分解できているのか？Matlabを疑うわけではないが、念のため確認
x-x2           %しつこく確認
            
% このファイルの使い方
% このファイルを、Matlabの入っているパソコンの適当なディレクトリにコピーする。
% Matlabを起動する。
% Current Directoryを、このファイルが入っているディレクトリに変更する。
% M-File Editorでこのファイルを読み込む。たとえば、「ファイル(F)→開く(O)→kiso2.mをダブルクリック」
% このファイルを走らせる。たとえば、カーソルをこのファイルの上において、F5キーを押す。
% 結果はCommand Windowに現れるはず。