I=8，Fs=25KHz时的8个调幅AM信号的信道化发射机仿真结果

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）