MATLAB solve求方程组所有的解并assume添加条件

qq_45111955

已于 2022-12-20 11:23:53 修改

阅读量3.2k

点赞数

分类专栏： matlab语法文章标签： matlab 开发语言

于 2022-12-20 11:10:09 首次发布

本文链接：https://blog.csdn.net/qq_45111955/article/details/128382149

版权

matlab语法专栏收录该内容

21 篇文章 3 订阅

订阅专栏

直接solve得到的是一个特解，而不是所有解

clc
clear
syms c1 c2 k real
assume(k~=1/2 & k>0)
eq1=k*sin(c1-c2)-sin(c1)==0
eq2=k*sin(c2-c1)-sin(c2)==0
disp('直接solve得到特解：')
[solx1,solx2]=solve(eq1,eq2,[c1,c2])
disp('全部解：')
[solx1,solx2, param, cond]=solve(eq1,eq2,[c1,c2], 'ReturnConditions', true)

结果：

eq1 =

k*sin(c1 - c2) - sin(c1) == 0

eq2 =

- sin(c2) - k*sin(c1 - c2) == 0

直接solve得到特解：
警告: Solutions are only valid under certain conditions. To include parameters
and conditions in the solution, specify the 'ReturnConditions' value as 'true'.

solx1 =

0
pi
0
2*atan(((2*k - 1)*(2*k + 1))^(1/2)/(2*k + 1))
-2*atan(((2*k - 1)*(2*k + 1))^(1/2)/(2*k + 1))

solx2 =

0
0
pi
-2*atan(((2*k - 1)*(2*k + 1))^(1/2)/(2*k + 1))
2*atan(((2*k - 1)*(2*k + 1))^(1/2)/(2*k + 1))

全部解：

solx1 =

z
-pi*(2*m - 1)
pi*l
2*atan(((2*k - 1)*(2*k + 1))^(1/2)/(2*k + 1)) + 2*pi*l
2*pi*l - 2*atan(((2*k - 1)*(2*k + 1))^(1/2)/(2*k + 1))

solx2 =

z1
pi*l
-pi*(2*m - 1)
2*pi*m - 2*atan(((2*k - 1)*(2*k + 1))^(1/2)/(2*k + 1))
2*atan(((2*k - 1)*(2*k + 1))^(1/2)/(2*k + 1)) + 2*pi*m

param =

[l, m, z, z1]

cond =

in(z/(2*pi), 'integer') & in(z1/(2*pi), 'integer') & 0 < 4*k^2 - 1 | in(z/(2*pi), 'integer') & in(z1/(2*pi), 'integer') & 4*k^2 - 1 < 0
in(l, 'integer') & in(m, 'integer') & 0 < 4*k^2 - 1 | in(l, 'integer') & in(m, 'integer') & 4*k^2 - 1 < 0
in(l, 'integer') & in(m, 'integer') & 0 < 4*k^2 - 1 | in(l, 'integer') & in(m, 'integer') & 4*k^2 - 1 < 0
in(l, 'integer') & in(m, 'integer') & 0 < 4*k^2 - 1
in(l, 'integer') & in(m, 'integer') & 0 < 4*k^2 - 1

>>