Math Help Forum

 Site Feed

 
Go Back   Math Help Forum > High School Math Help > Geometry
Reply
 
Thread Tools Display Modes
  #1  
Old 05-19-2008, 09:15 AM
Junior Member
 
Join Date: Apr 2008
Location: London
Posts: 12
Country:
Thanks: 1
Thanked 0 Times in 0 Posts
potticus is on a distinguished road
Default intersection of a circle and a cubic equation

Other people have given solutions for a circle and a parabola, but they mainly involve finding a quadratic in x^2, suitable as they want 2 solutions.

In this case, I need a single solution to a circle/cubic problem. I know there's only one solution because the radius of the circle has been set to the smallest distance between the line and the centre point.

I tried substituting the cubic equation into the equation for the circle, then after a lot of bracket expanding (which I've checked over and over), I still don't get any intersection points.

I end up with a 6th order equation (which I believe to be correct - just thinking logically, it seems there can be 6 intercepts between a cubic and a circle). However, finding the roots of this equation (using the roots function in matlab) doesn't work :s

This is where i'm at so far:

With (x-xp)^2 + (y-yp)^2 = dc^2
And y = ax^3 + bx^2 + cx + d

After substitution, expansion and simplification:
0=(a^2)x^6 + (2ab)x^5 + (2ac + b^2)x^4 + (2ad + 2bc - 2ayp)x^3 + (2bd + c^2 - 2byp + 1)x^2 + (2cd - 2cyp -2xp)x + (d^2 + yp^2 - 2dyp + xp^2 - dc^2)

I can't see why this doesn't give me a suitable value for x intercept point. Can someone please help me out?

Thanks
Reply With Quote
  #2  
Old 05-19-2008, 09:56 PM
CaptainBlack's Avatar
Grand Panjandrum


 
Join Date: Nov 2005
Location: Somewhere near the south coast
Posts: 9,364
Country:
Thanks: 432
Thanked 2,278 Times in 1,921 Posts
CaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond repute
Default

Quote:
Originally Posted by potticus View Post
Other people have given solutions for a circle and a parabola, but they mainly involve finding a quadratic in x^2, suitable as they want 2 solutions.

In this case, I need a single solution to a circle/cubic problem. I know there's only one solution because the radius of the circle has been set to the smallest distance between the line and the centre point.

I tried substituting the cubic equation into the equation for the circle, then after a lot of bracket expanding (which I've checked over and over), I still don't get any intersection points.

I end up with a 6th order equation (which I believe to be correct - just thinking logically, it seems there can be 6 intercepts between a cubic and a circle). However, finding the roots of this equation (using the roots function in matlab) doesn't work :s

This is where i'm at so far:

With (x-xp)^2 + (y-yp)^2 = dc^2
And y = ax^3 + bx^2 + cx + d

After substitution, expansion and simplification:
0=(a^2)x^6 + (2ab)x^5 + (2ac + b^2)x^4 + (2ad + 2bc - 2ayp)x^3 + (2bd + c^2 - 2byp + 1)x^2 + (2cd - 2cyp -2xp)x + (d^2 + yp^2 - 2dyp + xp^2 - dc^2)

I can't see why this doesn't give me a suitable value for x intercept point. Can someone please help me out?

Thanks
This has six roots altogether, of which 0, 2, 4 or 6 may be real. But this gets you not much further than finding the roots of a general sixth order polynomial.

But there is no simple general solution for this. With numeric coefficients all the real roots can be found, but with symbolic coefficients (without some helpfull structure to them - which I don't see here but there may be) you are probably out of luck.

RonL
__________________
"It is proof of a base and low mind for one to wish to think with the masses or majority, merely because the majority is the majority"
--Giordano Bruno
Reply With Quote
  #3  
Old 05-20-2008, 02:48 AM
Junior Member
 
Join Date: Apr 2008
Location: London
Posts: 12
Country:
Thanks: 1
Thanked 0 Times in 0 Posts
potticus is on a distinguished road
Default

The equation and centre/radius of the circle are random, in effect, they can be anything within a domain I specify.

I know I'm looking for 1 root, but surely solving this would give me that, I don't get why it's not possible to do it algebraic-ly, when all the letters represent are real numbers?

edit: When I come to use this, there will be real numbers in there, I need a case I can program in and run in a loop for a range of circle centres and a range of cubic co-efficients. So if I have a case with real numbers in it, for example the first set in my list - a=b=c=d=1, dc=0.9589, xp=0, yp=4. Substituting these in, and finding the roots of the 6th order eqn won't work
Reply With Quote
  #4  
Old 05-20-2008, 06:23 AM
CaptainBlack's Avatar
Grand Panjandrum


 
Join Date: Nov 2005
Location: Somewhere near the south coast
Posts: 9,364
Country:
Thanks: 432
Thanked 2,278 Times in 1,921 Posts
CaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond repute
Default

Quote:
Originally Posted by potticus View Post
The equation and centre/radius of the circle are random, in effect, they can be anything within a domain I specify.

I know I'm looking for 1 root, but surely solving this would give me that, I don't get why it's not possible to do it algebraic-ly, when all the letters represent are real numbers?

edit: When I come to use this, there will be real numbers in there, I need a case I can program in and run in a loop for a range of circle centres and a range of cubic co-efficients. So if I have a case with real numbers in it, for example the first set in my list - a=b=c=d=1, dc=0.9589, xp=0, yp=4. Substituting these in, and finding the roots of the 6th order eqn won't work
Indeed because that circle and that cubic do not intersect at a real point.

RonL
__________________
"It is proof of a base and low mind for one to wish to think with the masses or majority, merely because the majority is the majority"
--Giordano Bruno
Reply With Quote
  #5  
Old 05-20-2008, 06:26 AM
Junior Member
 
Join Date: Apr 2008
Location: London
Posts: 12
Country:
Thanks: 1
Thanked 0 Times in 0 Posts
potticus is on a distinguished road
Default

but they do, just at one point. that dc value is the smallest distance from the point to the cubic, so a circle with that radius will touch at that specific point.

edit: perhaps that is a matlab rounding error?
Reply With Quote
  #6  
Old 05-20-2008, 06:42 AM
CaptainBlack's Avatar
Grand Panjandrum


 
Join Date: Nov 2005
Location: Somewhere near the south coast
Posts: 9,364
Country:
Thanks: 432
Thanked 2,278 Times in 1,921 Posts
CaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond reputeCaptainBlack has a reputation beyond repute
Default

Quote:
Originally Posted by potticus View Post
but they do, just at one point. that dc value is the smallest distance from the point to the cubic, so a circle with that radius will touch at that specific point.

edit: perhaps that is a matlab rounding error?
My conclusion that they don't intersect is based on a plot which may be a bit crude. However XL solver thinks the minimum distance between the curves is 0.027.

RonL
__________________
"It is proof of a base and low mind for one to wish to think with the masses or majority, merely because the majority is the majority"
--Giordano Bruno
Reply With Quote
  #7  
Old 05-20-2008, 07:26 AM
Junior Member
 
Join Date: Apr 2008
Location: London
Posts: 12
Country:
Thanks: 1
Thanked 0 Times in 0 Posts
potticus is on a distinguished road
Default

hmm, my plots show nothing that small ... in the order of 1 appears to be correct.

i'll try just the plot in excel or something, and see what i get. i'm pretty sure that distance is right though :s
Reply With Quote
Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are Off
Refbacks are Off
Forum Jump


All times are GMT -7. The time now is 07:20 PM.


Powered by vBulletin® Version 3.7.3
Copyright ©2000 - 2008, Jelsoft Enterprises Ltd.
SEO by vBSEO 3.2.0 ©2008, Crawlability, Inc.
©2005 - 2008 Math Help Forum


Math Help Forum is a community of maths forums with an emphasis on maths help in all levels of mathematics.
Register to post your math questions or just hang out and try some of our math games or visit the arcade.