Subversion Repositories viirf

% MIT License

%

% Copyright (c) 2017 Mario Mauerer

%

% Permission is hereby granted, free of charge, to any person obtaining a copy

% of this software and associated documentation files (the "Software"), to deal

% in the Software without restriction, including without limitation the rights

% to use, copy, modify, merge, publish, distribute, sublicense, and/or sell

% copies of the Software, and to permit persons to whom the Software is

% furnished to do so, subject to the following conditions:

%

% The above copyright notice and this permission notice shall be included in all

% copies or substantial portions of the Software.

%

% THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR

% IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,

% FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE

% AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER

% LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,

% OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE

% SOFTWARE.

%

% VIIRF - Versatile IIR Filter

%

% This function writes matlab-generated SOS/G matrices into a file as python lists,

% Such that it can be utilized for the configuration-script IIR_Config.py.

%

% filename:

%   String of desired output-file. Don't forget ".txt"

% sos:

%   Matrix containing the SOS-coefficients. Size: [numsec 6]

% g:

%   Vector containing the section-gains and/or the final filter gain.

%   If the filter does not use any gains, this can be an empty vector.

% Example-usage:

%   export_sos_python_lists('LPFilt.txt', SOSDat, GDat);

function [] = export_sos_python_lists(filename, sos, g)

l = size(sos);

numsec = l(1);

fid = fopen(filename,'w'); % Create the file or overwrite an existing file

% First, the SOS-coefficients:

fprintf(fid, 'SOS = [');

for sec = 1:numsec

    fprintf(fid, '[');

    for c = 1:6

        fprintf(fid, '%.18f', sos(sec, c));

        if c < 6

            fprintf(fid, ', ');

        end

    end

    fprintf(fid, ']');

    if sec < numsec

        fprintf(fid, ',\n\t\t');

    end

end

fprintf(fid, ']');

% The G-Vector:

fprintf(fid, '\n\n');

fprintf(fid, 'G = [');

for sec = 1:length(g)

    fprintf(fid, '%.18f', g(sec));

    if sec < length(g)

        fprintf(fid, ', ');

    end

end

fprintf(fid, ']');

fclose(fid); % We're done

end