matlab 如何解非线性超静定符号方程组

wxguang1

一下方程的K u n为未知的变量 其余为以知变量
syms m l n K u k_1 k_4
%求解使得K_ZDJX  -  K_joint_surface_diag=0;


K_ZDJX =

                      1/18*((2*m^2*l^2+2*n^2*m^2+2*n^2*l^2)*u-m^2*l^2-n^2*m^2-2*n^2*l^2)*E/n/m/l/(1+u)/(-1+2*u)
                      1/18*((2*m^2*l^2+2*n^2*m^2+2*n^2*l^2)*u-m^2*l^2-2*n^2*m^2-n^2*l^2)*E/n/m/l/(1+u)/(-1+2*u)
                      1/18*((2*m^2*l^2+2*n^2*m^2+2*n^2*l^2)*u-2*m^2*l^2-n^2*m^2-n^2*l^2)/n*E/(1+u)/l/(-1+2*u)/m
                           1/6*((2*m^2*l^2+2*n^2*m^2+n^2*l^2)*u-n^2*l^2-m^2*l^2-n^2*m^2)/l/n*E/m/(1+u)/(-1+2*u)
                    1/12*((4*m^2*l^2+4*n^2*m^2+2*n^2*l^2)*u-2*m^2*l^2-4*n^2*m^2-n^2*l^2)/l/n*E/m/(1+u)/(-1+2*u)
                    1/12*((4*m^2*l^2+4*n^2*m^2+2*n^2*l^2)*u-4*m^2*l^2-2*n^2*m^2-n^2*l^2)*E/n/m/l/(1+u)/(-1+2*u)
                      1/18*((2*m^2*l^2+2*n^2*m^2+2*n^2*l^2)*u-m^2*l^2-n^2*m^2-2*n^2*l^2)*E/n/m/l/(1+u)/(-1+2*u)
                      1/18*((2*m^2*l^2+2*n^2*m^2+2*n^2*l^2)*u-m^2*l^2-2*n^2*m^2-n^2*l^2)*E/n/m/l/(1+u)/(-1+2*u)
                      1/18*((2*m^2*l^2+2*n^2*m^2+2*n^2*l^2)*u-2*m^2*l^2-n^2*m^2-n^2*l^2)/n*E/(1+u)/l/(-1+2*u)/m
1/6*((2*n^2*m^2*l^2+6*m^2-6*m^2*n*l+2*n^2*l^4)*u-n^2*m^2*l^2-2*n^2*l^4+3*m^2*n*l-3*m^2)*n*E/l/(1+u)/(-1+2*u)/m
   1/6*((2*n^2*m^2*l^2+6*m^2-6*m^2*n*l+2*n^2*l^4)*u-2*n^2*m^2*l^2+6*m^2*n*l-6*m^2-n^2*l^4)*n*E/l/m/(-1+u+2*u^2)
                                                        1/6*n*E*(n^2*m^2*l^2-3*m^2*n*l+3*m^2+n^2*l^4)/l/m/(1+u)
                                                                                    1/6*m*E*(n^2+l^2)/n/l/(1+u)
                                                         1/6*((2*n^2+2*l^2)*u-2*n^2-l^2)*m*E/n/(1+u)/l/(-1+2*u)
                                                         1/6*((2*n^2+2*l^2)*u-n^2-2*l^2)*m*E/n/(1+u)/l/(-1+2*u)
                      1/18*((2*m^2*l^2+2*n^2*m^2+2*n^2*l^2)*u-m^2*l^2-n^2*m^2-2*n^2*l^2)*E/n/m/l/(1+u)/(-1+2*u)
                      1/18*((2*m^2*l^2+2*n^2*m^2+2*n^2*l^2)*u-m^2*l^2-2*n^2*m^2-n^2*l^2)*E/n/m/l/(1+u)/(-1+2*u)
                      1/18*((2*m^2*l^2+2*n^2*m^2+2*n^2*l^2)*u-2*m^2*l^2-n^2*m^2-n^2*l^2)/n*E/(1+u)/l/(-1+2*u)/m
                      1/18*((2*m^2*l^2+2*n^2*m^2+2*n^2*l^2)*u-m^2*l^2-n^2*m^2-2*n^2*l^2)*E/n/m/l/(1+u)/(-1+2*u)
                      1/18*((2*m^2*l^2+2*n^2*m^2+2*n^2*l^2)*u-m^2*l^2-2*n^2*m^2-n^2*l^2)*E/n/m/l/(1+u)/(-1+2*u)
                      1/18*((2*m^2*l^2+2*n^2*m^2+2*n^2*l^2)*u-2*m^2*l^2-n^2*m^2-n^2*l^2)/n*E/(1+u)/l/(-1+2*u)/m
                      1/18*((2*m^2*l^2+2*n^2*m^2+2*n^2*l^2)*u-m^2*l^2-n^2*m^2-2*n^2*l^2)*E/n/m/l/(1+u)/(-1+2*u)
                      1/18*((2*m^2*l^2+2*n^2*m^2+2*n^2*l^2)*u-m^2*l^2-2*n^2*m^2-n^2*l^2)*E/n/m/l/(1+u)/(-1+2*u)
                      1/18*((2*m^2*l^2+2*n^2*m^2+2*n^2*l^2)*u-2*m^2*l^2-n^2*m^2-n^2*l^2)/n*E/(1+u)/l/(-1+2*u)/m


K_joint_surface_diag =

1/9*m*l*k_1
1/9*m*l*k_1
1/9*m*l*k_4
1/9*m*l*k_1
1/9*m*l*k_1
1/9*m*l*k_4
1/9*m*l*k_1
1/9*m*l*k_1
1/9*m*l*k_4
1/9*m*l*k_1
1/9*m*l*k_1
1/9*m*l*k_4
1/9*m*l*k_1
1/9*m*l*k_1
1/9*m*l*k_4
1/9*m*l*k_1
1/9*m*l*k_1
1/9*m*l*k_4
1/9*m*l*k_1
1/9*m*l*k_1
1/9*m*l*k_4
1/9*m*l*k_1
1/9*m*l*k_1
1/9*m*l*k_4

谢谢各位了 急用非常感谢

messenger

非线性超静定符号方程组，solve就能解。

wxguang1

怎么解 举个例子可以吗 谢谢