数学中国

用户名  找回密码
 注册
帖子
热搜: 活动 交友 discuz
查看: 7351|回复: 3

[求助]高斯曲线弥合中的参数的意义???

[复制链接]
发表于 2006-11-7 15:20 | 显示全部楼层 |阅读模式
    我们的具体问题是:x为处理年限,y为处理对象的变化。在处理的前几年,y呈增长趋势;几年后,y呈现稳定趋势;稳定几年后,y呈现下降趋势;下降几年后,y稳定在一个最低值。我们已经明白:在这个公式中,a代表y的最大值,x0代表y到达最大值时所需的年限(X)。
在著名数据作图软件SigmaPlot中,当对一组符合正态分布的数据弥合曲线时,我们找到了Peak模型中的“Modified Gaussian, 4 parameters”模型,曲线最为符合野外实际规律,其公式为:
f=a*exp(-0.5*abs((x-x0)/b)^c)。
我们已经清楚了a和X0的含义,完全清楚了它们在我们的回归公式中所代表的变量,然而,我们不清楚b和c的含义。
请高手指教。
发表于 2006-11-9 22:32 | 显示全部楼层

[求助]高斯曲线弥合中的参数的意义???

[这个贴子最后由luyuanhong在 2009/03/15 10:45pm 第 1 次编辑]

参数b控制曲线的集中分散程度,b越小,曲线越集中,峰的形状越窄,b越大,曲线越分散,峰的形状越宽。
参数c控制曲线的形状,c越小,峰的形状越接近三角形,c越大,峰的形状越接近方形。
参看附图,上面一幅图是其他参数不变、b变化时的3种曲线形状,下面一幅图是其他参数不变,c变化时的3种曲线形状。
 楼主| 发表于 2006-11-10 11:17 | 显示全部楼层

[求助]高斯曲线弥合中的参数的意义???

[这个贴子最后由peltigera在 2006/11/10 11:19am 第 1 次编辑]

感谢您的回复。[br][br][color=#990000]-=-=-=-=- 以下内容由 peltigera 时添加 -=-=-=-=-
我是做生态调查的,利用有限的数据来弥合更大时空尺度上的变化可能性。我所感兴趣的是,用b和c如何解决实际问题。
以下是我的具体问题:
分析物种 X 的产量随 年限 Y 的变化而变化的趋势。其中物种 X 的产量不可能超过100。根据高斯模型(Modified Gaussian, 4 parameters)弥合的曲线见附件。所用的公式为“f=a*exp(-0.5*abs((x-x0)/b)^c)”,弥合后,得出:
a = 91.8655; b= 5.9908; c= 4.0037; x0= 9.7982。
得出的公式为:
y = 91.8655 * exp((-0.5) * |(x-9.7982)/5.9908|^4.0037)
图片请看附件。
由此公式可知,物种 X 达到最高产量需要至少9.7982年,其最大产量可达91.8655。那么参数 b 是否与这个曲线的拐点有关系?能否说明在某些年限之间物种 X 的变化趋势的变异?参数 c 是否和曲线上方的平台长短有关系,也就是和物种 X 的最大产量持续时间有关系?
发表于 2006-11-10 22:25 | 显示全部楼层

[求助]高斯曲线弥合中的参数的意义???

修正高斯(Modified Gaussian)曲线方程为 y=a*exp{-0.5*abs[(x-x0)/b]^c} 。
其中,参数 b 的意义为:
在高度 y=a*exp(-0.5)(大约等于0.6a)处曲线峰的宽度的一半。(见附图)
因为,当 y=a*exp(-0.5) 时,必有 abs((x-x0)/b)=1,即 b=abs(x-x0) 。
参数 c 的意义为:
取一个接近 0 的值 k1(例如 k1=0.05),设在高度 y=k1*a 处,曲线峰的宽度为 w1,称 w1 为“底宽”(见附图)。
由于当 y=k1*a 时,必有 -0.5*abs((x-x0)/b)^c=ln(k1),即有 abs((x-x0)/b)=[-2*ln(k1)]^(1/c),所以底宽 w1=2b*[-2*ln(k1)]^(1/c)。
取一个接近 1 的值 k2(例如 k2=0.95),设在高度 y=k2*a 处,曲线峰的宽度为 w2,称 w2 为“顶宽”(见附图)。
由于当 y=k2*a 时,必有 -0.5*abs((x-x0)/b)^c=ln(k2),即有 abs((x-x0)/b)=[-2*ln(k2)]^(1/c),所以顶宽 w2=2b*[-2*ln(k2)]^(1/c)。
底宽与顶宽之比为 w1/w2=[ln(k1)/ln(k2)]^(1/c)。因此有 (w1/w2)^c=ln(k1)/ln(k2),即有c=ln[ln(k1)/ln(k2)]/ln[w1/w2],简单说来,参数 c 就是底宽与顶宽之比的对数的倒数。
您需要登录后才可以回帖 登录 | 注册

本版积分规则

LaTEX预览输入 教程 符号库 加行内标签 加行间标签 
对应的 LaTEX 效果:

Archiver|手机版|小黑屋|数学中国 ( 京ICP备05040119号 )

GMT+8, 2025-7-14 09:38 , Processed in 0.082874 second(s), 16 queries .

Powered by Discuz! X3.4

Copyright © 2001-2020, Tencent Cloud.

快速回复 返回顶部 返回列表
\frac{\square}{\square}\sqrt{\square}\square_{\baguet}^{\baguet}\overarc{\square}\ \dot{\baguet}\left(\square\right)\binom{\square}{\square}\begin{cases}\square\\\square\end{cases}\ \begin{bmatrix}\square&\square\\\square&\square\end{bmatrix}\to\Rightarrow\mapsto\alpha\ \theta\ \pi\times\div\pm\because\angle\ \infty
\frac{\square}{\square}\sqrt{\square}\sqrt[\baguet]{\square}\square_{\baguet}\square^{\baguet}\square_{\baguet}^{\baguet}\sum_{\baguet}^{\baguet}\prod_{\baguet}^{\baguet}\coprod_{\baguet}^{\baguet}\int_{\baguet}^{\baguet}\lim_{\baguet}\lim_{\baguet}^{\baguet}\bigcup_{\baguet}^{\baguet}\bigcap_{\baguet}^{\baguet}\bigwedge_{\baguet}^{\baguet}\bigvee_{\baguet}^{\baguet}
\underline{\square}\overline{\square}\overrightarrow{\square}\overleftarrow{\square}\overleftrightarrow{\square}\underrightarrow{\square}\underleftarrow{\square}\underleftrightarrow{\square}\dot{\baguet}\hat{\baguet}\vec{\baguet}\tilde{\baguet}
\left(\square\right)\left[\square\right]\left\{\square\right\}\left|\square\right|\left\langle\square\right\rangle\left\lVert\square\right\rVert\left\lfloor\square\right\rfloor\left\lceil\square\right\rceil\binom{\square}{\square}\boxed{\square}
\begin{cases}\square\\\square\end{cases}\begin{matrix}\square&\square\\\square&\square\end{matrix}\begin{pmatrix}\square&\square\\\square&\square\end{pmatrix}\begin{bmatrix}\square&\square\\\square&\square\end{bmatrix}\begin{Bmatrix}\square&\square\\\square&\square\end{Bmatrix}\begin{vmatrix}\square&\square\\\square&\square\end{vmatrix}\begin{Vmatrix}\square&\square\\\square&\square\end{Vmatrix}\begin{array}{l|l}\square&\square\\\hline\square&\square\end{array}
\to\gets\leftrightarrow\nearrow\searrow\downarrow\uparrow\updownarrow\swarrow\nwarrow\Leftarrow\Rightarrow\Leftrightarrow\rightharpoonup\rightharpoondown\impliedby\implies\Longleftrightarrow\leftharpoonup\leftharpoondown\longleftarrow\longrightarrow\longleftrightarrow\Uparrow\Downarrow\Updownarrow\hookleftarrow\hookrightarrow\mapsto
\alpha\beta\gamma\Gamma\delta\Delta\epsilon\varepsilon\zeta\eta\theta\Theta\iota\kappa\varkappa\lambda\Lambda\mu\nu\xi\Xi\pi\Pi\varpi\rho\varrho\sigma\Sigma\tau\upsilon\Upsilon\phi\Phi\varphi\chi\psi\Psi\omega\Omega\digamma\vartheta\varsigma\mathbb{C}\mathbb{H}\mathbb{N}\mathbb{P}\mathbb{Q}\mathbb{R}\mathbb{Z}\Re\Im\aleph\partial\nabla
\times\cdot\ast\div\pm\mp\circ\backslash\oplus\ominus\otimes\odot\bullet\varnothing\neq\equiv\not\equiv\sim\approx\simeq\cong\geq\leq\ll\gg\succ\prec\in\ni\cup\cap\subset\supset\not\subset\not\supset\notin\not\ni\subseteq\supseteq\nsubseteq\nsupseteq\sqsubset\sqsupset\sqsubseteq\sqsupseteq\sqcap\sqcup\wedge\vee\neg\forall\exists\nexists\uplus\bigsqcup\bigodot\bigotimes\bigoplus\biguplus\bigcap\bigcup\bigvee\bigwedge
\because\therefore\angle\parallel\perp\top\nparallel\measuredangle\sphericalangle\diamond\diamondsuit\doteq\propto\infty\bowtie\square\smile\frown\bigtriangledown\triangle\triangleleft\triangleright\bigcirc \wr\amalg\models\preceq\mid\nmid\vdash\dashv\nless\ngtr\ldots\cdots\vdots\ddots\surd\ell\flat\sharp\natural\wp\clubsuit\heartsuit\spadesuit\oint\lfloor\rfloor\lceil\rceil\lbrace\rbrace\lbrack\rbrack\vert\hbar\aleph\dagger\ddagger

MathQuill输入:

Latex代码输入: