应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 关于delta函数的高斯积分问题
51/1返回列表
查看: 3179  |  回复: 32
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

chemistryxrw

金虫 (小有名气)

[交流] 关于delta函数的高斯积分问题

如图，A是归一化系数，P是概率密度分布函数，现在如果我想求xn的概率分布密度函数，我的解法是对其他维的变量进行负无穷到正无穷的积分，delta函数用其的傅里叶变换代入式子积分，这样就可以得到xn的概率密度分布函数，然后我用经典的box-muller方法来做gauss分布，得到了各个维数的x值，这样得到的x数值我认为应该没有问题，但是这时候我将各个维数的x值如delta函数内的式子加起来却发现结果不等于0(我做的大概是10-4左右)，向大神们请教。
(a>0,数量级大概在10-5左右，b>0,数量级在10-6左右， 我取了50000个x值，我用的是fortran90)

1.JPG

[ Last edited by chemistryxrw on 2013-3-13 at 22:40 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2013-03-13 22:15:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leedobb

金虫 (正式写手)
引用回帖:
16楼: Originally posted by chemistryxrw at 2013-03-19 08:53:12
那delta函数怎么解析掉？不太明白你的意思，呵呵。。如果要化成exp(-ax^2+i k b x)的形式不也是先傅里叶变换积分后数值计算吗？不太明白望指教~
能不能讲详细一些，十分感谢...

这个积分对x积，直接就积出来了
结果是
sqrt(pi/a) exp ( -k^2 b^2/4/a)
对每个x都做一次，最后再对k做，
也是高斯积分。

多查查数学用表好了。
赞一下 回复此楼
17楼2013-03-19 12:12:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 33 个回答

chemistryxrw

金虫 (小有名气)
引用回帖:
2楼: Originally posted by sundapeng111 at 2013-03-15 09:56:34
看迷糊了 帮顶

呵呵。。谢谢
回复此楼
3楼2013-03-15 17:32:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

userhung

禁虫 (文学泰斗)
★
chemistryxrw(金币+1): 谢谢参与
δ函数用傅里叶变换代入式子积分没弄错把？
赞一下 回复此楼
8楼2013-03-15 18:48:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chemistryxrw

金虫 (小有名气)
引用回帖:
8楼: Originally posted by userhung at 2013-03-15 18:48:51
δ函数用傅里叶变换代入式子积分没弄错把？

应该没弄错，其实我觉得可能就是取分布的问题，毕竟取x值的时候还是得满足这个等式的，那应该怎么取分布？
赞一下 回复此楼
9楼2013-03-18 09:20:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
chemistryxrw12楼
2013-03-18 09:21   回复  
引用回帖:
6楼: Originally posted by wgatssd at 2013-03-15 18:06:18

chemistryxrw13楼
2013-03-18 09:21   回复  
引用回帖:
7楼: Originally posted by zhihaikou at 2013-03-15 18:46:41 祝福

查看全部 33 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭