2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




Начать новую тему Ответить на тему
 
 Большая СЛАУ в символьном виде. Mathematica.
Сообщение06.05.2010, 16:27 


01/12/06
463
МИНСК
Здравствуйте. Может быть, кто-нибудь подскажет, как лучше решить СЛАУ с такой матрицей в Mathematica:
$$
\begin{array}{lllllllllllllll}
 -i \text{c66}(1) \text{w1} \text{k3}(1,1) \sqrt{\text{w1}^2+\text{w2}^2} & i \text{c66}(1) \text{w1} \text{k3}(1,1) \sqrt{\text{w1}^2+\text{w2}^2} & -i
   \text{c66}(1) \text{w1} \text{k3}(2,1) \sqrt{\text{w1}^2+\text{w2}^2} & i \text{c66}(1) \text{w1} \text{k3}(2,1) \sqrt{\text{w1}^2+\text{w2}^2} & -\frac{i
   \text{c66}(1) \text{w2} \sqrt{\text{w1}^2+\text{w2}^2}}{s(3,1)} & \frac{i \text{c66}(1) \text{w2} \sqrt{\text{w1}^2+\text{w2}^2}}{s(3,1)} & 0 & 0 & 0 & 0 & 0 & 0
   & 0 & 0 & 0 \\
 -i \text{c66}(1) \text{w2} \text{k3}(1,1) \sqrt{\text{w1}^2+\text{w2}^2} & i \text{c66}(1) \text{w2} \text{k3}(1,1) \sqrt{\text{w1}^2+\text{w2}^2} & -i
   \text{c66}(1) \text{w2} \text{k3}(2,1) \sqrt{\text{w1}^2+\text{w2}^2} & i \text{c66}(1) \text{w2} \text{k3}(2,1) \sqrt{\text{w1}^2+\text{w2}^2} & \frac{i
   \text{c66}(1) \text{w1} \sqrt{\text{w1}^2+\text{w2}^2}}{s(3,1)} & -\frac{i \text{c66}(1) \text{w1} \sqrt{\text{w1}^2+\text{w2}^2}}{s(3,1)} & 0 & 0 & 0 & 0 & 0 &
   0 & 0 & 0 & 0 \\
 \text{c66}(1) \text{k2}(1,1) \left(\text{w1}^2+\text{w2}^2\right) & \text{c66}(1) \text{k2}(1,1) \left(\text{w1}^2+\text{w2}^2\right) & \text{c66}(1)
   \text{k2}(2,1) \left(\text{w1}^2+\text{w2}^2\right) & \text{c66}(1) \text{k2}(2,1) \left(\text{w1}^2+\text{w2}^2\right) & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 & 0 &
   0 \\
 -i \text{c66}(1) \text{w1} \text{k3}(1,1) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{c66}(1) \text{w1}
   \text{k3}(1,1) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(1) \text{w1} \text{k3}(2,1)
   \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{c66}(1) \text{w1} \text{k3}(2,1) \sqrt{\text{w1}^2+\text{w2}^2}
   e^{-\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -\frac{i \text{c66}(1) \text{w2} \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(3,1)
   \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,1)} & \frac{i \text{c66}(1) \text{w2} \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(3,1)
   \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,1)} & i \text{c66}(2) \text{w1} \text{k3}(1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(1,2)
   \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(2) \text{w1} \text{k3}(1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} &
   i \text{c66}(2) \text{w1} \text{k3}(2,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(2) \text{w1}
   \text{k3}(2,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & \frac{i \text{c66}(2) \text{w2}
   \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,2)} & -\frac{i \text{c66}(2) \text{w2} \sqrt{\text{w1}^2+\text{w2}^2}
   e^{-\text{h1} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,2)} & 0 & 0 & 0 \\
 -i \text{c66}(1) \text{w2} \text{k3}(1,1) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{c66}(1) \text{w2}
   \text{k3}(1,1) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(1) \text{w2} \text{k3}(2,1)
   \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{c66}(1) \text{w2} \text{k3}(2,1) \sqrt{\text{w1}^2+\text{w2}^2}
   e^{-\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & \frac{i \text{c66}(1) \text{w1} \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(3,1)
   \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,1)} & -\frac{i \text{c66}(1) \text{w1} \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(3,1)
   \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,1)} & i \text{c66}(2) \text{w2} \text{k3}(1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(1,2)
   \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(2) \text{w2} \text{k3}(1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} &
   i \text{c66}(2) \text{w2} \text{k3}(2,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(2) \text{w2}
   \text{k3}(2,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -\frac{i \text{c66}(2) \text{w1}
   \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,2)} & \frac{i \text{c66}(2) \text{w1} \sqrt{\text{w1}^2+\text{w2}^2}
   e^{-\text{h1} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,2)} & 0 & 0 & 0 \\
 \text{c66}(1) \text{k2}(1,1) \left(\text{w1}^2+\text{w2}^2\right) e^{\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & \text{c66}(1) \text{k2}(1,1)
   \left(\text{w1}^2+\text{w2}^2\right) e^{-\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & \text{c66}(1) \text{k2}(2,1) \left(\text{w1}^2+\text{w2}^2\right)
   e^{\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & \text{c66}(1) \text{k2}(2,1) \left(\text{w1}^2+\text{w2}^2\right) e^{-\text{h1} s(2,1)
   \sqrt{\text{w1}^2+\text{w2}^2}} & 0 & 0 & \text{c66}(2) \text{k2}(1,2) \left(\text{w1}^2+\text{w2}^2\right) \left(-e^{\text{h1} s(1,2)
   \sqrt{\text{w1}^2+\text{w2}^2}}\right) & \text{c66}(2) \text{k2}(1,2) \left(\text{w1}^2+\text{w2}^2\right) \left(-e^{-\text{h1} s(1,2)
   \sqrt{\text{w1}^2+\text{w2}^2}}\right) & \text{c66}(2) \text{k2}(2,2) \left(\text{w1}^2+\text{w2}^2\right) \left(-e^{\text{h1} s(2,2)
   \sqrt{\text{w1}^2+\text{w2}^2}}\right) & \text{c66}(2) \text{k2}(2,2) \left(\text{w1}^2+\text{w2}^2\right) \left(-e^{-\text{h1} s(2,2)
   \sqrt{\text{w1}^2+\text{w2}^2}}\right) & 0 & 0 & 0 & 0 & 0 \\
 -i \text{w1} e^{\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w1} e^{-\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w1} e^{\text{h1}
   s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w1} e^{-\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w2} e^{\text{h1} s(3,1)
   \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w2} e^{-\text{h1} s(3,1) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{\text{h1} s(1,2)
   \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{-\text{h1} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{\text{h1} s(2,2)
   \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{-\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w2} e^{\text{h1} s(3,2)
   \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w2} e^{-\text{h1} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}} & 0 & 0 & 0 \\
 -i \text{w2} e^{\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w2} e^{-\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w2} e^{\text{h1}
   s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w2} e^{-\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{\text{h1} s(3,1)
   \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{-\text{h1} s(3,1) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w2} e^{\text{h1} s(1,2)
   \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w2} e^{-\text{h1} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w2} e^{\text{h1} s(2,2)
   \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w2} e^{-\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w1} e^{\text{h1} s(3,2)
   \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w1} e^{-\text{h1} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}} & 0 & 0 & 0 \\
 \alpha (1,1) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}} & \alpha (1,1) \sqrt{\text{w1}^2+\text{w2}^2} \left(-e^{-\text{h1}
   s(1,1) \sqrt{\text{w1}^2+\text{w2}^2}}\right) & \alpha (2,1) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}} & \alpha (2,1)
   \sqrt{\text{w1}^2+\text{w2}^2} \left(-e^{-\text{h1} s(2,1) \sqrt{\text{w1}^2+\text{w2}^2}}\right) & 0 & 0 & \alpha (1,2) \sqrt{\text{w1}^2+\text{w2}^2}
   \left(-e^{\text{h1} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}}\right) & \alpha (1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(1,2)
   \sqrt{\text{w1}^2+\text{w2}^2}} & \alpha (2,2) \sqrt{\text{w1}^2+\text{w2}^2} \left(-e^{\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}}\right) & \alpha (2,2)
   \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h1} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & 0 & 0 & 0 & 0 & 0 \\
 0 & 0 & 0 & 0 & 0 & 0 & -i \text{c66}(2) \text{w1} \text{k3}(1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i
   \text{c66}(2) \text{w1} \text{k3}(1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(2) \text{w1}
   \text{k3}(2,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{c66}(2) \text{w1} \text{k3}(2,2)
   \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -\frac{i \text{c66}(2) \text{w2} \sqrt{\text{w1}^2+\text{w2}^2}
   e^{\text{h2} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,2)} & \frac{i \text{c66}(2) \text{w2} \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(3,2)
   \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,2)} & -i \text{c66}(3) \text{w1} \text{k3}(1,3) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(1,3)
   \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(3) \text{w1} \text{k3}(2,3) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(2,3) \sqrt{\text{w1}^2+\text{w2}^2}} &
   -\frac{i \text{c66}(3) \text{w2} \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(3,3) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,3)} \\
 0 & 0 & 0 & 0 & 0 & 0 & -i \text{c66}(2) \text{w2} \text{k3}(1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i
   \text{c66}(2) \text{w2} \text{k3}(1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(2) \text{w2}
   \text{k3}(2,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{c66}(2) \text{w2} \text{k3}(2,2)
   \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & \frac{i \text{c66}(2) \text{w1} \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h2}
   s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,2)} & -\frac{i \text{c66}(2) \text{w1} \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(3,2)
   \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,2)} & -i \text{c66}(3) \text{w2} \text{k3}(1,3) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(1,3)
   \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{c66}(3) \text{w2} \text{k3}(2,3) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(2,3) \sqrt{\text{w1}^2+\text{w2}^2}} &
   \frac{i \text{c66}(3) \text{w1} \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(3,3) \sqrt{\text{w1}^2+\text{w2}^2}}}{s(3,3)} \\
 0 & 0 & 0 & 0 & 0 & 0 & \text{c66}(2) \text{k2}(1,2) \left(\text{w1}^2+\text{w2}^2\right) e^{\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & \text{c66}(2)
   \text{k2}(1,2) \left(\text{w1}^2+\text{w2}^2\right) e^{-\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & \text{c66}(2) \text{k2}(2,2)
   \left(\text{w1}^2+\text{w2}^2\right) e^{\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & \text{c66}(2) \text{k2}(2,2) \left(\text{w1}^2+\text{w2}^2\right)
   e^{-\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & 0 & 0 & \text{c66}(3) \text{k2}(1,3) \left(\text{w1}^2+\text{w2}^2\right) \left(-e^{-\text{h2} s(1,3)
   \sqrt{\text{w1}^2+\text{w2}^2}}\right) & \text{c66}(3) \text{k2}(2,3) \left(\text{w1}^2+\text{w2}^2\right) \left(-e^{-\text{h2} s(2,3)
   \sqrt{\text{w1}^2+\text{w2}^2}}\right) & 0 \\
 0 & 0 & 0 & 0 & 0 & 0 & -i \text{w1} e^{\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w1} e^{-\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i
   \text{w1} e^{\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w1} e^{-\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w2} e^{\text{h2}
   s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w2} e^{-\text{h2} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{-\text{h2} s(1,3)
   \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{-\text{h2} s(2,3) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w2} e^{-\text{h2} s(3,3)
   \sqrt{\text{w1}^2+\text{w2}^2}} \\
 0 & 0 & 0 & 0 & 0 & 0 & -i \text{w2} e^{\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w2} e^{-\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i
   \text{w2} e^{\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w2} e^{-\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{\text{h2}
   s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w1} e^{-\text{h2} s(3,2) \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w2} e^{-\text{h2} s(1,3)
   \sqrt{\text{w1}^2+\text{w2}^2}} & i \text{w2} e^{-\text{h2} s(2,3) \sqrt{\text{w1}^2+\text{w2}^2}} & -i \text{w1} e^{-\text{h2} s(3,3)
   \sqrt{\text{w1}^2+\text{w2}^2}} \\
 0 & 0 & 0 & 0 & 0 & 0 & \alpha (1,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}} & \alpha (1,2)
   \sqrt{\text{w1}^2+\text{w2}^2} \left(-e^{-\text{h2} s(1,2) \sqrt{\text{w1}^2+\text{w2}^2}}\right) & \alpha (2,2) \sqrt{\text{w1}^2+\text{w2}^2} e^{\text{h2}
   s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}} & \alpha (2,2) \sqrt{\text{w1}^2+\text{w2}^2} \left(-e^{-\text{h2} s(2,2) \sqrt{\text{w1}^2+\text{w2}^2}}\right) & 0 & 0 &
   \alpha (1,3) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2} s(1,3) \sqrt{\text{w1}^2+\text{w2}^2}} & \alpha (2,3) \sqrt{\text{w1}^2+\text{w2}^2} e^{-\text{h2}
   s(2,3) \sqrt{\text{w1}^2+\text{w2}^2}} & 0
\end{array}$$

 Профиль  
                  
 
 Re: Большая СЛАУ в символьном виде. Mathematica.
Сообщение10.05.2010, 02:49 
Заслуженный участник


20/04/10
1999
А в чём проблема?? Если вам в символьном виде решить нужно, то здесь действительно математика задумается, возможно и навсегда)). Но с численным решением никаких проблем не будет, решали с матрицами в разы побольше. Решение можно оформить как функцию ваших параметров и тогда символьное решение ни к чему совсем. Поэтому, если уж здесь Вам было не лень ентого крокодила набирать, то теперь тоже в математике, ну и функция LinearSolve[m,b]. Посмотрите в хэлпе, у неё есть опция Method, она может быть Вам полезна.

 Профиль  
                  
 
 Re: Большая СЛАУ в символьном виде. Mathematica.
Сообщение10.05.2010, 13:26 


01/12/06
463
МИНСК
Решение нужно в символьном виде. После некоторых доработок размерность системы стала меньше. И ее в явном виде удается решить методом Краммера, а LinearSolve, почему-то, надолго задумывается и выдает сообщение, что не хватает памяти. Я думал, что, может быть, есть что-то эффективнее.

P.S. В Mathematica есть функция TeXForm, поэтому здесь систему мне набирать не пришлось.

 Профиль  
                  
 
 Re: Большая СЛАУ в символьном виде. Mathematica.
Сообщение10.05.2010, 13:36 
Заслуженный участник


20/04/10
1999
Андрей123 в сообщении #317582 писал(а):
Решение нужно в символьном виде.
Честно признаться, сложно представить для чего. Вы бы могли написать функцию которая бы возвращала столбец решений, а её аргументами были бы параметры матрицы. Но воля ваша, главное обязательно проверьте совпадение численного решения при определенных параметрах с решением аналитическим при тех же параметрах, так, на всякий случай.
Андрей123 в сообщении #317582 писал(а):
P.S. В Mathematica есть функция TeXForm, поэтому здесь систему мне набирать не пришлось.
Понятно)

 Профиль  
                  
 
 Re: Большая СЛАУ в символьном виде. Mathematica.
Сообщение11.05.2010, 16:45 


01/12/06
463
МИНСК
Комбинацию решений с некоторыми функциями надо потом обращать по Фурье. С численным решением, кстати, тоже не все хорошо. При некоторых заначениях параметров пишет, что матрица плохо обсловлена.

 Профиль  
                  
Показать сообщения за:  Поле сортировки  
Начать новую тему Ответить на тему  [ Сообщений: 5 ] 

Модераторы: Karan, Toucan, PAV, maxal, Супермодераторы



Кто сейчас на конференции

Сейчас этот форум просматривают: choocha


Вы не можете начинать темы
Вы не можете отвечать на сообщения
Вы не можете редактировать свои сообщения
Вы не можете удалять свои сообщения
Вы не можете добавлять вложения

Найти:
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group