最近一直在整粗糙表面计算机模拟,根据陈辉、胡元中的《粗糙表面计算机模拟》模拟出了高斯随机表面;根据宋俊杰、田爱玲的《非高斯随机粗糙表面的数字模拟》写了一点Matlab非高斯粗糙表面建模程序,但是在使用Johnson转换系统生成非高斯序列的时候一直不得要领,有没有做过类似工作的研友?求指点
目前遇到的难点:
1.Johnson转换系统的选择
2.Johnson转换系统中无边界系统以及边界系统的参数求解
Johnson转换系统的选择:
常用的Johnson转换系统有四种:Johnson SB、Johnson SU、Johnson SL和Johnson SN。其中,Johnson SB和Johnson SU适用于对称分布,Johnson SL和Johnson SN适用于非对称分布。
如果您需要模拟非对称分布的随机序列,建议选择Johnson SL或Johnson SN。如果您需要模拟对称分布的随机序列,可以选择Johnson SB或Johnson SU。另外,如果您不确定要选择哪种转换系统,可以尝试使用多个转换系统进行比较,选择效果最好的那个。
Johnson转换系统中无边界系统以及边界系统的参数求解:
对于无边界系统,Johnson转换系统中有两种方法可以实现。一种是使用标准正态分布来代替Johnson转换系统,这种方法可以用于对称分布和非对称分布。另一种方法是使用Johnson UB和Johnson UD转换系统,它们是专门用于无边界系统的。
对于边界系统的参数求解,可以使用最大似然估计方法或矩估计方法。最大似然估计方法是基于概率密度函数的最大似然估计,可以使用Matlab中的fminsearch函数进行求解。矩估计方法是基于概率密度函数的矩估计,可以使用Matlab中的mle函数进行求解。另外,也可以使用其他软件包,例如R中的fitdistrplus包。
希望以上信息对您有帮助。如果您有其他问题或需要进一步的帮助,请随时提问。