矩阵奇异问题

lianmei1984

w0=2*pi*2000; wb=1;
[z,p,k]=cheb1ap(3,3); [b,a]=zp2tf(z,p,k);
[num,den]=lp2bs(b,a,w0,wb); bode(num,den);
w=linspace(1,3000,300)*2*pi; H=freqs(num,den,w);
plot(w/(2*pi),abs(H));grid
%设计带阻滤波器。见附件1
%当用双线性变换法离散时，
[b1,b2]=bilinear(num,den,8000);
Warning: Matrix is close to singular or badly scaled.
         Results may be inaccurate. RCOND = 6.980122e-041.
> In bilinear at 89
Warning: Matrix is close to singular or badly scaled.
         Results may be inaccurate. RCOND = 6.980122e-041.
> In bilinear at 90
请高手帮忙看看到底是什么地方错了阿？自己实在是想不明白阿

[ 本帖最后由 zhangnan3509 于 2009-4-9 10:22 编辑 ]

ChaChing

LZ这个好像已经不是编程问题了!
LZ好像是要设计一个在2000Hz频宽仅1Hz的带阻滤波器, 会不会太窄些?
太久没玩这块了, bilinear怎回事忘了, 将之移动版块!

ChaChing

操作错误, 应该移至信号处理方法