数学建模社区-数学中国
标题:
我是新手，MATLAB接触不久。有一题不会，求教求教
[打印本页]
作者:
玖玖
时间:
2013-7-15 16:44
标题:
我是新手，MATLAB接触不久。有一题不会，求教求教
对于如下线性规划问题（有3n个决策变量（x,r,s）和2n个约束）
8 |+ {( h- g' K# I
main(-xn)
% M. R. E4 K# ]& |8 j
st
2 D' ^$ q @- |+ N
4x1-4r1=1;
# ^- B8 {$ C3 X. @) i
x1+s1=1;
) L& C% r. X/ K6 Z X, K5 Q
4xj-x(j-1)-4rj=0;(j=2,3,4......)
& D; w ]0 `) T8 J" @! X2 _$ u
4xj+x(j-1)+4sj=0;(j=2,3,4......)
c& H. V. [8 }! q
xj,rj,sj>=0;(j=2,3,4......)
* T' u/ D4 D2 V' U
0 r5 \0 b- P5 l4 C
请分别对n的不同取值（如n=2,10,50等），利用MATLAB软件求解上述规划。
( _ k+ D# j! O5 L8 O# r
其中j，n是下标 。
作者:
玖玖
时间:
2013-7-15 20:29
有大神在不？帮忙啊！！！
作者:
玖玖
时间:
2013-7-16 14:54
好吧，经过不断地努力我自己做出来了。
作者:
玖玖
时间:
2013-7-16 14:57
syms n
$ I; W5 s [! J5 o8 M% v( `
n=input('input n please.(press enter)n=')
2 i5 t5 \# V+ p' j3 ?9 n0 @: E
if(n==1)
, x6 x5 k, T# a( _, u( c) y, D$ N
a=[4,-4,0;1,0,1];
9 r5 i' R2 b0 D" K2 W# O
b=[1;1];
. O: Q2 R& |! z
v=zeros(1,3);
. R- [ _7 S. j% D5 q+ D5 E
c=[-1;0;0];
9 Z0 c! [! y& R& Y, f( }" _+ N; |
[x,f]=linprog(c,[],[],a,b,v)
& R. u5 o3 U4 D( D0 G( v' I0 z6 g
else
. `7 }! W* l9 a( H) P
a=zeros(2*n,3*n);
1 \: U" |% u y6 Y! t* F+ A
for j=2:n
9 a- X1 }3 q3 l& K1 ~" E+ Z6 n/ x
x=zeros(1,n);
, T* j) Y# Y5 J+ `* P
r=zeros(1,n);
7 W4 x; g$ @+ L- \1 T
s=zeros(1,n);
1 \. h" } [+ c* y1 O% d
x(j)=4;
: E3 v) F1 {4 m5 T% A9 V, x
x(j-1)=-1;
0 v- ~, q% S' a2 H+ J
r(j)=-4;
, [, w3 ]3 w7 a4 {8 w' ?- C
a((j+1),
=[x,r,s];
2 Q5 {, F* K4 J B& I m
x(j)=4;
' G0 i& l! y+ n5 r8 S$ h
x(j-1)=1;
& X" g2 w: ]2 z9 j. Y# l7 l
r(j)=0;
* S' [) }; k2 w
s(j)=4;
8 R$ x3 {" ?2 c
a((j+n),
=[x,r,s];
' d$ |( t4 r$ d, c# _
end
9 V" F. ~! r* P- a3 v3 F/ s- p) r
v=zeros(1,3*n);
9 J& I% o! A" \5 w D
c1=-ones(n,1);
2 |$ Q6 `9 u, M0 ?# A2 P2 `8 @
c2=zeros(2*n,1);
# L& U$ L1 e H S: m) f1 J
c=[c1;c2];
+ P% |( E/ o: `0 o
x=zeros(1,n);
, U* R- C& l1 f2 ~
r=zeros(1,n);
% g4 _0 c. X7 C+ ~" x6 Q0 _
s=zeros(1,n);
5 }' a, h) ]! y, F
x(1)=4;
: j8 k8 V% R3 f( F0 _3 }& u! m; y
r(1)=-4;
/ b# V* K2 p- p( D3 B9 ?0 x: p- w
a(1,
=[x,r,s];
. ~! R" y1 a& L- q, e3 s
x(1)=1;
' Y4 f; P W9 m/ g% z Q
r(1)=0;
v) G1 @/ C, B, l) A0 F6 w$ n0 Y
s(1)=1;
! ]4 v: g- @9 u
a(2,
=[x,r,s];
, }8 Q+ |( a( \; l: S8 d
b=[1;1;zeros(n-1,1);4*eye(n-1,1)];
+ g& A% y- ?! l/ I; ]+ V
[x,f]=linprog(c,[],[],a,b,v)
$ C j& R. L- Z
end
作者:
江边一粟
时间:
2013-8-15 11:37
楼主，笑脸是什么？求解。。。昨晚就在做这题，到现在还没做出来
作者:
qianchuan11
时间:
2014-8-11 20:30
/wx ssssssssssssssssssss
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5
=[x,r,s];=[x,r,s];=[x,r,s];=[x,r,s];, }8 Q+ |( a( \; l: S8 d