找回密码
 立即注册
查看: 2740|回复: 2

有关线性规划函数lp的使用,matlab转r,转的执行结果不对,哪位大侠帮忙看看哈

[复制链接]
发表于 2015-12-21 13:37:20 | 显示全部楼层 |阅读模式


matlab程序如下:

c=[-0.4 -0.28 -0.32 -0.72 -0.64 -0.6];
          A=[0.01 0.01 0.01 0.03 0.03 0.03;0.02 0 0 0.05 0 0;0 0.02 0 0 0.05 0;0 0 0.03 0 0 0.08];
          b=[850;700;100;900];
          Aeq=[]; beq=[];
          vlb=[0;0;0;0;0;0]; vub=[];
[x,fval]=linprog(c,A,b,Aeq,beq,vlb,vub)

R程序如下:

c<-c(-0.4, -0.28, -0.32, -0.72, -0.64, -0.6);
A<-matrix(c(0.01, 0.01 ,0.01 ,0.03, 0.03, 0.03,0.02 ,0 ,0, 0.05, 0, 0,
         0, 0.02, 0, 0, 0.05, 0, 0, 0, 0.03, 0, 0, 0.08
         ),nrow = 4,byrow="F");
         
         
         
A1<- matrix(c(1,0,0,0,0,0,
         0,1,0,0,0,0,
         0,0,1,0,0,0,
         0,0,0,1,0,0,
         0,0,0,0,1,0,
         0,0,0,0,0,1),nrow=6,byrow="F") ;

Anew=rbind(A,A1);
library(lpSolve);
eg.lp<-lp(direction="max",objective.in=c,const.mat=Anew,const.rhs=b,const.dir=c(rep("<=",4),rep(">=",6)));
eg.lp$solution

回复

使用道具 举报

 楼主| 发表于 2015-12-22 08:47:29 | 显示全部楼层
问题查到了,原来是A矩阵的问题,上下不一致,嘿嘿
回复

使用道具 举报

 楼主| 发表于 2015-12-22 08:47:51 | 显示全部楼层
问题查到了,原来是A矩阵的问题,上下不一致,嘿嘿
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Archiver|手机版|小黑屋|R语言中文网

GMT+8, 2024-11-23 15:13 , Processed in 0.023207 second(s), 17 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表