数学建模社区-数学中国
标题:
我是新手，MATLAB接触不久。有一题不会，求教求教
[打印本页]
作者:
玖玖
时间:
2013-7-15 16:44
标题:
我是新手，MATLAB接触不久。有一题不会，求教求教
对于如下线性规划问题（有3n个决策变量（x,r,s）和2n个约束）
|8 w/ B; B8 ?' k( r" j
main(-xn)
5 h2 Z: h. Z* i3 k
st
! Q) }! e& S4 ]+ g0 ]/ N0 H% f$ Y
4x1-4r1=1;
, R, A! H/ Q, ^& g
x1+s1=1;
$ K: K$ U, n; P/ ]# F
4xj-x(j-1)-4rj=0;(j=2,3,4......)
3 a* H% k9 J' j' _1 p6 L9 B
4xj+x(j-1)+4sj=0;(j=2,3,4......)
, v! {6 F& V' o6 n
xj,rj,sj>=0;(j=2,3,4......)
1 ]3 d$ P+ b5 ~) f6 u, x
1 q! s! i) n" t2 W# \5 J" Q
请分别对n的不同取值（如n=2,10,50等），利用MATLAB软件求解上述规划。
) ~0 c3 y7 \ @/ f! Y' I: i
其中j，n是下标 。
作者:
玖玖
时间:
2013-7-15 20:29
有大神在不？帮忙啊！！！
作者:
玖玖
时间:
2013-7-16 14:54
好吧，经过不断地努力我自己做出来了。
作者:
玖玖
时间:
2013-7-16 14:57
syms n
9 `& Y! H4 P2 D8 q8 `
n=input('input n please.(press enter)n=')
) i* ?2 I7 `% f+ H7 W
if(n==1)
( g- ]" | b+ R7 \7 n5 s
a=[4,-4,0;1,0,1];
1 I$ B! b2 e+ M$ x1 _
b=[1;1];
- j2 O) E7 ]9 x' c2 @
v=zeros(1,3);
8 N0 r0 c) W) s. w
c=[-1;0;0];
# O- E3 Z/ @% Y, F% L- r
[x,f]=linprog(c,[],[],a,b,v)
1 v* k" B* D3 O5 E
else
+ z7 q; i6 V* B9 n) m
a=zeros(2*n,3*n);
7 ~8 t* a }: U' ?+ o
for j=2:n
3 U6 K( r3 m- n2 E3 Z, n$ G. c: n4 n4 m
x=zeros(1,n);
- Q$ ~' i' x& z* A% d5 C9 A& b
r=zeros(1,n);
& j4 z! n! [" ?% |- H3 L4 P; ^# B- K9 H
s=zeros(1,n);
9 L5 B2 ^+ X- t* j
x(j)=4;
- Q+ C' e1 L/ p. K7 Z V# x
x(j-1)=-1;
. R! ^. Z" a3 w5 U) Z& n
r(j)=-4;
3 H) N9 q. }! x3 C6 \- o4 U
a((j+1),
=[x,r,s];
S9 z! B0 {$ z' H
x(j)=4;
' b8 v) |9 \" a; U4 L
x(j-1)=1;
) R0 J/ j7 T$ S, ~0 K9 c1 K; ~
r(j)=0;
% N. N2 ]( l$ V
s(j)=4;
2 b1 y. ~/ O9 ]" \/ a0 ~" O
a((j+n),
=[x,r,s];
7 e% j& y! j0 f6 g
end
( ?3 O+ p& e# {: W; P2 H4 y, S- |
v=zeros(1,3*n);
# @) ]2 i$ ]. ~
c1=-ones(n,1);
! q4 A3 d. H2 v6 y- l/ c
c2=zeros(2*n,1);
% Z( j6 g I4 S& ?' u, r
c=[c1;c2];
4 `$ S/ `0 h4 _3 E8 ?, v
x=zeros(1,n);
9 X5 |, `6 ~! _! h1 i2 g& {) Y, s
r=zeros(1,n);
7 v J% h: s3 e( |$ \
s=zeros(1,n);
9 W8 x" t) g( h7 y% m
x(1)=4;
/ R- E3 u* C2 P& ]( s" B" p
r(1)=-4;
2 o; x. p8 V# E
a(1,
=[x,r,s];
' ]) T1 t1 S+ M9 F4 S' @" `
x(1)=1;
8 c& ^' H9 i, S' u
r(1)=0;
! S1 r: C x: ]
s(1)=1;
: o# B8 z- y" c
a(2,
=[x,r,s];
: ~/ H8 N/ d( g% S: Y- e
b=[1;1;zeros(n-1,1);4*eye(n-1,1)];
1 e5 v# F: [* G9 n. o) W
[x,f]=linprog(c,[],[],a,b,v)
; h( l5 j7 z6 N( o5 K: Q
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];7 e% j& y! j0 f6 g=[x,r,s];' ]) T1 t1 S+ M9 F4 S' @" `=[x,r,s];: ~/ H8 N/ d( g% S: Y- e