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

登录注册

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

24小时热门版块排行榜

返回列表

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

» 猜你喜欢

PRB投稿，Acknowledgment sent to author状态十天了已经有3人回复
PVA溶解已经有2人回复
物理学I论文润色/翻译怎么收费? 已经有217人回复
基金申请已经有44人回复
CSC与新西兰维多利亚大学PhD奖学金项目已经有0人回复
新西兰Robinson研究所招聘CSC公派访问人员已经有0人回复
帮我的英语口语老师找学生已经有0人回复
什么时候开奖？已经有13人回复
散金币祈福已经有94人回复
青基已中已经有1人回复

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

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 的主题更新

返回列表

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

[基金申请] 投票: 有多少人是今天查系统知道结果的？ +15	爱看书的可乐 2026-08-26	17/850	2026-08-27 02:26 by xxniao123
[基金申请] 我不理解！ +15	Edward_pc 2026-08-26	23/1150	2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7	process2012 2026-08-23	10/500	2026-08-26 19:23 by hmhminy
[基金申请] 怎么查啊 +5	huang1991js 2026-08-26	5/250	2026-08-26 16:09 by Equinoxhua
[基金申请] 2026国自然函评费到账 +22	羊腰板 2026-08-21	25/1250	2026-08-26 15:00 by zuocuiping
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25	木水思豆 2026-08-25	28/1400	2026-08-26 14:53 by draco1987
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +4	yuleib84 2026-08-26	5/250	2026-08-26 13:40 by yuleib84
[基金申请] 科研孤儿太难了 +21	我4大白菜 2026-08-20	22/1100	2026-08-26 11:29 by nlgza
[基金申请] 国际合作可查了，中了面上 +18	Ldrop2023 2026-08-26	18/900	2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4	yanglien 2026-08-26	5/250	2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3	paperzjh 2026-08-26	3/150	2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了，明天出结果吗 +19	angus9576 2026-08-25	23/1150	2026-08-26 10:03 by zp519
[基金申请] 牛来！米来！面来！ +8	beefly 2026-08-26	8/400	2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6	ziyangfang 2026-08-25	9/450	2026-08-25 20:26 by huagongfeihu
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17	医学老男孩 2026-08-20	22/1100	2026-08-25 15:36 by 医学老男孩
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +3	苏东坡二世 2026-08-23	3/150	2026-08-25 10:35 by shisan1313
[基金申请] 我面上完蛋了 +13	且听虎啸 2026-08-20	14/700	2026-08-25 09:10 by mrkang
[教师之家] 跳槽后在研项目怎么办？ +5	简单化xn 2026-08-22	10/500	2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳今天，20号变了 +5	archvillain 2026-08-20	5/250	2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了？ +8	chengyan1220 2026-08-21	11/550	2026-08-21 17:52 by dcqxinyang

24小时热门版块排行榜

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

» 本帖@通知

» 猜你喜欢

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

【答案】应助回帖

【答案】应助回帖

【答案】应助回帖

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