Re: [問題] 用MATLAB求解聯立方程式已回收

看板MATLAB作者 (誰說誰是誰)時間14年前 (2010/04/14 11:29), 編輯推噓0(000)
留言0則, 0人參與, 最新討論串4/4 (看更多)
※ 引述《anan0218 (大餅妹)》之銘言: : ※ 引述《zupo (MATLAB程式編輯員)》之銘言: : : 五條方程式打上來看一下好了 : : 參數也給一下 : : 是只要數值解就好嗎? : : 建議五條的聯立方程式還是不要求解析解. : : 答案可能不是人類所能閱讀的 : p1*[q*s1*h1+q*s2*h2+q*(1-s1-s2)*h3+(1-q)*s1*h10+(1-q)*s2*h11+ : (1-q)*(1-s1-s2)*h12]--------------------------------------------------1 : p2*[q*s1*h4+q*s2*h5+q*(1-s1-s2)*h6+(1-q)*s1*h13+(1-q)*s2*h14+ : (1-q)*(1-s1-s2)*h12]--------------------------------------------------2 : (1-p1-p2)*[q*s1*h7+q*s2*h8+q*(1-s1-s2)*h9+(1-q)*s1*h16+(1-q)*s2*h11+ : (1-q)*(1-s1-s2)*h18]-------------------------------------------3 : s1*[q*p1*t7+q*p2*t4+q*(1-p1-p2)*t7+(1-q)*p1*t10+(1-q)*p2*t13+ : (1-q)*(1-p1-p2)*t16]--------------------------------------------------4 : s2*[q*p1*t2+q*p2*t5+q*(1-p1-p2)*t9+(1-q)*p1*t11+(1-q)*p2*t14+ : (1-q)*(1-p1-p2)*t17]--------------------------------------------------5 : (1-s1-s2)*[q*p1*t3+q*p2*t6+q*(1-p1-p2)*t9+(1-q)*p1*t12+(1-q)*p2*t15+ : (1-q)*(1-p1-p2)*t18]-------------------------------------------6 : q*[p1*s1*f1+p1*s2*f2+p1*(1-s1-s2)*f3+p2*s1*f4+p2*s2*f5+p2*(1-s1-s2)*f6+ : (1-p1-p2)*s1*f7+(1-p1-p2)*s2*f8+(1-p1-p2)*(1-s1-s2)*f9]----------------7 : (1-q)*[p1*s1*f10+p1*s2*f11+p1*(1-s1-s2)*f12+p2*s1*f13+p2*s2*f14+ : p2*(1-s1-s2)*f15+(1-p1-p2)*s1*f16+(1-p1-p2)*s2*f17+ : (1-p1-p2)*(1-s1-s2)*f18]-------------------------------------------8 : 其中h1~h18,t1~t18,f1~f18為事先給定 : 求未知數p1,p2,p3=(1-p1-p2),s1,s2,s3=(1-s1-s2),q,(1-q) : 1=2 : 2=3 : 4=5 : 5=6 : 7=8 : 初步的想法是以solve解... syms 所有代數 eq1=式1 eq2=式2 eq3=式3 . . . slove(eq1,eq2....,'欲求解的代數') 這樣應該就完成了 -- 我只能說 你的簽名檔太強大了= ..= -- moon16star:屈原表示: 03/17 23:56

03/17 23:57,
找我嗎?
03/17 23:57

03/17 23:59,
來汨羅江找我吧!
03/17 23:59

03/18 00:00,
汨羅江online
03/18 00:00
-- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.120.45.66 ※ 編輯: laserguy 來自: 140.120.45.66 (04/14 11:31)
文章代碼(AID): #1BnJQtol (MATLAB)
文章代碼(AID): #1BnJQtol (MATLAB)