软件无线电仿真MATLAB源码

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