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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 谁能帮我用 高版本的 1stopt拟合一下，需要涉及复数
81/1返回列表
查看: 1058  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

阿德

金虫 (正式写手)

[求助] 谁能帮我用 高版本的 1stopt拟合一下，需要涉及复数

我自己大概编了下代码，不知道对不对。自变量是q，因变量是Itotal，需要通过数据拟合得到6个parametersd, beta, fai, deta_p, deta_pu, N。
Fc和Fa是一个高斯函数fc(x)和fa(x)傅里叶变换结果， 我已经算出（不知道1stopt也能不能做傅里叶变换？）
原来的fc（x）和fa（x）也列出，方便大家帮我验证下：
fc（x)=(1/(2*pi*(fai*d*beta)^2))*exp(-((x-fai*d)^2/(2*(fai*d*beta)^2)))
fa(x)=(1/(2*pi*((1-fai)*beta*d)^2))*exp(-((x-(1-fai)*d)^2/(2*((1-fai)*beta*d)^2)))

//Parameters d, beta, fai, deta_p, deta_pu, N;
//Variable q;
//Function Fc=exp(-i*fai*d*q)*exp(-0.5*(fai*beta*d)^2*q^2);
Fa=exp(-i*(1-fai)*d*q)*exp(-0.5*((1-fai)*beta*d)^2*q^2);
R1=(2*det_p^2/q^2)*real(((1-Fc)*(1-Fa))/(1-Fc*Fa)+(1/N)*Fa*((1-Fc)^2/(1-Fa*Fc)^2)*(1-(Fa*Fc)^N));
R2=(2*det_p^2)/(q^2*N)*real((det_pu^2)/(det_p^2)*(1-(Fc*Fa)^N)-(det_pu/det_p)*((1+Fa)*((1-Fc)/(1-Fc*Fa)*(1-(Fc*Fa)^N))));
Itotal=(R1+R2)/(4*pi*q^2);
q=[0.01128
0.01292
0.01544
0.02042
0.02628
0.03208
0.04002
0.05125
0.05585
0.05957
0.06707
0.07123
0.07709
0.08295
0.09209
0.09998
]

Itotal=[0.05389
0.04665
0.03979
0.02857
0.02304
0.01848
0.016
0.01713
0.01943
0.02038
0.01695
0.01465
0.01142
0.00894
0.00629
0.00476
]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-11-08 10:19:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
dbb627: 金币+2, 谢谢指导 2013-11-23 13:12:26
CODE:
Parameters d, beta1, fai, det_p, det_pu, N;
ComplexStr = i;
ConstStr Fc=exp(-i*fai*d*q)*exp(-0.5*(fai*beta1*d)^2*q^2),
         Fa=exp(-i*(1-fai)*d*q)*exp(-0.5*((1-fai)*beta1*d)^2*q^2),
         R1=(2*det_p^2/q^2)*(((1-Fc)*(1-Fa))/(1-Fc*Fa)+(1/N)*Fa*((1-Fc)^2/(1-Fa*Fc)^2)*(1-(Fa*Fc)^N)),
         R2=(2*det_p^2)/(q^2*N)*((det_pu^2)/(det_p^2)*(1-(Fc*Fa)^N)-(det_pu/det_p)*((1+Fa)*((1-Fc)/(1-Fc*Fa)*(1-(Fc*Fa)^N))));
Variable q, Itotal;
Function Itotal=(R1+R2)/(4*pi*q^2);
Data;
q        Itotal
0.01128        0.05389
0.01292        0.04665
0.01544        0.03979
0.02042        0.02857
0.02628        0.02304
0.03208        0.01848
0.04002        0.016
0.05125        0.01713
0.05585        0.01943
0.05957        0.02038
0.06707        0.01695
0.07123        0.01465
0.07709        0.01142
0.08295        0.00894
0.09209        0.00629
0.09998        0.00476

均方差(RMSE):0.00183697334929586
残差平方和(SSE):5.39915373763721E-5
相关系数(R): 0.990875148759682
相关系数之平方(R^2): 0.981833560429521
确定系数(DC): 0.981825482681031
F统计(F-Statistic): 110.0441879638

参数                  最佳估算
--------------------        -------------
d        -0.00249508280554617
beta1        27.8068118307281
fai        1.47620332846665
det_p        0.404674691653248
det_pu        0.665500525801728
n        -39159.4491242722
谁能帮我用 高版本的 1stopt拟合一下，需要涉及复数
f1.jpg
赞一下 回复此楼
2楼2013-11-08 14:38:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿德

金虫 (正式写手)
引用回帖:
2楼: Originally posted by dingd at 2013-11-08 14:38:30
Parameters d, beta1, fai, det_p, det_pu, N;
ComplexStr = i;
ConstStr Fc=exp(-i*fai*d*q)*exp(-0.5*(fai*beta1*d)^2*q^2),
         Fa=exp(-i*(1-fai)*d*q)*exp(-0.5*((1-fai)*beta1*d)^2*q^2),
        ...

非常感谢，但是参数不太对噢
我忘了设定限制范围了，d大于0，beta（0-1）, det_p(0-1),det_pu(0-1), N为大于零的自然数，谢谢
赞一下 回复此楼
3楼2013-11-08 15:59:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿德

金虫 (正式写手)
引用回帖:
3楼: Originally posted by 阿德 at 2013-11-08 15:59:02
非常感谢，但是参数不太对噢
我忘了设定限制范围了，d大于0，beta（0-1）, det_p(0-1),det_pu(0-1), N为大于零的自然数，谢谢...

fai也是0到1范围
赞一下 回复此楼
4楼2013-11-08 16:00:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

均方差(RMSE):0.0025341163115579
残差平方和(SSE):0.000102747927688061
相关系数(R): 0.982569988533571
相关系数之平方(R^2): 0.965443782366861
确定系数(DC): 0.965413209513977
F统计(F-Statistic): 57.8255447208437

参数                  最佳估算
--------------------        -------------
d        17.4646969649632
beta1        2.905559798969E-8
fai        0.482202770714691
det_p        0.00929539851544782
det_pu        0.00417818750227686
n        2

No        目标 itotal.实部        计算 itotal.实部
1        0.05389        0.0544772112041915
2        0.04665        0.0460209586008805
3        0.03979        0.0378363387352564
4        0.02857        0.0294667125736607
5        0.02304        0.0247399787304851
6        0.01848        0.0220435583582139
7        0.016        0.0195685453014309
8        0.01713        0.0169229616389267
9        0.01943        0.0159459032114299
10        0.02038        0.0151753667190481
11        0.01695        0.0136525716887068
12        0.01465        0.0128199810244441
13        0.01142        0.01166022312157
14        0.00894        0.0105185901266382
15        0.00629        0.00878837236948551
16        0.00476        0.00736376960386167
赞一下 回复此楼
5楼2013-11-08 20:24:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿德

金虫 (正式写手)
引用回帖:
5楼: Originally posted by dingd at 2013-11-08 20:24:04
均方差(RMSE):0.0025341163115579
残差平方和(SSE):0.000102747927688061
相关系数(R): 0.982569988533571
相关系数之平方(R^2): 0.965443782366861
确定系数(DC): 0.965413209513977
F统计(F-Statistic): 57. ...

奇怪了，没法按照限定范围拟合么？如果fai在0.5到1之间，beta在0到0.5之间呢

[ 发自小木虫客户端 ]
赞一下 回复此楼
6楼2013-11-10 08:56:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

仔细回看你自己前面的帖子，fai和beta的范围都是在0到1之间。
赞一下 回复此楼
7楼2013-11-10 11:23:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿德

金虫 (正式写手)
引用回帖:
7楼: Originally posted by dingd at 2013-11-10 11:23:31
仔细回看你自己前面的帖子，fai和beta的范围都是在0到1之间。

是的，我是想再缩小下范围

[ 发自小木虫客户端 ]
赞一下 回复此楼
8楼2013-11-10 17:43:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 阿德 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭