Question

我正在尝试编写代码以绘制负阶跃函数（该函数在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)))

倍频程-负阶跃功能

1 个答案: