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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 谁能帮我用 高版本的 1stopt拟合一下，需要涉及复数
81/1返回列表
查看: 1123  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 326求调剂 +3 诺贝尔化学奖觊� 2026-03-15 3/150 2026-03-16 13:51 by houyaoxu
[考研] 0703化学调剂 290分有科研经历，论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 304求调剂 +4 素年祭语 2026-03-15 4/200 2026-03-16 09:42 by 闲人终南山
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 材料与化工（0856）304求B区调剂 +7 邱gl 2026-03-10 11/550 2026-03-14 12:18 by 邱gl
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 一志愿湖师大化学289求调剂 +6 XMCMM3.14159 2026-03-10 6/300 2026-03-14 00:28 by JourneyLucky
[考研] b区环境工程求调剂 +4 Maps1 2026-03-10 6/300 2026-03-14 00:23 by JourneyLucky
[考研] 材料工程,326分，求调剂 +6 KRSLSR 2026-03-10 6/300 2026-03-13 23:47 by JourneyLucky
[考研] 0703，333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 336求调剂 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 304求调剂（085602一志愿985） +12 化工人999 2026-03-09 12/600 2026-03-13 12:02 by JourneyLucky
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 0856化工原理 +6 z2839474511 2026-03-10 6/300 2026-03-13 10:41 by houyaoxu
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生，2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
信息提示
关闭
请填处理意见
关闭