<div dir="ltr">Hi Rafael<div><br></div><div>Thank you, I shall print and study. :-)</div><div><br></div><div>Cheers,</div><div>Claus</div></div><br><div class="gmail_quote"><div dir="ltr">On Fri, May 25, 2018 at 8:42 PM Rafael Guerra <<a href="mailto:jrafaelbguerra@hotmail.com">jrafaelbguerra@hotmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">







<div lang="EN-US" link="blue" vlink="purple">
<div class="m_-860483954252869115WordSection1">
<p class="MsoNormal"><font size="2" color="black" face="Calibri"><span style="font-size:11.0pt;color:black">Hi Claus,<u></u><u></u></span></font></p>
<p class="MsoNormal"><font size="2" color="black" face="Calibri"><span style="font-size:11.0pt;color:black"><u></u> <u></u></span></font></p>
<p class="MsoNormal"><font size="2" color="black" face="Calibri"><span style="font-size:11.0pt;color:black">I am not aware of such function. However, you can find simple code here below for both phase modulation and demodulation, which is straightforward to
 translate in Scilab:<u></u><u></u></span></font></p>
<p class="MsoNormal"><font size="2" color="black" face="Calibri"><span style="font-size:11.0pt;color:black"><a href="https://www.gaussianwaves.com/2017/06/phase-demodulation-using-hilbert-transform-application-of-analytic-signal/" target="_blank">https://www.gaussianwaves.com/2017/06/phase-demodulation-using-hilbert-transform-application-of-analytic-signal/</a><u></u><u></u></span></font></p>
<p class="MsoNormal"><font size="2" color="black" face="Calibri"><span style="font-size:11.0pt;color:black"><u></u> <u></u></span></font></p>
<p class="MsoNormal"><font size="2" color="black" face="Calibri"><span style="font-size:11.0pt;color:black">Note that the phase modulation is coded differently from you snippet below.<u></u><u></u></span></font></p>
<p class="MsoNormal"><font size="2" color="black" face="Calibri"><span style="font-size:11.0pt;color:black"><u></u> <u></u></span></font></p>
<p class="MsoNormal"><font size="2" color="black" face="Calibri"><span style="font-size:11.0pt;color:black">Regards,<u></u><u></u></span></font></p>
<p class="MsoNormal"><font size="2" color="black" face="Calibri"><span style="font-size:11.0pt;color:black">Rafael<u></u><u></u></span></font></p>
<p class="MsoNormal"><font size="2" color="black" face="Calibri"><span style="font-size:11.0pt;color:black"><u></u> <u></u></span></font></p>
<p class="MsoNormal"><b><font size="2" face="Calibri"><span style="font-size:11.0pt;font-weight:bold">From:</span></font></b><span> users [mailto:<a href="mailto:users-bounces@lists.scilab.org" target="_blank">users-bounces@lists.scilab.org</a>]
<b><span style="font-weight:bold">On Behalf Of </span></b>Claus Futtrup<br>
<b><span style="font-weight:bold">Sent:</span></b> Friday, May 25, 2018 7:17 PM<br>
<b><span style="font-weight:bold">To:</span></b> International users mailing list for Scilab. <<a href="mailto:users@lists.scilab.org" target="_blank">users@lists.scilab.org</a>><br>
<b><span style="font-weight:bold">Subject:</span></b> [Scilab-users] Simulating phase modulation<u></u><u></u></span></p>
<p class="MsoNormal"><font size="2" face="Calibri"><span style="font-size:11.0pt"><u></u> <u></u></span></font></p>
<div>
<p class="MsoNormal"><font size="2" face="Calibri"><span style="font-size:11.0pt">Hi there<u></u><u></u></span></font></p>
<div>
<p class="MsoNormal"><font size="2" face="Calibri"><span style="font-size:11.0pt"><u></u> <u></u></span></font></p>
</div>
<div>
<p class="MsoNormal"><font size="2" face="Calibri"><span style="font-size:11.0pt">In a loudspeaker the driver can move several
<span class="m_-860483954252869115gmail-ginger-module-highlighter-mistake-anim">millimeter</span> in an attempt to reproduce a low-frequency note. If the speaker also at the same time produce a higher tone, this second tone is phase modulated by the first one. This is a distortion
 of the original signal which I'd like to simulate / illustrate with some simple Scilab code, if possible.<u></u><u></u></span></font></p>
</div>
<div>
<p class="MsoNormal"><font size="2" face="Calibri"><span style="font-size:11.0pt"><u></u> <u></u></span></font></p>
</div>
<div>
<p class="MsoNormal"><font size="2" face="Calibri"><span style="font-size:11.0pt">In Matlab this can be simulated with
<span class="m_-860483954252869115gmail-ginger-module-highlighter-mistake-anim">pmmod</span>.<u></u><u></u></span></font></p>
</div>
<div>
<p class="MsoNormal"><font size="2" face="Calibri"><span style="font-size:11.0pt"><a href="https://matlabandsimulink.wordpress.com/2013/03/12/phase-modulation/" target="_blank">https://matlabandsimulink.wordpress.com/2013/03/12/phase-modulation/</a><u></u><u></u></span></font></p>
</div>
<div>
<p class="MsoNormal"><font size="2" face="Calibri"><span style="font-size:11.0pt"><u></u> <u></u></span></font></p>
</div>
<div>
<p class="MsoNormal"><font size="2" face="Calibri"><span style="font-size:11.0pt">Is there a similar function in Scilab? (name - please ?)<u></u><u></u></span></font></p>
</div>
<div>
<p class="MsoNormal"><font size="2" face="Calibri"><span style="font-size:11.0pt"><u></u> <u></u></span></font></p>
</div>
<div>
<p class="MsoNormal"><font size="2" face="Calibri"><span style="font-size:11.0pt">Here's the code I have written so far - this is the part that shows the input signal (the un-distorted signal):<u></u><u></u></span></font></p>
</div>
<div>
<p class="MsoNormal"><font size="2" face="Calibri"><span style="font-size:11.0pt"><u></u> <u></u></span></font></p>
</div>
<div>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">sample_rate</span></font><font color="#5c5c5c" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#5c5c5c">=</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">20000</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">;</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"><u></u><u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">t</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"> <font color="#5c5c5c"><span style="color:#5c5c5c">=</span></font> <font color="rosybrown"><span style="color:rosybrown">0</span></font><font color="#ffaa00"><span style="color:#ffaa00">:</span></font><font color="rosybrown"><span style="color:rosybrown">1</span></font><font color="#5c5c5c"><span style="color:#5c5c5c">/</span></font><font color="black"><span style="color:black">sample_rate</span></font><font color="#ffaa00"><span style="color:#ffaa00">:</span></font><font color="rosybrown"><span style="color:rosybrown">0.6</span></font><font color="black"><span style="color:black">;</span></font><u></u><u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">N</span></font><font color="#5c5c5c" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#5c5c5c">=</span></font><font color="#32b9b9" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#32b9b9">size</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">t,</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">'*'</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">)</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">;</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"> <i><font color="#64ae64"><span style="color:#64ae64;font-style:italic">//number of samples</span></font></i><u></u><u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">y1</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"> <font color="#5c5c5c"><span style="color:#5c5c5c">=</span></font> <font color="#32b9b9"><span style="color:#32b9b9">sin</span></font><font color="#4a55db"><span style="color:#4a55db">(</span></font><font color="rosybrown"><span style="color:rosybrown">2</span></font><font color="#5c5c5c"><span style="color:#5c5c5c">*</span></font><font color="orchid"><span style="color:orchid">%pi</span></font><font color="#5c5c5c"><span style="color:#5c5c5c">*</span></font><font color="rosybrown"><span style="color:rosybrown">50</span></font><font color="#5c5c5c"><span style="color:#5c5c5c">*</span></font><font color="black"><span style="color:black">t</span></font><font color="#4a55db"><span style="color:#4a55db">)</span></font><font color="black"><span style="color:black">;</span></font><u></u><u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">y2</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"> <font color="#5c5c5c"><span style="color:#5c5c5c">=</span></font> <font color="rosybrown"><span style="color:rosybrown">0.5</span></font><font color="#5c5c5c"><span style="color:#5c5c5c">*</span></font><font color="#32b9b9"><span style="color:#32b9b9">sin</span></font><font color="#4a55db"><span style="color:#4a55db">(</span></font><font color="rosybrown"><span style="color:rosybrown">2</span></font><font color="#5c5c5c"><span style="color:#5c5c5c">*</span></font><font color="orchid"><span style="color:orchid">%pi</span></font><font color="#5c5c5c"><span style="color:#5c5c5c">*</span></font><font color="rosybrown"><span style="color:rosybrown">500</span></font><font color="#5c5c5c"><span style="color:#5c5c5c">*</span></font><font color="black"><span style="color:black">t</span></font><font color="#4a55db"><span style="color:#4a55db">)</span></font><font color="black"><span style="color:black">;</span></font><u></u><u></u></span></font></pre>
<pre><i><font size="2" color="#64ae64" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:#64ae64;font-style:italic">// y2 = 0.5*sin(2*%pi*500*t+%pi/4);</span></font></i><font face="Monospaced"><span style="font-family:"Monospaced",serif"><u></u><u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">s</span></font><font color="#5c5c5c" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#5c5c5c">=</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">y1</span></font><font color="#5c5c5c" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#5c5c5c">+</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">y2</span></font><font color="#5c5c5c" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#5c5c5c">+</span></font><font color="#32b9b9" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#32b9b9">grand</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">1</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">,N,</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">'nor'</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">,</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">0</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">,</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">1</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">)</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">;</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"><u></u><u></u></span></font></pre>
<pre><font size="2" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif"><u></u> <u></u></span></font></pre>
<pre><i><font size="2" color="#64ae64" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:#64ae64;font-style:italic">// Plot time-domain</span></font></i><font face="Monospaced"><span style="font-family:"Monospaced",serif"><u></u><u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">endplot</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"> <font color="#5c5c5c"><span style="color:#5c5c5c">=</span></font> <font color="#32b9b9"><span style="color:#32b9b9">round</span></font><font color="#4a55db"><span style="color:#4a55db">(</span></font><font color="black"><span style="color:black">N</span></font><font color="#5c5c5c"><span style="color:#5c5c5c">/</span></font><font color="rosybrown"><span style="color:rosybrown">15</span></font><font color="#4a55db"><span style="color:#4a55db">)</span></font><font color="black"><span style="color:black">;</span></font><u></u><u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">twoplots</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"> <font color="#5c5c5c"><span style="color:#5c5c5c">=</span></font> <font color="black"><span style="color:black">scf</span></font><font color="#4a55db"><span style="color:#4a55db">()</span></font><font color="black"><span style="color:black">;</span></font> <i><font color="#64ae64"><span style="color:#64ae64;font-style:italic">// Set Current Figure (Graphics Window)</span></font></i><u></u><u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">subplot</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">211</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">)</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">;</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"><u></u><u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">plot</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">t</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">1</span></font><font color="#ffaa00" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#ffaa00">:</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">endplot</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">)</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">,y1</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">1</span></font><font color="#ffaa00" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#ffaa00">:</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">endplot</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">)</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">,t</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">1</span></font><font color="#ffaa00" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#ffaa00">:</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">endplot</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">)</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">,y2</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">1</span></font><font color="#ffaa00" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#ffaa00">:</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">endplot</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">))</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">;</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"><u></u><u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">subplot</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">212</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">)</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">;</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"><u></u><u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">plot</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">t</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">1</span></font><font color="#ffaa00" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#ffaa00">:</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">endplot</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">)</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">,y1</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">1</span></font><font color="#ffaa00" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#ffaa00">:</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">endplot</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">)</span></font><font color="#5c5c5c" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#5c5c5c">+</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">y2</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">1</span></font><font color="#ffaa00" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#ffaa00">:</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">endplot</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">))</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">;</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"><u></u><u></u></span></font></pre>
<pre><font size="2" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif"><u></u> <u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">y</span></font><font color="#5c5c5c" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#5c5c5c">=</span></font><font color="#32b9b9" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#32b9b9">fft</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">s</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">)</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">;</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"><u></u><u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">ymax</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"> <font color="#5c5c5c"><span style="color:#5c5c5c">=</span></font> <font color="#32b9b9"><span style="color:#32b9b9">max</span></font><font color="#4a55db"><span style="color:#4a55db">(</span></font><font color="#32b9b9"><span style="color:#32b9b9">abs</span></font><font color="#4a55db"><span style="color:#4a55db">(</span></font><font color="black"><span style="color:black">y</span></font><font color="#4a55db"><span style="color:#4a55db">))</span></font><font color="black"><span style="color:black">;</span></font><u></u><u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">y</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"> <font color="#5c5c5c"><span style="color:#5c5c5c">=</span></font> <font color="black"><span style="color:black">y</span></font> <font color="#5c5c5c"><span style="color:#5c5c5c">./</span></font> <font color="black"><span style="color:black">ymax;</span></font> <i><font color="#64ae64"><span style="color:#64ae64;font-style:italic">// Normalize</span></font></i><u></u><u></u></span></font></pre>
<pre><font size="2" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif"><u></u> <u></u></span></font></pre>
<pre><i><font size="2" color="#64ae64" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:#64ae64;font-style:italic">// s is real so the fft response is conjugate symmetric</span></font></i><font face="Monospaced"><span style="font-family:"Monospaced",serif"><u></u><u></u></span></font></pre>
<pre><i><font size="2" color="#64ae64" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:#64ae64;font-style:italic">// and we retain only the first N/2 points</span></font></i><font face="Monospaced"><span style="font-family:"Monospaced",serif"><u></u><u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">f</span></font><font color="#5c5c5c" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#5c5c5c">=</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">sample_rate</span></font><font color="#5c5c5c" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#5c5c5c">*</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">0</span></font><font color="#ffaa00" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#ffaa00">:</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">N</span></font><font color="#5c5c5c" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#5c5c5c">/</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">2</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">))</span></font><font color="#5c5c5c" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#5c5c5c">/</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">N;</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"> <i><font color="#64ae64"><span style="color:#64ae64;font-style:italic">//associated frequency vector</span></font></i><u></u><u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">n</span></font><font color="#5c5c5c" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#5c5c5c">=</span></font><font color="#32b9b9" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#32b9b9">size</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">f,</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">'*'</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">)</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">;</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"><u></u><u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">fftplots</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"> <font color="#5c5c5c"><span style="color:#5c5c5c">=</span></font> <font color="black"><span style="color:black">scf</span></font><font color="#4a55db"><span style="color:#4a55db">()</span></font><font color="black"><span style="color:black">;</span></font><u></u><u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">plot</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">f</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">2</span></font><font color="#ffaa00" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#ffaa00">:$</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">)</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">,</span></font><font color="#32b9b9" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#32b9b9">abs</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">y</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">(</span></font><font color="rosybrown" face="Monospaced"><span style="font-family:"Monospaced",serif;color:rosybrown">2</span></font><font color="#ffaa00" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#ffaa00">:</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">n</span></font><font color="#4a55db" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#4a55db">)))</span></font><font color="black" face="Monospaced"><span style="font-family:"Monospaced",serif;color:black">;</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"> <i><font color="#64ae64"><span style="color:#64ae64;font-style:italic">// drop first datapoint, f = 0 (it prevents log-plot)</span></font></i><u></u><u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">a</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"> <font color="#5c5c5c"><span style="color:#5c5c5c">=</span></font> <font color="black"><span style="color:black">gca</span></font><font color="#4a55db"><span style="color:#4a55db">()</span></font><font color="black"><span style="color:black">;</span></font><u></u><u></u></span></font></pre>
<pre><font size="2" color="black" face="Monospaced"><span style="font-size:10.0pt;font-family:"Monospaced",serif;color:black">a</span></font><font color="#5c5c5c" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#5c5c5c">.</span></font><font color="#aaaaaa" face="Monospaced"><span style="font-family:"Monospaced",serif;color:#aaaaaa">log_flags</span></font><font face="Monospaced"><span style="font-family:"Monospaced",serif"> <font color="#5c5c5c"><span style="color:#5c5c5c">=</span></font> <font color="rosybrown"><span style="color:rosybrown">"lnn"</span></font><font color="black"><span style="color:black">;</span></font><u></u><u></u></span></font></pre>
<p class="MsoNormal"><font size="2" face="Calibri"><span style="font-size:11.0pt"><u></u> <u></u></span></font></p>
</div>
<div>
<p class="MsoNormal"><font size="2" face="Calibri"><span style="font-size:11.0pt">Best regards,<u></u><u></u></span></font></p>
</div>
<div>
<p class="MsoNormal"><font size="2" face="Calibri"><span style="font-size:11.0pt">Claus<u></u><u></u></span></font></p>
</div>
</div>
</div>
</div>

_______________________________________________<br>
users mailing list<br>
<a href="mailto:users@lists.scilab.org" target="_blank">users@lists.scilab.org</a><br>
<a href="http://lists.scilab.org/mailman/listinfo/users" rel="noreferrer" target="_blank">http://lists.scilab.org/mailman/listinfo/users</a><br>
</blockquote></div>