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

zxcvbvbnmn

solve函数是解方程经常使用的函数，有很多坛友在使用该函数时经常遇到问题，现简单汇总如下 一、使用用法老旧 3 K+ i/ `5 U4 r4 f( t  n  n/ O. P+ N例 syms x a=1; equ='x^2==a'; xs=solve(equ,x) / J, U. N: P) u4 }' k( M ) {( |/ `8 O3 L, V( l, N; R[color=rgb(51, 102, 153) !important]复制代码 ; e2 f0 v" i! _4 t" f; Y6 a/ [! A% } 会报错 6 w+ |- s8 t: Y9 ^# p8 _8 N4 PError using solve>getEqns (line 418) List of equations must not be empty. 错误在于使用用法过于陈旧，现使用的MATLAB多为高版本 解决办法 使用用法改变，具体参见solve帮助文档 例 5 F5 N; Q% N, a4 _' r: _# g6 ^ syms x a=1; equ=x^2==a; xs=solve(equ,x) ) |; K  p: L6 h" ~+ R0 y; i4 }( j  J* y[color=rgb(51, 102, 153) !important]复制代码 xs = -1 9 s! U' U4 P  f/ C& _  1 二、循环求解 1 \3 I) J1 x1 W! X解方程经常遇到方程的一个或多个参数是变化的 5 Q% p, Z8 d, q$ T1 l例如例子的变量a分别等于1、4、9、16 ( N1 m  B1 }8 b经常会出现以下代码 / o, N0 b1 W9 R! }8 N syms x a=[1,4,9,16]; equ=x^2==a; xs=solve(equ,x) [color=rgb(51, 102, 153) !important]复制代码 2 u9 A' @' _. r( I* m0 A  }MATLAB运算无解 + A5 I7 T; _- wxs = - J/ c9 L! b2 M: g* e: b; F# S; B Empty sym: 0-by-1 这种情况应使用循环求解 4 d( K0 q9 }9 X但有时会遇到以下情况 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]复制代码 $ n  H- m4 i* v/ K, d, ^$ ~$ c报错如下 / Y4 e$ r! y( v7 d) D8 Dx = , I& s/ c, s9 s4 t-1 + m4 A7 s$ F' s+ A3 u6 C# D/ L  1 / L" d0 e; k( y% p) b( J" j; l: ^Error using sym.getEqnsVars>checkVariables (line 92) Second argument must be a vector of symbolic variables. 原因在于经过第一次循环后，变量x以转变为具体数值（-1，1），即导致在第二次循环中方程里无未知量 * Z- i3 F5 C# H% N  Q& A+ t/ T- ~1 b解决办法保存解得变量名改变 + k, b% A6 W& M9 y9 h: W: c: d8 T5 X循环求解方程并保存变量的代码如下 5 ^' n- l6 Z1 [% {& ?& v! d syms x a=[1,4,9,16]; for i=1:length(a) equ=x.^2==a(i); xs(i,=solve(equ,x)； end [color=rgb(51, 102, 153) !important]复制代码 运行后在命令窗口输入xs即可看到全部解如下 1 A( j3 X5 j; l7 R& U) x>> xs 5 D* j  k1 D  c6 q ! V  E" ^( n( k8 M$ b6 zxs = [ -1, 1] & D( |, z* Y. w- A% K[ -2, 2] [ -3, 3] [ -4, 4] # }0 W4 D6 h9 s% Z, n9 @% B+ x$ _希望对需要的坛友有所帮助，也希望大家多提宝贵意见，共同进步

fantasyqqq

楼主很细心