全国大学生数学建模竞赛2021A1

[fj1,title1]=xlsread("附件1.csv");             %节点
fj1=fj1.*ones(2226,3);
fj1=fj1'
f=inf;
R=300.4;
F=0.466*R;
B2=0;
for B=290:0.01:310
    e=0;
    a=fj1(1,:).^2+fj1(2,:).^2;
    b=4*(B-R+F);
    c=-B-fj1(3,:);
    d=abs(a./b+c);
    e=sum(d);
    
    if e<f
        f=e;
        B2=B;
    end
end

B2

i=0:0.01:250;
f=inf;
d=0;
j=[];
k2=0;
G=0;
for B=300:0.01:301
    a=-(R^2-i.^2).^0.5;
    b=i.^2./(4.*(B-R+F))-B;
    c=abs(a-b);
    d=sum(c);
    j=[j d];
    p=0;
    k=0;
    
    for o=c
        if o<0.6
            k=k+1;
        end
    end
   if k>k2
       k2=k;
       G=B;
   end
    if d<f
        f=d;
        g=B;
    end
end
g