------------------------------------------------------------------------ TITLE请着名软体类别 按ctrl+v可预览 发文前请把灰色的注解删除 可使用ctrl+y删除一整行 [软体程式类别]: 请填入软体程式类别 例如SAS、SPSS、R、EVIEWS...等 gauss [程式问题]: 资料处理、回归、叙述统计、logistic、probit...等 线性规划 [软体熟悉度]: 请把以下不需要的部份删除 新手(不到1个月) 根据不同熟悉度 板上高手可以给予你不同指导 新手或许需要的还是语法的熟悉 而其他部分则可能往往是盲点上的点破 [问题叙述]: 请详尽叙述遭遇到的问题，可能的话，分点叙述你要处理的流程 能够分点叙述处理流程，基本上在撰写程式或者跑统计时是相当有帮助的 甚至有经验的人可以藉此告知你步骤是否错误 [程式范例]: 虽然张贴程式很可怕，但基本上有些程式还是要张贴才能解决 x1 = ones(n1,1)~dat1[.,4:24]; @ x matrix = constraint coefs @ x2 = ones(n2,1)~dat2[.,4:24]; x3 = ones(n3,1)~dat3[.,4:24]; x4 = ones(n4,1)~dat4[.,4:24]; x5 = ones(n5,1)~dat5[.,4:24]; x = x1|x2|x3|x4|x5; yhat1 = x1*parm[.,1]; @ right-hand-sides of constraints@ yhat2 = x2*parm[.,2]; yhat3 = x3*parm[.,3]; yhat4 = x4*parm[.,4]; yhat5 = x5*parm[.,5]; yhat = yhat1|yhat2|yhat3|yhat4|yhat5; n = rows(x); @ no. of observations @ q = zeros(22,22); r = -x'yhat; @ LP objective function @ @q = x'x; r = x'yhat;@ @ QP objective function @ {b,u1,u2,u3,u4,ret} = QProg(parm[.,6],q,r,0,0,x,yhat,0); @ minimise objective function @ 请问如果要改他的限制式方向要怎麽改呢?就是此篇code的限制式是 X'(b*) >= X'(b) <不等式右边已知> 现在我想改成X x'(b*) <= X '(b) 也就是限制式理面的不等式 符号方向相反.... 感谢....... ----------------------------------------------------------------------------- -- And as I draw near The scene becomes clear Like watching my life on a screen --

※ 发信站: 批踢踢实业坊(ptt.cc) ◆ From: 61.59.238.241