search.m

用matlab编写的线性规划之单纯型方法用来求解一些简单的线性规划问题

%ra =[43.5249 2.3400 18.1473 14.7889 8.2488 2.7854 1.0404 2.2395 0.3840
%1.7995 41.0597 13.4133 7.4060 3.2977 10.6992 6.2435 15.4154 6.0546 4.2460 1.6254 3.1913]
clear
ra=[43.5249 41.0597 15.4154]/100;

r=[2.3400
   18.1473
   14.7889
  %  8.2488
    2.7854    
    1.0404 
    2.2395 
    0.3840
    1.7995
   13.4133
    7.4060
    3.2977
   10.6992
    6.2435
    6.0546
    4.2460
    1.6254
    3.1913
    ];
r=r/100;
d=100;
a=[1 9 14];%;[1 5 10]
b=[8 13 17];%;[4 9 13]
for i=50:1000
    x=i*r;
    y=round(x);
    N=sum(y);
    
    e=abs((y/N-r))./r;
    %er=abs(y-r*N);
    for j=1:3    
        e2(j)=abs(sum(y(a(j):b(j)))/N-ra(j))/ra(j);
    end

    if d>max(e) && max(e2)<=0.01
        d=max(e);
        e1=e;
        N1=N;
        y1=y;
        x1=x;
        e4=abs((y1/N1-r))./r;
        er=abs(y1-r*N1);
        for j=1:3    
            e3(j)=abs(sum(y1(a(j):b(j)))/N1-ra(j))/ra(j);
            nn(j)=sum(y1(a(j):b(j)));
        end
    end
    if max(e)<=0.05 && max(e2)<=0.01
        break
    end
end
max_e=d
N=N1
n=y1
%x=x1
e=e1
E=[nn;e3]