2

計算で 50 の個々の振幅のそれぞれを使用できるように、連続Gaussianパルスを分割する必要があります。50 discrete partsこれが私が試したことです：

% Gauss pulse discretisation
tg = 20*10^(-3);         % pulse duration [sec]
B1 = 1;                  % max amplitude [muT]
t  = -tg/2:tg/50:tg/2;   % sampling times
sd = 0.25;               % pulse standard deviation
% pulse shape
p  = B1*exp(-((t-tg/2).^2)/(2*sd.^2));
plot(t,p);

ただし、プロットは20ms in duration!のガウス パルスのようには見えません。sampling timeの定義方法に問題はありますか? たとえば、が次のsampling timeように定義されている場合

t  = -1:tg/50:1

パルスはガウスのように見えますが、5001 の部分に分解されます。誰かが私を正しい方向に向けてもらえますか?

4

1 に答える 1

0

ガウスをプロットしたときにガウスのように見えるようにするには、次のことを確認する必要があります。(1) 中心付近でサンプリングする、(2) サンプリング間隔が標準偏差 (SD) よりもはるかに小さい、および(3) 各側に少なくとも 2 つまたは 3 つの SD をサンプリングして、減衰を確認します。したがって、あなたの例では、ガウシアンは を中心tg\2にしているため、SD はsd = 0.25(ところで、SD は ではなくパルス持続時間を設定しますtg)、SD を測定値として (および ではなく) 使用してサンプリング間隔を延長し、tgそれを移動します。平均を中心にしています。を使用すると、これらすべてを簡単に実行できますlinspace。

t = linspace(-3*sd, 3*sd, 50) + tg\2;

さらに 20msec のパルス幅が必要な場合はsd、 ではなく 20msec のオーダーで作成してくださいtg。また、「持続時間」はガウス分布の定義の問題であり、+- 無限大に及ぶことに注意してください。「パルス持続時間は -2 SD から +2 SD まで」のように定義する必要があります。これは、有効なパルス持続時間がテールの減衰量に応じて定義されることを意味します。

% Gauss pulse discretisation
tg = 0;         % pulse center [sec]
B1 = 1;         % max amplitude [muT]
sd = .5*20*10^(-3);% half the pulse duration (msec)
t  = tg/2 + linspace(-3*sd,3*sd,50);  
% pulse shape
p  = B1*exp(-((t-tg/2).^2)/(2*sd.^2));
plot(t,p,'.-');

ガウス

于 2015-11-29T20:00:28.830 に答える