<div dir="auto">Dear Steven<div dir="auto"><br></div><div dir="auto">I guess it is because the spectra is mirrored at the nyquist frequency and you stay with half the spectra....</div><div dir="auto"><br></div><div dir="auto">Best</div><div dir="auto">Stephan</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El mar., 26 nov. 2019 22:03, Hoffman, Steven <<a href="mailto:steven.hoffman5@montana.edu">steven.hoffman5@montana.edu</a>> escribió:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">




<div dir="ltr">
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
 <br>
</div>
<div>
<div style="font-family:Calibri,Arial,Helvetica,sans-serif;font-size:12pt;color:rgb(0,0,0)">
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
To whom may be able to help me out,</div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
<br>
</div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
We were trying to recover the actual amplitude units for our power spectra, and I noticed something interesting in the output from FieldTrip. I ran a simulation where I take a 20Hz sine wave of +/- 5 "Volts", and calculate the power on it (Using the mtfft method).
 What we expected the amplitude to be was 25 (since power is the V^2 /Hz). However FieldTrip gave a value closer to about 12 V for the amplitude. </div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
<br>
</div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
I then calculated the power spectra of the simulated +/- 5 "Volt" sine wave using the Matlab built in "periodogram.m" function. This too yielded an amplitude not of 25, but of 12.5 . </div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
<br>
</div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
After some digging I found some documentation in the matlab that said that the units of the power spectra were actually (Amp ^2 / 2) . So in this case that makes sense as (5^2/2 = 12.5). </div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
<br>
</div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
My question is why the divide by two? I don't think I've ever seen amplitude expressed this way, its always V^2 / Hz ...</div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
<br>
</div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
Attached are figures of my matlab simulations (One for Field Trip, and for matlab's periodogram.m), and the snippet from the Matlab documentation for periodogram.m where it talks about amplitude being Amp^2 / 2. </div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
<br>
</div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
<br>
</div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
Thanks much for any help </div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
<br>
</div>
<div style="margin:0px;font-size:12pt;font-family:Calibri,Helvetica,sans-serif;background-color:rgb(255,255,255)">
-Steven Hoffman</div>
</div>
</div>
</div>

_______________________________________________<br>
fieldtrip mailing list<br>
<a href="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip" rel="noreferrer noreferrer" target="_blank">https://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
<a href="https://doi.org/10.1371/journal.pcbi.1002202" rel="noreferrer noreferrer" target="_blank">https://doi.org/10.1371/journal.pcbi.1002202</a><br>
</blockquote></div>