Path: news.mathworks.com!not-for-mail From: <HIDDEN> Newsgroups: comp.soft-sys.matlab Subject: Re: solving a system of multi-variable polynomial equations in mupad Date: Tue, 9 Oct 2012 08:02:08 +0000 (UTC) Organization: The MathWorks, Inc. Lines: 41 Message-ID: <k50lm0$qrj$1@newscl01ah.mathworks.com> References: <k4ib2c$eel$1@newscl01ah.mathworks.com> <k4l3c0$2ed$1@newscl01ah.mathworks.com> <k4n7o4$1hf$1@newscl01ah.mathworks.com> <k4nprj$9rg$1@newscl01ah.mathworks.com> <k4vg5t$nss$1@newscl01ah.mathworks.com> Reply-To: <HIDDEN> NNTP-Posting-Host: www-02-blr.mathworks.com Content-Type: text/plain; charset=UTF-8; format=flowed Content-Transfer-Encoding: 8bit X-Trace: newscl01ah.mathworks.com 1349769728 27507 172.30.248.47 (9 Oct 2012 08:02:08 GMT) X-Complaints-To: news@mathworks.com NNTP-Posting-Date: Tue, 9 Oct 2012 08:02:08 +0000 (UTC) X-Newsreader: MATLAB Central Newsreader 1187260 Xref: news.mathworks.com comp.soft-sys.matlab:780205 "Daniel" wrote in message <k4vg5t$nss$1@newscl01ah.mathworks.com>... > ....... I am equally interested in understanding when the system has no solutions as when it does have solutions. ....... - - - - - - - - - - - Daniel, you have said, "... I am equally interested in understanding when the system has no solutions as when it does have solutions." I will give you a summary here of my conclusions on that question. To begin with, there can be no solutions unless the two quantities which I have called F are zero. That is, s-t_h and t-s_h must both be zero. In fact none of the sets of just six equations can have a solution unless the associated F is zero. This can be demonstrated using the rotation method which I described earlier which with nonzero F's leads to contradictions in each of the three separate cases B^2-A*C < 0, B^2-A*C = 0, and B^2-A*C > 0. Assuming zero F's, I find that seeing the situation from the point of view of six x-y points in two-dimensional space and two A-B-C points in three dimensional space to be useful. The six points are the two sets, (a1,b1), (c1,d1), (e1,f1) and (a2,b2), (c2,d2), (e2,f2) which I have called (x1,y1), (x2,y2), (x3,y3), and (x4,y4), (x5,y5), and (x6,y6). The two A-B-C points are (A1,B1,C1) = (alpha1-beta1_h,alpha2-beta2_h,alpha3-beta3_h), and (A2,B2,C2) = (beta1-alpha1_h,beta2-alpha2_h,beta3-alpha3_h) . First, if the six x-y points are all located at the origin, that is, all zero, then any two A-B-C points anywhere in their 3D space will constitute a solution. That is quite obvious. Second, if all six x-y points lie along some one line through the origin with not all points at the origin and if a is the angle between that line and the x-axis, then any two points (A1,B1,C1) and (A2,B2,C2) will constitute a solution to the 24 equations provided that each is orthogonal to the 3D vector (cos(a)^2,2*sin(a)*cos(a),sin(a)^2), that is if A1*cos(a)^2 + 2*B1*sin(a)*cos(a) + C1*sin(a)^2 = 0 and A2*cos(a)^2 + 2*B2*sin(a)*cos(a) + C2*sin(a)^2 = 0 . Thus the A-B-C points must lie in a certain plane containing the origin. Showing that to be true for the equations of the form of A*x2*x3 + B*(x2*y3+x3*y2) + C*y2*y3 = 0 is also accomplished using the aforementioned rotation method. Finally, if the three points of the first x-y set lie along one line through the origin, not all at the origin, and if the second set of three points lie along a different line through the origin, again not all at the origin, and if a and b are the respective angles between the x-axis and these lines, then the set of A-B-C point pairs that provide solutions are all pairs that are orthogonal to both (cos(a)^2,2*sin(a)*cos(a),sin(a)^2) and (cos(b)^2,2*sin(b)*cos(b),sin(b)^2) . Stated equivalently, it is those A-B-C points which lie along the cross product of these vectors. However we can remove the common nonzero factor, sin(a-b), from this cross product and arrive at the simpler vector (-2*sin(a)*sin(b),sin(a+b),-2*cos(a)*cos(b)) . Thus any two A-B-C points along this direction will constitute a solution to the 24 equations. The A-B-C points must therefore lie in a certain one-dimensional line through the origin for such a set of six x-y points. All other combinations of points in these respective x-y and A-B-C spaces are non-solutions. Again I point out that MuPAD would have had a very difficult time explaining to humans, not to mention arriving at, such a complicated situation. Roger Stafford