八度-负步进功能

Question

我正在尝试编写代码来绘制负阶跃函数（当 t*< t 时等于零，当 t>=t* 时等于特定值的函数）。我想要一个代码，我可以自己更改阶跃的幅度和阶跃发生的时间。我尝试了以下功能，但它不适合我要找的东西：

  x = zeros (size (t));
  ind2= t<1000;
  x(ind2)=0;
  ind = t >= 1000 & t < inf;
  x(ind) = -20*10^-5;
  endfunction

有人可以建议我并回答吗？谢谢

Answer 1

这行得通吗？

clear;clf;

## usage: Pristep = Heaviside (t,tstar,amplitude)
##
##
function Pristep = Heaviside (t,tstar,amplitude)
  Pristep=zeros(size(t));Heavisideon=find(t>tstar)
  Pristep(Heavisideon)=amplitude*ones(size(Heavisideon));
endfunction

%% example:

nt=101;t=linspace(-10,10,nt);
plot(t,Heaviside(t,pi,-exp(1)))

八度-负步进功能

Octave- Negative step function

octave