DSPRelated.com
Forums

orthogonality of complex (or trigonometric) sinusoids...

Started by fisico32 May 3, 2010
hello forum,

Two sinusoids of different frequency w1 and w2 are orthogonal if integral
of their product (one time the complex conj of the other) gives zero.
That is only true if the integral of integration is T, and and the period
of w1 and the period of w2 fit an integer number of times in T...

That means that two sinusoids are orthogonal only if their frequencies are
different and one an integer multiple of the other....
Correct?

In all other cases where the two frequencies are different and not multiple
of each other the integral is nonzero....or is there some interval of
integration over which they are actually orthogonal?

For instance, 
w1=(4pi) rad/s 
w2=(3.2pi) rad/s 

y(t)=sin(w1*t+phi) 
x(t)=sin(w2*t+theta)
Is the integral zero for certain limits of integration? How do I find the
limits? Clearly the phases theta and phi matter....

thanks
fisico32  
On May 3, 5:20&#4294967295;pm, "fisico32" <marcoscipioni1@n_o_s_p_a_m.gmail.com>
wrote:
> hello forum, > > Two sinusoids of different frequency w1 and w2 are orthogonal if integral > of their product (one time the complex conj of the other) gives zero. > That is only true if the integral of integration is T, and and the period > of w1 and the period of w2 fit an integer number of times in T... > > That means that two sinusoids are orthogonal only if their frequencies are > different and one an integer multiple of the other.... > Correct?
No. Both periods dividing into T does not mean one freq is an integer multiple of the other. Ex. f1=2Hz, f2=3Hz, T=1.
> > In all other cases where the two frequencies are different and not multiple > of each other the integral is nonzero....or is there some interval of > integration over which they are actually orthogonal? > > For instance, > w1=(4pi) rad/s > w2=(3.2pi) rad/s > > y(t)=sin(w1*t+phi) > x(t)=sin(w2*t+theta) > Is the integral zero for certain limits of integration? How do I find the > limits? Clearly the phases theta and phi matter.... > > thanks > fisico32 &#4294967295;
>On May 3, 5:20=A0pm, "fisico32" <marcoscipioni1@n_o_s_p_a_m.gmail.com> >wrote: >> hello forum, >> >> Two sinusoids of different frequency w1 and w2 are orthogonal if
integral
>> of their product (one time the complex conj of the other) gives zero. >> That is only true if the integral of integration is T, and and the
period
>> of w1 and the period of w2 fit an integer number of times in T... >> >> That means that two sinusoids are orthogonal only if their frequencies
ar=
>e >> different and one an integer multiple of the other.... >> Correct? > >No. Both periods dividing into T does not mean one freq is an integer >multiple of the other. > >Ex. f1=3D2Hz, f2=3D3Hz, T=3D1. > >> >> In all other cases where the two frequencies are different and not
multip=
>le >> of each other the integral is nonzero....or is there some interval of >> integration over which they are actually orthogonal? >> >> For instance, >> w1=3D(4pi) rad/s >> w2=3D(3.2pi) rad/s >> >> y(t)=3Dsin(w1*t+phi) >> x(t)=3Dsin(w2*t+theta) >> Is the integral zero for certain limits of integration? How do I find
the
>> limits? Clearly the phases theta and phi matter.... >> >> thanks >> fisico32 =A0 >
Hello Dirk, thanks for the reply. you are right. I guess I meant that if both frequencies have periods T1 and T2 that fit a different integer number or times in the intervalof integration T, then their frequencies w1 and w2 are automatically multiples of each other and the signals are orthogonal over T. My question is about orthogonality between two frequencies that are not multiple of each other. Over which interval are they orthogonal?
>
On 5/3/2010 5:40 PM, fisico32 wrote:
>> On May 3, 5:20=A0pm, "fisico32"<marcoscipioni1@n_o_s_p_a_m.gmail.com> >> wrote: >>> hello forum, >>> >>> Two sinusoids of different frequency w1 and w2 are orthogonal if > integral >>> of their product (one time the complex conj of the other) gives zero. >>> That is only true if the integral of integration is T, and and the > period >>> of w1 and the period of w2 fit an integer number of times in T... >>> >>> That means that two sinusoids are orthogonal only if their frequencies > ar= >> e >>> different and one an integer multiple of the other.... >>> Correct? >> >> No. Both periods dividing into T does not mean one freq is an integer >> multiple of the other. >> >> Ex. f1=3D2Hz, f2=3D3Hz, T=3D1. >> >>> >>> In all other cases where the two frequencies are different and not > multip= >> le >>> of each other the integral is nonzero....or is there some interval of >>> integration over which they are actually orthogonal? >>> >>> For instance, >>> w1=3D(4pi) rad/s >>> w2=3D(3.2pi) rad/s >>> >>> y(t)=3Dsin(w1*t+phi) >>> x(t)=3Dsin(w2*t+theta) >>> Is the integral zero for certain limits of integration? How do I find > the >>> limits? Clearly the phases theta and phi matter.... >>> >>> thanks >>> fisico32 =A0 >> > > Hello Dirk, > thanks for the reply. you are right. > I guess I meant that if both frequencies have periods T1 and T2 that fit a > different integer number or times in the intervalof integration T, then > their frequencies w1 and w2 are automatically multiples of each other and > the signals are orthogonal over T. > > My question is about orthogonality between two frequencies that are not > multiple of each other. Over which interval are they orthogonal?
Aside from Dirk's correction, still no. A sine and cosine of the same frequency are orthogonal. All sinusoids of different frequencies are orthogonal in that the time-average of their products goes asymptotically to zero as the period of integration becomes progressively longer. The orthogonality is the basis for the Fourier transform. The residual product sum for short integration times gives rise to what we call leakage. Jerry -- "I view the progress of science as ... the slow erosion of the tendency to dichotomize." --Barbara Smuts, U. Mich. &#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;&#4294967295;
On Mon, 03 May 2010 16:40:46 -0500, "fisico32"
<marcoscipioni1@n_o_s_p_a_m.gmail.com> wrote:

>I guess I meant that if both frequencies have periods T1 and T2 that fit a >different integer number or times in the intervalof integration T, then >their frequencies w1 and w2 are automatically multiples of each other
Nope. As Dirk mentioned, if one signal is 2 Hz and the other is 3 Hz, then the interval of integration T is 1 second. But 2 Hz is not a multiple of 3 Hz. They are both multiples of 1 Hz. That's the key: the integration period T is the reciprocal of the greatest common factor of the two frequencies. Example: f1 = 14 Hz, f2 = 16 Hz. 14 = 2 * 7 16 = 2 * 8 The largest common factor is 2. The interval T is 1/2 second. One-half second is 7 periods of f1 and 8 periods of f2. If the two frequencies have an irrational relationship, then T is infinite. Example: f1 = 2 Hz, f2 = pi Hz. There is no common factor.
>My question is about orthogonality between two frequencies that are not >multiple of each other. Over which interval are they orthogonal?
See above. Greg
On May 3, 5:20&#4294967295;pm, "fisico32" <marcoscipioni1@n_o_s_p_a_m.gmail.com>
wrote:
> hello forum, > > Two sinusoids of different frequency w1 and w2 are orthogonal if integral > of their product (one time the complex conj of the other) gives zero. > That is only true if the integral of integration is T, and and the period > of w1 and the period of w2 fit an integer number of times in T... > > That means that two sinusoids are orthogonal only if their frequencies are > different and one an integer multiple of the other.... > Correct? > > In all other cases where the two frequencies are different and not multiple > of each other the integral is nonzero....or is there some interval of > integration over which they are actually orthogonal? > > For instance, > w1=(4pi) rad/s > w2=(3.2pi) rad/s > > y(t)=sin(w1*t+phi) > x(t)=sin(w2*t+theta) > Is the integral zero for certain limits of integration? How do I find the > limits? Clearly the phases theta and phi matter.... > > thanks > fisico32 &#4294967295;
The phases don't matter. Basically for orthogonality find the ratio of w1 to w2. If it is irrational you need infinite limits. Otherwise express as ratio of integers in the form of a/b. Then for an arbitrary lower limit of "p", your upper limit of "q" is simply q=p+n*a*2*pi/w1 where n is any integer. you may also write q=p+m*b*2*pi/w2 where m is any integer For your example w1=4pi and w2=3.2pi thus a/b=1.25=5/4 q=p+n*5*(0.5)=p+2.5*n IHTH, Clay
fisico32 wrote:
> hello forum, > > Two sinusoids of different frequency w1 and w2 are orthogonal if integral > of their product (one time the complex conj of the other) gives zero.
The definition I have seen uses the _average_ of the integral of their product. This makes no difference for a finite interval -- more on that later.
> That is only true if the integral of integration is T, and and the period > of w1 and the period of w2 fit an integer number of times in T... > > That means that two sinusoids are orthogonal only if their frequencies are > different and one an integer multiple of the other.... > Correct?
Incorrect, even by your definition of orthogonality.
> In all other cases where the two frequencies are different and not multiple > of each other the integral is nonzero....or is there some interval of > integration over which they are actually orthogonal? > > For instance, > w1=(4pi) rad/s > w2=(3.2pi) rad/s > > y(t)=sin(w1*t+phi) > x(t)=sin(w2*t+theta) > Is the integral zero for certain limits of integration? How do I find the > limits? Clearly the phases theta and phi matter....
Take the integral of the _average_ of your two sinusoids over the interval -T/2 to +T/2, and take the limit of the average as T goes to infinity. What conditions must be met for the two sinusoids to _not_ be orthogonal? -- Tim Wescott Control system and signal processing consulting www.wescottdesign.com
fisico32 <marcoscipioni1@n_o_s_p_a_m.gmail.com> wrote:
 
> Two sinusoids of different frequency w1 and w2 are orthogonal if integral > of their product (one time the complex conj of the other) gives zero. > That is only true if the integral of integration is T, and and the period > of w1 and the period of w2 fit an integer number of times in T...
Look at the Wikipedia pages for Orthogonal_functions and Inner_product.
> That means that two sinusoids are orthogonal only if their > frequencies are different and one an integer multiple of the other.... > Correct?
Well, you can generalize a little more than that. You can have finite or infinite limits on the integral. Also, some inner products have a weight function in addition to the two functions contributing to the product. With infinite limits, sinusoids are orthogonal when the frequency is different. With finite limits, it is usual that only sinusoids with an integer number of periods over the range of integration are used. (Depending on boundary conditions, sometimes a half integer also works.) Another possibility is a Gaussian weighting function of a specified width. The integration limits would be infinite, though the main contribution comes from a finite width region. -- glen
Clay <clay@claysturner.com> wrote:
> On May 3, 5:20?pm, "fisico32" wrote:
>> Two sinusoids of different frequency w1 and w2 are orthogonal if integral >> of their product (one time the complex conj of the other) gives zero. >> That is only true if the integral of integration is T, and and the period >> of w1 and the period of w2 fit an integer number of times in T...
(snip)
> Basically for orthogonality find the ratio of w1 to w2. If it is > irrational you need infinite limits. Otherwise express as ratio of > integers in the form of a/b. > Then for an arbitrary lower limit of "p", your upper limit of "q" is > simply
> q=p+n*a*2*pi/w1 where n is any integer.
I suppose it could work that way, but it is usual for the limits to come first, along with the boundary conditions and any weighting function, and then find which basis functions satisfy the boundary conditions. Two common boundary conditions are for the function or its derivative to go to zero at the boundary, which allows for a whole number of periods over the region of interest. A little less common, is to have the function go to zero at one end, and the derivative at the other, which defines the modes for a cylindrical tube closed at one end, such that an odd number of half periods of sine or cosine fit. (And not counting end effects for real tubes.) Even more rare are mixed boundary conditions, allowing for some linear combination of the function and its derivative. And all that with the assumption that the basis functions are sinusoids. For a circular drum head with uniform tension the radial basis functions are Bessel functions, which have their own orthogonality relationship. -- glen
On May 3, 10:10&#4294967295;pm, glen herrmannsfeldt <g...@ugcs.caltech.edu> wrote:
> Clay <c...@claysturner.com> wrote: > > On May 3, 5:20?pm, "fisico32" wrote: > >> Two sinusoids of different frequency w1 and w2 are orthogonal if integral > >> of their product (one time the complex conj of the other) gives zero. > >> That is only true if the integral of integration is T, and and the period > >> of w1 and the period of w2 fit an integer number of times in T... > > (snip) > > > Basically for orthogonality find the ratio of w1 to w2. If it is > > irrational you need infinite limits. Otherwise express as ratio of > > integers in the form of a/b. > > Then for an arbitrary lower limit of "p", your upper limit of "q" is > > simply > > q=p+n*a*2*pi/w1 &#4294967295;where n is any integer. > > I suppose it could work that way, but it is usual for the limits > to come first, along with the boundary conditions and any > weighting function, and then find which basis functions satisfy > the boundary conditions. >
Glen, I constrained my response to the OP's particular example. Since the OP was having trouble with it, I didn't even want to venture into weight functions, Hilbert spaces, and Sturm-Liouville systems. 'cause we know one can delve deeply into orthogonality and its various generalizations. Clay .
> > -- glen