07-38.txt

这是《MATLAB 7简明教程》一书配套的实例源码。

例7-38 使用solve函数求解一般代数方程组。
解：在命令窗口中输入如下命令，并按Enter键确认。
>>  solve('p*sin(x) = r')
ans = 
asin(r/p)
>>
这里x是未知量，继续在命令窗口中输入以下程序，并按Enter键确认。
>> [x,y] = solve('x^2 + x*y + y = 3','x^2 - 4*x + 3 = 0')
x = 
[ 1]
[ 3]
y =
[    1]
[ -3/2]
>> 
继续在命令窗口中输入以下程序，并按Enter键确认。
>>S = solve('x^2*y^2 - 2*x - 1 = 0','x^2 - y^2 - 1 = 0')
S = 
    x: [8x1 sym]
    y: [8x1 sym]
>>
这里S是一个结构体。继续在命令窗口中输入以下程序，并按Enter键确认。
>> [u,v] = solve('a*u^2 + v^2 = 0','u - v = 1')
u = 
[ 1/2/(a+1)*(-2*a+2*(-a)^(1/2))+1]
[ 1/2/(a+1)*(-2*a-2*(-a)^(1/2))+1]
v =
[ 1/2/(a+1)*(-2*a+2*(-a)^(1/2))]
[ 1/2/(a+1)*(-2*a-2*(-a)^(1/2))]
>> 
这里a被作为参数求解关于u和v的方程组。继续在命令窗口中输入以下程序，并按Enter键确认。
>>S = solve('a*u^2 + v^2','u - v = 1','a,u')
S = 
    a: [1x1 sym]
    u: [1x1 sym]
>>
这里v被作为参数来求解方程组，并且返回S.a 和S.u。继续在命令窗口中输入以下程序，并按Enter键确认。
>> [a,u,v] = solve('a*u^2 + v^2','u - v = 1','a^2 - 5*a + 6')
a =
[ 2]
[ 2]
[ 3]
[ 3]
u =
[ 1/3+1/3*i*2^(1/2)]
[ 1/3-1/3*i*2^(1/2)]
[ 1/4+1/4*i*3^(1/2)]
[ 1/4-1/4*i*3^(1/2)]
v =
 [ -2/3+1/3*i*2^(1/2)]
[ -2/3-1/3*i*2^(1/2)]
[ -3/4+1/4*i*3^(1/2)]
[ -3/4-1/4*i*3^(1/2)]
>>
上边的程序段求解了关于a、u和v 的方程组。继续在命令窗口中输入以下程序，并按Enter键确认。
>> [x,y] = solve('sin(x+y)-exp(x)*y = 0','x^2-y = 2')
x =
-6.0173272500593065641097297117905
y =
34.208227234306296508646214438330
>>