matrixbase::GetAmplitude

Description

Get the amplitude(modulus) of the Complex matrix.

Syntax

BOOL GetAmplitude( matrix & mAmplitude )

Parameters

mAmplitude
[output] matrix containing amplitude(modulus) matrix for this Complex matrix

Return

Returns TRUE on successful exit and FALSE on failure.

Examples

EX1

void matrixbase_GetAmplitude_ex1()
{
    int ii,jj;
    matrix mAmp;
    matrix<complex> mComplex = {
        { 1,  1+1i,  1i, -1+1i},
        {-2, -2-2i, -2i,  2-1i},
        { 0,  99,    99,  99  }
    };
    for(ii=0; ii<3; ii++)
        for(jj=0; jj<4; jj++) 
            if(mComplex[ii][jj]==99)
                mComplex[ii][jj]=NANUM;
    int rc = mComplex.GetAmplitude( mAmp );
    if(!rc) 
        printf("Error: GetAmplitude failed. rc=%d\n", rc); 
    else{
        printf("The matrix is:\n");
        for(int ii=0; ii< mAmp.GetNumRows(); ii++){
                        for(int jj=0; jj< mAmp.GetNumCols(); jj++) 
                                printf("%g  ", mAmp[ii][jj]);
                        printf("\r\n");
        }
    }
}

EX2

// Get the phase and amplitude matrix of a matrix
void matrixbase_GetAmplitude_ex2()
{
    int ii,jj;

    matrix<complex> mComplex = {
        { 1,  1+1i,  1i, -1+1i},
        {-2, -2-2i, -2i,  2-1i},
        { 0,  99,    99,  99  }
    };
    for(ii=0; ii<3; ii++)
        for(jj=0; jj<4; jj++) 
            if(mComplex[ii][jj]==99)
                mComplex[ii][jj]=NANUM;
//  Input matrix is:
//        { 1,  1+1i,  1i, -1+1i}
//        {-2, -2-2i, -2i,  2-2i}
//        { 0,  --,    --,  --  }
//  Output phase matrix is:
//        {0        0.7854   1.5708  2.35619 }
//        {3.14159 -2.35619 -1.5708 -0.463648}
//        {0        --       --      --      }
//    that is,
//        {0   PI/4    PI/2  PI*3/4}
//        {PI -PI*3/4 -PI/2 -PI/4  }
//        {0   --      --    --    }

//  Output amplitude matrix is:
//        {1  1.41421  1    1.41421}
//        {2  2.82843  2    2.82843}
//        {0  --       --   --     }
//    that is,
//        {1  sqrt(2)    1    sqrt(2)  }
//        {2  2*sqrt(2)  2    2*sqrt(2)}
//        {0  --         --   --       }

    matrix mAmp1, mAng1;
    PageBase pgbase;

    MatrixPage MatPg1;
    MatPg1.Create("Origin");
    MatrixLayer MatLy1 = MatPg1.Layers(0);
    MatLy1.SetInternalData(FSI_COMPLEX);   // Set the internal data type to complex
    Matrix<complex> Mat1(MatLy1);
    pgbase = Project.Pages(); // Get the active page
    pgbase.Rename("Original");
    Mat1 = mComplex;
    printf("  Original matrix is %s.\n",Mat1.GetName());
    printf("\n");
    
    MatrixPage MatPg2;
    MatPg2.Create("Origin");
    MatrixLayer MatLy2 = MatPg2.Layers(0);
    Matrix Mat2(MatLy2);
    pgbase = Project.Pages(); // Get the active page
    pgbase.Rename("Phase");

    int rc=mComplex.GetPhase(mAng1);      // Get phase values

    if(!rc) 
        printf("  Error: GetPhase failed.\n");
    else {
        Mat2 = mAng1;
        printf("  Phase matrix is %s.\n",Mat2.GetName());
        printf("  Observe that -PI<angle<=PI.\n"); 
        printf("  Observe that GetPhase(0)=0, and GetPhase(--)=--.\n"); 
    }
    printf("\n");
    
    MatrixPage MatPg3;
    MatPg3.Create("Origin");
    MatrixLayer MatLy3 = MatPg3.Layers(0);
    Matrix Mat3(MatLy3);
    pgbase = Project.Pages(); // Get the active page
    pgbase.Rename("Amplitude");

    rc=mComplex.GetAmplitude(mAmp1);      // Get amplitude values

    if(!rc) 
        printf("  Error: GetAmplitude failed.\n");
    else {
        Mat3 = mAmp1;
        printf("  Amplitude amplitude matrix is %s.\n",Mat3.GetName());
        printf("  Observe that 0<=Amplitude.\n"); 
        printf("  Observe that GetAmplitude(--)=--.\n"); 
    }
}

Remark

Get the amplitude(modulus) of the Complex matrix. Causes a runtime error if the underlying base type of the matrix is not Complex.

See Also

matrixbase::GetReal, matrixbase::GetImaginary, matrixbase::GetPhase, matrixbase::Conjugate, matrixbase::Cross, matrixbase::Inverse, matrixbase::MakeComplex

header to Include

origin.h