【问题】解矩阵但解须符合特定条件

楼主: fyq (推销员)   2014-07-01 23:47:43
各位好,
小弟有一个问题困扰许久
现在有一个连续方程组
S1X1+S2X2+error=S
T1X1+T2X2+error=T
PV1X1+PV2X2+error=PV
X1+X2+error=1
这是我拿来解水团成分的
其中S1 S2 T1 T2 PV1 PV2 和S T PV都是已知数
要解出X1 X2,
分别表示1水团和2水团分别在我的实际观测数据之中占了多少百分比。
所以我就写成
A=[S1 S2 1;T1 T2 1; PV1 PV2 1; 1 1 1];
U=[S; T; PV; 1];
A/U之后就是X1 X2 和error了
但是
这边我遇到一个问题
这样解出来的X1和X2有可能超过1或者是负值
而这是不可能在现实世界发生的
请问
该如何设入条件让0<X1<1 以及 0<X2<1呢?
谢谢大家的帮忙!!!
作者: jack10313 (彩色笔)   2014-07-02 00:11:00
关键字 Non-negative matrix factorization
楼主: fyq (推销员)   2014-07-02 15:11:00
J<大,我查了你说的关键字,看完还是不知道如何使用...可以请您多说明一点吗?
作者: jack10313 (彩色笔)   2014-07-02 22:59:00
没有办法 这属于硕班的能力 要自己研究才行
楼主: fyq (推销员)   2014-07-03 10:01:00
好的,谢谢J大!不过后来想了一个方法,就是将姐出来的解normalize一样可以得到0<=X<=1的解答J大提供的程式码我会在自己研究,谢谢!
作者: kurt28   2014-07-03 16:57:00
满好奇normalize得到的结果为何会跟另外加入constraint一样

Links booklink

Contact Us: admin [ a t ] ucptt.com