数学建模社区-数学中国
标题:
我是新手，MATLAB接触不久。有一题不会，求教求教
[打印本页]
作者:
玖玖
时间:
2013-7-15 16:44
标题:
我是新手，MATLAB接触不久。有一题不会，求教求教
对于如下线性规划问题（有3n个决策变量（x,r,s）和2n个约束）
$ W8 O4 r! k7 @: X1 V4 O( I9 m. s
main(-xn)
. g" p+ ^, J9 e: F; O% P
st
7 L: y4 a# O! y$ I' ?' n a
4x1-4r1=1;
/ ~& T% y; O- s& k/ b, M
x1+s1=1;
! s. R' x9 l% Q+ E1 s- Z0 }9 r
4xj-x(j-1)-4rj=0;(j=2,3,4......)
9 @" w* S- |. K, T
4xj+x(j-1)+4sj=0;(j=2,3,4......)
! N$ G/ V1 P5 z& s) V5 b2 D) I
xj,rj,sj>=0;(j=2,3,4......)
" @3 K; z6 P0 o" m' c( S; @1 C
7 u0 h8 b: N$ N. Q; U4 `' v1 K
请分别对n的不同取值（如n=2,10,50等），利用MATLAB软件求解上述规划。
0 h! B. @, V' V* w
其中j，n是下标 。
作者:
玖玖
时间:
2013-7-15 20:29
有大神在不？帮忙啊！！！
作者:
玖玖
时间:
2013-7-16 14:54
好吧，经过不断地努力我自己做出来了。
作者:
玖玖
时间:
2013-7-16 14:57
syms n
' R! a/ [7 ^. V6 h% e. t2 [
n=input('input n please.(press enter)n=')
7 h5 Q9 A7 R# j- q( I
if(n==1)
) f6 y5 O+ m( S" z3 p4 c- Z
a=[4,-4,0;1,0,1];
. O* k2 K0 T; |% B/ {6 `8 Q
b=[1;1];
7 `* h& ^% s2 k9 Q; h" C* X
v=zeros(1,3);
, G* A5 H- v' ] ^
c=[-1;0;0];
{1 x2 B u" {9 u! q- H
[x,f]=linprog(c,[],[],a,b,v)
' U$ l7 T/ h, b: G3 [2 o% }
else
. z0 W1 V4 p* X# b
a=zeros(2*n,3*n);
/ h+ x3 f4 g6 F1 J8 a
for j=2:n
" V- R. }: J* \- _
x=zeros(1,n);
5 O: T4 u6 v) l% l( J* o
r=zeros(1,n);
- Z5 B' e2 E& c2 _1 k" O
s=zeros(1,n);
0 J* r0 v# g, ]/ `1 |
x(j)=4;
) b/ i$ a& \, J# I
x(j-1)=-1;
- S1 N9 s% i w1 L
r(j)=-4;
) P6 F, X: I8 m, V' \
a((j+1),
=[x,r,s];
+ \3 t& i; K# |
x(j)=4;
! E) b& d2 a1 V0 D
x(j-1)=1;
; D. Q* l' a8 G
r(j)=0;
- m% g4 h" M7 {: m% m. E
s(j)=4;
! l1 Q8 j9 U! f, _5 o" A) u
a((j+n),
=[x,r,s];
9 L, Y- R; H1 i* ^4 g7 q
end
9 O& y) }+ y& d: h
v=zeros(1,3*n);
4 g/ a/ ?% H) M
c1=-ones(n,1);
% N# X( a) Q+ I, R3 ?
c2=zeros(2*n,1);
1 b3 P) Z* d& l, [
c=[c1;c2];
: S: s" z# `- x% e
x=zeros(1,n);
" i) A$ K9 v/ J2 ^( R
r=zeros(1,n);
" I( @$ n$ H# h$ I
s=zeros(1,n);
* ^# W+ m1 S8 u6 L- z
x(1)=4;
; f8 u5 j" l4 k/ s
r(1)=-4;
/ d+ a, p9 O" ~. | O# }
a(1,
=[x,r,s];
& B3 R m7 Z5 X- O! r
x(1)=1;
6 t7 K) \# |% N
r(1)=0;
- Q6 y' ?* y) B& G6 M7 {
s(1)=1;
5 x+ @3 s! Y9 D6 F/ C
a(2,
=[x,r,s];
A" m2 m3 b* M- c9 j. j
b=[1;1;zeros(n-1,1);4*eye(n-1,1)];
, v) V. e7 n M% ^7 N. V
[x,f]=linprog(c,[],[],a,b,v)
6 \0 b' Z$ h( ^. Z. ~3 h
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];+ \3 t& i; K# |=[x,r,s];=[x,r,s];& B3 R m7 Z5 X- O! r=[x,r,s];