24小时热门版块排行榜    

查看: 1601  |  回复: 4

oscar1860

金虫 (小有名气)

[交流] 【求助】推求“稳健变异函数”可通过什么软件实现?

我知道GS+、AICGIS是可以计算变异函数的,但是用中位数中位偏差法或中位数调节法这些公式来推求稳健变异函数,与传统求变异函数公式不一样了,用什么软件可以实现呢?软件源代码我也改不了呀!急求!

[ Last edited by oscar1860 on 2010-10-6 at 13:02 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Ky_20060924

银虫 (正式写手)

★
oscar1860(金币+1):谢谢参与
oscar1860(金币+5): 2010-10-12 01:04:02
变异函数
地质统计学中用以研究区域化变量空间变化特征和强度的手段和工具,它被定义为区域化变量增量平方的数学期望,即区域化变量增量的方差。在以向量h相隔两点x、x+h处的两个区域化变量值Z(x)与Z(x+h)之间的变异,用变异函数2γ(x,h)表征:2γ(x,h)=E{〔Z(x+h) Z(x)〕2}=∫〔Z(x+h)-Z(x)〕2dx实际工作中只有有限个观测点,难以上式计算变异函数。根据内蕴假设,增量〔Z(x+h) Z(x)〕只与h有关,而与x无关。因此,实际计算的实验变异函数2γ*(h)是在以向量h相隔的N对点的两个观测值间增量平方的平均值,即2γ*(h)=1N(h)∑N(h)i=1〔Z(xi+h) Z(xi)〕22γ*(h)为增量方差之半,又叫半变异函数,简称变异函数。利用上式对一系列间距h计算出一系列γ*(h),用以编制出经验半变异函数曲线,亦称半变差图(见图)。该图中的一些特征值反映了矿体变化性。变程a:经验半变异函数曲线及相应理论曲线图两点间距小于变程a时,两点变量有自相关性,大于变程a自相关性消失。可以根据变程a确定工程或样品间距。基台值C(o)反映区域化变量变异性大小,在满足平衡条件时等于数据的先验方差。块金常数C0表示原点处变异函数的不连续性,代表了观测误差、矿化微观变化等导致的随机变化。拱高C表示总变异中的空间变化。典型变异函数曲线分为抛物线型(连续型)、线性型、间断型(块金型)、随机型(纯块金型)、转变型,它们代表了具有不同连续性和随机性的地质体参数的变化性特点

用Matlab可以实现!
2楼2010-10-11 19:09:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oscar1860

金虫 (小有名气)

老兄,不知你是否有编好的程序的代码?  我觉得计算稳健变异函数值的程序比较容易实现,但是模型拟合及拟合出来的各参数这部分的程序代码似乎有点难度。
3楼2010-10-12 01:12:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Ky_20060924

银虫 (正式写手)

oscar1860(金币+5): 2010-10-12 09:12:42
不好意思,理论上的计算我懂,计算简单,模拟是有一定困难,介绍几篇文章,相信对你有一定帮助:
应用GEO-EAS拟合煤矿床估值中变异函数
变异函数的自动拟合研究
基于粒子群优化的理论变异函数拟合方法研究
遗传算法拟合变异函数的参数http://www.doc88.com/p-5750304721.html
4楼2010-10-12 07:56:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yan2688

新虫 (初入文坛)

★
oscar1860(金币+1): 谢谢参与
楼主您好,我最近也在做这个,请教下您是怎么做的
5楼2018-07-12 18:37:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 oscar1860 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-11 4/200 2026-11-12 02:56 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 Z7qfCsdIiD1T 2026-11-11 4/200 2026-11-12 02:53 by 9fOTqj074tIg
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 Z7qfCsdIiD1T 2026-11-10 5/250 2026-11-12 02:45 by 9fOTqj074tIg
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Z7qfCsdIiD1T 2026-11-10 5/250 2026-11-12 02:36 by 9fOTqj074tIg
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 02:31 by 9fOTqj074tIg
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 czAiolzR0IJC 2026-11-10 4/200 2026-11-12 02:12 by 9fOTqj074tIg
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 4/200 2026-11-12 02:12 by 9fOTqj074tIg
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 5/250 2026-11-12 02:09 by 9fOTqj074tIg
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 G2rhSvqvhNzl 2026-11-10 3/150 2026-11-12 01:50 by 9fOTqj074tIg
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 G2rhSvqvhNzl 2026-11-10 4/200 2026-11-12 01:50 by 9fOTqj074tIg
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-12 01:50 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 YHiw1gou1Eda 2026-11-10 5/250 2026-11-12 01:17 by 9fOTqj074tIg
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 YHiw1gou1Eda 2026-11-10 6/300 2026-11-12 01:08 by 9fOTqj074tIg
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 6/300 2026-11-12 01:06 by 9fOTqj074tIg
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 0TEMgxdT1eIw 2026-11-09 8/400 2026-11-12 00:33 by 9fOTqj074tIg
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-11 17:58 by TgvkbZZl9tqB
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 zSbvVX9JFl9J 2026-11-09 7/350 2026-11-11 16:41 by TgvkbZZl9tqB
[考博] 27博士 +3 锦鲤送福至 2026-11-10 3/150 2026-11-11 15:53 by 可淡不可忘
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
信息提示
请填处理意见