2016-04-28 8 views
0

Когда я пытаюсь сделать это:Решить уравнение с комплексно сопряженное

(%i1) declare (z, complex); 
(%o1)        done 
(%i2) eq1: z^3 + 3 * %i * conjugate(z) = 0; 
               3 
(%o2)      3 %i conjugate(z) + z = 0 
(%i3) solve(eq1, z); 
        1/6 5/6  1/3    1/3 
      (- 1) (3 %i - 3 ) conjugate(z) 
(%o3) [z = - -----------------------------------------, 
           2 
     1/6 5/6  1/3    1/3 
    (- 1) (3 %i + 3 ) conjugate(z) 
z = -----------------------------------------, 
         2 
      1/6 1/3    1/3 
z = - (- 1) 3 conjugate(z) ] 

конъюгаты не упрощена. И решение для z в терминах z не очень полезно. Есть ли способ упростить его?

Также, как я могу упростить часть (-1)^(1/6)?

Кроме того, это уравнение явно имеет 0 в качестве его корня, но это не в наборе решений, почему?

+1

Для записи я вижу, что 'to_poly_solve' (который видит) может решить некоторые уравнения, содержащие' conjugate'. Для примера, приведенного здесь, он просто создает большой беспорядок, но, возможно, он более полезен для других проблем. –

+0

Спасибо @RobertDodier, что полезно знать. – wvxvw

ответ

2

не знаю solve ничего не знает о conjugate. Попробуйте это, чтобы решить эту проблему с реальными и мнимыми частями z как две переменные. Как это:

(%i2) declare ([zr, zi], real) $ 
(%i3) z : zr + %i*zi $ 
(%i4) eq1: z^3 + 3 * %i * conjugate(z) = 0; 
(%o4) (zr+%i*zi)^3+3*%i*(zr-%i*zi) = 0 
(%i5) solve (eq1, [zr, zi]); 
(%o5) [[zr = %r1, 
     zi = (sqrt(9*%r1^2-%i)+3*%r1)^(1/3)-%i/(sqrt(9*%r1^2-%i)+3*%r1)^(1/3) 
              +%i*%r1], 
     [zr = %r2, 
     zi = ((sqrt(3)*%i)/2-1/2)*(sqrt(9*%r2^2-%i)+3*%r2)^(1/3) 
      -(%i*((-(sqrt(3)*%i)/2)-1/2))/(sqrt(9*%r2^2-%i)+3*%r2)^(1/3) 
      +%i*%r2], 
     [zr = %r3, 
     zi = ((-(sqrt(3)*%i)/2)-1/2)*(sqrt(9*%r3^2-%i)+3*%r3)^(1/3) 
      -(%i*((sqrt(3)*%i)/2-1/2))/(sqrt(9*%r3^2-%i)+3*%r3)^(1/3)+%i*%r3]] 

Примечание переменные %r1, %r2 и %r3 в растворе. Они представляют собой произвольные значения.

+0

Спасибо за ответ! Но, честно говоря, я не очень доволен этим решением. Это дает мнимую часть в терминах реальной части, если я правильно понимаю это ... Возможно, есть способ сделать это в полярной форме? Я ожидаю, что это приведет к результату, близкому к '-3^(1/3) cis (% pi/2) + k *% pi'. – wvxvw

+0

Ну, вы можете поместить его обратно через 'ev (z: zr +% i * zi, zr =% r3, zi = );", а затем попробуйте упростить это с помощью 'radcan' и' polarform'. Но когда я пытаюсь, я получаю довольно большой беспорядок. –

+0

ОК, я вижу. Спасибо за усилия, хотя! – wvxvw