Re: Finding inverse of non-linear transformation
- To: mathgroup at smc.vnet.net
- Subject: [mg116420] Re: Finding inverse of non-linear transformation
- From: "Christopher O. Young" <cy56 at comcast.net>
- Date: Mon, 14 Feb 2011 04:27:51 -0500 (EST)
- References: <ij8d29$396$1@smc.vnet.net>
On 2/13/11 5:50 AM, in article ij8d29$396$1 at smc.vnet.net, "Christopher O. Young" <cy56 at comcast.net> wrote: > I'm trying to find the inverse of the simple non-linear transformation > > (a x + b x + c) ({ > {x}, > {y}}) > > > I'm trying to use this as a sort of counter-example to the usual projective > transformation, which has the polynomial in the denominator. Does > Mathematica have a way to do this? InverseFunction doesn't seem to work > here. > > In[1]:= F[x_, y_, a, b, c] := (a x + b x + c) ( { {x}, {y} } ) > > In[2]:= InverseFunction[[x_, y_, a, b, c]] > >> > During evaluation of In[2]:= Part::pspec: Part specification x_ is neither > an integer nor a list of integers. >> > > > > Out[2]= InverseFunction[[x_, y_, a, b, c]] > I didn't have underscores for the a, b, and c the way I should have, but that wasn't the main problem. InverseFunction doesn't seem designed to find the inverse of vector transformations. But I think using "Solve" with added variables on the other sides of the equations gets me the result I'm looking for: Solve[{ (a x + b x + c) x == X, (a x + b x + c) y == Y }, {x, y} ] {{y -> (-((c Y)/(a + b)) - (Sqrt[c^2 + 4 a X + 4 b X] Y)/(a + b))/( 2 X), x -> (-c - Sqrt[c^2 + 4 a X + 4 b X])/( 2 (a + b))}, {y -> (-((c Y)/(a + b)) + ( Sqrt[c^2 + 4 a X + 4 b X] Y)/(a + b))/(2 X), x -> (-c + Sqrt[c^2 + 4 a X + 4 b X])/(2 (a + b))}}