create_fis_noise_set_with_velocity_double.m

input=zeros(21*15,2);
output= zeros(21*15,2);
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),r(18000)];
    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(:,1),opt);
[in,out,rule] = getTunableSettings(fis);
fisout = tunefis(fis,[in;out],input_full,output_full(:,1),tunefisOptions("Method","anfis"));



    fis = genfis(input_full,output_full(:,2),opt);
[in,out,rule] = getTunableSettings(fis);
fisout_main = tunefis(fis,[in;out],input_full,output_full(:,2),tunefisOptions("Method","anfis"));

fis_set={fisout_main,fisout};