声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1209|回复: 3

[FFT] 傅里叶变换编程求助

[复制链接]
发表于 2012-3-19 03:31 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
各位朋友,这段时间为一个拟合的问题,而弄得是寝食难安。
   我的问题如下:
我所采用的模型为复数形式,其表达式子为:
A* = A' +iA"
A' 为实部,A''为虚部
在时间域上表达为:
A*(t)=B*exp[-(t/C)^D]
t为时间
我自己通过傅里叶变换到频率域上的表达式为:
A*(w)=积分{exp(-iwt)*(-dA*/dt)}dt  (因为是自己变换的,希望是正确的)
其中w为角频率,这个值我是可以通过实验测得的值
积分下限为0,积分上限为正无穷;
我测得数据都是在频率域上测得的
我想要的值是B,C和D。因为这个方程没有解析解,只有数值解。B的初值大概为1, C是1.5, D 大概为0.5左右。

我的测量数据如下:
角频率omiga
[12.56        
10.048        
8.0384        
6.28        
5.1496        
4.11566        
3.29253        
2.63402        
2.10722        
1.68577
1.34862        
1.0789        
0.86312        
0.69049        
0.55239        
0.44192        
0.35353        
0.28283        
0.22626        
0.18101        
0.14481        
0.11585        
0.09268        
0.07414        
]

A'的测量值为:
[0.76538
0.74637
0.72135
0.69147
0.66726
0.63551
0.59859
0.56865
0.52808
0.48369
0.44231
0.39634
0.35425
0.31168
0.2748
0.23469
0.20617
0.17367
0.15272
0.13405
0.11922
0.10928
0.10173
0.09586
]

A''的实验测量值为:
[0.1607
0.17198
0.18464
0.19813
0.20762
0.2193
0.23032
0.23912
0.24767
0.25374
0.25661
0.25599
0.25135
0.24307
0.23132
0.21352
0.19722
0.17348
0.15463
0.13355
0.11465
0.09804
0.0863
0.07535
]
我是每间隔30秒取样一次,随着频率的降低,其系统本身取样的时间要变长的。

请求大家帮忙,非常感谢!!
回复
分享到:

使用道具 举报

 楼主| 发表于 2012-3-20 02:10 | 显示全部楼层
继续求助,非常感谢
 楼主| 发表于 2012-3-20 14:53 | 显示全部楼层
请大侠指教,非常感谢
 楼主| 发表于 2012-3-26 00:13 | 显示全部楼层
希望能得到大家的帮助,再次感谢
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-26 20:09 , Processed in 0.077651 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表