generatePRRFile.m

Go to the documentation of this file.

function generatePRRFile( filePath, targetName )
    s2p = sparameters( filePath );
 
    freqSamples = s2p.Frequencies * 1e-9;
 
    fitData = struct("poles", {}, "residues", {}, "remainder", {}, "proportional", {});
    S = repmat( fitData, s2p.NumPorts, s2p.NumPorts );
 
    prrFile = fopen( targetName, 'w' );
    fprintf( prrFile, "%d", 50 );
        
    % Find a complex rational model for freq domain data
    for a = 1 : 1 : s2p.NumPorts
        for b = 1 : 1 : s2p.NumPorts
            [ S(a, b).poles, S(a, b).residues, S(a, b).remainder, S(a, b).proportional ] = findComplexRationalApproximation(freqSamples, squeeze(s2p.Parameters( a, b, : )) );     
            
            fprintf( prrFile, "\n%g %g\n", real( S(a, b).remainder ), imag( S(a, b).remainder ) );      
            fprintf( prrFile, "%g %g", real( S(a, b).poles( 1 ) ), imag( S(a, b).poles( 1 ) ) );
            for i = 2 : 1 : length( S(a, b).poles )
                fprintf( prrFile, " %g %g", real( S(a, b).poles( i ) ), imag( S(a, b).poles( i ) ) );
            end
 
            fprintf( prrFile, "\n%g %g", real( S(a, b).residues( 1 ) ), imag( S(a, b).residues( 1 ) ) );
            for i = 2 : 1 : length( S(a, b).residues )
                fprintf( prrFile, " %g %g", real( S(a, b).residues( i ) ), imag( S(a, b).residues( i ) ) );
            end
        end
    end
end