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

Question

你好，我想用Octave 3.8.1创建一个双向脉冲对，或者Lilly波(Lilly wave)，就像matlab一样。

这就是我发现的莉莉波礼来波

​

​

以前用于神经生理学和神经外科的波形是在单向电流通过大脑时损伤神经元的。礼来博士开发了一种新的电波形式来平衡电流，首先是在一个方向，然后在一个短暂的间隔后，在另一个方向。因此，在神经元中运动的离子首先会以一种方式被推进，然后以另一种方式快速推进，刺激神经元，并将离子留在原来的位置上。这种新的波形被称为平衡双向脉冲对，或Lilly波。这种平衡脉冲对刺激大脑的微观研究表明，这种刺激电流波形对神经元网络没有损伤:由矩形脉冲的准微分和无源电学元件引起的电流脉冲对。按峰值的2%测量，正脉冲(向上)的持续时间为34秒，负脉冲(向下)的持续时间为28秒。刺激。

我有两个问题: 1)这是创建“礼来波”输出到音频源的最佳方法吗? 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)

产额

​