Re: [問題] 用MATLAB求解聯立方程式已回收
※ 引述《zupo (MATLAB程式編輯員)》之銘言：
: ※ 引述《anan0218 (大餅妹)》之銘言：
: : 如題~
: : 不才小妹在下我因為論文的關係來請教各位前輩
: : 首先以第一篇論文的概念來解說
: : 機率p,(1-p),q,(1-q),2個未知數(不同假設狀況對應的機率)
: : a,b,c,d,e,f,g,h為不同假設狀況的值,事先給定
: : 以
: : p[q*a+(1-q)*b]=(1-p)[q*c+(1-q)*d]----1
: : q[p*e+(1-p)*f]=(1-q)[p*g+(1-p)*h]----2
: : 兩條聯立方程式來求解(平面的概念)
: : 當初的作法是先手算求出p,q公式解,再代回原式
: : 以MATLAB算出均衡解
: : but..第二篇論文要以5條聯立方程式求解5個未知數(立體的概念)
: : 機率分別為p1,p2,p3=(1-p1-p2),q,(1-q),s1,s2,s3=(1-s1-s2)
: : 不同假設狀況的值為54個,事先給定
: : 是否可以直接在MATLAB裡列式,直接跑5條聯立方程式呢?
: : (因為手算展開實在是太多了...)
: : 已爬過文,但還是看不懂...
: 五條方程式打上來看一下好了
: 參數也給一下
: 是只要數值解就好嗎?
: 建議五條的聯立方程式還是不要求解析解.
: 答案可能不是人類所能閱讀的
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解...
--
推倒學長,比推導公式簡單...
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 134.208.33.244
討論串 (同標題文章)