24小时热门版块排行榜    

查看: 952  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

tj郭楠

新虫 (初入文坛)

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

已知隐函数,形状是一个椭球,x,y,z有9组观测值,参数有8个,α1-α8;求如何拟合出这8个参数
最近一直在思考这个,无奈本人编程基础太差,无从下手,求各位大神给个思路;
如果需要,我把发出来(实在太长了)。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tj郭楠

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by smitest at 2018-05-31 20:48:18
这也太复杂了,啥模型啊
...

是一个材料屈服准则
6楼2018-05-31 22:24:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

smitest

木虫 (小有名气)

2楼2018-05-31 12:59:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tj郭楠

新虫 (初入文坛)

引用回帖:
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;
3楼2018-05-31 20:25:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tj郭楠

新虫 (初入文坛)

是不是发在计算模拟板块会比较好啊?
4楼2018-05-31 20:32:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见