傅里叶变换编程求助
各位朋友,这段时间为一个拟合的问题,而弄得是寝食难安。我的问题如下:
我所采用的模型为复数形式,其表达式子为:
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秒取样一次,随着频率的降低,其系统本身取样的时间要变长的。
请求大家帮忙,非常感谢!! 继续求助,非常感谢
请大侠指教,非常感谢 希望能得到大家的帮助,再次感谢
页:
[1]