xiaomosan 发表于 2012-3-19 03:31

傅里叶变换编程求助

各位朋友,这段时间为一个拟合的问题,而弄得是寝食难安。
   我的问题如下:
我所采用的模型为复数形式,其表达式子为:
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秒取样一次,随着频率的降低,其系统本身取样的时间要变长的。

请求大家帮忙,非常感谢!!

xiaomosan 发表于 2012-3-20 02:10

继续求助,非常感谢

xiaomosan 发表于 2012-3-20 14:53

请大侠指教,非常感谢

xiaomosan 发表于 2012-3-26 00:13

希望能得到大家的帮助,再次感谢
页: [1]
查看完整版本: 傅里叶变换编程求助