Sat 01 Jun 2013 01:35:35 PM UTC, original submission:
grpdelay fails if the coefficient vectors are columns.
This script:
version
warning ("error", "Octave:broadcast");
DR=[1.00000 -0.00000 -3.37219 0.00000 ...
5.45710 -0.00000 -5.24394 0.00000 ...
3.12049 -0.00000 -1.08770 0.00000 0.17404];
N=[-0.0139469 -0.0222376 0.0178631 0.0451737 ...
0.0013962 -0.0259712 0.0016338 0.0165189 ...
0.0115098 0.0095051 0.0043874];
[T,w]=grpdelay(N,DR,1024);
plot(w/pi,T)
T=[];
warning ("off", "Octave:broadcast");
[T,w]=grpdelay(N',DR',1024);
plot(w/pi,T)
produces:
octave:1> version
ans = 3.6.4
octave:2> warning ("error", "Octave:broadcast");
octave:3> DR=[1.00000 -0.00000 -3.37219 0.00000 ...
> 5.45710 -0.00000 -5.24394 0.00000 ...
> 3.12049 -0.00000 -1.08770 0.00000 0.17404];
octave:4> N=[-0.0139469 -0.0222376 0.0178631 0.0451737 ...
> 0.0013962 -0.0259712 0.0016338 0.0165189 ...
> 0.0115098 0.0095051 0.0043874];
octave:5> [T,w]=grpdelay(N,DR,1024);
octave:6> plot(w/pi,T)
octave:7> T=[];
octave:8> warning ("off", "Octave:broadcast");
octave:9> [T,w]=grpdelay(N',DR',1024);
warning: grpdelay: setting group delay to 0 at singularity
octave:10> plot(w/pi,T)
octave:11>
|