24小时热门版块排行榜

返回列表

tj郭楠

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 20
帖子: 17
在线: 5.3小时
虫号: 5944925
注册: 2017-03-12
性别: GG
专业: 金属材料的力学行为

[求助] 新人求助，matlab拟合隐函数参数

已知隐函数 $f(x,y,z)$ ，形状是一个椭球，x，y，z有9组观测值，参数有8个，α1-α8；求如何拟合出这8个参数
最近一直在思考这个，无奈本人编程基础太差，无从下手，求各位大神给个思路；
如果需要，我把 $f(x,y,z)$ 发出来（实在太长了）。

回复此楼

» 猜你喜欢

售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急已经有8人回复
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有9人回复
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急已经有8人回复
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急已经有8人回复
各位基础老师们，看完中央《从严治党》的决定后，对学术界更加有信心了吧已经有5人回复
2027博士申请-临床检验诊断学、病原生物学、微生物方向已经有5人回复

1楼 2018-05-30 22:59:35

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

smitest

木虫 (小有名气)

应助: 4 (幼儿园)
金币: 2825.7
红花: 4
帖子: 191
在线: 20小时
虫号: 992136
注册: 2010-04-08
专业: 原子和分子物理

发出来看看

发自小木虫Android客户端

回复此楼

2楼2018-05-31 12:59:03

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

tj郭楠

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 20
帖子: 17
在线: 5.3小时
虫号: 5944925
注册: 2017-03-12
性别: GG
专业: 金属材料的力学行为

引用回帖:

2楼: Originally posted by smitest at 2018-05-31 12:59:03
发出来看看

F=@(p,x)(0.5...
*(abs(...
0.5*(p(1)*(2*x(:,1)-x(:,2))/3.0+p(2)*(2*x(:,2)-x(:,1))/3.0)+0.5*sqrt((p(1)*(2*x(:,1)-x(:,2))/3.0-p(2)*(2*x(:,2)-x(:,1))/3.0).^2+4*p(7)*x(:,3).^2)...%a1
-0.5*(p(1)*(2*x(:,1)-x(:,2))/3.0+p(2)*(2*x(:,2)-x(:,1))/3.0)-0.5*sqrt((p(1)*(2*x(:,1)-x(:,2))/3.0-p(2)*(2*x(:,2)-x(:,1))/3.0).^2+4*p(7)*x(:,3).^2)...%a2
).^8.0...%a
+abs(...
((p(3)-4*(p(4)+p(5)-p(6)))*x(:,2)+2*(4*p(5)-p(6)+p(4)-p(3))*x(:,1))/9.0+((p(6)-4*(p(4)+p(5)-p(3)))*x(:,1)+2*(4*p(4)-p(6)+p(5)-p(3))*x(:,2))/9.0+sqrt(( ((p(3)-4*(p(4)+p(5)-p(6)))*x(:,2)+2*(4*p(5)-p(6)+p(4)-p(3))*x(:,1))/9.0-((p(6)-4*(p(4)+p(5)-p(3)))*x(:,1)+2*(4*p(4)-p(6)+p(5)-p(3))*x(:,2))/9.0).^2+4*p(8)*x(:,3).^2)...%2*b1
+0.5*(((p(3)-4*(p(4)+p(5)-p(6)))*x(:,2)+2*(4*p(5)-p(6)+p(4)-p(3))*x(:,1))/9.0+((p(6)-4*(p(4)+p(5)-p(3)))*x(:,1)+2*(4*p(4)-p(6)+p(5)-p(3))*x(:,2))/9.0-sqrt(( ((p(3)-4*(p(4)+p(5)-p(6)))*x(:,2)+2*(4*p(5)-p(6)+p(4)-p(3))*x(:,1))/9.0-((p(6)-4*(p(4)+p(5)-p(3)))*x(:,1)+2*(4*p(4)-p(6)+p(5)-p(3))*x(:,2))/9.0).^2+4*p(8)*x(:,3).^2))...%b2
).^8.0...
+abs(...
0.5*(((p(3)-4*(p(4)+p(5)-p(6)))*x(:,2)+2*(4*p(5)-p(6)+p(4)-p(3))*x(:,1))/9.0+((p(6)-4*(p(4)+p(5)-p(3)))*x(:,1)+2*(4*p(4)-p(6)+p(5)-p(3))*x(:,2))/9.0+sqrt(( ((p(3)-4*(p(4)+p(5)-p(6)))*x(:,2)+2*(4*p(5)-p(6)+p(4)-p(3))*x(:,1))/9.0-((p(6)-4*(p(4)+p(5)-p(3)))*x(:,1)+2*(4*p(4)-p(6)+p(5)-p(3))*x(:,2))/9.0).^2+4*p(8)*x(:,3).^2))...%b1
+((p(3)-4*(p(4)+p(5)-p(6)))*x(:,2)+2*(4*p(5)-p(6)+p(4)-p(3))*x(:,1))/9.0+((p(6)-4*(p(4)+p(5)-p(3)))*x(:,1)+2*(4*p(4)-p(6)+p(5)-p(3))*x(:,2))/9.0-sqrt(( ((p(3)-4*(p(4)+p(5)-p(6)))*x(:,2)+2*(4*p(5)-p(6)+p(4)-p(3))*x(:,1))/9.0-((p(6)-4*(p(4)+p(5)-p(3)))*x(:,1)+2*(4*p(4)-p(6)+p(5)-p(3))*x(:,2))/9.0).^2+4*p(8)*x(:,3).^2)...%2*b2
).^8.0...
)...%b
).^1.0/8.0;