混合像元分解中的完全约束的最小二乘怎么解

快马加鞭

现有两个矩阵：Sample矩阵(112*120)，Test矩阵(112*1)，满足  a*Sample(:,1:30)+b*Sample(:,31:35)+c*Sample(:,36:90)+d*Sample(:,91:120) = Test;
同时要求满足约束条件：a>0,b>0.c>0,d>0;a+b+c+d=1;
想请教如何在matlab中求解最优的a、b、c、d参数呢？

这个问题的道理极其类似于混合像元分解中的有约束条件的最小二乘，如果要自己写程序求解，应该如何解呢

luxinloss

用c#程序来写