create_fis_noise_set_with_velocity.m

input=zeros(21*15,2);
output= zeros(21*15,1);
k=0;
for dist=0:0.5:10
    for vel=0:0.1:1.5
        k=k+1;
        current_var=[0.1;1.3]*0.4/(1+dist/5)*vel/1.5;
        %current_var=[0.1;1.3]*0.4;
        max_res=0;
        rests=zeros(1,20000);
        for i=1:20000
            res=normrnd([0;0],current_var);
            rests(i)=res(2);
        end
        r=sort(rests);
        input(k,:)=[dist+10,vel];
        output(k)=r(19800);
    end
end
output_full=[output;output];
input_full=[input;flip(input(:,1)-10),input(:,2)];
opt = genfisOptions('GridPartition');
opt.NumMembershipFunctions = [3,3];
opt.InputMembershipFunctionType = "trimf";
    fis = genfis(input_full,output_full,opt);
[in,out,rule] = getTunableSettings(fis);
fisout = tunefis(fis,[in;out],input_full,output_full,tunefisOptions("Method","anfis"));
fis_set={fisout};