minimax_gauss.m

来自「无线移动信道的matlab源码 用于对无线通信信道的编程实现」· M 代码 · 共 34 行

clear all;
close all;
clc;
N=100;f_max=91;sigma_0_2=1;N_i=25;PLOT=1;
f_c=sqrt(log(2))*f_max;

kappa_c=f_max/f_c;

F_list=[];
save F_list F_list

tau_max=N_i/(2*kappa_c*f_c);
tau=linspace(0,0.1,N);
r_mm=sigma_0_2*exp(-(pi*f_c/sqrt(log(2))*tau).^2);
f_i_n=zeros(N_i,1);
c_i_n=sqrt(sigma_0_2)*sqrt(2/N_i)*ones(N_i,1);
 for ii=1:N_i
     f_i_n(ii)=f_c/sqrt(log(2))*erfinv((2*ii-1)/(2*N_i));
 end
   

save data r_mm tau N_i c_i_n tau_max PLOT
tic

xo=f_i_n;
   options = optimset('MaxIter',10);
   [x,fval]=fminimax('fun_Gauss1',xo,[],[],[],[],[],[],[],options);

load x

f_i_n=sort(abs(x));
toc