《量化投资:以MATLAB为工具》

MATLAB技术论坛

 找回密码
 注册账号
查看: 697|回复: 1
收起左侧

[提问] 1stopt软件常微分非线性方程组参数拟合问题!急啊!毕不了业了!!求各位大佬帮助

[复制链接]
发表于 2017-6-15 14:31:02 | 显示全部楼层 |阅读模式
之前从来没用过1stopt,下面是自己写的程序。要拟合的参数是Cr1,Cr2,Cu1,Cu2,mu,自变量是t,因变量是xi1,xi2,lambda1,lambda2,s,现在实验数据是lambda1和s的。
Title "lashennihe";
Parameters Cr1,Cr2,Cu1,Cu2,mu;
Variable t,s,lambda1,xi1,xi2,lambda2;
ODEFunction s=Cr1*(2*lambda1-2*lambda1^(-3)*lambda2^(-2))+Cr2*(-2*lambda1^(-3)+2*lambda1*lambda2^2)+(Cu1-Cr1)*(2*lambda1*xi1^(-2)-2*lambda1^(-3)*lambda2^(-2)*xi1^2*xi2^2)+(Cu2-Cr2)*(-2*lambda1^(-3)*xi1^2+2*lambda1*lambda2^2*xi1^(-2)*xi2^(-2));
              Cr1*(2*lambda2-2*lambda2^(-3)*lambda1^(-2))+Cr2*(-2*lambda2^(-3)+2*lambda2*lambda1^2)+(Cu1-Cr1)*(2*lambda2*xi2^(-2)-2*lambda2^(-3)*lambda1^(-2)*xi1^2*xi2^2)+(Cu2-Cr2)*(-2*lambda2^(-3)*xi2^2+2*lambda2*lambda1^2*xi1^(-2)*xi2^(-2))=0;
              lambda1=1.176471*t;
              xi1'=-1/mu*((Cu1-Cr1)*(-2*lambda1^2*xi1^(-3)+2*lambda1^(-2)*lambda2^(-2)*xi1*xi2^2))+(Cu2-Cr2)*(2*lambda1^(-2)*xi1-2*lambda1^2*lambda2^2*xi1^(-3)*xi2^(-2)));
              xi2'=-1/mu*((Cu1-Cr1)*(-2*lambda2^2*xi2^(-3)+2*lambda1^(-2)*lambda2^(-2)*xi2*xi1^2))+(Cu2-Cr2)*(2*lambda2^(-2)*xi2-2*lambda1^2*lambda2^2*xi2^(-3)*xi1^(-2)));
Data;
// s lambda1
实验数据完整的在excel表格里,有6638组
6606.6        1
6840.116667        1
6883.366667        1.000044219
6778.4        1
6443.183333        1.000002212
6995.716667        1.000044219
7890.883333        1.000044219

stressstrain.xls

391 KB, 下载次数: 0, 下载积分: 贝壳 -1

数据

发表于 2017-6-15 21:31:51 | 显示全部楼层
你要拟合的公式是什么形式,还有拟合参数的变量取值范围是多少
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 注册账号

本版积分规则

QQ|网站地图|MATLAB技术论坛|Simulink仿真论坛 ( 蜀ICP备19014457号 

GMT+8, 2019-8-23 21:23 , Processed in 0.085664 second(s), 13 queries , Gzip On, MemCache On.

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表