一、实验内容
P206第六题
function f=wuyan2(c)
y=[3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.41 38.6 50.2 62.9 76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.0 226.5 251.4 281.4]
t=[0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 170 180 190 200 210] f=y-c(1)/(1+c(1)/3.9-1)*exp^(-c(2)*t)
c0=[1 1]
c=lsqnonlin('wuyan2',c0)
P206第七题
function f=wuyan1(c)
q=[0.4518 0.4862 0.5295 0.5934 0.7171 0.8964 1.0202 1.1963 1.4928 1.6909 1.8548 2.1618 2.6638 3.4634 4.6759 5.8478 6.7885 7.4463 7.8345 8.2068 8.9468 9.7315 10.5172 11.7390 13.6876 ];
k=[0.0911 0.0961 0.1230 0.1430 0.1860 0.2543 0.3121 0.3792 0.4754 0.4410 0.4517 0.5595 0.8080 1.3072 1.7042 2.0019 2.2914 2.4941 2.8406 2.9855 3.2918 3.7214 4.3500 5.5567 7.0477];
l=[4.2361 4.3725 4.5295 4.6436 4.8179 4.9873 5.1282 5.2783 5.4334 5.5329 6.4749
6.5491 6.6152 6.6808 6.7455 6.8065 6.8950 6.9820
7.0637 7.1394 7.2085 7.3025 7.3470
7.4432 7.5200];
f=q-c(1)*k.^c(2).*l.^c(3)
c0=[1 1 1]
c=lsqnonlin('wuyan1',c0)
c =
0.4091 0.6401 1.1446
a=0.4091 α=0.6401 β=1.1446
P239第五题
c=[-20 -30];
A=[1 2;5 4];
b=[20 70];
v1=[0 0];
[x,f,ef,out,lag]=linprog(c,A,b,[],[],v1)
z=-f
x =
10.0000
5.0000
z =
350.0000
生产甲10件生产乙5件收入最大350元
A原料+1
c=[-20 -30];
A=[1 2;5 4];
b=[21 70];
v1=[0 0];
[x,f,ef,out,lag]=linprog(c,A,b,[],[],v1)
z=-f
x =
9.3333
5.8333
z =
361.6667
生产甲 9.3333件生产乙5.8333件收入最大361.6667元多11.6667元B原料+1
c=[-20 -30];
A=[1 2;5 4];
b=[20 71];
v1=[0 0];
[x,f,ef,out,lag]=linprog(c,A,b,[],[],v1)
z=-f
x =
10.3333
4.8333
z =
351.6667
生产甲10.3333件生产乙4.8333件收入最大 351.6667元多1.6667元
若以6元买A ,2元买B,当全买A时利润最大
P239第十题
c=[-21,-11,-1,-31,-21,-11,-41,-31,-21];
a=[1,1,1,0,0,0,0,0,0;
0,0,0,1,1,1,0,0,0;
0,0,0,0,0,0,1,1,1;
-2,0,0,-4,0,0,2,0,0;
0,-4,0,0,2,0,0,0,0;
0,0,-6,0,0,0,0,0,-2;
-0.5,0,0,1,0,0,2,0,0;
0,-1.5,0,0,0,0,0,1,0;
0,0,-2.5,0,0,-1,0,0,0;
1,1,1,1,1,1,1,1,1]
b=[5000;5000;5000;0;0;0;0;0;0;14000];
l=[0 0 0 0 0 0 0 0 0]
aeq=[1,0,0,1,0,0,1,0,0;
0,1,0,0,1,0,0,1,0;
0,0,1,0,0,1,0,0,1;];
beq=[3000;2000;1000];
[x,r]=linprog(c,a,b,aeq,beq,l,[])
x =
1.0e+003 *
2.4000
0.8000
0.0000
0.0000
0.0000
0.0000
0.6000
1.2000
1.0000
r =
-1.4200e+005
所以在用A 原材料 2400 ，800, 0 分别来加工‘甲’，‘乙’‘丙’产品；
用B 原材料 0，0 ， 0 分别来加工‘甲’，‘乙’‘丙’产品
用C 原材料 600，1200， 1000分别来加工‘甲’，‘乙’‘丙’产品时该工厂的利润有最大值，最大值是1.4200e+005
加广告费后
c=[-21,-11,-1,-31,-21, -11,-41,-31,-21,1,1,1];
a=[1,1,1,0,0,0,0,0,0,0,0,0;
0,0,0,1,1,1,0,0,0,0,0,0;
0,0,0,0,0,0,1,1,1,0,0,0;
-2,0,0,-4,0,0,2,0,0,0,0,0;
0,-4,0,0,2,0,0,0,0,0,0,0;
0,0,-6,0,0,0,0,0,-2,0,0,0;
-0.5,0,0,1,0,0,2,0,0,0,0,0;
0,-1.5,0,0,0,0,0,1,0,0,0,0;
0,0,-2.5,0,0,-1,0,0,0,0,0,0;
1,1,1,1,1,1,1,1,1,0,0,0]
b=[5000;5000;5000;0;0;0;0;0;0;14000]
l=[0 0 0 0 0 0 0 0 0]
aeq=[1,0,0,1,0,0,1,0,0,-10,0,0;
0,1,0,0,1,0,0,1,0,0,-10,0;
0,0,1,0,0,1,0,0,1,0,0,-10;];
beq=[3000;2000;1000];
[x,r]=linprog(c,a,b,aeq,beq,l,[])
x =
1.0e+003 *
2.9308
2.0692
0.0000
1.1430
3.3479
0.0000
0.1612
3.1038
1.2441
0.1235
0.6521
0.0244
r =
-3.1820e+005
既在用A 原材料 2930.8，2069．2 ，0.0000 分别来加工‘甲’，‘乙’‘丙’产品；
用B 原材料 1143.0, 3347．9, 0.0000 分别来加工‘甲’，‘乙’‘丙’产品
用C 原材料 161.2 ，3103.8 ， 1244.1 分别来加工‘甲’，‘乙’‘丙’产品
且对与甲，乙，丙三种汽油的广告投入分别 123.5 ，652.1 ，24.4 时，该工厂的利润有最大值，最大值是3.1820e+005。