创建一个平衡的双向脉冲对，也称为Lilly Wave，并在octave / matlab中更改它的样本量

Question

问候我想用Octave 3.8.1创建双向脉冲对或Lilly Wave（Lilly wave），就像matlab一样。

这是我发现的Lilly浪潮 Lilly wave

  

以前在神经生理学和神经外科中使用的波形   当单向电流通过时，神经元受损   脑。莉莉博士开发了一种新的电波形式来平衡   当前，首​​先在一个方向，然后，在一个短暂的间隔后，在   另一个。因此，在神经元中移动的离子将首先被推动一个   方式，然后快速反过来，刺激神经元和   将离子留在神经元内的先前位置。这个   新波形被称为平衡双向脉冲对，或者   莉莉波用这种方法刺激大脑的显微镜研究   平衡脉冲对表明没有神经元损伤   这种网络的刺激电流波形：脉冲对   准分化产生的电流，被动的   电子元件，矩形脉冲。测量值为2％   峰值，正脉冲的持续时间（向上）是34秒，并且   负脉冲的持续时间（向下）是28秒。刺激

我有两个问题： 1）这是创造一个＆＃34;莉莉波的最佳方式吗？导出到音频源？ 2）如何将信号变为44100个样本而不是154350个样本。

以下代码：

clear all
graphics_toolkit gnuplot %use this for now it's older but allows zoom
figure
clf
% SCRIPT BEGINS
t=linspace(0,1,44100);
freq=1; %how many in 1 sec
A = 1; % amplitude
T = 1/freq; % period of the signal

%sine wave
ysin=sin(2*pi*freq*t);

square=0*t;

lilly=[ysin(1:length(t)/2),square(1:length(t)/2),-ysin(1:length(t)/2),square,square];
figure;
plot(lilly)

Answer 1

这是第2点的答案。

脉冲有尾巴，就像高斯人一样

# a simple gaussian
gauss = @(t, t0, g) exp(-((t-t0)/g).^2);

# sampling
t=0:1:44100;

# pulses peak positions (s)
t1 = 10000; 
t2 = 30000; 
# pulses width (at 1/e^2) (s)
g = 2000;

lilly = gauss(t, t1, g) - gauss(t, t2, g);

plot(t, lilly)

产生