24小时热门版块排行榜    

CyRhmU.jpeg
查看: 563  |  回复: 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 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见