| 查看: 951 | 回复: 2 | |||
wtisunday木虫 (正式写手)
|
[求助]
MATLAB 非线性拟合 fminsearch中包含非线性方程的求解问题
|
|
format short global a1 a2 a3 a4 a5 a6 a7 a8 b1 b2 b3 b4 b5 b6 b7 b8 k0=[[3 3 3 2 2 2 1 1];[1 1 1 1 1 1 1 1]]'; %参数的初始值 x0=[1,0]'; %拉伸应力初始值 t=[0.00423,0.01621,0.06216,0.23841,0.91437,3.50692,14.6856,101.343]'; G=[336370,64246.4,44849.7,10649.7,5151.99,1224.43,63.4524,2.67362]'; %relaxation time,relaxation modulus实验数据[t,G] xexp=[0.03162 0.1 0.31623 1 3.16228 10.5635 20.395]'; yexp=[84834.3 84834.3 79215.4 62322.2 52509.5 32501.6 25570.4]'; %xexp,yexp实验数据[频率x、拉伸粘度y] [k,fval]=fminsearch(@ObjFunc,k0,x0,xexp,yexp);%非线性拟合 function f=ObjFunc(x0,xexp,yexp) G=[336370,64246.4,44849.7,10649.7,5151.99,1224.43,63.4524,2.67362]'; t=[0.00423,0.01621,0.06216,0.23841,0.91437,3.50692,14.6856,101.343]'; [x]=fsolve(@kineticseqs103,x0); ff=sum(x(:,2)-x(:,1))./xexp; f=(sum(ff'-yexp))^2; function pp=rheolomodel(k,x,g,xexp) global a1 a2 a3 a4 a5 a6 a7 a8 b1 b2 b3 b4 b5 b6 b7 b8 k=[[a1,a2,a3,a4,a5,a6,a7,a8];[b1,b2,b3,b4,b5,b6,b7,b8]]'; q=1+(x(:,1)+2*x(:,2))/3./G.^0.5; u=2.*k(:,1)*exp(2/k(:,2)-1)*(1-1/q.^2)+1/q.^2*(1-0.3/k(:,2).*... ((x(:,1)).^2+2*(x(:,2)).^2)/3/g.^2); pp1=u.*x(:,1)-2*t.*xexp.*(x(:,1)+G)+G.*(u-1)+... 0.3./k(:,2)./G.*(x(:,1))^2; pp2=u.*x(:,2)+t.*xexp.*x(:,2)+G.*(u-1)+... 0.3./k(:,2)./G.*(x(:,2))^2+t.*G.*xexp; pp=[pp1,pp2]'; 以上是我依照前人经验编的程序，由于是新手，在程序出错和修改上解决不了，还望各位虫友帮帮忙啊！ |
回复此楼» 猜你喜欢
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
物理学I论文润色/翻译怎么收费?
已经有67人回复
-
桂林理工大学物理学专业招收调剂
已经有18人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
-
【方案分享】单环磁场+轴心控制+偏转导出电子束约束系统（可行性实验）
已经有6人回复
-
【修正版】单环用磁约束低速电子实验方案（简化版）
已经有0人回复
-
桂林理工大学物理学专业招收调剂，还有三个名额！！！
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 合同和编制的啥区别 求解
已经有10人回复
- 关于一级A标准执行问题的一点看法
已经有8人回复
- 2012标书格式问题
已经有3人回复
- 本科毕业生求解
已经有9人回复
- 关于测ee值的问题
已经有20人回复
- matlab非线性参数拟合问题
已经有7人回复
- 请教一个数据非线性拟合的问题
已经有11人回复
- 求助：非线性方程系数回归问题？
已经有10人回复
- 关于origin里面的非线性拟合的问题
已经有6人回复
- 求助origin高手 非线性拟合问题
已经有6人回复
- 【求助】求解非线性方程
已经有8人回复
- 请教关于origin中自定义函数非线性回归拟合问题！
已经有10人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2012-02-16 08:23:48
wtisunday
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1742.3
- 散金: 30
- 红花: 1
- 帖子: 937
- 在线: 568.8小时
- 虫号: 444484
- 注册: 2007-10-27
- 性别: GG
- 专业: 高分子材料的加工与成型
|
我的问题是多模态问题 就是已知t的i个分量和G的i个分量，通过两个方程求出x的i个分量和y的i个分量，再通过x和y分量和建立优化函数，通过优化实现a的2i个分量的确定。 比如 i=3 u(i)=a(i)*exp(a(i+4)-t(i))+1/(x(i)+y(i))/G(i).....中间变量表达式 f1=u(i)*x(i)-t(i)*(x(i)+y(i))+g(i)*((u(i)-1)+expdata; f2=u(i)*y(i)+t(i)*y(i)+y(i)^(2*xexpdata); f=[f1,f2]';......分量方程组（每一分量均有此关系） Y=sum(x(i)-y(i))/xexpdata.....计算yexpdata Op=sum(Y-yexpdata)^2......优化函数 |
3楼2012-02-16 08:51:17
20