| 查看: 594 | 回复: 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 |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
国自然面上复盘~欢迎讨论
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
基金不中，共勉
已经有5人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有3人回复
-
我不理解！
已经有23人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有10人回复
-
系统查不到
已经有10人回复
-
怎么查啊
已经有5人回复
2楼2019-08-11 15:37:20
