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

登录注册

应《网络安全法》要求，自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

» 猜你喜欢

反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有272人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复

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

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

返回列表

24小时热门版块排行榜

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

» 本帖@通知

» 猜你喜欢

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

【答案】应助回帖

【答案】应助回帖

【答案】应助回帖

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