24小时热门版块排行榜    

查看: 849  |  回复: 0

小朱w

铁虫 (初入文坛)

[求助] 关于用MATLAB拟合求解多元方程组未知参数

各位前辈好:现有如下问题,x、y、z为已知量,怎么拟合求未知参数a11,a12,a22,b11,b12,b22,c11,c12,c22,d11,d12,d22和f(i)、p(i)、m(i).

已知数据如下:
x=[ -0.0008  -0.0010  -0.0013   -0.0014   -0.0016 -0.0017   -0.0018  -0.0018   -0.0017   -0.0016  -0.0014  -0.0012];
y=[ -0.0288-0.0336 -0.0378 -0.0413 -0.04399 -0.0458  -0.0467 -0.0466      -0.0455 -0.04354  -0.0407 -0.0371];
z=[-0.0593   -0.0687   -0.0769  -0.0836   -0.0887 -0.0920   -0.0934  -0.0929   -0.0905   -0.0862 -0.0803   -0.0728];

方程组如下:
N=12;
for i=1:N   % 3N个方程
  0  = a11*y(i)+a12*z(i) + p(i)*(b11*y(i)+b12*z(i))-f(i);
  0= a12*y(i)+a22*z(i) +p(i)*(b12*y(i)+b22*z(i))-m(i);
  0= t^2*p(i)/12 - (c11*y(i)^2+2*c12*y(i)*z(i)+c22*z(i)^2)-p(i)*(d11*y(i)^2+2*d12*y(i)*z(i)+d22*z(i)^2)-x(i);
end
thet(1)= 0; xx(1) = 1+x(1); yy(1) = y(1);
for i=2:N
  thet(i) = sum(z(1:i-1));
  xx(i)  = cos(thet(i))*(1+x(i))-sin(thet(i))*y(i);
  yy(i)  = sin(thet(i))*(1+x(i))+cos(thet(i))*y(i);%
  %3(N-1)个方程
  0= cos(thet(i))*p(i)-sin(thet(i))*f(i) - p(1);
  0= cos(thet(i))*f(i)+sin(thet(i))*p(i)- f(1);
  0= m(i)+ (1+x(i))*f(i) - y(i)*p(i) - m(i-1);
end
希望各位学哥、学姐,大侠们能给小弟指点一下!!跪谢!!!不胜感激!!!!!1真的把小弟难住了。。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小朱w 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 07:41 by OEbVnUOu01ol
[教师之家] 跳槽后在研项目怎么办? +4 简单化xn 2026-08-22 8/400 2026-08-23 07:20 by ilovexiaobin
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:43 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见