📄 minimax_gauss.m
字号:
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
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -