在Matlab中矢量化函数数组

Question

以下示例说明了我的一般问题：

f=@(x,y) cos(x.*y);
Yvalues = linspace(0,1,50);
W = @(x) f(x,Yvalues);

如果我只想一次评估一个点，那么

可以正常工作。例如：

norm(W(pi/3)-f(pi/3,Yvalues))
ans =

      0

但我如何在任意数量的点上评估W？

提前致谢。

Answer 1

如果你改变了

f=@(x,y) cos(x.*y);

到

f=@(x,y) cos(x'*y);

你可以执行 W（[1 2 3]）

例如，

>> f = @(x,y) cos(x'*y);
>> yv = linspace(0,1,5);
>> W = @(x) f(x,yv);
>> W(1)
ans =
    1.0000    0.9689    0.8776    0.7317    0.5403
>> W(2)
ans =
    1.0000    0.8776    0.5403    0.0707   -0.4161
>> W(3)
ans =
    1.0000    0.7317    0.0707   -0.6282   -0.9900
>> W([1 2 3])
ans =
    1.0000    0.9689    0.8776    0.7317    0.5403
    1.0000    0.8776    0.5403    0.0707   -0.4161
    1.0000    0.7317    0.0707   -0.6282   -0.9900