H(z) = yulewalk(N,frq,mag) - finds the N-th order iir filter:<br><br><img title="%5Cnormalsize%5C%21H%28z%29%3D%5Cfrac%7BB%28z%29%7D%7BA%28z%29%7D%3D%5Cfrac%7Bb%281%29%20%2B%20b%282%29z%5E%7B-1%7D%2B...%2Bb%28n%29z%5E%7B-%28n-1%29%7D%7D%7B1%2Ba%281%29z%5E%7B-1%7D%2B...%2Ba%28n%29z%5E%7B-%28n-1%29%7D%7D.gif" alt="%5Cnormalsize%5C%21H%28z%29%3D%5Cfrac%7BB%28z%29%7D%7BA%28z%29%7D%3D%5Cfrac%7Bb%281%29%20%2B%20b%282%29z%5E%7B-1%7D%2B...%2Bb%28n%29z%5E%7B-%28n-1%29%7D%7D%7B1%2Ba%281%29z%5E%7B-1%7D%2B...%2Ba%28n%29z%5E%7B-%28n-1%29%7D%7D.gif" src="http://www.texify.com/img/%5Cnormalsize%5C%21H%28z%29%3D%5Cfrac%7BB%28z%29%7D%7BA%28z%29%7D%3D%5Cfrac%7Bb%281%29%20%2B%20b%282%29z%5E%7B-1%7D%2B...%2Bb%28n%29z%5E%7B-%28n-1%29%7D%7D%7B1%2Ba%281%29z%5E%7B-1%7D%2B...%2Ba%28n%29z%5E%7B-%28n-1%29%7D%7D.gif"><br>
<br><span style="font-family:arial,helvetica,sans-serif">H(z) : filter B(z)/A(z)</span><br style="font-family:arial,helvetica,sans-serif"><span style="font-family:arial,helvetica,sans-serif">N    : integer (order of desired filter)</span><br style="font-family:arial,helvetica,sans-serif">
<span style="font-family:arial,helvetica,sans-serif">frq  : real row vector (non-decreasing order), frequencies.</span><br style="font-family:arial,helvetica,sans-serif"><span style="font-family:arial,helvetica,sans-serif">mag  : non negative real row vector (same size as frq), desired magnitudes.</span><br style="font-family:courier new,monospace">
<br>----------------------------------------<br><br>yulewalk() function in &quot;/scipy/signal/filter_design.py&quot; is empty. It perhaps can be calculated using scipy.optimize.leastsq() but I lack skills to figure out how?<br>
<br>All I can do right now is port it quick and dirty from Matlab to Python, but I feel it&#39;s bad idea and waste of time<br><br>Can someone provide tip how can I calculate 10-th order yulewalk with scipy?<br><br>