Question

向您请教matlab中solve函数的解方程使用方法

如图所示，当使用0.6时可以算出结果，但是当使用一个变量k代替0.6时会报错，应该怎么用k代替呢？
我要在一个for循环中使用，k值会不断变化

Answer 1

Matlab中solve函数主要是用来求解线性方程组的解析解或者精确解。对于得出的结果是符号变量，可以通过vpa()得出任意位数的数值解！
solve函数的语法定义主要有以下四种：
solve(eq)
solve(eq, var)
solve(eq1, eq2, …, eqn)
g = solve(eq1, eq2, …, eqn, var1, var2, …, varn)
eq代表方程，var代表的是变量。
例1：
syms a b c x; solve(‘a*x^2 + b*x + c’)
当没有指定变量的时候matlab默认求解的是关于x的一元二次方程的解，求解的结果为：
ans = -(b + (b^2 – 4*a*c)^(1/2))/(2*a) -(b – (b^2 – 4*a*c)^(1/2))/(2*a)
当指定变量为b的时候：
syms a b c x; solve(‘a*x^2 + b*x + c’,'b’)
求解的结果为：
ans = -(a*x^2 + c)/x
从上面的例子很容易理解语法1,2。
例2：对于方程组的情况
syms x; S = solve(‘x + y = 1′,’x – 11*y = 5′); S = [S.x S.y]
求解的结果为：
S = [ 4/3, -1/3]
例3：
syms a u v; A = solve(‘a*u^2 + v^2′, ‘u – v = 1′, ‘a^2 – 5*a + 6′)
的求解结果为
A = a: [4x1 sym] u: [4x1 sym] v: [4x1 sym] 对于查看具体的数值可以通过
Aa = A.a Au = A.u Av = A.v
命令来查看。
PS：对于solve求解的方程，默认的为eq=0，eq1=0，eq2=0….eqn=0;

追问

答非所问！！！

Answer 2

clear all;close all;clc
syms x y
k=0.6;
assume(x>0)
y=psi(x)-k;
solve(y,x)
ans =
 2.2999264943186564750273431892452