01 /*
02 * $Id: ComplexRationalPolynomialMatrix2.java,v 1.5 2008/02/02 05:53:02 koga Exp $
03 *
04 * Copyright (C) 2004 Koga Laboratory. All rights reserved.
05 *
06 */
07 package matxbook.chap17;
08
09 import org.mklab.nfc.matrix.RationalPolynomialMatrix;
10 import org.mklab.nfc.scalar.DoubleComplexNumber;
11 import org.mklab.nfc.scalar.Polynomial;
12 import org.mklab.nfc.scalar.RationalPolynomial;
13 import org.mklab.tool.matrix.Simplify;
14
15
16 /**
17 * @author koga
18 * @version $Revision: 1.5 $, 2004/04/19
19 */
20 public class ComplexRationalPolynomialMatrix2 {
21
22 /**
23 * メインメソッド
24 *
25 * @param args コマンドライン引数
26 */
27 @SuppressWarnings("nls")
28 public static void main(String[] args) {
29 Polynomial s = new Polynomial("s");
30 RationalPolynomial a11 = s.add(1).divide(s.add(2));
31 RationalPolynomial a12 = s.add(3).divide(s.add(4));
32 RationalPolynomial a21 = s.add(5).divide(s.add(6));
33 RationalPolynomial a22 = s.add(7).divide(s.add(8));
34 RationalPolynomialMatrix re = new RationalPolynomialMatrix(new RationalPolynomial[][] { {a11, a12}, {a21, a22}});
35 RationalPolynomial b11 = s.add(3).divide(s.add(2));
36 RationalPolynomial b12 = s.add(5).divide(s.add(4));
37 RationalPolynomial b21 = s.add(7).divide(s.add(6));
38 RationalPolynomial b22 = s.add(1).divide(s.add(8));
39 RationalPolynomialMatrix im = new RationalPolynomialMatrix(new RationalPolynomial[][] { {b11, b12}, {b21, b22}});
40 DoubleComplexNumber i = new DoubleComplexNumber(0, 1);
41 RationalPolynomialMatrix a = (RationalPolynomialMatrix)re.add(im.multiply(i));
42 Simplify.simplify(a).roundToZeroElementWise().print("a");
43 }
44 }
|