成分の計算

  • こちらのつづき
  • さらに付け加える条件
    • a_{i,j}a_{4-i,4-j}=0
  • Maxima で各マスへの移動を計算してみる
x:0.1;
y:0.2;
nl:0.3;
nt:0.4;
solve([-a11-a21-a31+a13+a33+a23=x,a11+a12+a13-a31-a32-a33=y, (a11-x)^2+(a21-x)^2+(a31-x)^2+(a13-x)^2+(a23-x)^2+(a33-x)^2=nl,(a11-y)^2+(a12-y)^2+(a13-y)^2+(a31-y)^2+(a32-y)^2+(a33-y)^2=nt,a11*a33=0,a12*a32=0,a13*a31=0,a21*a23=0],[a11,a12,a13,a21,a23,a31,a32,a33]);
  • たくさんの解がでたが a_{i,j}a_{4-i,4-j}=0 をみたしているものは一つしかなかった
    •  {\it a_{11}}=0.0 ,\\ {\it a_{12}}=- 0.0389828727540191 ,\\ {\it a_{13}}=0.0 ,\\ {\it a_{21}}=- 0.293585307426138 ,\\ {\it a_{23}}=0.0 ,\\ {\it a_{31}}=- 0.0226987884695725 ,\\ {\it a_{32}}=0.0 ,\\ {\it a_{33}}=- 0.216284074605452
  • 全て負になってしまっているので正の条件をいれたい