使用Octave绘制函数

Question

我是GNU Octave的新手，我想在[{1}}范围内绘制函数psi：

[0 : 2000]

怎么做？

Answer 1

首先，您需要修改psi和H函数才能执行element-wise operations，以便传递s个值的数组它并收到一组y值

function y = H(x)
    y = double(x > 0);
endfunction

function y = psi(s)
  t = 200
  phiabs = 500
  K = 1000
  n0 = 1000
  y = -n0 * exp(-(s - phiabs) ./ t) .* H(s - phiabs) - K .* H(s) .* H(phiabs - s)
endfunction

另请注意，您需要使用exp而不是e.^()。

然后，您将要使用plot在图中绘制结果

s = 0:2000;
y = psi(s);

plot(s, y)