| 查看: 974 | 回复: 5 | |||
[求助]
matlab求解非线性方程组
|
|
clear syms OX OY AX AY BX BY CX CY DX DY EX EY FX FY GX GY HX HY IX TY JX JY en1=sym('((OX-AX)^2+(OY-AY)^2)^0.5=30'); en2=sym('((BX-AX)^2+(BY-AY)^2)^0.5=62'); en3=sym('((BX-DX)^2+(DX-BY)^2)^0.5=16.5'); en4=sym('((BX-CX)^2+(BY-CY)^2)^0.5=48.5'); en5=sym('((EX-BX)^2+(EY-BY)^2)^0.5=30'); en6=sym('((GX-DX)^2+(GY-DY)^2)^0.5=68'); en7=sym('((FX-GX)^2+(FY-GY)^2)^0.5=64.882'); en8=sym('((HX-GX)^2+(HY-GY)^2)^0.5=52.647'); en9=sym('((FX-HX)^2+(FY-HY)^2)^0.5=114'); en10=sym('((IX-EX)^2+(IY-EY)^2)^0.5=114'); en11=sym('((HX-IX)^2+(HY-IY)^2)^0.5=32'); en12=sym('(AX-DX)^2+(AY-DY)^2=62^2+16.5^2'); en13=sym('((JX-IX)^2+(JY-IY)^2)^0.5=90'); en14=sym('((EX-CX)^2+(EY-CY)^2)^0.5=48.5+30'); en15=sym('((JX-HX)^2+(JY-HY)^2)^0.5=32+90'); en16=sym('FX=111'); en17=sym('abs(FY)=5'); en18=sym('CX=70.5'); en19=sym('abs(CY)=15'); en20=sym('OX=0'); en21=sym('OY=0'); en22=sym('XA=0'); %en23=sym('YA=-30'); [JX,JY]=solve(en1,en2,en3,en4,en5,en6,en7,en8,en9,en10,en11,en12,en13,en14,en15,en16,en17,en18,en19,en20,en21,en22,OX, OY, AX ,AY ,BX, BY,CX, CY, DX, DY ,EX ,EY, FX ,FY ,GX, GY, HX, HY ,IX, TY, JX ,JY); 这个方程组求出来的是无解的状态，为什么？有没有其他办法？ |
回复此楼» 猜你喜欢
290调剂生物0860
已经有32人回复
-
271求调剂
已经有28人回复
-
105500药学求调剂
已经有4人回复
-
274求调剂
已经有10人回复
-
335求调剂
已经有20人回复
-
化工学硕294分，求导师收留
已经有33人回复
-
290求调剂
已经有22人回复
-
291 求调剂
已经有36人回复
-
食品与营养（0955）271求调剂
已经有16人回复
-
22408 312求调剂
已经有10人回复
WashingGong
金虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 16794.7
- 散金: 20
- 红花: 2
- 帖子: 1856
- 在线: 170.1小时
- 虫号: 2307033
- 注册: 2013-02-28
- 性别: GG
- 专业: 传动机械学
2楼2015-11-26 16:51:42
|
function f=fun1(x) OX=0; OY=0; CX=70.5; CY=15; FX=111; FY=5; XA=0; en1=((OX-XA)^2+(OY-x(15))^2)^0.5-30; en2=((x(1)-XA)^2+(x(2)-x(15))^2)^0.5-62; en3=((x(1)-x(3))^2+(x(3)-x(2))^2)^0.5-16.5; en4=((x(1)-CX)^2+(x(2)-CY)^2)^0.5-48.5; en5=((x(5)-x(1))^2+(x(6)-x(2))^2)^0.5-30; en6=((x(7)-x(3))^2+(x(8)-x(4))^2)^0.5-68; en7=((FX-x(7))^2+(FY-x(8))^2)^0.5-64.882; en8=((x(9)-x(7))^2+(x(10)-x(8))^2)^0.5-52.647; en9=((FX-x(9))^2+(FY-x(10))^2)^0.5-114; en10=((x(11)-x(5))^2+(x(12)-x(6))^2)^0.5-114; en11=((x(9)-x(11))^2+(x(10)-x(12))^2)^0.5-32; en12=(XA-x(3))^2+(x(15)-x(4))^2-62^2-16.5^2; en13=((x(13)-x(11))^2+(x(14)-x(12))^2)^0.5-90; en14=((x(5)-CX)^2+(x(6)-CY)^2)^0.5-48.5-30; en15=((x(13)-x(9))^2+(x(14)-x(10))^2)^0.5-32-90; f=[en1;en2;en3;en4;en5;en6;en7;en8;en9;en10;en11;en12;en13;en14;en15]; |
3楼2015-11-27 19:48:25
4楼2015-11-27 19:52:54
5楼2015-11-29 10:34:39
6楼2015-11-30 10:02:13
5