版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (5374)
>考研 (2355)
>导师招生 (1689)
>虫友互识 (459)
>文献求助 (229)
>硕博家园 (152)
>休闲灌水 (131)
>招聘信息布告栏 (118)
>考博 (108)
>论文投稿 (46)
>博后之家 (41)
>基金申请 (37)
>教师之家 (32)
>公派出国 (25)
>海外博后 (23)
>绿色求助(高悬赏) (22)

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

yj112358

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 82.4
散金: 16
帖子: 96
在线: 103.5小时
虫号: 2837253
注册: 2013-11-29
专业: 胶体与界面化学

[求助] 求大神帮我用1stopt拟合带有复数的函数已有1人参与

求大神用1stopt帮我拟合一下这个带有复数的函数啊，拜托了，ConstStr里是一系列的关系式，函数RA最终是一个实数，
ComplexStr = i;
Parameters kmax, d, sild;
constant nord=1.46;next=1.57;n1=1;c=3*(10^8);alpha=pi/2;fa1=pi*36/180;n2=1.402;k2=0.0157;v0=2923;fw=16.50098;k0=2*pi*v/1000;
Variable v, RA ;
ConstStr nx=next*(sin(sild))^2+nord*(cos(sild))^2,f=0.5*(3*(cos(sild))^2-1),kx_max=kmax*(f*(sin(alpha))^2/2+(1-f)/3),V=2*pi*c*fw;,delta=2*pi*c*(v-v0),kx=kx_max*V^2/(4*delta^2+V^2),Nx=nx+1i*kx,N2=n2+1i*k2,I1=d*(Nx^2-N2^2),I2=k0/N2*sin(fa1)*I1,I3=1i*I1,fa2=asin(n1/N2*sin(fa1)),b=sin(fa2-fa1)+I3,b1=sin(fa1+fa2)-I3,rs=b/b1,rf=sin(fa2-fa1)/sin(fa2+fa1),Rsf=rf*conj(rf),Rs=rs*conj(rs);
Function RA=-log10(Rs/Rsf) ;
Data;
v RA
2890.43 -0.000144037
2891.77 -0.000148991
2893.12 -0.000158899
2894.81 -0.000163853
2896.49 -0.000173761
2898.52 -0.00018367
2900.54 -0.000198532
2902.9 -0.000223303
2904.92 -0.000243119
2906.94 -0.00026789
2908.29 -0.000292661
2909.3 -0.000307523
2909.97 -0.000317431
2910.64 -0.000327339
2911.32 -0.000342202
2912.33 -0.000366972
2911.65 -0.00035211
2912.66 -0.000376881
2913.33 -0.000391743
2913.67 -0.000406606
2914.01 -0.000421468
2914.34 -0.000431376
2914.68 -0.000451193
2915.35 -0.000466055
2915.68 -0.000485872
2916.02 -0.000505688
2916.69 -0.000530459
2917.03 -0.000550275
2917.36 -0.000565138
2918.03 -0.00058
2918.37 -0.000594862
2918.7 -0.000609725
2919.38 -0.000629541
2919.71 -0.000644404
2920.39 -0.000659266
2921.06 -0.000679083
2921.73 -0.000688991
2922.4 -0.000698899
2923.42 -0.000703853
2924.43 -0.000693945
2924.77 -0.000688991
2925.44 -0.000679083
2925.78 -0.00066422
2926.12 -0.000654312
2926.46 -0.000644404
2927.14 -0.000624587
2927.48 -0.000614679
2927.81 -0.000604771
2928.15 -0.000584954
2928.83 -0.000570092
2929.51 -0.000550275
2930.18 -0.000535413
2930.52 -0.000515596
2931.2 -0.00049578
2931.54 -0.000480917
2932.22 -0.000466055
2932.56 -0.000446239
2932.9 -0.000431376
2933.91 -0.000401651
2934.93 -0.000371927
2935.94 -0.000342202
2936.96 -0.000322385
2937.97 -0.000297615
2938.65 -0.000277798
2939.33 -0.000257982
2940 -0.000243119
2941.02 -0.000223303
2942.37 -0.000198532
2943.05 -0.000173761
2944.4 -0.000153945
2945.41 -0.000148991
2946.76 -0.000129174
2947.44 -0.000114312
2948.79 -9.94495E-05
2950.48 -8.95413E-05
2951.49 -8.45872E-05
2952.84 -8.45872E-05
2953.85 -8.45872E-05

回复此楼

» 本帖@通知

@dbb627

» 猜你喜欢

［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生已经有0人回复
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有97人回复
0702一志愿吉大B区求调剂本科期间发表一篇Sci 已经有2人回复
070200求调剂，一志愿某211，288分已经有18人回复
法国博士后职位已经有0人回复
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入机器视觉与3D光学成像课题组。已经有0人回复
广州大学光电信息工程专业调剂，招收物理学专业学生已经有1人回复
0702一志愿吉大B区求调剂本科期间发表一篇Sci 已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1stOpt 求助！！！已经有4人回复
用1stopt非线性回归，拟合度较差已经有3人回复
求大哥大姐帮跑一下1stopt~新手啊，版本低啊已经有6人回复
1stOpt函数中含有间接变量怎么办已经有3人回复
求大神帮忙用1stOpt 拟合复数数据已经有10人回复
求大神跑下1stOpt，一个sharemodel，谢谢！已经有7人回复
求大神帮忙拟合一个非线性方程，求出模型参数已经有15人回复
用1stOpt进行复数曲线拟合时，结果出错。急求教！！！已经有5人回复
求大神帮忙用1stOpt 拟合复数数据已经有13人回复
求助一个逻辑斯蒂函数的拟合已经有8人回复
大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？已经有8人回复
matlab高手帮忙做个自定义复杂函数拟合，含有复数，已经有25人回复
用MATLAB或者ORIGIN怎么拟合复数数据已经有10人回复
跪求1stOpt2.5以上版本，或者请高人帮忙拟合包含复数的函数的参数已经有8人回复

1楼 2015-07-01 11:13:52

已阅关注TA 给TA发消息送TA红花 TA的回帖

回帖支持 ( 显示支持度最高的前 50 名 )

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
yj112358(月只蓝代发): 金币+5, 感谢耐心应助 2016-11-30 13:18:22
月只蓝: 金币+5, 尊重应助者的劳动和热情。 2016-11-30 13:23:36

CODE:

Title "Type your title here";

ComplexStr = i;

Parameters kmax, d, sild,nmax;

constant n1=1,

c=3*(10^8),

alpha=pi/2,

fa1=pi*36/180,n2=1.402,k2=0.0157,v0=2923,fw=16.50098,y=2*pi*c*fw;

ConstStr f=0.5*(3*(cos(sild))^2-1),

delta=2*pi*c*(v-v0),

k0=2*pi*v/1000,

kx_max=kmax*(f*(sin(alpha))^2/2+(1-f)/3),

kz_max=kmax*(f*(cos(alpha))^2+(1-f)/3),

nx=nmax-2*y*kx_max*delta/(4*delta^2+y^2),

nz=nmax-2*y*kz_max*delta/(4*delta^2+y^2),

kx=kx_max*y^2/(4*delta^2+y^2),

kz=kz_max*y^2/(4*delta^2+y^2),

n_x=nx+kx*i,

n_z=nz+kz*i,

n_2=n2+k2*i,

I1=d*((n_x)^2-(n_2)^2),

I2=((n_z)^2-(n_2)^2)*d/(n_z)^2,

I3=i*k0*sin(fa1)*I1/(n_2),

fa2=arcsin((n_z)/(n_2)*sin(fa1)),

b=I3-sin(fa2-fa1),

b1=sin(fa1+fa2)-I3,

rs=b/b1,

rf=sin(fa2-fa1)/sin(fa2+fa1),

rf_1=rf*Conjugate(rf),

rs_1=rs*Conjugate(rs);

Variable v, ra ;

Function ra=-log10(rs_1/rf_1) ;

Data;

v         RA

2890.43        -0.000144037

2891.77        -0.000148991

2893.12        -0.000158899

....

均方差(RMSE):2.38604321991028E-5
残差平方和(SSE):4.44069775287826E-8
相关系数(R): 0.993149879702572
相关系数之平方(R^2): 0.986346683553234
确定系数(DC): 0.984753822545759
F统计(F-Statistic): 1617.89187753453

参数最佳估算
-------------------- -------------
kmax -0.974602060214981
d 4.8823738972729E-5
sild -1.10099552620844
nmax 2.0621648149797

c164.jpg

赞一下(1人)

回复此楼

7楼2015-07-01 21:20:35

已阅关注TA 给TA发消息送TA红花 TA的回帖

普通回帖

yj112358

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 82.4
散金: 16
帖子: 96
在线: 103.5小时
虫号: 2837253
注册: 2013-11-29
专业: 胶体与界面化学

求求各位大神了

回复此楼

2楼2015-07-01 12:34:23

已阅关注TA 给TA发消息送TA红花 TA的回帖

yj112358

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 82.4
散金: 16
帖子: 96
在线: 103.5小时
虫号: 2837253
注册: 2013-11-29
专业: 胶体与界面化学

发现函数写错了，现改好放在下面，求各位大侠帮忙跑跑？？？
Title "Type your title here";
ComplexStr = i;
Parameters kmax, d, sild;
constant n1=1;
c=3* (10^8);
alpha=pi/2;
fa1=pi*36/180;n2=1.402;k2=0.0157;v0=2923;fw=16.50098;k0=2*pi*v/1000;
Variable v, RA ;
ConstStr f=0.5*(3*(cos(sild))^2-1);
delta=2*pi*c*(v-v0);
V=2*pi*c*fw;
kx_max=kmax*(f*(sin(alpha))^2/2+(1-f)/3);
kz_max=kmax*(f*(cos(alpha))^2+(1-f)/3);
nx=nmax-2*V*kx_max*delta/(4*delta^2+V^2);
nz=nmax-2*V*kz_max*delta/(4*delta^2+V^2);
kx=kx_max*V^2/(4*delta^2+V^2);
kz=kz_max*V^2/(4*delta^2+V^2);
Nx=nx+1i*kx;
Nz=nz+kz*1i;
N2=n2+1i*k2;
I1=d*(Nx^2-N2^2);
I2=(Nz^2-N2^2)*d/Nz^2;
I3=1i*k0/N2*sin(fa1)I1;
fa2=asin(Nz/N2*sin(fa1));
b=-sin(fa2-fa1)+I3;
b1=sin(fa1+fa2)-I3;
rs=b/b1;
rf=sin(fa2-fa1)/sin(fa2+fa1);
Rsf=rf*conj(rf);
Rs=rs*conj(rs);
Function RA=-log10(Rs/Rsf) ;
Data;
v RA
2890.43 -0.000144037
2891.77 -0.000148991
2893.12 -0.000158899
2894.81 -0.000163853
2896.49 -0.000173761
2898.52 -0.00018367
2900.54 -0.000198532
2902.9 -0.000223303
2904.92 -0.000243119
2906.94 -0.00026789
2908.29 -0.000292661
2909.3 -0.000307523
2909.97 -0.000317431
2910.64 -0.000327339
2911.32 -0.000342202
2912.33 -0.000366972
2911.65 -0.00035211
2912.66 -0.000376881
2913.33 -0.000391743
2913.67 -0.000406606
2914.01 -0.000421468
2914.34 -0.000431376
2914.68 -0.000451193
2915.35 -0.000466055
2915.68 -0.000485872
2916.02 -0.000505688
2916.69 -0.000530459
2917.03 -0.000550275
2917.36 -0.000565138
2918.03 -0.00058
2918.37 -0.000594862
2918.7 -0.000609725
2919.38 -0.000629541
2919.71 -0.000644404
2920.39 -0.000659266
2921.06 -0.000679083
2921.73 -0.000688991
2922.4 -0.000698899
2923.42 -0.000703853
2924.43 -0.000693945
2924.77 -0.000688991
2925.44 -0.000679083
2925.78 -0.00066422
2926.12 -0.000654312
2926.46 -0.000644404
2927.14 -0.000624587
2927.48 -0.000614679
2927.81 -0.000604771
2928.15 -0.000584954
2928.83 -0.000570092
2929.51 -0.000550275
2930.18 -0.000535413
2930.52 -0.000515596
2931.2 -0.00049578
2931.54 -0.000480917
2932.22 -0.000466055
2932.56 -0.000446239
2932.9 -0.000431376
2933.91 -0.000401651
2934.93 -0.000371927
2935.94 -0.000342202
2936.96 -0.000322385
2937.97 -0.000297615
2938.65 -0.000277798
2939.33 -0.000257982
2940 -0.000243119
2941.02 -0.000223303
2942.37 -0.000198532
2943.05 -0.000173761
2944.4 -0.000153945
2945.41 -0.000148991
2946.76 -0.000129174
2947.44 -0.000114312
2948.79 -9.94495E-05
2950.48 -8.95413E-05
2951.49 -8.45872E-05
2952.84 -8.45872E-05
2953.85 -8.45872E-05

赞一下

回复此楼

3楼2015-07-01 14:34:44

已阅关注TA 给TA发消息送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

感谢参与，应助指数 +1

好多错误啊：
1：
constant n1=1;
c=3* (10^8);
alpha=pi/2;
fa1=pi*36/180;n2=1.402;k2=0.0157;v0=2923;fw=16.50098;k0=2*pi*v/1000;
中间用“,”号，别用“;”号，最后的“k0=2*pi*v/1000”对吗？v可不是常数；
2：
1stOpt不区分大小写，代码里混写的太多了；
3：I3=1i*k0/N2*sin(fa1)I1,有低级错误吧
其它自己慢慢改。

赞一下

回复此楼

4楼2015-07-01 16:15:34

已阅关注TA 给TA发消息送TA红花 TA的回帖

yj112358

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 82.4
散金: 16
帖子: 96
在线: 103.5小时
虫号: 2837253
注册: 2013-11-29
专业: 胶体与界面化学

引用回帖:

4楼: Originally posted by dingd at 2015-07-01 16:15:34
好多错误啊：
1：
constant n1=1;
c=3* (10^8);
alpha=pi/2;
fa1=pi*36/180;n2=1.402;k2=0.0157;v0=2923;fw=16.50098;k0=2*pi*v/1000;
中间用“,”号，别用“;”号，最后的“k0=2*pi*v/1000”对吗？v可不是常 ...

实在是不好意思，我慢慢改去。

赞一下

回复此楼

5楼2015-07-01 18:18:49

已阅关注TA 给TA发消息送TA红花 TA的回帖

yj112358

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 82.4
散金: 16
帖子: 96
在线: 103.5小时
虫号: 2837253
注册: 2013-11-29
专业: 胶体与界面化学

引用回帖:

大侠你看这下改的怎么样了，想通过拟合求参数 kmax, d, sild,nmax。

Title "Type your title here";
ComplexStr = i;
Parameters kmax, d, sild,nmax;
constant n1=1,
c=3*(10^8),
alpha=pi/2,
fa1=pi*36/180,n2=1.402,k2=0.0157,v0=2923,fw=16.50098,y=2*pi*c*fw;
Variable v, ra ;
ConstStr f=0.5*(3*(cos(sild))^2-1);
delta=2*pi*c*(v-v0);
k0=2*pi*v/1000;
kx_max=kmax*(f*(sin(alpha))^2/2+(1-f)/3);
kz_max=kmax*(f*(cos(alpha))^2+(1-f)/3);
nx=nmax-2*y*kx_max*delta/(4*delta^2+y^2);
nz=nmax-2*y*kz_max*delta/(4*delta^2+y^2);
kx=kx_max*y^2/(4*delta^2+y^2);
kz=kz_max*y^2/(4*delta^2+y^2);
n_x=nx+kx*i;
n_z=nz+kz*i;
n_2=n2+k2*i;
I1=d*((n_x)^2-(n_2)^2);
I2=((n_z)^2-(n_2)^2)*d/(n_z)^2;
I3=i*k0*sin(fa1)*I1/(n_2);
fa2=asin((n_z)/(n_2)*sin(fa1));
b=I3-sin(fa2-fa1);
b1=sin(fa1+fa2)-I3;
rs=b/b1;
rf=sin(fa2-fa1)/sin(fa2+fa1);
(rf_1)=rf*conj(rf);
(rs_1)=rs*conj(rs);
Function ra=-log10(rs_1/rf_1) ;

赞一下

回复此楼

6楼2015-07-01 19:07:40

已阅关注TA 给TA发消息送TA红花 TA的回帖

yj112358

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 82.4
散金: 16
帖子: 96
在线: 103.5小时
虫号: 2837253
注册: 2013-11-29
专业: 胶体与界面化学

引用回帖:

7楼: Originally posted by dingd at 2015-07-01 21:20:35
Title "Type your title here";
ComplexStr = i;
Parameters kmax, d, sild,nmax;
constant n1=1,
c=3*(10^8),
alpha=pi/2,
fa1=pi*36/180,n2=1.402,k2=0.0157,v0=2923,fw=16.50098,y=2*pi*c*fw; ...

太感谢你了！不过这里要求参数都正数？不知道能不能限定？
拜托了！

赞一下

回复此楼

8楼2015-07-01 21:33:50

已阅关注TA 给TA发消息送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

问问题是提前说清。有时间自己慢慢算，结果也许会更好。

均方差(RMSE):2.28871243744847E-5
残差平方和(SSE):4.08579960463844E-8
相关系数(R): 0.995627394072845
相关系数之平方(R^2): 0.991273907828284
确定系数(DC): 0.985972288752503
F统计(F-Statistic): 1758.42418135516

参数最佳估算
-------------------- -------------
kmax 42.1993762794446
d 6.1139360393898E-9
sild 0.0544860772094876
nmax 1.40553526345549

赞一下

回复此楼

9楼2015-07-01 22:21:59

已阅关注TA 给TA发消息送TA红花 TA的回帖

yj112358

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 82.4
散金: 16
帖子: 96
在线: 103.5小时
虫号: 2837253
注册: 2013-11-29
专业: 胶体与界面化学

引用回帖:

9楼: Originally posted by dingd at 2015-07-01 22:21:59
问问题是提前说清。有时间自己慢慢算，结果也许会更好。

均方差(RMSE):2.28871243744847E-5
残差平方和(SSE):4.08579960463844E-8
相关系数(R): 0.995627394072845
相关系数之平方(R^2): 0.991273907828284
...

太感谢大侠了！我想问一下这个带有复数的函数对1stopt版本有多大的要求啊？1.5行不？谢谢

赞一下

回复此楼

10楼2015-07-01 22:41:45

已阅关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 yj112358 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 301求调剂 +13	121. 2026-04-04	13/650	2026-04-06 14:20 by lqwchd
[考研] 346分的生物与医药08600求调剂 +6	常雨阳上岸 2026-04-05	7/350	2026-04-06 12:36 by lys0704
[考研] 材料与化工371求调剂 +14	陪琳看海 2026-04-04	15/750	2026-04-06 06:59 by houyaoxu
[考研] 找调剂 +10	楚乔乔 2026-04-01	10/500	2026-04-05 22:19 by syh9288
[考研] 材料调剂 +18	一样YWY 2026-04-02	19/950	2026-04-04 22:14 by hemengdong
[考研] 0703求调剂 +6	zizimo 2026-03-31	6/300	2026-04-04 14:16 by 无际的草原
[考研] 土木304求调剂 +4	兔突突突， 2026-03-31	4/200	2026-04-04 13:34 by 1753564080
[考研] 265求调剂 +17	林深温澜 2026-04-01	20/1000	2026-04-04 01:09 by userper
[考研] 总分328生物与医药考数学求调剂 +7	aaadim 2026-04-02	9/450	2026-04-03 22:53 by syh9288
[考研] 数一英一285求调剂 +7	AZMK 2026-04-03	9/450	2026-04-03 13:03 by ms629
[考研] 315分 085602 求调剂 +15	26考研上岸版26 2026-04-02	15/750	2026-04-03 12:45 by xingguangj
[考研] 273求调剂 +20	李芷新1 2026-03-31	20/1000	2026-04-03 09:58 by linyelide
[考研] 生物学308分求调剂（一志愿华东师大） +6	相信必会光芒万� 2026-03-31	7/350	2026-04-02 23:16 by JourneyLucky
[考研] 交通运输考试264分求工科调剂 +4	jike777 2026-04-02	4/200	2026-04-02 21:53 by zllcz
[考研] 考研调剂 +3	李木子0120 2026-04-02	5/250	2026-04-02 21:45 by dongzh2009
[考研] 环境工程297分求调剂一志愿杭高院 +15	GENJIOW 2026-03-31	16/800	2026-04-02 17:56 by cyh—315
[考研] 材料专业求调剂 +10	月月鸟木 2026-04-01	10/500	2026-04-02 12:57 by wxiongid
[考研] 354求调剂 +4	lxb598 2026-03-31	5/250	2026-04-02 09:55 by Jaylen.
[考研] 土木304求调剂 +6	兔突突突， 2026-03-31	7/350	2026-04-02 09:06 by coolminer
[考研] 085602化学工程268分蹲调剂 +8	月照花林。 2026-04-01	8/400	2026-04-01 22:08 by 无际的草原

24小时热门版块排行榜

[求助] 求大神帮我用1stopt拟合带有复数的函数 已有1人参与

» 本帖@通知

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

【答案】应助回帖

【答案】应助回帖

【答案】应助回帖

[求助] 求大神帮我用1stopt拟合带有复数的函数已有1人参与