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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 物理 » 其他交流 » 色散关系 一维傅里叶变换 1D-fft 数据解读
121/2返回列表12下一页
查看: 2828  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

hqyzt

木虫 (小有名气)

[求助] 色散关系 一维傅里叶变换 1D-fft 数据解读 已有2人参与

第一幅图中的任意一条直线代表一个前沿波。第二幅图表示色散关系。我的最终目的是获得每一条波对应的色散关系。
我的做法是：我采用getdata，对图像建立坐标，然后获取了波（x，y）的值对图中任意一个前沿波进行一维傅里叶变换，得到若干个频率，振幅，相位等数据点。
而问题是，频率，振幅和相位所对应的数据点如何能够转换成该前沿波的色散关系，也就是波数与角频率之间的关系？若不能，那么问题出在哪？需要补充哪些条件后，才能知晓。

附上参考文献的做法：
After a background subtraction, the front line is tracked down by an intensity threshold criterion and subjected to a one-dimensional Fourier transform. To resolve the continuous spectrum correctly, a Hamming window is applied before the fast Fourier transform. For the determination of the dispersion relation, an exponential function A(k,t)= A0(k)exp{n(k)t} is fitted to the time evolution of the lowest 30 modes, yielding a growth rate n(k) for each of these modes.

色散关系 一维傅里叶变换 1D-fft 数据解读
2.jpg


色散关系 一维傅里叶变换 1D-fft 数据解读-1
3.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2014-04-25 08:29:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

孙?�730

至尊木虫 (知名作家)

优秀版主
引用回帖:
3楼: Originally posted by hqyzt at 2014-04-25 10:28:05
1. 您让我分清的是我的波是在时域还是频域的吗？ 我第一幅图中的波，是沿着纵坐标整体向下，时间间隔是6秒。貌似既不是时间，也不是频域，而是空间上的。
2. 频域参数指的是哪些？
3. 抱歉，我看不懂。...

频域参数是f或w（差2pi）
与时间t是一对
k与x是一对 上面那个写错了
f频域 t时域
x实空间 k波矢空间
赞一下(3人) 回复此楼
sunyu730_emuch@126.com
6楼2014-04-25 16:25:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

孙?�730

至尊木虫 (知名作家)

优秀版主

【答案】应助回帖

感谢参与，应助指数 +1
1.分清时域和频域
2.色散关系 频域参数是横坐标
3.k频域 t时域
赞一下(1人) 回复此楼
sunyu730_emuch@126.com
2楼2014-04-25 08:39:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhengyongyb

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
hqyzt: 金币+50, ★★★很有帮助, 感谢 2014-04-25 15:41:25
你对图中任意一个前沿波（比如ti 时刻的）进行一维傅里叶变换，相当于把t时刻的A(x)--->A(k), 得到的不是频率，而是波数，当然了，是t时刻的, 记作A(k，ti);
你把图中每一个前时刻的前沿波都进行这样的一维傅里叶变换，得到每个时刻的A(k，ti);

每个k, 对A(k，ti) 作时间傅里叶变换，就得到
A(k，w) ,
w是频率；
每个k, A(k，w) 显然是w的函数，但是一般很尖锐，认为其中使 A(k，w)最大的那个w附近A(k，w)不为0， 其他地方为0；这个w我们记为wk,
每个k对应的wk显然是不一样的，对应画出来就是色散图了
赞一下(1人) 回复此楼
4楼2014-04-25 13:01:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhengyongyb

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
hqyzt: 金币+10, ★★★很有帮助, 感谢 2014-04-26 07:54:58
引用回帖:
5楼: Originally posted by hqyzt at 2014-04-25 15:38:22
1. 第一个一维傅里叶变换，是傅里叶正变换还是负变换?用什么软件可以实现呢？...

？
你不是做过吗？
等距取2N个点，然后快速傅里叶变换啊，这个用matlab可以，也可以自己编一个fortran程序啊，代码可以直接在网上下
赞一下(1人) 回复此楼
8楼2014-04-25 20:47:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhengyongyb

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
hqyzt: 金币+10, ★★★很有帮助, 感谢 2014-04-26 07:55:04
引用回帖:
5楼: Originally posted by hqyzt at 2014-04-25 15:38:22
1. 第一个一维傅里叶变换，是傅里叶正变换还是负变换?用什么软件可以实现呢？...

注意每条线取点时采用的间距都要相等，这样出来的k单位才一样啊
赞一下(1人) 回复此楼
9楼2014-04-25 20:50:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhengyongyb

金虫 (正式写手)
引用回帖:
7楼: Originally posted by hqyzt at 2014-04-25 17:06:10
每个k, 对A(k，ti) 作时间傅里叶变换. 这是什么意思？抱歉，我才疏学浅，实在没看懂。...

你每条线作fft不是都会出来一组{A(k)}嘛，其中 k=1,2,.....,2N-1, 2N，
对时间的fft这样做：
每个k都做一次, 比如，k=8情况，我把每条线fft变换出的A(8)都拿出来，如果有2M条线，不就有2M个A(8)值嘛，相当于t=0s的A(8)值，t=6s的A(8)值，t=12s的A(8)值，。。。。。t=6*2M秒的A(8)值，2M个A(8)值分别对应2M个时刻点，直接fft就行啦啊
赞一下(1人) 回复此楼
10楼2014-04-25 21:05:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

hqyzt

木虫 (小有名气)
???????:
2?: Originally posted by ????730 at 2014-04-25 08:39:26
1.???????????
2.????? ?????????????
3.k??? t???

1. ??????????????????????????????? ????????е????????????????????????????????6??ò??????????????????????????
2. ?????????????Щ??
3. ??????????????
赞一下 回复此楼
3楼2014-04-25 10:28:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hqyzt

木虫 (小有名气)
引用回帖:
4楼: Originally posted by zhengyongyb at 2014-04-25 13:01:55
你对图中任意一个前沿波（比如ti 时刻的）进行一维傅里叶变换，相当于把t时刻的A(x)--->A(k), 得到的不是频率，而是波数，当然了，是t时刻的, 记作A(k，ti);
你把图中每一个前时刻的前沿波都进行这样的一维傅里 ...

1. 第一个一维傅里叶变换，是傅里叶正变换还是负变换?用什么软件可以实现呢？
赞一下 回复此楼
5楼2014-04-25 15:38:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hqyzt

木虫 (小有名气)
引用回帖:
4楼: Originally posted by zhengyongyb at 2014-04-25 13:01:55
你对图中任意一个前沿波（比如ti 时刻的）进行一维傅里叶变换，相当于把t时刻的A(x)--->A(k), 得到的不是频率，而是波数，当然了，是t时刻的, 记作A(k，ti);
你把图中每一个前时刻的前沿波都进行这样的一维傅里 ...

每个k, 对A(k，ti) 作时间傅里叶变换. 这是什么意思？抱歉，我才疏学浅，实在没看懂。
赞一下 回复此楼
7楼2014-04-25 17:06:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hqyzt 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 15:18 by KM5EcsNQRBPn
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 14:57 by KM5EcsNQRBPn
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 6/300 2026-08-23 14:45 by KM5EcsNQRBPn
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-21 5/250 2026-08-23 14:33 by KM5EcsNQRBPn
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 6/300 2026-08-23 14:21 by KM5EcsNQRBPn
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 2026国自然函评费到账 +15 羊腰板 2026-08-21 16/800 2026-08-23 10:45 by process2012
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 10:44 by IXZuIJ2Q7OVy
[基金申请] 让我中一个面上吧！ +12 大萍1987 2026-08-20 14/700 2026-08-23 10:30 by wrm
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 今日不放榜？网传国自然预计 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 紫杉醇
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 估计是周四 +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
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
关闭
请填处理意见
关闭