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（[8?12.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（x00khk）