24小时热门版块排行榜    

CyRhmU.jpeg
查看: 743  |  回复: 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 的主题更新
信息提示
请填处理意见