search.asv

用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.297710.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
    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 5 10];
b=[4 9 13];
for i=30:200
    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 max(e)<=0.04 && max(e2)<=0.01
        break
    end
    if d>max(e) && max(e2)<=0.01
        d=max(e);
        e1=e;
        N1=N;
        y1=y;
        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);
            n(j)=sum(y1(a(j):b(j)));
        end
    end
end
max_e=d;
N=N1
y1
%e=e1
E=[n;e3]