new22解代数方程组.pdf

解代数方程组和复杂方程 解代数方程组 求解方程组 2 2 ⁄ = 0; − 2 = syms x y a [solx,soly] = solve(x^2*y^2 == 0, x-y/2 == a) solx = 0, a; soly = -2*a, 0 2 2 将第一个方程改变为 = 1 ，新方程组具有更多的解 [solx,soly] = solve(x^2*y^2 == 1, x-y/2 == a) solx = a/2 - (a^2 - 2)^(1/2)/2 soly = - a - (a^2 - 2)^(1/2) a/2 - (a^2 + 2)^(1/2)/2 - a - (a^2 + 2)^(1/2) a/2 + (a^2 - 2)^(1/2)/2 (a^2 - 2)^(1/2) - a a/2 + (a^2 + 2)^(1/2)/2 (a^2 + 2)^(1/2) - a 解方程组 2 2 2 2 − = ; + = 1; −2 ∗= 3 syms u v a 用 S = solve(u^2 - v^2 == a^2, u + v == 1, a^2 - 2*a == 3) 结 S = a: [2x1 sym] 构 u: [2x1 sym] v : [2x1 sym] 数 S.a %查看a 的内容 组 ans = -1, 3 输 M = [S.a, S.u, S.v] %查看a,u和v M = [ -1, 1, 0] 出 [ 3, 5, -4] 解 s2 = [S.a(2), S.u(2), S.v(2)] %查看第2个解 s2 = [ 3, 5, -4] 用左除解线性方程组 clear u v x y syms u v x y eqns = [x + 2*y == u, 4*x + 5*y == v]; S = solve(eqns); sol = [S.x; S.y] %输出solve 的解 sol = (2*v)/3 - (5*u)/3 (4*u)/3 - v/3 [A,b] = equationsToMatrix(eqns,x,y); % 求方程的矩阵形式 z = A\b z = (2*v)/3 - (5*u)/3 (4*u)/3 - v/3 复杂问题求简单解 1.求实数解 syms x , solve(x^5 - 1 == 0, x) ans = 1 - (2^(1/2)*(5 - 5^(1/2))^(1/2)*1i)/4 - 5^(1/2)/4 - 1/4 (2^(1/2)*(5 - 5^(1/2))^(1/2)*1i)/4 - 5^(1/2)/4 - 1/4 5^(1/2)/4 - (2^(1/2)*(5^(1/2) + 5)^(1/2)*1i)/4 - 1/4 5^(1/2)/4 + (2^(1/2)*(5^(1/2