Java Code Examples for org.apache.commons.math3.complex.ComplexUtils#polar2Complex()

The following examples show how to use org.apache.commons.math3.complex.ComplexUtils#polar2Complex() . You can vote up the ones you like or vote down the ones you don't like, and go to the original project or source file by following the links above each example. You may check out the related API usage on the sidebar.

Example 1

Source File: Cascade.java From chart-fx with Apache License 2.0

6 votes

public Complex response(final double normalizedFrequency) {
    final double w = 2 * Math.PI * normalizedFrequency;
    final Complex czn1 = ComplexUtils.polar2Complex(1., -w);
    final Complex czn2 = ComplexUtils.polar2Complex(1., -2 * w);
    Complex ch = new Complex(1);
    Complex cbot = new Complex(1);

    for (int i = 0; i < mNumBiquads; i++) {
        final Biquad stage = mBiquads[i];

        Complex ct = new Complex(stage.getB0() / stage.getA0()); // NOPMD
        ct = ct.add(czn1.multiply(stage.getB1() / stage.getA0()));
        ct = ct.add(czn2.multiply(stage.getB2() / stage.getA0()));

        Complex cb = new Complex(1); // NOPMD
        cb = cb.add(czn1.multiply(stage.getA1() / stage.getA0()));
        cb = cb.add(czn2.multiply(stage.getA2() / stage.getA0()));

        ch = ch.multiply(ct);
        cbot = cbot.multiply(cb);
    }

    return ch.divide(cbot);
}

Example 2

Source File: LinkData.java From powsybl-core with Mozilla Public License 2.0

6 votes

public static BranchAdmittanceMatrix calculateBranchAdmittance(double r, double x, double ratio1, double alpha1,
    double ratio2, double alpha2, Complex ysh1, Complex ysh2) {

    Complex a1 = ComplexUtils.polar2Complex(ratio1, alpha1);
    Complex a2 = ComplexUtils.polar2Complex(ratio2, alpha2);

    Complex ytr;
    if (r == 0.0 && x == 0.0) {
        ytr = Complex.ZERO;
    } else {
        ytr = new Complex(r, x).reciprocal();
    }

    BranchAdmittanceMatrix branchAdmittance = new BranchAdmittanceMatrix();

    branchAdmittance.y11 = ytr.add(ysh1).divide(a1.conjugate().multiply(a1));
    branchAdmittance.y12 = ytr.negate().divide(a1.conjugate().multiply(a2));
    branchAdmittance.y21 = ytr.negate().divide(a2.conjugate().multiply(a1));
    branchAdmittance.y22 = ytr.add(ysh2).divide(a2.conjugate().multiply(a2));

    return branchAdmittance;
}

Example 3

Source File: SV.java From powsybl-core with Mozilla Public License 2.0

6 votes

public SV otherSide(double r, double x, double g, double b, double ratio) {
    Complex z = new Complex(r, x); // z=r+jx
    Complex y = new Complex(g, b); // y=g+jb
    Complex s1 = new Complex(p, q); // s1=p1+jq1
    Complex u1 = ComplexUtils.polar2Complex(u, Math.toRadians(a));
    Complex v1 = u1.divide(Math.sqrt(3f)); // v1=u1/sqrt(3)

    Complex v1p = v1.multiply(ratio); // v1p=v1*rho
    Complex i1 = s1.divide(v1.multiply(3)).conjugate(); // i1=conj(s1/(3*v1))
    Complex i1p = i1.divide(ratio); // i1p=i1/rho
    Complex i2 = i1p.subtract(y.multiply(v1p)); // i2=i1p-y*v1p
    Complex v2 = v1p.subtract(z.multiply(i2)); // v2=v1p-z*i2
    Complex s2 = v2.multiply(3).multiply(i2.conjugate()); // s2=3*v2*conj(i2)

    Complex u2 = v2.multiply(Math.sqrt(3f));
    return new SV(-s2.getReal(), -s2.getImaginary(), u2.abs(), Math.toDegrees(u2.getArgument()));
}

Example 4

Source File: SV.java From powsybl-core with Mozilla Public License 2.0

6 votes

public SV otherSide(double r, double x, double g1, double b1, double g2, double b2, double ratio) {
    Complex z = new Complex(r, x); // z=r+jx
    Complex y1 = new Complex(g1, b1); // y1=g1+jb1
    Complex y2 = new Complex(g2, b2); // y2=g2+jb2
    Complex s1 = new Complex(p, q); // s1=p1+jq1
    Complex u1 = ComplexUtils.polar2Complex(u, Math.toRadians(a));
    Complex v1 = u1.divide(Math.sqrt(3f)); // v1=u1/sqrt(3)

    Complex v1p = v1.multiply(ratio); // v1p=v1*rho
    Complex i1 = s1.divide(v1.multiply(3)).conjugate(); // i1=conj(s1/(3*v1))
    Complex i1p = i1.divide(ratio); // i1p=i1/rho
    Complex i2p = i1p.subtract(y1.multiply(v1p)); // i2p=i1p-y1*v1p
    Complex v2 = v1p.subtract(z.multiply(i2p)); // v2p=v1p-z*i2
    Complex i2 = i2p.subtract(y2.multiply(v2)); // i2=i2p-y2*v2
    Complex s2 = v2.multiply(3).multiply(i2.conjugate()); // s2=3*v2*conj(i2)

    Complex u2 = v2.multiply(Math.sqrt(3f));
    return new SV(-s2.getReal(), -s2.getImaginary(), u2.abs(), Math.toDegrees(u2.getArgument()));
}

Example 5

Source File: Cascade.java From iirj with Apache License 2.0

6 votes

public Complex response(double normalizedFrequency) {
	double w = 2 * Math.PI * normalizedFrequency;
	Complex czn1 = ComplexUtils.polar2Complex(1., -w);
	Complex czn2 = ComplexUtils.polar2Complex(1., -2 * w);
	Complex ch = new Complex(1);
	Complex cbot = new Complex(1);

	for (int i = 0; i < m_numBiquads; i++) {
		Biquad stage = m_biquads[i];
		Complex cb = new Complex(1);
		Complex ct = new Complex(stage.getB0() / stage.getA0());
		ct = MathSupplement.addmul(ct, stage.getB1() / stage.getA0(), czn1);
		ct = MathSupplement.addmul(ct, stage.getB2() / stage.getA0(), czn2);
		cb = MathSupplement.addmul(cb, stage.getA1() / stage.getA0(), czn1);
		cb = MathSupplement.addmul(cb, stage.getA2() / stage.getA0(), czn2);
		ch = ch.multiply(ct);
		cbot = cbot.multiply(cb);
	}

	return ch.divide(cbot);
}

Example 6

Source File: Biquad.java From iirj with Apache License 2.0

6 votes

public Complex response(double normalizedFrequency) {
    double a0 = getA0();
    double a1 = getA1();
    double a2 = getA2();
    double b0 = getB0();
    double b1 = getB1();
    double b2 = getB2();

    double w = 2 * Math.PI * normalizedFrequency;
    Complex czn1 = ComplexUtils.polar2Complex(1., -w);
    Complex czn2 = ComplexUtils.polar2Complex(1., -2 * w);
    Complex ch = new Complex(1);
    Complex cbot = new Complex(1);

    Complex ct = new Complex(b0 / a0);
    Complex cb = new Complex(1);
    ct = MathSupplement.addmul(ct, b1 / a0, czn1);
    ct = MathSupplement.addmul(ct, b2 / a0, czn2);
    cb = MathSupplement.addmul(cb, a1 / a0, czn1);
    cb = MathSupplement.addmul(cb, a2 / a0, czn2);
    ch = ch.multiply(ct);
    cbot = cbot.multiply(cb);

    return ch.divide(cbot);
}

Example 7

Source File: Biquad.java From chart-fx with Apache License 2.0

5 votes

public Complex response(final double normalizedFrequency) {
    final double a0 = getA0();
    final double a1 = getA1();
    final double a2 = getA2();
    final double b0 = getB0();
    final double b1 = getB1();
    final double b2 = getB2();

    final double w = 2 * Math.PI * normalizedFrequency;
    final Complex czn1 = ComplexUtils.polar2Complex(1., -w);
    final Complex czn2 = ComplexUtils.polar2Complex(1., -2 * w);
    Complex ch = new Complex(1);
    Complex cbot = new Complex(1);

    Complex ct = new Complex(b0 / a0);

    ct = ct.add(czn1.multiply(b1 / a0));
    ct = ct.add(czn2.multiply(b2 / a0));

    Complex cb = new Complex(1);
    cb = cb.add(czn1.multiply(a1 / a0));
    cb = cb.add(czn2.multiply(a2 / a0));

    ch = ch.multiply(ct);
    cbot = cbot.multiply(cb);

    return ch.divide(cbot);
}

Example 8

Source File: Butterworth.java From chart-fx with Apache License 2.0

5 votes

public void design() {
    reset();
    final double n2 = 2.0 * nPoles;
    final int pairs = nPoles / 2;
    for (int i = 0; i < pairs; ++i) {
        final Complex c = ComplexUtils.polar2Complex(1F, Math.PI / 2.0 + (2 * i + 1) * Math.PI / n2);
        addPoleZeroConjugatePairs(c, Complex.INF);
    }

    if ((nPoles & 1) == 1) {
        add(new Complex(-1), Complex.INF);
    }
}

Example 9

Source File: StateVariablesAdder.java From powsybl-core with Mozilla Public License 2.0

5 votes

private static Complex complexVoltage(double r, double x, double g, double b,
    double v, double angle, double p, double q) {
    BranchAdmittanceMatrix adm = LinkData.calculateBranchAdmittance(r, x, 1.0, 0.0, 1.0, 0.0,
        new Complex(g * 0.5, b * 0.5), new Complex(g * 0.5, b * 0.5));
    Complex v1 = ComplexUtils.polar2Complex(v, Math.toRadians(angle));
    Complex s1 = new Complex(p, q);
    return (s1.conjugate().divide(v1.conjugate()).subtract(adm.y11().multiply(v1))).divide(adm.y12());
}

Example 10

Source File: LinkData.java From powsybl-core with Mozilla Public License 2.0

5 votes

static Flow flowBothEnds(Complex y11, Complex y12, Complex y21, Complex y22,
    double u1, double theta1, double u2, double theta2) {

    Complex v1 = ComplexUtils.polar2Complex(u1, theta1);
    Complex v2 = ComplexUtils.polar2Complex(u2, theta2);

    return flowBothEnds(y11, y12, y21, y22, v1, v2);
}

Example 11

Source File: BranchDataTest.java From powsybl-core with Mozilla Public License 2.0

5 votes

Bus calcStarBusV1V2V3Y(BranchTestCase w1, BranchTestCase w2, BranchTestCase w3) {
    Complex v1 = ComplexUtils.polar2Complex(w1.bus1.u, w1.bus1.theta);
    Complex v2 = ComplexUtils.polar2Complex(w2.bus1.u, w2.bus1.theta);
    Complex v3 = ComplexUtils.polar2Complex(w3.bus1.u, w3.bus1.theta);
    Complex ytr1 = new Complex(w1.branch.end1.r, w1.branch.end1.x).reciprocal();
    Complex ytr2 = new Complex(w2.branch.end1.r, w2.branch.end1.x).reciprocal();
    Complex ytr3 = new Complex(w3.branch.end1.r, w3.branch.end1.x).reciprocal();

    // FIXME consider tap.rho and tap.alpha
    Complex a01 = new Complex(w1.branch.end2.ratedU / w1.branch.end1.ratedU, 0);
    Complex a1 = new Complex(1, 0);
    Complex a02 = new Complex(w2.branch.end2.ratedU / w2.branch.end1.ratedU, 0);
    Complex a2 = new Complex(1, 0);
    Complex a03 = new Complex(w3.branch.end2.ratedU / w3.branch.end1.ratedU, 0);
    Complex a3 = new Complex(1, 0);

    Complex ysh01 = new Complex(w1.branch.end2.g, w1.branch.end2.b);
    Complex ysh02 = new Complex(w2.branch.end2.g, w2.branch.end2.b);
    Complex ysh03 = new Complex(w3.branch.end2.g, w3.branch.end2.b);
    Complex y01 = ytr1.negate().divide(a01.conjugate().multiply(a1));
    Complex y02 = ytr2.negate().divide(a02.conjugate().multiply(a2));
    Complex y03 = ytr3.negate().divide(a03.conjugate().multiply(a3));
    Complex y0101 = ytr1.add(ysh01).divide(a01.conjugate().multiply(a01));
    Complex y0202 = ytr2.add(ysh02).divide(a02.conjugate().multiply(a02));
    Complex y0303 = ytr3.add(ysh03).divide(a03.conjugate().multiply(a03));

    Complex v0 = y01.multiply(v1).add(y02.multiply(v2)).add(y03.multiply(v3)).negate()
            .divide(y0101.add(y0202).add(y0303));

    Bus starBus = new Bus();
    starBus.u = v0.abs();
    starBus.theta = v0.getArgument();
    return starBus;
}

Example 12

Source File: Butterworth.java From iirj with Apache License 2.0

5 votes

public void design() {
	reset();
	double n2 = 2 * nPoles;
	int pairs = nPoles / 2;
	for (int i = 0; i < pairs; ++i) {
		Complex c = ComplexUtils.polar2Complex(1F, Math.PI/2.0
				+ (2 * i + 1) * Math.PI / n2);
		addPoleZeroConjugatePairs(c, Complex.INF);
	}

	if ((nPoles & 1) == 1)
		add(new Complex(-1), Complex.INF);
}

Example 13

Source File: TransformerModel.java From ipst with Mozilla Public License 2.0

5 votes

public StateVariable toSv2(StateVariable sv1) {
    Complex s1 = new Complex(-sv1.p, -sv1.q); // s1=p1+jq1
    Complex u1 = ComplexUtils.polar2Complex(sv1.u, Math.toRadians(sv1.theta));
    Complex v1 = u1.divide(SQUARE_3); // v1=u1/sqrt(3)
    Complex v1p = v1.multiply(ratio); // v1p=v1*rho
    Complex i1 = s1.divide(v1.multiply(3)).conjugate(); // i1=conj(s1/(3*v1))
    Complex i1p = i1.divide(ratio); // i1p=i1/rho
    Complex i2 = i1p.subtract(y.multiply(v1p)).negate(); // i2=-(i1p-y*v1p)
    Complex v2 = v1p.subtract(z.multiply(i2)); // v2=v1p-z*i2
    Complex s2 = v2.multiply(3).multiply(i2.conjugate()); // s2=3*v2*conj(i2)
    Complex u2 = v2.multiply(SQUARE_3);
    return new StateVariable(-s2.getReal(), -s2.getImaginary(), u2.abs(), Math.toDegrees(u2.getArgument()));
}

Example 14

Source File: TransformerModel.java From ipst with Mozilla Public License 2.0

5 votes

public StateVariable toSv1(StateVariable sv2) {
    Complex s2 = new Complex(-sv2.p, -sv2.q); // s2=p2+jq2
    Complex u2 = ComplexUtils.polar2Complex(sv2.u, Math.toRadians(sv2.theta));
    Complex v2 = u2.divide(SQUARE_3); // v2=u2/sqrt(3)
    Complex i2 = s2.divide(v2.multiply(3)).conjugate(); // i2=conj(s2/(3*v2))
    Complex v1p = v2.add(z.multiply(i2)); // v1'=v2+z*i2
    Complex i1p = i2.negate().add(y.multiply(v1p)); // i1'=-i2+v1'*y
    Complex i1 = i1p.multiply(ratio); // i1=i1p*ration
    Complex v1 = v1p.divide(ratio); // v1=v1p/ration
    Complex s1 = v1.multiply(3).multiply(i1.conjugate()); // s1=3*v1*conj(i1)
    Complex u1 = v1.multiply(SQUARE_3);
    return new StateVariable(-s1.getReal(), -s1.getImaginary(), u1.abs(), Math.toDegrees(u1.getArgument()));
}

Example 15

Source File: LinkData.java From powsybl-core with Mozilla Public License 2.0

2 votes

static Complex flowYshunt(Complex ysh, double u, double theta) {

        Complex v = ComplexUtils.polar2Complex(u, theta);

        return ysh.conjugate().multiply(v.conjugate().multiply(v));
    }