| 查看: 944 | 回复: 4 | ||||
[求助]
含有参变量的非线性方程组拟合，1stOpt无法运行，求各位大神帮我看一下程序 已有1人参与
|
原公式为 需要拟合得到K11, K12, K21（t为参变量） 数据为： x1 x2 y 0.10 5.3111E-05 0.0383 0.20 2.10486E-05 0.0205 0.30 1.16083E-05 0.0153 0.40 7.1922E-06 0.0124 0.50 4.93786E-06 0.0114 0.60 3.78207E-06 0.0106 在1stOpt 5.5版本中程序如下，但是无法运行，不知道是什么原因，还请大神帮我修改一下，谢谢哈： parameter k11,k12,k21; parameter t[0,1]; variable x1,x2,y; t=x1-x2*(t*k11+k21*t*x2+2*k12*t^2); y=x2/x1*(k11*t+2*k21*t*x2+k12*t^2); data; 0.10 5.3111E-05 0.0383 0.20 2.10486E-05 0.0205 0.30 1.16083E-05 0.0153 0.40 7.1922E-06 0.0124 0.50 4.93786E-06 0.0114 0.60 3.78207E-06 0.0106  搜狗截图15年08月05日1739_1.png |
回复此楼» 猜你喜欢
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有130人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有252人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有10人回复
-
散金币祈福
已经有94人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 申报生产常见问答--国家局
已经有37人回复
- CDE共性问题解答解答20120725
已经有26人回复
2楼2015-08-05 17:46:28
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
Charles201: 金币+50, ★★★★★最佳答案, 谢谢版主，非常有帮助 2015-08-06 19:17:26
感谢参与，应助指数 +1
Charles201: 金币+50, ★★★★★最佳答案, 谢谢版主，非常有帮助 2015-08-06 19:17:26
|
t=x1-x2*(t*k11+k21*t*x2+2*k12*t^2)其中t是隐式的，不过这是关于t的二次方程组，可求得 t=(1/4)*(-k21*x2^2-k11*x2-1+sqrt(k21^2*x2^4+2*k11*k21*x2^3+k11^2*x2^2+8*k12*x1*x2+2*k21*x2^2+2*k11*x2+1))/(k12*x2)； 记为（1） t=-(1/4)*(k21*x2^2+k11*x2+sqrt(k21^2*x2^4+2*k11*k21*x2^3+k11^2*x2^2+8*k12*x1*x2+2*k21*x2^2+2*k11*x2+1)+1)/(k12*x2)； 记为（2） 对于式（1） 计算结果为： 均方差(RMSE): 0.00035314293391921 残差平方和(SSE): 7.48259590662405E-7 相关系数(R): 0.999330443819298 相关系数之平方(R^2): 0.998661335944076 决定系数(DC): 0.998652436046244 卡方系数(Chi-Square): 2.83851550549827E-5 F统计(F-Statistic): 2984.05363623455 参数 最佳估算 ---------- ------------- k11 -78.3805841562086 k12 4785.45080564761 k21 3627201.14766568 ====== 结果输出 ===== No 实测值y 计算值y 1 0.0383 0.0383270 2 0.0205 0.0201765 3 0.0153 0.0157170 4 0.0124 0.0128832 5 0.0114 0.0110885 6 0.0106 0.0102280 对于式（2）则没有好的拟合结果。 |
3楼2015-08-06 09:04:53
|
请求帮忙跑一下代码，拜托了！ Title "math"; Parameters w,j; Variable t,d; Function d=(cos(pi/180*10.6305)*sin((t+(j-116.391389)/15-12)*15*pi/180))/(sin(pi/180*w)*cos(pi/180*10.6305)*cos((t+(j-116.391389)/15-12)*15*pi/180)-cos(pi/180*w)*(pi/180*10.6305)); data; //t d 14.7 0.479788 14.75 0.470044 14.8 0.460681 14.85 0.451726 14.9 0.443062 14.95 0.434765 15 0.426655 15.05 0.418884 15.1 0.411403 15.15 0.404113 15.2 0.397062 15.25 0.390132 15.3 0.383499 15.35 0.376979 15.4 0.370675 15.45 0.364494 15.5 0.358438 15.55 0.352565 15.6 0.34683 15.65 0.341217 15.7 0.335668 ; |
4楼2015-09-12 14:28:34
5楼2015-09-13 21:43:09
40