MathGroup Archive 2013

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Trig identity for linear combination of sines

  • To: mathgroup at smc.vnet.net
  • Subject: [mg131514] Re: Trig identity for linear combination of sines
  • From: "Kevin J. McCann" <kjm at KevinMcCann.com>
  • Date: Fri, 16 Aug 2013 07:07:20 -0400 (EDT)
  • Delivered-to: l-mathgroup@mail-archive0.wolfram.com
  • Delivered-to: l-mathgroup@wolfram.com
  • Delivered-to: mathgroup-outx@smc.vnet.net
  • Delivered-to: mathgroup-newsendx@smc.vnet.net
  • References: <h13vr9$8qh$1@smc.vnet.net> <kuho3l$l2d$1@smc.vnet.net>

Expand the sin(x+alpha)=sin(x)cos(alpha)+cos(x)sin(alpha) terms in the 
first expression. Find the coefficients of sin(x) and cos(x). Then 
equate these to the corresponding expressions in the second expansion.
Solve for A and z. It is not a simple result.

Here is the code:

rule1 = Sin[m_ + n_] -> Sin[m] Cos[n] + Cos[m] Sin[n]
expr = a Sin[x + q] + b Sin[x + r] + c Sin[x + t] /. rule1 // Expand
eq1 = Coefficient[expr, Sin[x]] == A Cos[z]
eq2 = Coefficient[expr, Cos[x]] == A Sin[z]
result=Solve[{eq1, eq2}, {A, z}] // Simplify

Note that there is a ConditionalExpression that results from the 
periodicity of the trig functions.

You can check it out by choosing some numbers. In addition, I set the 
constant C[1] in the result to zero:

a = 12.7; b = 2.2; c = 4.7;
q = 85.4 .1; r = 27.9; t = 16.14;
x = 75.031;
{a Sin[x + q] + b Sin[x + r] + c Sin[x + t],A Sin[x + z] /. result /. 
C[1] -> 0}

A Sin[x + z] /. result /. C[1] -> 0

These give the same answers.

Kevin

On 8/15/2013 1:14 AM, orenpedatzur at gmail.com wrote:
> On Monday, June 15, 2009 2:07:21 AM UTC+3, Chelly wrote:
>> Hi:
>>
>> I have a function which is a linear combination of sines, of the form
>>
>> y = a*sin(x+q) + b*sin(x+r) + c*sin(x+t)
>>
>> and I am looking for a trignometric identity that would reduce the above to the form
>>
>> y = A*sin(x+Z)
>>
>> I would also like to find out how to implement this rule in Mathematica.
>>
>> Thanks
>> Chelly
>
> Hi Chelly
>
> This kind of linear combination of sines is covered in wikipedia article on trigonometric identities:
>
> http://en.wikipedia.org/wiki/List_of_trigonometric_identities
>
> I would also like to know how to implement this in mathematica.
>
> bye!
>
> Oren P
>



  • Prev by Date: Re: question related to (-1)^(1/3)
  • Next by Date: Re: fyi, making animation into PDF file from images genenerated by Mathematica
  • Previous by thread: Re: Trig identity for linear combination of sines
  • Next by thread: Application of decision trees to waveform recognition and GDP per