24小时热门版块排行榜    

查看: 765  |  回复: 6
当前主题已经存档。

positron

木虫 (职业作家)

[交流] 求助傅里叶变换相关问题

有限元法求到立方体空间三维势能矩阵Vr(30,30,30)

三维实数矩阵

如何进行傅里叶变换到另一个空间

是这种形式么?V(h)=∫exp(-i*h*r)*Vr*dr

只是每个格点的数值不是r的函数,而是特定的一个值,不知道怎么变换了

谢谢回答帮助
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

senlia

木虫 (小有名气)


positron(金币+1):谢谢参与
您这个题目貌似有点没讲清..第一句没看懂
Vr是一个数阵,不是函数,怎么用V(h)=∫exp(-i*h*r)*Vr*dr这个公式呢?
h是网格步长还是大概看得懂,但是r是什么?
你要用的是不是离散fourier变换呀?这个在fortran和matlab里面都有相关函数的,C有没有我不清楚
书到用时方恨少, 天涯何处无芳草.
2楼2009-12-19 17:03:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

positron

木虫 (职业作家)

引用回帖:
Originally posted by senlia at 2009-12-19 17:03:
您这个题目貌似有点没讲清..第一句没看懂
Vr是一个数阵,不是函数,怎么用V(h)=∫exp(-i*h*r)*Vr*dr这个公式呢?
h是网格步长还是大概看得懂,但是r是什么?
你要用的是不是离散fourier变换呀?这个在fortran和 ...

多谢版主回答,问题没有表述清楚,请见谅

问题是这样的:

已经求得坐标空间的一个势能矩阵Vr,是30*30*30的三维实数矩阵

想把它变换到动量空间去,上边我用了h表示动量空间,又没有说明,引起误解,请见谅
应该是离散傅里叶变换

但是不知道是怎么变换,我以为是矩阵Vr点乘exp(-i*h*r)积分

现在看来可能不是的

不知道matlab求离散傅里叶变换是啥函数,知道的fourier傅里叶变换是对一个函数
求傅里叶变换的,不知道对矩阵如何处理

还有离散傅里叶变换公式是否也是乘以exp(-i*h*r)然后积分?


谢谢帮助,数学太差,也不太清楚如何更好表述问题
3楼2009-12-19 17:43:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

senlia

木虫 (小有名气)

★ ★ ★ ★ ★
positron(金币+5,VIP+0):谢谢版主,学习啦 12-19 17:59
离散fourier变换就是在 不知道原函数的具体形式,但是知道在fourier点[均匀节点]上的函数值 的情况下 对fourier变换的离散化
无非就是把积分∫(...)dx换成求和形式Segma(...)/N 其中N是一维数组使用的Fourier点数
fourier变换和离散fourier变换都适用于周期函数的形式
matlab里的离散fourier变换函数是fft 逆变换是ifft 请在matlab中help查一下 这两个函数是关于一维数组的
至于多维的情况 貌似可以一维一维的变换

[ Last edited by senlia on 2009-12-19 at 17:52 ]
书到用时方恨少, 天涯何处无芳草.
4楼2009-12-19 17:51:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zap65535

木虫 (著名写手)

两院大学士后


positron(金币+1):谢谢参与
不能沉,我也要学习
5楼2009-12-19 19:48:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

波不动

木虫 (正式写手)

Wave No Move


positron(金币+1):谢谢参与
注意到,在matlab运用fft做多维傅立叶变换的时候,中间要进行矩阵的转置才行。
端好自己的碗,吃好自己的饭。
6楼2009-12-20 15:17:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

波不动

木虫 (正式写手)

Wave No Move

★ ★
positron(金币+2,VIP+0):多谢波不动版主,非常感谢,我再试试 12-20 20:42
刚又查了一下,matlab还有一个fftn的函数。就是N维傅立叶变换,就不用矩阵转置了。

比如你输入了一个三维的Vr(30,30,30),
然后打Vh=fftn(Vr)就可以得到结果Vh了。
端好自己的碗,吃好自己的饭。
7楼2009-12-20 15:25:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 positron 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 14/700 2026-08-21 16:51 by applepetty
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 人气不行了 +7 fansofjerry 2026-08-21 7/350 2026-08-21 14:13 by Vivilian
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 时间戳又变了8-15 +15 archvillain 2026-08-15 27/1350 2026-08-21 11:32 by tim76
[基金申请] 今天基金会出结果吗?20260819 +15 kkkl_v 2026-08-19 16/800 2026-08-21 11:22 by 365372687
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 13/650 2026-08-21 08:30 by 北京莱茵编辑
[基金申请] 时间戳今天,20号变了 +4 archvillain 2026-08-20 4/200 2026-08-21 08:10 by gatelove
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见