从索引列表构建向量

Question

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

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 矩阵，其中包含 d 和 s 上 P 中的所有索引排列，并且只找到最小行（可能通过范数，但可能只是通过行和）。是否有一个表达式的行为类似于函数？

Answer 1

对于编辑 2，有一个元素除法运算符和一个向量除法运算符。b = d（P）./（d（P） - s（P））（带点）是我需要的。

Answer 2

正如

Divakar指出的那样，解决方案是w（P）。谢谢！