毕业设计matlab程序里面的一点问题！急求大哥大姐帮忙！

求助大神，拜托

第一种是N=41,W=0.40-0.48.这两个参数代入下面的程序，得到mes的结果图，并绘制图像。（程序正常运行了)第二种是取W=0.45时，N取分别取31 33 35 37 41 97 99 101 127时，两个参数代入下面的程序，得到mes的结果图，并绘制图像。
运行不了，得不出结果图。求帮帮我，毕设。

1. 第一种
   
2. w1=zeros(10);  %%%定义一个0数组
   
3. y=zeros(10);   %%%定义一个0数组
   
4. w1=[0.40 0.41 0.42 0.43 0.44 0.45 0.46 0.47 0.48 ]; %%赋值
   
5. N=41;
   
6. for i=1:1:9   
   
7.   w=w1(i);   %%每次循环都会将数组w1中的值赋给w，参与运算
   
8.   %下面的程序没动
   
9.   
   
10. [h0,h1,g0,g1]=firpr2chfb(N,w); [H1z,w]=freqz(h0,1,512);
    
11. H1_abs=abs(H1z);H1_db=20*log10(H1_abs); [H2z,w]=freqz(h1,1,512);
    
12. H2_abs=abs(H2z);H2_db=20*log10(H2_abs);
    
13. n=1:500; T=0.2;  x=sin(n*T);
    
14. hlp=mfilt.firdecim(2,h0); hhp=mfilt.firdecim(2,h1);
    
15. glp=mfilt.firinterp(2,g0); ghp=mfilt.firinterp(2,g1);
    
16. x0=filter(hlp,x); x0=filter(glp,x0);
    
17. x1=filter(hhp,x); x1=filter(ghp,x1);
    
18. xidle=x0+x1;  xshift=[zeros(1,N) x(1:end-N)];
    
19. e=xidle-xshift;
    
20. mes=sum(abs(e).^2)/length(e);
    
21. y(i)=mes; %每次运算的mes值都赋给相应的y（i）
    
22. end
    
23. figure(7);
    
24. i=1:1:9;
    
25. sprintf('%f',w1(1))
    
26. plot(w1(i),y(i),'--b');
    
27. 第二种
    
28. w1=zeros(10);  %%%定义一个0数组
    
29. y=zeros(10);   %%%定义一个0数组
    
30. w1=[31 33 35 37 41 97 99 101 127]; %%赋值
    
31. w=0.45;
    
32. for i=1:1:9   
    
33.   N=w1(i);   %%每次循环都会将数组w1中的值赋给w，参与运算
    
34.   %下面的程序没动
    
35.   
    
36. [h0,h1,g0,g1]=firpr2chfb(N,w); [H1z,w]=freqz(h0,1,512);
    
37. H1_abs=abs(H1z);H1_db=20*log10(H1_abs); [H2z,w]=freqz(h1,1,512);
    
38. H2_abs=abs(H2z);H2_db=20*log10(H2_abs);
    
39. n=1:500; T=0.2;  x=sin(n*T);
    
40. hlp=mfilt.firdecim(2,h0); hhp=mfilt.firdecim(2,h1);
    
41. glp=mfilt.firinterp(2,g0); ghp=mfilt.firinterp(2,g1);
    
42. x0=filter(hlp,x); x0=filter(glp,x0);
    
43. x1=filter(hhp,x); x1=filter(ghp,x1);
    
44. xidle=x0+x1;  xshift=[zeros(1,N) x(1:end-N)];
    
45. e=xidle-xshift;
    
46. mes=sum(abs(e).^2)/length(e);
    
47. y(i)=mes; %每次运算的mes值都赋给相应的y（i）
    
48. end
    
49. figure(7);
    
50. i=1:1:9;
    
51. sprintf('%f',w1(1))
    
52. plot(w1(i),y(i),'--b');
    

复制代码

求助大神，拜托

大哥大姐帮帮忙，关乎小弟能否顺利毕业!忙完还要准备二战呢，呜呜呜

thx1234

The length of the window specified must exceed the filter order by one.