24小时热门版块排行榜

返回列表

763443895

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 155.4
散金: 35
帖子: 49
在线: 30.5小时
虫号: 2683946
注册: 2013-09-26
性别: GG
专业: 传动机械学

[求助] 变步长的非线性方程组求解

function y=t1(x)
%%%%%%%中间变量%%%%%%%%%%%%%%%%%%%
M=[cosd(x(1)-x(2))+cosd(x(1))];
N=[sind(x(1)-x(2))+sind(x(1))];
k=x(1)/sind(34.32);
q=[cosd(x(4)+x(5)-63.77)+cosd(x(4)+24.81)];
w=[sind(x(4)+x(5)-63.77)+sind(x(4)+24.81)];
e=x(4)/sind(48.743);
g=x(1)-x(2)+66.793;
h=x(4)+x(5)-63.77;
%%%%%%% 6个非线性方程组  %%%%%%%%%%%
y(1)=M+cosd(x(6))*[-sind(e)*q+cosd(e)*w]-sind(x(6))*[cosd(e)*q+sind(e)*w];
y(2)=cosd(x(3))*[sind(k)*M-cosd(k)*N]-sind(x(3))*[-cosd(k)*M-sind(k)*N]-q;
y(3)=sind(x(3))*[sind(k)*M-cosd(k)*N]+cosd(x(3))*[-cosd(k)*M-sind(k)*N]-sind(x(6))*[-sind(e)*q+cosd(e)*w]-cosd(x(6))*[cosd(e)*q+sind(e)*w];
y(4)=cosd(g)+cosd(x(6))*[sind(h)*sind(e)+sind(h)*cosd(e)]-sind(x(6))*[sind(e)*sind(h)+cosd(e)*cosd(h)];
y(5)=cosd(x(3))*(cosd(k)*sind(g)-sind(k)*cosd(g))+sind(x(3))*(sind(g)*sind(k)+cosd(g)*cosd(k))-cosd(h);
y(6)=cosd(x(3))*(sind(k)*sind(g)+cosd(k)*cosd(g))-sind(x(3))*(sind(g)*cosd(k)-sind(g)*cosd(k))-sind(x(6))*(sind(e)*cosd(h)+sind(h)*cosd(e))
      -cosd(x(6))*(sind(h)*sind(e)+cosd(e)*cosd(h));

6个未知数，假设已知6个未知数的一组解为[3,4,5,6,3,4],
想让x(3)按照步长为1求解其他5个参数！！！
主程序该怎么写，新手，求大神帮忙，
已经知道用fsolve进行求解，可是结果不是想要的
结果格式应该是:
x(1) 3
x(2) 4
x(3) 5 4 3 2 1
x(4) 6
x(5) 3
x(6) 4

回复此楼