1
$\begingroup$

The following is an error correcting equation for a sidereal astrophotography tracking mount I'm building and $t$ represents the amount of time before the tracking is off by a quarter of a stepper motor step. I need to solve for $t$, but I've having trouble solving this.

I can plot it with a value of $n$ and find an ok approximation, but the problem is that I need it for several values of $n$. Any good ways? $n$ is the integer number of corrections applied thus far. If it can't be solved in a general sense, what's a good way to generate approximate solutions for a few hundred values of $n$, starting at 0?

$$-0.25=\frac{\sqrt{2\cdot 150^2-2\cdot 150^2\cdot \cos(t\cdot 0.000072733)} - (t\cdot 0.0109170306+(n\cdot -0.25\cdot 0.005))}{0.005}$$

  • 2
    Please evaluate the constants, like 2*150^2=45000, then clear the denominator (multiply by 200) before posting. If you do that, take the second term on the right to the left and square, you get a quadratic. Do you know how to solve those?2011-03-28
  • 2
    @Ross: Am I missing something? I see $t$ both inside and outside the cosine, so I would not expect an algebraic solution...2011-03-28
  • 0
    Is the cosine inside the radical supposed to be squared? Just wondering, as that would simplify the radical.2011-03-28
  • 0
    You are right, $t$ is inside the cosine, but you could still simplify as far as you can. With $t$ inside the cosine you are probably into a numeric solution instead of algebraic. Seeing the constants clearly can help-sometimes an iterative solution is appropriate, sometimes a basic root find. Also, if $t$ is small compared with 1/.00007=1.4E+4, you can use a small angle approximation.2011-03-28
  • 1
    Using Mathematica to numerically approximate the solution to the original equation and numerically solve the polynomial equation resulting from using the 4th degree power series approximation for cosine, it appears that $\cos x\approx 1-\frac{x^2}{2}+\frac{x^4}{24}$ is close enough that the solution is within 0.1% for $0\le n\le 169$, 0.2% for $0\le n\le 429$.2011-03-28
  • 0
    @Isaac: And I'll bet that if you discard the $x^4$ term, it is still pretty good.2011-03-28
  • 0
    @Ross: Not so much, because the solution values for $t$ are in the thousands and with a $t$ inside the cosine. The 2nd degree series just isn't good enough (the error seems to be greater than 10% for $n\ge 3$, greater than 100% for $n\ge 19$, greater than 1000% for $n\ge 338$).2011-03-28
  • 0
    @Isaac: Good point. And an example of how we (I) tend to think of variables as being of order of magnitude 1, even without reason to do so. Maybe OP will give us a range of interest.2011-03-28
  • 0
    We can use a half-angle formula to simplify the radical expression. If $0.0109...$ is supposed to be exactly $150 \cdot 0.0000727...$, we can simplify the entire equation rather nicely: $n + 1 = 120000 ( A t - 2 \sin( A t /2 ) )$, where $A := 0.000072733...$.2011-03-28
  • 0
    @DayLateDon - It would be if exact construction were possible. In an exact construction, 150 would actually be 150.0973588306.2011-03-28

1 Answers 1

2

Here is a table of values of $t$ for integers $0\le n\le 500$, generated using Mathematica's FindRoot.

Here are the first few values:

n=  0:  t= 174.7 n=  1:  t= 339.8 n=  2:  t= 489.7 n=  3:  t= 623.7 n=  4:  t= 743.3 n=  5:  t= 850.4 n=  6:  t= 947.2 n=  7:  t=1035.4 n=  8:  t=1116.4 n=  9:  t=1191.3 n= 10:  t=1261.0 

Here's the rough source code:

Table[{n,    t /. FindRoot[-0.25 == (Sqrt[          2*150^2 -            2*150^2*Cos[             t*0.000072733]] - (t*0.0109170306 + (n*-0.25*.005)))/       0.005, {t, 200 + 150 n}]}, {n, 0, 500}] 
  • 0
    @rfusca: You're welcome. It's probably worth some checking to make sure that the numbers are sensible, as I really did just blindly beat the problem over the head with Mathematica.2011-03-28
  • 0
    I had computed the first 12, and they line up, so I'm willing to run with it2011-03-28
  • 0
    Just an FYI for anybody looking at this, what's interesting is that from an engineering perspective this means that starting out it will need to be error corrected every few minutes but after an hour its diverging fast enough at it needs error correction every few seconds. This is an isosceles Haig mount for those interested.2011-03-28