24小时热门版块排行榜    

查看: 2916  |  回复: 11

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
回复此楼
已阅   回复此楼   关注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波矢空间
sunyu730_emuch@126.com
6楼2014-04-25 16:25:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

孙?�730

至尊木虫 (知名作家)

优秀版主

【答案】应助回帖

感谢参与,应助指数 +1
1.分清时域和频域
2.色散关系 频域参数是横坐标
3.k频域 t时域
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显然是不一样的,对应画出来就是色散图了
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程序啊,代码可以直接在网上下
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单位才一样啊
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就行啦啊
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 的主题更新
信息提示
请填处理意见