Solving quadratic trig equations

A few days ago I wrote about how to systematically solve trig equations. That post was abstract and general. This post will be concrete and specific, looking at the special case of quadratic equations in sines and cosines, i.e. any equation of the form

As outlined earlier, we turn the equation into a system of equations in s and c.

The resultant of

and

as a function of sis

where

Example 1

Let's look at a particular example. Suppose we want to solve

Then the possible sine values are the roots of

This equation as four real roots: s = -0.993462, -0.300859, -0.0996236, or 0.966329.

So any solution to our original equation must have sine equal to one of these values. Now sine takes on each value twice during each period, so we have a little work left to find the values of . Take the last root for example. If we take the arcsine of 0.966329 we get 1.31056, and = 1.31056 is not a solution to our equation. But arcsin(y) returns only one possible solution to the equation sin(x) = y. In this case, = - 1.31056 is the solution we're looking for.

The full set of solutions for 0 < 2 are

In the example above our polynomial in s had four real roots in [-1, 1]. In general we could have roots outside this interval, including complex roots. If we're looking for solutions with real values of then we discard these.

Example 2

Now suppose we want to solve

Our resultant is

and the roots are 0.119029, 0.987302, and -0.766973 0.319513i.

If we're only interested in real values of then the two solutions are arcsin(0.119029) = 0.119312 and arcsin(0.987302) = 1.41127. But there are two complex solutions, = 3.91711 0.433731i.

The post Solving quadratic trig equations first appeared on John D. Cook.