| 查看: 1058 | 回复: 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]'; 以上是我依照前人经验编的程序，由于是新手，在程序出错和修改上解决不了，还望各位虫友帮帮忙啊！ |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有290人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 合同和编制的啥区别 求解
已经有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
30