solve求解方程的若干常见问题

zxcvbvbnmn

solve函数是解方程经常使用的函数，有很多坛友在使用该函数时经常遇到问题，现简单汇总如下 * G! F' j0 q# q( F$ M, r7 g/ \% D一、使用用法老旧 9 N2 V0 h5 o, P4 I例 syms x a=1; equ='x^2==a'; xs=solve(equ,x) ( f) A. A  P# W1 I 5 G. U2 _: k  R/ W6 |2 M0 h[color=rgb(51, 102, 153) !important]复制代码 # T7 r) X# f7 U1 ^* { 会报错 ( M& G; R; Y- J/ l. D9 X: EError using solve>getEqns (line 418) List of equations must not be empty. 错误在于使用用法过于陈旧，现使用的MATLAB多为高版本   G: \! m  H  s解决办法 使用用法改变，具体参见solve帮助文档 例 + q- |( F0 o! C% z syms x a=1; equ=x^2==a; xs=solve(equ,x) 2 j. T5 d* e0 ?) V1 j$ z  [[color=rgb(51, 102, 153) !important]复制代码 0 o$ H* G9 J+ a/ c. A xs = -1   1 - S1 @2 M$ P' f' J- }% K0 q" r$ ~二、循环求解 1 p/ D3 C! R* h; I0 F解方程经常遇到方程的一个或多个参数是变化的 + M1 e6 d9 u3 H7 E. e+ p: P例如例子的变量a分别等于1、4、9、16 9 @2 E- I" o8 n' \经常会出现以下代码 - Q! p& R# l9 D. L' _9 p% F syms x a=[1,4,9,16]; equ=x^2==a; xs=solve(equ,x) [color=rgb(51, 102, 153) !important]复制代码 0 ]' [5 [( [( ?( c/ R5 x/ U9 ] MATLAB运算无解 xs = ; _- J7 N1 K# o* d+ b; c5 j' }Empty sym: 0-by-1 & F5 b2 _( u  y0 ^6 f这种情况应使用循环求解 2 G, [. B0 L7 F. p" P0 m' H但有时会遇到以下情况 syms x a=[1,4,9,16]; for i=1:length(a) equ=x.^2==a(i); x=solve(equ,x) end [color=rgb(51, 102, 153) !important]复制代码 % y2 t/ R3 O1 ~" ^! ^3 ?, _报错如下 x = ) S% ?% O' }3 J  o" ^( S: B7 N-1   1 % n& c* I; Y0 B, t# T 2 i0 _! j! \# d4 `8 s+ NError using sym.getEqnsVars>checkVariables (line 92) - g2 N4 S, x+ b8 G3 QSecond argument must be a vector of symbolic variables. $ w' H& G9 x, l1 B. `; |- E2 i. ?3 H2 I0 I原因在于经过第一次循环后，变量x以转变为具体数值（-1，1），即导致在第二次循环中方程里无未知量 解决办法保存解得变量名改变 循环求解方程并保存变量的代码如下 syms x a=[1,4,9,16]; for i=1:length(a) equ=x.^2==a(i); xs(i,=solve(equ,x)； end ; k# W6 f+ G" z4 r3 M [color=rgb(51, 102, 153) !important]复制代码 运行后在命令窗口输入xs即可看到全部解如下 >> xs % J0 `& Y( u6 Q# b, P- L, ^" Xxs = 6 m5 N" C; q. R5 T2 G* z0 ] 8 @: D. r" h' O[ -1, 1] [ -2, 2] ' ~, F1 X. ?; B8 e# N! A9 X' f[ -3, 3] [ -4, 4] ' p5 U8 c/ E! x% b6 r' A7 _希望对需要的坛友有所帮助，也希望大家多提宝贵意见，共同进步

fantasyqqq

楼主很细心