[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Octave-bug-tracker] [bug #52604] Matlab invfreqz(h, w, 4, 5, [], 30) ex
From: |
Stephen |
Subject: |
[Octave-bug-tracker] [bug #52604] Matlab invfreqz(h, w, 4, 5, [], 30) example doesn't give sample output as octave |
Date: |
Wed, 6 Dec 2017 12:41:43 -0500 (EST) |
User-agent: |
Mozilla/5.0 (Windows NT 6.1; Win64; x64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/62.0.3202.94 Safari/537.36 |
Follow-up Comment #1, bug #52604 (project octave):
Also, invfreqz(h,w,3,3,wt,30) shown below doesn't return the same result as
matlab (see https://edoras.sdsu.edu/doc/matlab/toolbox/signal/spectop9.html)
>> [b,a] = butter(4,0.4)
b =
0.046583 0.186332 0.279497 0.186332 0.046583
a =
1.000000 -0.782095 0.679979 -0.182676 0.030119
>> [h,w] = freqz(b,a,64);
>> [b4,a4] = invfreqz(h,w,4,4)
b4 =
0.046583 0.186332 0.279497 0.186332 0.046583
a4 =
1.000000 -0.782095 0.679979 -0.182676 0.030119
>> wt = ones(size(w));
>> [b30,a30] = invfreqz(h,w,3,3,wt,30)
b30 =
0.046389 0.178546 0.244578 0.127586
a30 =
1.00000 -0.95016 0.73823 -0.20056
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?52604>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/
- [Octave-bug-tracker] [bug #52604] Matlab invfreqz(h, w, 4, 5, [], 30) example doesn't give sample output as octave, Stephen, 2017/12/06
- [Octave-bug-tracker] [bug #52604] Matlab invfreqz(h, w, 4, 5, [], 30) example doesn't give sample output as octave,
Stephen <=