求助：使用matlab里的solve函数求解十个方程联立的方程组的显式解求解不出。

C JJ

已于 2024-06-22 15:36:30 修改

阅读量279

点赞数 7

分类专栏：求助文章标签： matlab

于 2024-06-22 15:36:05 首次发布

本文链接：https://blog.csdn.net/DUCK_cjj/article/details/139883267

版权

求助专栏收录该内容

1 篇文章 0 订阅

订阅专栏

求助：使用matlab里的solve函数求解十个方程联立的方程组的显式解求解不出。

应该是solve求解不了十个方程的原因。请教各位大佬解决办法。

equ = sym(zeros(10,1));

equ(1) = - z(645)*x_Helicopterpp - z(650)*phi_Helicopterpp - z(646)*theta_Helicopterpp - z(647)*psi_Helicopterpp - z(649)*alpha_xpp - z(651)*alpha_ypp + 0.5*z(653)*beta_xpp + 0.5*z(654)*beta_ypp - z(739);

equ(2) = - z(645)*y_Helicopterpp - z(659)*phi_Helicopterpp - z(656)*theta_Helicopterpp - z(657)*psi_Helicopterpp - z(658)*alpha_xpp - z(660)*alpha_ypp + 0.5*z(661)*beta_xpp + 0.5*z(662)*beta_ypp - z(740);

equ(3) = - z(645)*z_Helicopterpp - z(666)*phi_Helicopterpp - z(667)*theta_Helicopterpp - z(664)*psi_Helicopterpp - z(665)*alpha_xpp - z(668)*alpha_ypp + 0.5*z(669)*beta_xpp + 0.5*z(670)*beta_ypp - z(741);

equ(4) = - z(650)*x_Helicopterpp - z(659)*y_Helicopterpp - z(666)*z_Helicopterpp - z(673)*phi_Helicopterpp - z(680)*theta_Helicopterpp - z(674)*psi_Helicopterpp - z(678)*alpha_xpp - z(679)*alpha_ypp - z(683)*beta_xpp - z(676)*beta_ypp - z(742);

equ(5) = - z(646)*x_Helicopterpp - z(656)*y_Helicopterpp - z(667)*z_Helicopterpp - z(690)*phi_Helicopterpp - z(685)*theta_Helicopterpp - z(686)*psi_Helicopterpp - z(688)*alpha_xpp - z(689)*alpha_ypp - z(687)*beta_xpp - z(691)*beta_ypp - z(743);

equ(6) = - z(647)*x_Helicopterpp - z(657)*y_Helicopterpp - z(664)*z_Helicopterpp - z(695)*phi_Helicopterpp - z(693)*theta_Helicopterpp - z(694)*psi_Helicopterpp - z(697)*alpha_xpp - z(698)*alpha_ypp - z(696)*beta_xpp - z(699)*beta_ypp - z(744);

equ(7) = - z(649)*x_Helicopterpp - z(658)*y_Helicopterpp - z(665)*z_Helicopterpp - z(708)*phi_Helicopterpp - z(706)*theta_Helicopterpp - z(707)*psi_Helicopterpp - z(704)*alpha_xpp - z(705)*alpha_ypp - z(711)*beta_xpp - z(701)*beta_ypp - z(745);

equ(8) = - z(651)*x_Helicopterpp - z(660)*y_Helicopterpp - z(668)*z_Helicopterpp - z(718)*phi_Helicopterpp - z(716)*theta_Helicopterpp - z(717)*psi_Helicopterpp - z(714)*alpha_xpp - z(715)*alpha_ypp - z(719)*beta_xpp - z(713)*beta_ypp - z(746);

equ(9) = 0.5*z(653)*x_Helicopterpp + 0.5*z(661)*y_Helicopterpp + 0.5*z(669)*z_Helicopterpp - z(726)*phi_Helicopterpp - z(724)*theta_Helicopterpp - z(725)*psi_Helicopterpp - z(722)*alpha_xpp - z(723)*alpha_ypp - z(721)*beta_xpp - 0.5*z(729)*beta_ypp - z(747);

equ(10) = 0.5*z(654)*x_Helicopterpp + 0.5*z(662)*y_Helicopterpp + 0.5*z(670)*z_Helicopterpp - z(737)*phi_Helicopterpp - z(735)*theta_Helicopterpp - z(736)*psi_Helicopterpp - z(733)*alpha_xpp - z(734)*alpha_ypp - 0.5*z(729)*beta_xpp - z(732)*beta_ypp - z(748);

% 解方程组
vars = [x_Helicopterpp, y_Helicopterpp, z_Helicopterpp, phi_Helicopterpp, theta_Helicopterpp, psi_Helicopterpp, alpha_xpp, alpha_ypp, beta_xpp, beta_ypp];
sol = solve(equ, vars);

% 显示解
x_Helicopterpp = sol.x_Helicopterpp;
y_Helicopterpp = sol.y_Helicopterpp;
z_Helicopterpp = sol.z_Helicopterpp;
phi_Helicopterpp = sol.phi_Helicopterpp;
theta_Helicopterpp = sol.theta_Helicopterpp;
psi_Helicopterpp = sol.psi_Helicopterpp;
alpha_xpp = sol.alpha_xpp;
alpha_ypp = sol.alpha_ypp;
beta_xpp = sol.beta_xpp;
beta_ypp = sol.beta_ypp;