Octave 倍频程半高宽计算

Octave 倍频程半高宽计算,octave,distribution,gaussian,Octave,Distribution,Gaussian,我在计算数据的半高宽时遇到了一些问题。因为信号包中的“fwhm”函数产生的值比我预期的要大100倍 我所做的是 根据高斯分布函数(你可以在维基百科上找到),我生成了一些数据。在该函数中,您可以给出一个特定的σ(RMS)值(FWHM=σ*2.355)。这是我为了解情况而写的剧本 x=10:0.01:40; x0=25; sigma=0.25; y=(1/(sigma*sqrt(2*pi)))*exp(-((x-x0).^2)/(2*sigma^2)); z=fwhm(y)/2.355; plot(

我在计算数据的半高宽时遇到了一些问题。因为信号包中的“fwhm”函数产生的值比我预期的要大100倍

我所做的是

根据高斯分布函数(你可以在维基百科上找到),我生成了一些数据。在该函数中,您可以给出一个特定的σ(RMS)值(FWHM=σ*2.355)。这是我为了解情况而写的剧本

x=10:0.01:40;
x0=25;
sigma=0.25;
y=(1/(sigma*sqrt(2*pi)))*exp(-((x-x0).^2)/(2*sigma^2));
z=fwhm(y)/2.355;
plot(x,y)
当我比较结果时,“fwhm”函数(24.999)的输出比我在函数中使用的(0.25)大100倍

如果你有任何想法,这将是非常有帮助的


提前谢谢

你的z是100倍大,因为你在x中的步长是1/100(0.01)。如果使用半高宽(y),则x中的步长应为1。如果没有,您必须指定

在您的情况下,您应该:

z=fwhm(x, y)/2.355
z =  0.24999

哪个与你的西格玛相匹配?你的z是100倍大,因为你在x中的步长是1/100(0.01)。如果使用半高宽(y),则x中的步长应为1。如果没有,您必须指定

在您的情况下,您应该:

z=fwhm(x, y)/2.355
z =  0.24999
哪个符合你的标准差