Question

我正在为一个类的最终项目实现NNLS algorithm，我正在尝试从现有的向量w和索引列表中找到一种构建新向量w'的优雅方法P.我笨重的解决方案如下所示。

function [w_prime] = vector_from_indices[w, P]
w_prime = zeros(numel(P));
for i = 1:numel(P)
    w_prime(i) = w(P(i));
end
end

有没有办法优雅地在一条线上或更高效地完成这项工作？我是CS专业的高级数学课程，所以我想我会将程序美学的严谨性应用到我的最终项目中:)

新问题：

function [a] = minimizer(P, d, s)
a = d(P(1))/(d(P(1)) - s(P(1)));
for i =1:numel(P)
    if a > d(P(i))/(d(P(i)) - s(P(i)))
        a = d(P(i))/(d(P(i)) - s(P(i)));
    end
end
end

>> P 
P =
     1     3     4
>> d 
d =
    0.4387
    0.3816
    0.7655
    0.7952
    0.1869
>> s
s =
    0.4898
    0.4456
    0.6463
    0.7094
    0.7547
>> a = minimizer(P,d,s)
a =
   -8.5995
>> d(1)/(d(1) - s(1))
d(3)/(d(3) - s(3))
d(4)/(d(4) - s(4))
ans =
   -8.5995
ans =
    6.4219
ans =
    9.2643
>> b = min(d(P)/(d(P) - s(P)))
b =
         0    3.6806         0
>> d(P)/(d(P) - s(P))

ans =

         0    3.6806         0
         0    6.4219         0
         0    6.6709         0

我创建了一个函数，它最小化了P中索引的d /（d-s）值。有没有办法使用MATLAB的内置函数找到最小化？显然，单行表达式（我指定给b）不起作用，我想我理解为什么;它生成一个3x3矩阵，其中所有索引的排列都在P和d上，并且只是找到最小行（可能是通过一个标准，但可能只是采用行和）。是否有一个表达式就像函数一样？

Answer 1

正如迪瓦卡指出的那样，解决方案是w（P）。谢谢！

Answer 2

对于编辑2，有一个元素分割算子和一个矢量分割算子。 b = d（P）./（d（P） - s（P））（带点）是我所需要的。

从索引列表构建向量

2 个答案: