技術士(総合技術監理・機械部門)のブログ

技術士がお届けする,技術,工学,技術者倫理などの話題

scilab 19 '

2017-10-10 20:27:45 | 工学

さて、前回のmatlabとscilabの極配置の結果の違いは
アルゴリズムの違いと言うしかない
以下matlabで

>> ke=place(Abb',Aab',L)';
>> ke
ke =

14.4000 0.6000
0.3000 15.7000

と確かに、0.3だ.
さて

>> eig(Abb'-Aab'*ke)

ans =

-15.1000
-15.9000

あれ、L=[-15 -16]
だったので、誤差がある.0.1
一方、scilabの結果では 0.3を0にして
>> ke1=[14.4 0.6;0.0 15.7];
>> eig(Abb'-Aab'*ke1)

ans =

-15.2354
-15.7646

まあ、この誤差の方が大きい

これはmatlabに軍配か