| 查看: 610 | 回复: 1 | |||
763443895银虫 (初入文坛)
|
[求助]
变步长的非线性方程组求解
|
|
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 |
回复此楼» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有12人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有9人回复
-
国社科系统bug了，是不是要放榜了？
已经有6人回复
-
申博发邮件
已经有6人回复
-
各位大神，目前国内有哪些比较好用的逆合成软件？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
2楼2019-08-11 15:37:20
20