matlab解方程组解析解

问题描述

对于一个如下形式的方程组
$$\{\begin{array}{l}{y}_1=f_1(x_1,x_2,...,x_n)\\ y_2=f_2(x_1,x_2,...,x_n)\\ ...\\ y_m=f_m(x_1,x_2,...,x_n)\end{array}$$
需要求解其反函数，即
$$\{\begin{array}{l}{x}_1=g_1(y_1,y_2,...,y_m)\\ x_2=g_2(y_1,y_2,...,y_m)\\ ...\\ x_n=g_m(y_1,y_2,...,y_m)\end{array}$$

matlab求解

可以通过matlab的slove函数可以实现上述功能

syms x1 x2 ....  xn, y1 y2 ... ym
eq1 = y1 == f_1(x_1, x_2, ..., x_n)
eq2 = y2 == f_2(x_1, x_2, ..., x_n)
...
eqm = ym == f_m(x_1, x_2, ..., x_n)
[x1_, x2_, ..., xn_] = solve(eq1, eq2, eq3, ..., eqm, x_1, x_2, ..., x_n)

具体实例

syms kappa d phi varepsilon_0 varepsilon_a varepsilon_b varepsilon_c
equ1 = varepsilon_a == kappa*d*cos(-phi) + varepsilon_0;
equ2 = varepsilon_b == kappa*d*cos(-phi+2/3*pi) + varepsilon_0;
equ3 = varepsilon_c == kappa*d*cos(-phi+4/3*pi) + varepsilon_0;

result = solve(equ1,equ2,equ3, kappa, phi, varepsilon_0)