高阶欧拉法离散滤波器程序

lianmei1984

我做了一个模拟陷波器，要求用欧拉法和龙格库塔法进行离散为数字陷波器，高阶欧拉法离散后为什么成了窄带带通滤波器了?高阶欧拉法程序 在附件4中。
程序如下：
模拟陷波器
f=2000; num=[1 0 (2*pi*f)^2];
den=[1 pi*f/100 (2*pi*f)^2];
bode(num,den)
w=linspace(1,3000,300)*2*pi;
H=freqs(num,den,w);
plot(w/(2*pi),abs(H));grid
双线性变换法离散
[numd,dend]=bilinear(num,den,44100);
fs=44100;
[H,W]=freqz(numd,dend);
%绘制频率响应曲线
figure(1)
plot(W*fs/(2*pi),abs(H))
欧拉法离散
G=tf(num,den) ; %假设一个连续的传递函数
fs=44100;h=1/fs;
F1=highEulerfun(G,h);
G1=sym2tf(F1);% (get(G1);[nn1,dd1]=tfdata(G1,'v');)获取分子分母系数
num1=G1.num{1};den1=G1.den{1};
[H1,W1]=freqz(num1,den1);
%绘制频率响应曲线
figure(2)
plot(W1*fs/(2*pi),abs(H1))
请高手帮忙看看程序是哪里错了阿？万分感谢阿