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

zxcvbvbnmn

solve函数是解方程经常使用的函数，有很多坛友在使用该函数时经常遇到问题，现简单汇总如下 % M* X; l1 G, G5 {* M, D/ }, T一、使用用法老旧 例 syms x a=1; equ='x^2==a'; xs=solve(equ,x) ; Y2 {5 z" s. J0 L, k [color=rgb(51, 102, 153) !important]复制代码 ) d# B2 ?; w4 H8 k会报错 Error using solve>getEqns (line 418) 1 C9 s; ?. J4 }$ Z$ T- f0 l) JList of equations must not be empty. 错误在于使用用法过于陈旧，现使用的MATLAB多为高版本 解决办法 $ m: |1 Q9 \$ k2 J3 J使用用法改变，具体参见solve帮助文档 5 q9 A5 q: L/ P2 W& k例 syms x a=1; equ=x^2==a; xs=solve(equ,x)   L1 e( q2 R( A, r4 |+ ` [color=rgb(51, 102, 153) !important]复制代码 " S) y2 e- q: K/ L- o/ ]$ Dxs = 5 ~* a# j, Z. @% i-1   1 二、循环求解 解方程经常遇到方程的一个或多个参数是变化的 例如例子的变量a分别等于1、4、9、16 : j4 V+ v* I1 W8 L( I8 E( ^6 H经常会出现以下代码 syms x a=[1,4,9,16]; equ=x^2==a; xs=solve(equ,x) [color=rgb(51, 102, 153) !important]复制代码 MATLAB运算无解 xs = ( u: \+ I" }( `7 F8 P" s Empty sym: 0-by-1 这种情况应使用循环求解 - C! K$ h0 p0 M% N5 C+ Q. R- Z+ N$ V但有时会遇到以下情况 9 K* N/ k/ v4 k7 K4 z% m7 N syms x a=[1,4,9,16]; for i=1:length(a) equ=x.^2==a(i); x=solve(equ,x) end : ^) @+ Z" T- x* P/ M2 k# Y3 ^0 l[color=rgb(51, 102, 153) !important]复制代码 1 ?; _( l4 R+ z报错如下 x = 0 K4 c! C6 D* ]" o  Y! k5 W- ]-1 1 ^. t2 t, z0 O3 a  1 9 f" T! V- q) }& ~" T$ r$ L; k6 m Error using sym.getEqnsVars>checkVariables (line 92) # d0 N1 C: _" |0 ?- Y- vSecond argument must be a vector of symbolic variables. 原因在于经过第一次循环后，变量x以转变为具体数值（-1，1），即导致在第二次循环中方程里无未知量 解决办法保存解得变量名改变 / ^/ y8 ~9 n) k" s循环求解方程并保存变量的代码如下 : v% O3 @" H7 `4 v( B5 b8 |1 N syms x a=[1,4,9,16]; for i=1:length(a) equ=x.^2==a(i); xs(i,=solve(equ,x)； end 8 `* T$ @: L$ |, Z. t5 ` [color=rgb(51, 102, 153) !important]复制代码 运行后在命令窗口输入xs即可看到全部解如下 >> xs * x; U! j  D, {- y/ q" ]0 p2 O- s xs = " K: I! N( S* s' w  r [ -1, 1] [ -2, 2] [ -3, 3] $ s: P! O0 Z; u) i4 j[ -4, 4] 希望对需要的坛友有所帮助，也希望大家多提宝贵意见，共同进步

% L: y: G3 c& [0 L' J

fantasyqqq

楼主很细心