数学建模社区-数学中国
标题:
我是新手，MATLAB接触不久。有一题不会，求教求教
[打印本页]
作者:
玖玖
时间:
2013-7-15 16:44
标题:
我是新手，MATLAB接触不久。有一题不会，求教求教
对于如下线性规划问题（有3n个决策变量（x,r,s）和2n个约束）
$ t6 e/ p5 Y' y' I5 l
main(-xn)
% D' [; O/ Y* X) c% e( P. {
st
& O, T/ L% I2 M9 k
4x1-4r1=1;
: |% `" T: ]8 R( J
x1+s1=1;
& ]: H: J. G1 j+ J
4xj-x(j-1)-4rj=0;(j=2,3,4......)
# i0 _1 V+ O/ M+ g0 P2 j
4xj+x(j-1)+4sj=0;(j=2,3,4......)
6 W, n1 N# w4 }6 L4 J
xj,rj,sj>=0;(j=2,3,4......)
, |5 u4 [( o- N6 r5 E# S% Y
. w- H/ |; u! Y- N# I
请分别对n的不同取值（如n=2,10,50等），利用MATLAB软件求解上述规划。
' J+ W% j. p) f! S/ V
其中j，n是下标 。
作者:
玖玖
时间:
2013-7-15 20:29
有大神在不？帮忙啊！！！
作者:
玖玖
时间:
2013-7-16 14:54
好吧，经过不断地努力我自己做出来了。
作者:
玖玖
时间:
2013-7-16 14:57
syms n
: v$ q; D7 Y: Y, l
n=input('input n please.(press enter)n=')
! d+ `" h1 ~: \7 H) D j
if(n==1)
; q7 d* s; m. _3 x" i
a=[4,-4,0;1,0,1];
! G+ [) c( Q2 g! o) k1 E
b=[1;1];
; p( `6 P( j8 L* P9 a. Y8 T2 x. d
v=zeros(1,3);
b' u2 [2 S1 o/ k& u
c=[-1;0;0];
% `4 m9 {) K& k# W! }8 W8 k
[x,f]=linprog(c,[],[],a,b,v)
( c2 ^ b$ Z1 d; n- N
else
[1 @! o: }5 ^: t' Z' U& C" d
a=zeros(2*n,3*n);
( I& }& k( w& x4 f
for j=2:n
0 N# B6 }! Z! J, ]' B) R) c
x=zeros(1,n);
: X& ]! [: b, s$ \7 u. p
r=zeros(1,n);
$ K# _1 |4 |( O% i; {. d
s=zeros(1,n);
! Z/ ], Z- L9 N5 u1 z7 E
x(j)=4;
1 k& [1 w0 X/ [
x(j-1)=-1;
! J: q+ H4 u) |! n! ]
r(j)=-4;
$ m! N) K6 T' G5 u& E# o5 O
a((j+1),
=[x,r,s];
) ~* y Y8 E0 V
x(j)=4;
* Q4 {$ m, I1 k9 p! _. e/ v
x(j-1)=1;
! w2 Y3 x' O4 Z, ^: L
r(j)=0;
. I D* k, w& L& N
s(j)=4;
9 z2 a. J9 m1 h$ Y7 d( x( O8 t! b
a((j+n),
=[x,r,s];
7 o% Z% ?9 t3 h) c% I
end
4 d- w6 A! F/ f% W6 i7 f6 U
v=zeros(1,3*n);
. {4 i5 _- z+ e* X( P5 x' n
c1=-ones(n,1);
* N& ^) j4 V2 t
c2=zeros(2*n,1);
" z6 r( k1 x+ q% g
c=[c1;c2];
8 S8 f6 `6 C# ~) g9 L6 q Z! y4 E
x=zeros(1,n);
1 @: e9 h5 U8 B/ J0 i, f
r=zeros(1,n);
' A( \- k6 i' R; |, ?
s=zeros(1,n);
4 ?# x* H. f- K
x(1)=4;
$ I! a: n4 a4 {0 x. B( E4 A+ O
r(1)=-4;
$ u( O% \% P% G
a(1,
=[x,r,s];
' h1 I; I! a% n5 @$ I- W
x(1)=1;
7 A$ x$ T5 }. M# p ~- ?
r(1)=0;
" G0 E/ H* n* w8 V( a9 q
s(1)=1;
. _' ]# |+ `+ y) |, x9 Z
a(2,
=[x,r,s];
9 D7 a* ~1 k! A" J* p7 j' s# v
b=[1;1;zeros(n-1,1);4*eye(n-1,1)];
2 T$ [6 m4 z% L1 z% Y! t/ V: d
[x,f]=linprog(c,[],[],a,b,v)
- y8 S' V" l7 u) M; s
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];) ~* y Y8 E0 V=[x,r,s];7 o% Z% ?9 t3 h) c% I=[x,r,s];' h1 I; I! a% n5 @$ I- W=[x,r,s];