a=[1.0?1.0?1.0?1.0?1.0?1.0?1.0?1.0];f=[1.0?2.0?3.0?4.0?5.0?6.0?7.0?8.0];
I=8;n1=200;Fs=25.0;fs=Fs*2*I;kf=12.5;
for?k=1:I
??for?r=1:（n1+I）
????m（kr）=a（k）*（1.0+0.5*cos（2*pi*f（k）/Fs*（r-1）））;
?????end
??end
[n0f0m0w]=remezord（[812.5][1?0][0.001?0.001]fs）;
b=remez（287f0m0w）;figure（1）plot（20*log10（abs（fft（b））））;grid?on;
for?r=1:36?
???for?k=1:I
?????h（kr）=b（（r-1）*I+k）;
??end?
end
??for?r=1:（n1+I）
????for?k=1:I
???????mk（k）=?m（kr）;
end
????mfft=fft（mk）;
for?k=1:I
?????x0（kr）=mfft（k）*exp（j*pi/（2*I）*（k-1））;
??end
end
?for?r=1:（n1+I）?
?for?k=1:I
??x00（k（2*r-1）=x0（kr）;
??x00（k2*r）=0;
??end
end


?for?k=1:I???
????for?r=1:（n1+I）?
??????x00k（r）=x00（kr）;
end
for?r=1:36
hk（r）=h（kr）;
end
y0=conv