Hello, I'm trying to solve vector polygon question. It's quite complicated,so i wrote it to matlab, and i'm getting error. Can anybody help me?
3 次查看(过去 30 天)
显示 更早的评论
I'm trying to solve vector polygon question. It's quite complicated,so i wrote it to matlab, and i'm getting error.
Program:
"
x = sym('x');
fi1 = sym('fi1');
fi3 = sym('fi3');
i3 = sym('i3');
i1 = sym('33');
i2 = sym('30');
i7=sym('18');
y = sym('y');
fi3=2*pi+2*(atan((i3-sqrt(2)* sqrt(-162+i3^2))/(18+i3)))
%fi3=2*pi+2*(atan((sqrt(2)* sqrt(-162+i3^2))/(18+i3)))
fi1=2*pi+2*atan((((2178-y^2)^0.5)-33)/(x-33))
%fi1=2*pi+2*atan(((-(2178-y^2)^0.5)-33)/(x-33))
d=solve(i1^2+i3^2-i2^2+i1*i3*cos(fi1-fi3)+x*[i1*cos(fi1)+i3*cos(fi3)]+x^2,i3)
d2=((20718*x^2*y^2 - 817638624*x - 324*x^3*y^2 - 3*x^2*y^4 + 6*x^4*y^2 + 883872*x^2*(2178 - y^2)^(1/2) - 15048*x^3*(2178 - y^2)^(1/2) + 264*x^4*(2178 - y^2)^(1/2) + 266112*y^2*(2178 - y^2)^(1/2) + x*(319519672582252032*x + 15908494561344*x^2*y^2 - 222521904000*x^3*y^2 - 4667000976*x^2*y^4 + 5889852144*x^4*y^2 + 41675040*x^3*y^4 - 41884128*x^5*y^2 + 463068*x^2*y^6 - 1044666*x^4*y^4 + 697500*x^6*y^2 - 1032*x^3*y^6 + 2340*x^5*y^4 - 2088*x^7*y^2 - 7*x^2*y^8 + 28*x^4*y^6 - 42*x^6*y^4 + 28*x^8*y^2 + 372213278908416*x^2*(2178 - y^2)^(1/2) - 6535809865728*x^3*(2178 - y^2)^(1/2) + 178286619456*x^4*(2178 - y^2)^(1/2) - 1809343008*x^5*(2178 - y^2)^(1/2) + 32997888*x^6*(2178 - y^2)^(1/2) - 220176*x^7*(2178 - y^2)^(1/2) + 2640*x^8*(2178 - y^2)^(1/2) - 7574107419648*y^2*(2178 - y^2)^(1/2) - 2877259968*y^4*(2178 - y^2)^(1/2) + 247104*y^6*(2178 - y^2)^(1/2) - 274110554234880*x*y^2 + 75823045056*x*y^4 - 7688736*x*y^6 + 126*x*y^8 - 6849698145725952*x*(2178 - y^2)^(1/2) - 17341459944956928*x^2 + 305128327477248*x^3 - 8410551732000*x^4 + 86586895296*x^5 - 1487713392*x^6 + 7897824*x^7 - 116559*x^8 + 654*x^9 - 7*x^10 + 788518683152640*y^2 + 106220572128*y^4 - 38733552*y^6 + 657*y^8 + 41241014899983360*(2178 - y^2)^(1/2) - 255558644352*x^2*y^2*(2178 - y^2)^(1/2) + 3232291392*x^3*y^2*(2178 - y^2)^(1/2) + 50903424*x^2*y^4*(2178 - y^2)^(1/2) - 84148416*x^4*y^2*(2178 - y^2)^(1/2) - 312048*x^3*y^4*(2178 - y^2)^(1/2) + 486288*x^5*y^2*(2178 - y^2)^(1/2) - 2640*x^2*y^6*(2178 - y^2)^(1/2) + 7920*x^4*y^4*(2178 - y^2)^(1/2) - 7920*x^6*y^2*(2178 - y^2)^(1/2) + 4317131934720*x*y^2*(2178 - y^2)^(1/2) - 815756832*x*y^4*(2178 - y^2)^(1/2) + 45936*x*y^6*(2178 - y^2)^(1/2) - 1928017446574222080)^(1/2) - 33*(319519672582252032*x + 15908494561344*x^2*y^2 - 222521904000*x^3*y^2 - 4667000976*x^2*y^4 + 5889852144*x^4*y^2 + 41675040*x^3*y^4 - 41884128*x^5*y^2 + 463068*x^2*y^6 - 1044666*x^4*y^4 + 697500*x^6*y^2 - 1032*x^3*y^6 + 2340*x^5*y^4 - 2088*x^7*y^2 - 7*x^2*y^8 + 28*x^4*y^6 - 42*x^6*y^4 + 28*x^8*y^2 + 372213278908416*x^2*(2178 - y^2)^(1/2) - 6535809865728*x^3*(2178 - y^2)^(1/2) + 178286619456*x^4*(2178 - y^2)^(1/2) - 1809343008*x^5*(2178 - y^2)^(1/2) + 32997888*x^6*(2178 - y^2)^(1/2) - 220176*x^7*(2178 - y^2)^(1/2) + 2640*x^8*(2178 - y^2)^(1/2) - 7574107419648*y^2*(2178 - y^2)^(1/2) - 2877259968*y^4*(2178 - y^2)^(1/2) + 247104*y^6*(2178 - y^2)^(1/2) - 274110554234880*x*y^2 + 75823045056*x*y^4 - 7688736*x*y^6 + 126*x*y^8 - 6849698145725952*x*(2178 - y^2)^(1/2) - 17341459944956928*x^2 + 305128327477248*x^3 - 8410551732000*x^4 + 86586895296*x^5 - 1487713392*x^6 + 7897824*x^7 - 116559*x^8 + 654*x^9 - 7*x^10 + 788518683152640*y^2 + 106220572128*y^4 - 38733552*y^6 + 657*y^8 + 41241014899983360*(2178 - y^2)^(1/2) - 255558644352*x^2*y^2*(2178 - y^2)^(1/2) + 3232291392*x^3*y^2*(2178 - y^2)^(1/2) + 50903424*x^2*y^4*(2178 - y^2)^(1/2) - 84148416*x^4*y^2*(2178 - y^2)^(1/2) - 312048*x^3*y^4*(2178 - y^2)^(1/2) + 486288*x^5*y^2*(2178 - y^2)^(1/2) - 2640*x^2*y^6*(2178 - y^2)^(1/2) + 7920*x^4*y^4*(2178 - y^2)^(1/2) - 7920*x^6*y^2*(2178 - y^2)^(1/2) + 4317131934720*x*y^2*(2178 - y^2)^(1/2) - 815756832*x*y^4*(2178 - y^2)^(1/2) + 45936*x*y^6*(2178 - y^2)^(1/2) - 1928017446574222080)^(1/2) + 422928*x*y^2 + 30*x*y^4 + 12388464*x*(2178 - y^2)^(1/2) - 33166584*x^2 + 727056*x^3 - 22239*x^4 + 294*x^5 - 3*x^6 - 23757624*y^2 + 1521*y^4 - 817638624*(2178 - y^2)^(1/2) - 264*x^2*y^2*(2178 - y^2)^(1/2) - 2376*x*y^2*(2178 - y^2)^(1/2) + 40473111888)/(4*(132*y^2*(2178 - y^2)^(1/2) - 2*x^2*y^2 - 132*x^2*(2178 - y^2)^(1/2) - 574992*x + 132*x*y^2 + 8712*x*(2178 - y^2)^(1/2) + 13068*x^2 - 132*x^3 + x^4 - 13068*y^2 + y^4 - 574992*(2178 - y^2)^(1/2) + 28462104)))^(1/2)
%result y(x)
d3=solve(y^2+(x-(i2^2-y^2)^0.5)^2-i1^2-d2^2-i1*d2*cos(fi1-fi3),y)
%{d3 should be my y(x) function, but matlab shows error:
d3
Warning: Explicit solution could not be found.
> In solve at 179
In Untitled at 36
d3 =
[ empty sym ]
}
1 个评论
pfb
2015-4-16
can you format the code (there is a button for that, in the editor) and provide more information?
回答(1 个)
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!