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

zxcvbvbnmn

solve函数是解方程经常使用的函数，有很多坛友在使用该函数时经常遇到问题，现简单汇总如下 / O+ F5 Y! C% M+ `8 E' ]+ d5 P一、使用用法老旧 ; [( e" g: ]& A% B5 M0 R例 syms x a=1; equ='x^2==a'; xs=solve(equ,x) $ J/ p2 \7 C# O [color=rgb(51, 102, 153) !important]复制代码 会报错 Error using solve>getEqns (line 418) % |- y( p$ E+ RList of equations must not be empty. 错误在于使用用法过于陈旧，现使用的MATLAB多为高版本 解决办法 使用用法改变，具体参见solve帮助文档 . G8 f0 j& H! C2 T. y, [例 ! T  a* c# Q0 w* s5 i: H) T+ [ syms x a=1; equ=x^2==a; xs=solve(equ,x) / n, U! e8 g# T( X% _ [color=rgb(51, 102, 153) !important]复制代码 ; f1 [+ {% Y) c" k4 Y7 g( |/ h xs = & ~' M( E, r+ @! h" o1 Q -1   1 二、循环求解 1 N3 e! T  V+ Q( e; ]0 z- w0 K8 m解方程经常遇到方程的一个或多个参数是变化的 例如例子的变量a分别等于1、4、9、16 ! o  k9 B' @4 C! o3 z9 c3 @经常会出现以下代码 syms x a=[1,4,9,16]; equ=x^2==a; xs=solve(equ,x) [color=rgb(51, 102, 153) !important]复制代码 0 S8 J0 s. H2 H. ^$ ]) O ' b) ?& Q5 V0 r2 f/ C' ?! ?( _" N+ AMATLAB运算无解 xs = Empty sym: 0-by-1 1 @- J: o( p/ |这种情况应使用循环求解 但有时会遇到以下情况 syms x a=[1,4,9,16]; for i=1:length(a) equ=x.^2==a(i); x=solve(equ,x) end ' _, q; w9 `* Y% B [color=rgb(51, 102, 153) !important]复制代码 . N( v; H6 J& h+ }3 b$ \! R7 ~ # x# V- e; ^* \( L& w7 i报错如下 9 T" I, t9 d# D) @x = 7 V1 H& [1 a; X4 y4 p) C1 @# J0 _ 3 a" F8 N7 v. N-1   1 Error using sym.getEqnsVars>checkVariables (line 92) 2 S( a+ J7 |; c% NSecond argument must be a vector of symbolic variables. 4 o$ H9 \3 G1 `, r" w" X5 L+ C原因在于经过第一次循环后，变量x以转变为具体数值（-1，1），即导致在第二次循环中方程里无未知量 解决办法保存解得变量名改变 1 M- h- Q% R% F循环求解方程并保存变量的代码如下 9 u5 e6 L1 y5 j8 O syms x a=[1,4,9,16]; for i=1:length(a) equ=x.^2==a(i); xs(i,=solve(equ,x)； end ' Y: B) y8 M! C/ r$ S [color=rgb(51, 102, 153) !important]复制代码 ) o+ s- ~5 l2 @6 t 8 F' I) x6 o0 J* N, r2 _/ \运行后在命令窗口输入xs即可看到全部解如下 >> xs ) H/ M9 _2 ]! ^& t  A+ `3 { xs = " P% c* |: m: l* M 8 d& P  _" G+ B1 W[ -1, 1] [ -2, 2]   p2 M6 `  C" g9 `1 }. t; x, o, `[ -3, 3] 4 b& W* H7 B: `! T" _! R) }0 W[ -4, 4] 希望对需要的坛友有所帮助，也希望大家多提宝贵意见，共同进步

, z# e: p4 k1 A+ C/ G3 ]

fantasyqqq

楼主很细心