[其他] 用LibreOffice解非線性規劃問題

[bugmens (2013新年快樂)]

這是剛才有網友所提問的問題
 (1.72*10^4)*x[1]+x[9]*((8.4*10^-5)-x[1])=1.702
 (1.72*10^4)*x[2]+x[9]*((8.4*10^-5)-x[2])=1.705
 (1.72*10^4)*x[3]+x[9]*((8.4*10^-5)-x[3])=1.673
 (1.72*10^4)*x[4]+x[9]*((8.4*10^-5)-x[4])=1.629
 (1.72*10^4)*x[5]+x[9]*((8.4*10^-5)-x[5])=1.608
 (1.72*10^4)*x[6]+x[9]*((8.4*10^-5)-x[6])=1.634
 (1.72*10^4)*x[7]+x[9]*((8.4*10^-5)-x[7])=1.617
 (1.72*10^4)*x[8]+x[9]*((8.4*10^-5)-x[8])=1.609
 (1.15*10^4)*x[1]+x[10]*((8.4*10^-5)-x[1])=1.271
 (1.15*10^4)*x[2]+x[10]*((8.4*10^-5)-x[2])=1.286
 (1.15*10^4)*x[3]+x[10]*((8.4*10^-5)-x[3])=1.265
 (1.15*10^4)*x[4]+x[10]*((8.4*10^-5)-x[4])=1.239
 (1.15*10^4)*x[5]+x[10]*((8.4*10^-5)-x[5])=1.222
 (1.15*10^4)*x[6]+x[10]*((8.4*10^-5)-x[6])=1.241
 (1.15*10^4)*x[7]+x[10]*((8.4*10^-5)-x[7])=1.226
 (1.15*10^4)*x[8]+x[10]*((8.4*10^-5)-x[8])=1.22

找出x[1]~x[10]滿足上面16個式子
------------------------------------------------
結果
x1      -0.0000718352
x2      -0.0000794992
x3      -0.0000687698
x4      -0.0000554851
x5      -0.0000467990
x6      -0.0000565097
x7      -0.0000488421
x8      -0.0000457684
x9      18546.4610124886
x10     13457.1978330652



軟體執行畫面https://dl.dropboxusercontent.com/u/23455489/06.18.gif
檔案下載https://dl.dropboxusercontent.com/u/23455489/06.18.ods

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 1.164.220.129
※ 編輯: bugmens         來自: 1.164.220.129        (06/18 23:52)
推 i19890813   :先推 我研究一下!!  感謝!   單看結果不如預期QQ
→ bugmens     :等一下,我的式子還有錯,我再檢驗看看
※ 編輯: bugmens         來自: 1.164.220.129        (06/19 00:01)
※ 編輯: bugmens         來自: 1.164.220.129        (06/19 00:06)
→ bugmens     :我改好了,只是後面8個式子的誤差都在0.3~0.4
→ bugmens     :我跑了好幾次也沒辦法再讓誤差減少
→ bugmens     :你看看這樣的結果能不能接受
→ i19890813   :抱歉我...後面八條的第一個係數打錯....
→ i19890813   :1.72*10^4  -->  1.15*10^4   非常抱歉QQ
※ 編輯: bugmens         來自: 1.164.220.129        (06/19 00:24)
→ bugmens     :我改好了,誤差總和最小值為0.14
推 i19890813   :謝謝bug大