| 查看: 1632 | 回复: 5 | |||
[交流]
拟合非线性方程系数的精度问题
|
|||
|
本人想拟合一个非线性方程的系数(4个未知数),精度要求很高,期望要到0.01mm。用MATLAB中优化函数lsqnonlin、fmincon、fminunc等函数和遗传算法计算时还得不到满意的精度.不知道上述提到的计算方法大家在实际工程应用中都能达到多少精度啊!并且为了能达到我想要的精度还有什么好的方法不? 目标函数为求和sum(i=1,6,z(i)),z(i)=m(i)^2+a*m(i)+b*m(i)*sin(w(i))+c*sin(w(i))+d*m(i)*cos(w(i))+e*cos(w(i))+f,求它的最小值,接近与0.其中a=-2*x(1); b=2*(x(4)-x(2));c=2*(x(3)*x(2)-x(4)*x(1));d=2*(x(1)-x(3)); e=2*(x(3)*x(1)+x(2)*x(4)-x(1)^2-x(2)^2); f=2*(x(1)^2+x(2)^2)+x(3)^2+x(4)^2-2*x(2)*x(4)-2*x(1)*x(3)-(R+r)^2;R=50;r=1.5 w=n*pi/180 m=[41.76554 21.16515 12.83882 10.98004 14.71395 24.88945] n=[0 2.5 5 7.5 10 12.5] x(1)~x(4)为所求未知数,理想解为x(1)=422.443 ;x(2)=20 ;x(3)=62.918 ;x(4)=46.95,结果精度要求0.01 麻烦各位大侠了! [ Last edited by myb19880127 on 2012-3-31 at 10:43 ] |
» 猜你喜欢
070200求调剂,一志愿某211,288分
已经有20人回复
法国博士后职位
已经有0人回复
物理学I论文润色/翻译怎么收费?
已经有111人回复
重庆交大26年硕士生招生拟调剂通知已出!欢迎加入机器视觉与3D光学成像课题组。
已经有0人回复
**
已经有1人回复
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
基底STO,薄膜SRO,XRD里面的振荡,是laue震荡还是kiessig振荡? 怎么判断?
已经有2人回复
B区学生调剂-兰州交通大学材料科学与工程学院
已经有8人回复
山西大同大学物理学专业还有调剂名额,欢迎调剂!
已经有9人回复
» 本主题相关价值贴推荐,对您同样有帮助:
在matlab中拟合曲面的表达式
已经有29人回复
求助:非线性方程系数回归问题?
已经有10人回复
求非线性方程中的参数问题
已经有13人回复
求助关于线性拟合的
已经有1人回复
» 抢金币啦!回帖就可以得到:
邀您投稿 Polymers 特刊-高分子材料加工主题
+1/87
上海应用技术大学刘敏课题组招生---欢迎材料,物理和化学相关背景
+1/46
B区学生调剂-兰州交通大学材料科学与工程学院
+1/42
大连工业大学纺织科学与工程接收硕士调剂
+1/40
华南理工大学 基础医学学硕和临床医学学硕调剂
+1/20
江苏师范大学招收化学及材料专业研究生
+1/19
河北大学燕麦基因组学与分子育种课题组招收2026级考研调剂生
+1/18
复试调剂-招收调剂生(一志愿07开头专业),海洋生物专业4个名额
+1/13
大连大学-贵州省煤炭洁净利用重点实验室联合培养研究生 化学5人+环境工程7人
+1/10
~【西昌学院资源与环境学院硕士招生调剂】~
+1/9
湖南理工大学调剂系统已开。 化学、化工、材料化工 调剂指标充足,录取率高,欢迎调剂
+1/8
河南师大全国重点实验室—0780药学、08600生物与医药仍有调剂指标!
+1/8
江西科技师范大学申亮组招研究生
+1/6
贵州财经大学新一代电子信息技术(含量子技术等)硕士专业(085401)接收调剂
+1/5
福建理工大学材料学院招收专业代码08开头学术型硕士研究生
+1/4
河南理工大学化学化工学院招收材料、化工类专业调剂研究生
+1/4
招收化学、海洋科学、材料、环境专业调剂
+2/4
招07及08开头化学、材料、环境、食品、生物与医药等调剂生
+2/2
青岛农业大学化学与药学院招收调剂,欢迎报名!
+1/2
招收荧光探针/有机合成/生物学博士生
+1/1
|
目标函数为z=m(i)^2+a*m(i)+b*m(i)*sin(w(i))+c*sin(w(i))+d*m(i)*cos(w(i))+e*cos(w(i))+f,求它的最小值,接近与0.其中a=-2*x(1); b=2*(x(4)-x(2));c=2*(x(3)*x(2)-x(4)*x(1));d=2*(x(1)-x(3)); e=2*(x(3)*x(1)+x(2)*x(4)-x(1)^2-x(2)^2); f=2*(x(1)^2+x(2)^2)+x(3)^2+x(4)^2-2*x(2)*x(4)-2*x(1)*x(3)-(R+r)^2;R=50;r=1.5 w=n*pi/180 m=[41.76554 21.16515 12.83882 10.98004 14.71395 24.88945] n=[0 2.5 5 7.5 10 12.5] x(1)~x(4)为所求未知数,理想解为x(1)=422.443 ;x(2)=20 ;x(3)=62.918 ;x(4)=46.95,结果精度要求0.01 麻烦各位大侠了! [ Last edited by myb19880127 on 2012-3-31 at 10:38 ] |
3楼2012-03-31 09:57:53
2楼2012-03-30 16:53:41
4楼2012-03-31 10:23:11
5楼2012-03-31 10:36:37
★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
dbb627: 金币+2, 谢谢应助! 2012-04-01 21:19:08
小木虫: 金币+0.5, 给个红包,谢谢回帖
dbb627: 金币+2, 谢谢应助! 2012-04-01 21:19:08
|
你的目标函数z应该取绝对值吧?否则最小值远远小于0。 用1stOpt求解: 目标函数值(最小): 0.288024770736229 x1: 422.733761186019 x2: 19.3667780059412 x3: 62.9959407161319 x4: 46.9203589652975 将你的理想解(x(1)=422.443 ;x(2)=20 ;x(3)=62.918 ;x(4)=46.95)代入验证,目标函数值(最小)=5.506768407 从理论上看,结果已经好于你的理想值了。 |
6楼2012-03-31 10:49:49













回复此楼