用Matlab编写一个12元的方程组,并求符号解,程序也没调好,麻烦各位了
各位大侠能帮我编一个12元的方程组,求这12个方程组的符号解,系数都认为是符号。我刚学编程,编了一个,只是一直没有调试好,万分感谢。下面的是我推导的公式:
(r2-r1)*x1+i*g1*x2-i*g1*x3+r2=0
i*g1*x1+(-i*w2+i*w1)*x2+i*g2*x4+i*g2*x5-i*g1*x10=0
-i*g1*x1+(-i*w1+i*w2)*x3-i*g2*x7-i*g2*x9+i*g1*x10=0
(-i*w2+i*w1+i*w3)*x4+i*g1*x8=0
(-i*w2+i*w1-i*w3)*x5+i*g1*x6=0
i*g1*x5+(-i*w3-r3/2-r1+r2)*x6-i*g1*x7+r2*x11=0 (1)
-i*g1*x6+(i*w2-i*w1-i*w3)*x7=0
i*g1*x4+(i*w3-r3/2-r1+r2)*x8-i*g1*x9+r2*x12=0
-i*g1*x8+(i*w2-i*w1+i*w3)*x9=0
-i*g1*x2+i*g1*x3=0
(-i*w3-r3/2)*x11+i*g2*x1=0
(i*w3-r3/2)*x12-i*g2*x1=0)@wuming524 返回小木虫查看更多
今日热帖
syms x1 x2 x3 x4 x5 x6 x7 x8 x9 x10 x11 x12 r1 r2 r3 g1 g2 w1 w2 w3
eqns=[(r2-r1)*x1+i*g1*x2-i*g1*x3+r2==0,...
i*g1*x1+(-i*w2+i*w1)*x2+i*g2*x4+i*g2*x5-i*g1*x10==0,...
-i*g1*x1+(-i*w1+i*w2)*x3-i*g2*x7-i*g2*x9+i*g1*x10==0,...
(-i*w2+i*w1+i*w3)*x4+i*g1*x8==0,...
(-i*w2+i*w1-i*w3)*x5+i*g1*x6==0,...
i*g1*x5+(-i*w3-r3/2-r1+r2)*x6-i*g1*x7+r2*x11==0,...
-i*g1*x6+(i*w2-i*w1-i*w3)*x7==0,...
i*g1*x4+(i*w3-r3/2-r1+r2)*x8-i*g1*x9+r2*x12==0,...
-i*g1*x8+(i*w2-i*w1+i*w3)*x9==0,...
-i*g1*x2+i*g1*x3==0,...
(-i*w3-r3/2)*x11+i*g2*x1==0,...
(i*w3-r3/2)*x12-i*g2*x1==0]
S=solve(eqns,x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12)
sol=[S.x1,S.x2,S.x3,S.x4,S.x5,S.x6,S.x7,S.x8,S.x9,S.x10,S.x11,S.x12]
但是好像没有解哎,
膜拜一下!
看得头晕。。。大体看起来还行。
我老师解出来了,他让我自己调试,我解不出来,我再试试吧。
符号解,还是mathmatica更合适吧