public static void Test()
{
{
/*
* 部分分数分解のテスト。
* s + 3
* X(s) = -------------
* (s+1)s(s-2)
*
* 部分分数分解すると
*
* 2/3 -3/2 5/6
* X(s) = ------- + ------ + -------
* s + 1 s s - 2
*/
var numerCoeffs = new WWComplex [] {
new WWComplex(3, 0),
new WWComplex(1, 0)
};
var dRoots = new WWComplex [] {
new WWComplex(-1, 0),
WWComplex.Zero(),
new WWComplex(2, 0)
};
var p = WWPolynomial.PartialFractionDecomposition(numerCoeffs, dRoots);
for (int i = 0; i < p.Count; ++i)
{
Console.WriteLine(p[i].ToString("s"));
if (i != p.Count - 1)
{
Console.WriteLine(" + ");
}
}
Console.WriteLine("");
}
{
// 約分のテスト
// p-1 (p+1) -(p+1) + (p-1) -2
// ───── ⇒ ────────────────────── ⇒ 1 + ─────
// p+1 p+1 p+1
var numerC = new List <WWComplex>();
numerC.Add(new WWComplex(-1, 0)); // 定数項。
numerC.Add(WWComplex.Unity()); // 1乗の項。
var denomR = new List <WWComplex>();
denomR.Add(new WWComplex(-1, 0));
var r = Reduction(numerC.ToArray(), denomR.ToArray());
r.Print("p");
}
{
// 約分のテスト
// p^2+3x+2 (p^2+3x+2) -(p^2+7x+12) -4x-10
// ──────────── ⇒ 1+ ───────────────────────── ⇒ 1 + ────────────
// (p+3)(p+4) p^2+7x+12 (p+3)(p+4)
var numerC = new WWComplex [] {
new WWComplex(2, 0), // 定数項。
new WWComplex(3, 0), // 1乗の項。
WWComplex.Unity()
}; // 2乗の項。
var denomR = new WWComplex [] {
new WWComplex(-3, 0),
new WWComplex(-4, 0)
};
var r = Reduction(numerC, denomR);
r.Print("p");
}
{
// 部分分数分解。
// -4x-10 2 -6
// ──────────── ⇒ ───── + ─────
// (p+3)(p+4) p+3 p+4
var numerC = new WWComplex [] {
new WWComplex(-10, 0),
new WWComplex(-4, 0)
};
var denomR = new WWComplex [] {
new WWComplex(-3, 0),
new WWComplex(-4, 0)
};
var p = WWPolynomial.PartialFractionDecomposition(numerC, denomR);
for (int i = 0; i < p.Count; ++i)
{
Console.WriteLine(p[i].ToString("s"));
if (i != p.Count - 1)
{
Console.WriteLine(" + ");
}
}
Console.WriteLine("");
}
{
var deriv = new RealPolynomial(new double[] { 1, 1, 1, 1 }).Derivative();
Console.WriteLine("derivative of p^3+p^2+p+1={0}",
deriv.ToString("p"));
}
{
var r2 = AlgebraicLongDivision(new RealPolynomial(new double[] { 6, 3, 1 }),
new RealPolynomial(new double[] { 1, 1 }));
Console.WriteLine("(p^2+3x+6)/(p+1) = {0} r {1}",
r2.quotient.ToString(), r2.remainder.ToString());
}
}
