半値幅からσが決まるのであれば、標準正規分布を作ればいい気がします
データ点数が決まれば、それに対応する最大1のガウス分布を掛ければいいだけかなと
関数にするなら
function y = mygausswin(FWHM, Dataspan, Data_length)
みたいな感じですかね
---------------------------------------------------------------------------------------------------
以下例です
必要なパラメータを記入します
FWHM =25*10^-9;
Data_Span = 100*10^-9;
Data_length = 10^3;
標準ガウス分布を作成します
Norm = 1/Data_Span;
x = linspace(-Data_Span/2,Data_Span/2,Data_length)*Norm;
mu = mean(x);
σは半値幅の式を使います
sigma = FWHM/(2*sqrt(2*log(2)))*Norm;
f = (1/sigma/sqrt(2*pi))*exp(-(x-mu).^2/(2*sigma^2));
y = f./max(f);
プロットします
plot(x,y)
xlim([-1/2 1/2])