| 查看: 1071 | 回复: 3 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
1stOpt函数中含有间接变量怎么办
|
|||
刚接触1stOpt就遇到了很大的困难：由于函数表达式过于复杂，不得不在其中引入间接变量（如代码中的F11,F12等等），而运行结果却并不是设定的parameter。请高手指点迷津，不胜感激 结果如下： 参数 最佳估算 ---------- ------------- x1 0 x2 0 x3 0 rate1 0 rate2 0 rate3 0 r1 0 r2 0 r3 0 rs 0 代码如下： Constant K=23.05,u=21.58,c1=4.6948,cs=0.00001,phi=0.1412,pi=3.1415926; Parameters x1,a1,x2,a2,x3,a3; Variable p,c; Function x1+x3+x3=phi; Ia1=2*pi*a1*(arccos(a1)-a1*(1-a1^2)^0.5)/(1-a1^2)^1.5; F11=6*u*Ia1/(2*pi*(3*K+4*u)); F12=6*u*(3*Ia1-4*pi)/(4*pi*(3*K+4*u)); F13=a1^2*(3-9*Ia1/(4*pi))*(6*K+2*u)/(2*(1-a1^2)*(3*K+4*u))+(6*u*Ia)/(8*pi*(3*K+4*u)); F14=0.5*((3-9*Ia1/(4*pi))*(6*K+2*u)/(2*(1-a1^2)*(3*K+4*u))-3*u*(1-Ia1/pi)/(3*K+4*u)); rate1=-p/(F11-F12*F13/(F13+F14))/K; //if rate1 <-1 // rate1=-1 //end Ia2=2*pi*a2*(arccos(a2)-a2*(1-a2^2)^0.5)/(1-a2^2)^1.5; F21=6*u*Ia2/(2*pi*(3*K+4*u)); F22=6*u*(3*Ia2-4*pi)/(4*pi*(3*K+4*u)); F23=a2^2*(3-9*Ia2/(4*pi))*(6*K+2*u)/(2*(1-a2^2)*(3*K+4*u))+(6*u*Ia)/(8*pi*(3*K+4*u)); F24=0.5*((3-9*Ia2/(4*pi))*(6*K+2*u)/(2*(1-a2^2)*(3*K+4*u))-3*u*(1-Ia2/pi)/(3*K+4*u)); rate2=-p/(F21-F22*F23/(F23+F24))/K; Ia3=2*pi*a3*(arccos(a3)-a3*(1-a3^2)^0.5)/(1-a3^2)^1.5; F31=6*u*Ia3/(2*pi*(3*K+4*u)); F32=6*u*(3*Ia3-4*pi)/(4*pi*(3*K+4*u)); F33=a3^2*(3-9*Ia3/(4*pi))*(6*K+2*u)/(2*(1-a3^2)*(3*K+4*u))+(6*u*Ia)/(8*pi*(3*K+4*u)); F34=0.5*((3-9*Ia3/(4*pi))*(6*K+2*u)/(2*(1-a3^2)*(3*K+4*u))-3*u*(1-Ia3/pi)/(3*K+4*u)); rate3=-p/(F31-F32*F33/(F33+F34))/K; L1z=1/(1-(a1*(1+rate1))^2)-(a1*(1+rate1))*arccos((a1*(1+rate1)))/(1-(a1*(1+rate1))^2)^1.5; L1x=(1-L1z)/2; L1y=L1x; R1=1/9*(1/(L1x*c1+(1-L1x)*c)+1/(L1y*c1+(1-L1y)*c)+1/(L1z*c1+(1-L1z)*c)); L2z=1/(1-(a2*(1+rate2))^2)-(a2*(1+rate2))*arccos((a2*(1+rate2)))/(1-(a2*(1+rate2))^2)^1.5; L2x=(1-L2z)/2; L2y=L2x; R2=1/9*(1/(L2x*c1+(1-L2x)*c)+1/(L2y*c1+(1-L2y)*c)+1/(L2z*c1+(1-L2z)*c)); L3z=1/(1-(a3*(1+rate3))^2)-(a3*(1+rate3))*arccos((a3*(1+rate3)))/(1-(a3*(1+rate3))^2)^1.5; L3x=(1-L3z)/2; L3y=L3x; R3=1/9*(1/(L3x*c1+(1-L3x)*c)+1/(L3y*c1+(1-L3y)*c)+1/(L3z*c1+(1-L3z)*c)); Rs=1/9*(1/(1/3*cs+(1-1/3)*c)+1/(1/3*cs+(1-1/3)*c)+1/(1/3*cs+(1-1/3)*c)); (1-(x1*(1+rate1)+x2*(1+rate2)+x3*(1+rate3)))*(cs-c)*Rs+x1*(1+rate1)*(c1-c)*R1+x2*(1+rate2)*(c1-c)*R2+x3*(1+rate3)*(c1-c)*R3=0; data; 0.008 0.1325 0.115 0.1309 0.020 0.1299 0.026 0.1290 0.040 0.1277 0.060 0.1269 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
人气不行了
已经有10人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 1stopt为何结果是 计算值全为0
已经有8人回复
- 1stopt自变量超多时是否有简单语法？AND 自变量多 数据量大的话1stopt都会算的特慢么
已经有7人回复
- 求大神帮忙用1stOpt 拟合复数数据
已经有10人回复
- 求高手给做一道题，用MATLAB或1stOpt，请把结果直接给出来，我的软件就不行，谢谢
已经有5人回复
- 1stopt提示函数字符串有误！真心求教！
已经有8人回复
- 求助--用1stopt拟合蓝宝石衬底的GaN薄膜红外反射光谱
已经有6人回复
- 1stopt解方程组
已经有12人回复
- 用1stopt软件做微分方程组的参数估计
已经有13人回复
- 用1stOpt进行复数曲线拟合时，结果出错。急求教！！！
已经有5人回复
- 这个编号后为什么1stopt没有反应，不能执行
已经有5人回复
- 哪位朋友有注册版1stOpt软件？急需
已经有5人回复
- 1stOPT1.5版每次仿真结果不一样，求助！急！
已经有3人回复
- 运算过程有复数出现，求高版本1stOpt
已经有6人回复
- 1stopt运行程序，总是出现字符函数错误的提示，请高手帮助
已经有5人回复
- 大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？
已经有8人回复
- 1stopt4阶非线性常微分方程，帮小弟看一下，跪求。
已经有7人回复
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
- 麻烦高人指点一下1stopt和spss的非线性回归中weibull的两参数的公式应该怎么表达啊？
已经有9人回复
3楼2014-08-07 21:53:42
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3881
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
2楼2014-07-31 16:20:25
feixiaolin
荣誉版主 (文坛精英)
-
专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3881
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
4楼2014-08-07 21:56:15
