matlab 递归限制,matlab中超过最大递归数，哥哥姐姐帮忙看一下吧。

wp=0.4*pi;

ws=0.6*pi;

Rp=0.45;

As=80;

%给定指标

delta1=(10^(Rp/20)-1)/(10^(Rp/20)+1);

delta2=(1+delta1)*(10^( -As/20));

%求波动指数

weights=[delta2/delta1 1];

deltaf=(ws-wp)/(2*pi);

%给定权函数和△f=wp-ws

N=ceil((-20*log10(sqrt(delta1*delta2))-13)/(14.6*deltaf)+1);

N=N+mod(N-1,2);

%估算阶数N

f=[0 wp/pi ws/pi 1]; A=[1 1 0 0];

%给定频率点和希望幅度值

h=remez(N-1,f,A,weights);

%求冲击响应

[db,mag,pha,grd,w]=freqz_m(h,[1]);

delta_w=2* pi/1000; wsi= ws/delta_w+1;

wpi=wp/delta_w+1;

Asd=-max(db(wsi:1:500));

%求阻带衰减

sbplot(2,2,1); n=0:1:N-1;stem(n,h);

axis([0,52,-0.1,0.3;title('脉冲响应');

xlabel('n');

ylabel('hd(n)')

grid on;