Octave 最小y的倍频程x值

我试图在输入到倍频程的简单函数中辨别最小y的x值。我有最小的y，我可以回去搜索x给y的值，但我试图让它一次显示坐标。我肯定有一个简单的方法去做，但我肯定错过了一些东西。下面是函数：

%W1
%function W=kopt(R,k);
R=3;
k=[0:0.001:3];
F=(-0.5*k.^2)+(k.^2-k-(R.^(-1))*(1+k.*R).*(e.^(-2*k.*R)).+k.*(k.-2).* 
(1+k.*R).*e.^(-k.*R))./(1+(e.^(-k.*R)).*(1+k*R+(k.^2)*(R.^2)/3));
W=F+1/R;
plot(k,W)
min(W)

谢谢你的帮助

可以返回最小值的索引：

> [minval, iminval] = min(W)
minval = -0.31700
iminval =  981
> k(iminval)
ans =  0.98000

所以在

绘制（k，W）

之后，您可以

[minval, iminval] = min(W);
hold()
plot(k(iminval), minval)

可以返回最小值的索引：

> [minval, iminval] = min(W)
minval = -0.31700
iminval =  981
> k(iminval)
ans =  0.98000

所以在

绘制（k，W）

之后，您可以

[minval, iminval] = min(W);
hold()
plot(k(iminval), minval)