Re: [運算] 請問十元一次聯立方程式的解法已回收

看板MATLAB作者vicoju (啾)時間15年前 (2011/03/23 13:11)推噓3(3推 0噓 5→)

留言8則, 4人參與討論串2/2 (看更多)

※ 引述《vicoju (啾)》之銘言： : 小弟我在算一顆子彈的參數不過十元一次聯立方程式真的難倒我了 : 請各位高手幫我解答一下a.b.c.d.e.f.g.h.i.j是我要求的未知數 : eq1=sym('0=a/11.5*cos(8.5*180/pi)+b/32*cos(60*180/pi)+c/54*cos(68*180/pi)+d/84*cos(72*180/pi)+e/111*cos(74*180/pi)+f/138*cos(75.5*180/pi)+g/164*cos(76.7*180/pi)+h/190*cos(77.1*180/pi)+i/217*cos(77.4*180/pi)+j/244*cos(77.6*180/pi)') : eq2=sym('0=a/29*(-cos(65*180/pi))+b/13*cos(7*180/pi)+c/32*cos(66*180/pi)+d/58*cos(71*180/pi)+e/84*cos(80*180/pi)+f/110*cos(82*180/pi)+g/137*cos(83*180/pi)+h/164*cos(85*180/pi)+i/190*cos(86*180/pi)+j/216*cos(87*180/pi)') : eq3=sym('0=a/54*(-cos(63*180/pi))+b/29*(-cos(77*180/pi))+c/13.5*cos(5*180/pi)+d/31*cos(62*180/pi)+e/56*cos(76*180/pi)+f/81.5*cos(80*180/pi)+g/108.4*cos(83*180/pi)+h/134.6*cos(85.3*180/pi)+i/161*cos(85.9*180/pi)+j/187*cos(86.1*180/pi)') : eq4=sym('0=a/83*(-cos(113*180/pi))+b/57*(-cos(108*180/pi))+c/31*(-cos(95*180/pi))+d/19*(-cos(37*180/pi))+e/29.5*cos(28*180/pi)+f/54*cos(43*180/pi)+g/80*cos(47*180/pi)+h/106*cos(50*180/pi)+i/132.5*cos(52*180/pi)+j/159*cos(54*180/pi)') : eq5=sym('0=a/111.5*(-cos(88.5*180/pi))+b/85.5*(-cos(75*180/pi))+c/61*(-cos(68*180/pi))+d/37.5*(-cos(50.5*180/pi))+e/24.5*(-cos(3*180/pi))+f/34.8*cos(44*180/pi)+g/57.5*cos(62*180/pi)+h/83*cos(75*180/pi)+i/103.5*cos(78.5*180/pi)+j/134*cos(80.4*180/pi)') : eq6=sym('0=a/137*(-cos(80*180/pi))+b/113*(-cos(78*180/pi))+c/86.9*(-cos(73*180/pi))+j/150*cos(78*180/pi)+d/61.5*(-cos(68*180/pi))+e/38*(-cos(51*180/pi))+f/23*(-cos(10*180/pi))+g/32*cos(45*180/pi)+h/54*cos(65*180/pi)+i/80*cos(71*180/pi)') : eq7=sym('0=a/168*(-cos(83*180/pi))+b/141*(-cos(80*180/pi))+c/115*(-cos(78*180/pi))+d/89*(-cos(74*180/pi))+e/64*(-cos(69*180/pi))+f/40*(-cos(55*180/pi))+g/23*(-cos(15*180/pi))+h/30.5*cos(43*180/pi)+i/52*cos(65*180/pi)+j/77*cos(72*180/pi)') : eq8=sym('0=a/196.5*(-cos(83.5*180/pi))+b/170*(-cos(82*180/pi))+c/144*(-cos(81*180/pi))+d/116.5*(-cos(79*180/pi))+e/91*(-cos(77.5*180/pi))+f/66*(-cos(70*180/pi))+g/41.5*(-cos(58*180/pi))+h/24*(-cos(20*180/pi))+i/30*cos(40.5*180/pi)+j/51.5*cos(60*180/pi)') : eq9=sym('0=a/224*(-cos(81*180/pi))+b/198*(-cos(80*180/pi))+c/171*(-cos(79*180/pi))+d/144*(-cos(77.5*180/pi))+e/119*(-cos(75*180/pi))+f/92*(-cos(71*180/pi))+g/66*(-cos(69*180/pi))+h/42*(-cos(58*180/pi))+i/24.5*(-cos(21*180/pi))+j/28*cos(38*180/pi)') : eq10=sym('0=a/247*(-cos(85*180/pi))+b/221*(-cos(83*180/pi))+c/193*(cos(82*180/pi))+d/167*(-cos(81*180/pi))+e/142*(-cos(79*180/pi))+f/116*(-cos(78*180/pi))+g/89*(-cos(76*180/pi))+h/64*(-cos(69*180/pi))+i/40*(-cos(55*180/pi))+j/24*(-cos(18*180/pi))') : sol=fsolve(eq1,eq2,eq3,eq4,eq5,eq6,eq7,eq8,eq9,eq10) : 請問是哪邊出錯了~~~!!!!!!!! : a .b .c. d.e.f.g.h.i.j為我要求的未知數其他都是已知 10個未知數 10個方程式 : 應該解的出來才對請各位高手指點一下我改成矩陣 A=[a/11.5*cos(8.5) b/32*cos(60) c/54*cos(68) d/84*cos(72) e/111*cos(74) f/138*cos(75.5) g/164*cos(76.7) h/190*cos(77.1) i/217*cos(77.4) j/244*cos(77.6) a/29*(-cos(65)) b/13*cos(7) c/32*cos(66) d/58*cos(71) e/84*cos(80) f/110*cos(82) g/137*cos(83) h/164*cos(85) i/190*cos(86) j/216*cos(87) a/54*(-cos(63)) b/29*(-cos(77)) c/13.5*cos(5) d/31*cos(62) e/56*cos(76) f/81.5*cos(80) g/108.4*cos(83) h/134.6*cos(85.3) i/161*cos(85.9) j/187*cos(86.1) a/83*(-cos(113)) b/57*(-cos(108)) c/31*(-cos(95)) d/19*(-cos(37)) e/29.5*cos(28) f/54*cos(43) g/80*cos(47) h/106*cos(50) i/132.5*cos(52) j/159*cos(54) a/111.5*(-cos(88.5)) b/85.5*(-cos(75)) c/61*(-cos(68)) d/37.5*(-cos(50.5)) e/24.5*(-cos(3)) f/34.8*cos(44) g/57.5*cos(62) h/83*cos(75) i/103.5*cos(78.5) j/134*cos(80.4) a/137*(-cos(80)) b/113*(-cos(78)) c/86.9*(-cos(73)) j/150*cos(78) d/61.5*(-cos(68)) e/38*(-cos(51)) f/23*(-cos(10)) g/32*cos(45) h/54*cos(65) i/80*cos(71) a/168*(-cos(83)) b/141*(-cos(80)) c/115*(-cos(78)) d/89*(-cos(74)) e/64*(-cos(69)) f/40*(-cos(55)) g/23*(-cos(15)) h/30.5*cos(43) i/52*cos(65) j/77*cos(72) a/196.5*(-cos(83.5)) b/170*(-cos(82)) c/144*(-cos(81)) d/116.5*(-cos(79)) e/91*(-cos(77.5)) f/66*(-cos(70)) g/41.5*(-cos(58)) h/24*(-cos(20)) i/30*cos(40.5) j/51.5*cos(60) a/224*(-cos(81)) b/198*(-cos(80)) c/171*(-cos(79)) d/144*(-cos(77.5)) e/119*(-cos(75)) f/92*(-cos(71)) g/66*(-cos(69)) h/42*(-cos(58)) i/24.5*(-cos(21)) j/28*cos(38) a/247*(-cos(85)) b/221*(-cos(83)) c/193*(cos(82)) d/167*(-cos(81)) e/142*(-cos(79)) f/116*(-cos(78)) g/89*(-cos(76)) h/64*(-cos(69)) i/40*(-cos(55)) j/24*(-cos(18))] B=[0 0 0 0 0 0 0 0 0 0] ??? Undefined function or variable 'a'. 我的未知數在分子要如何寫才能A\B呢 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 192.83.195.240

推

Gwaewluin

03/23 13:47, , 1^F

03/23 13:47, 1^F

→

Gwaewluin

03/23 13:48, , 2^F

03/23 13:48, 2^F

→

Gwaewluin

03/23 13:49, , 3^F

03/23 13:49, 3^F