Matlab project on signal operation in signal and system

Author: m-umarjan11

ssProject/StarWars3.wav

@@ -0,0 +1,13 @@
+function extracted_watermark = audio_watermark_extraction(watermarked_audio_file, watermark_length, host)
+
+    [w_sig, ~] = audioread(watermarked_audio_file);
+    
+    [a_sig, ~] = audioread(host);
+    
+    extracted_watermark = w_sig(1:watermark_length) - a_sig(1:watermark_length);
+
+    extracted_watermark = extracted_watermark > 0.5;
+    
+end
+
+

ssProject/StarWars60.wav

@@ -0,0 +1,11 @@
+function extracted_watermark = audio_watermark_extraction_emd(processed, host, watermark_length)
+    [sig1, ~] = audio_read(processed);
+    residual_1 = sig1;
+    [sig2, ~] = audio_read(host);
+    residual_2 = sig2;
+    
+    imf1 = residual_1 - mean(residual_1);
+    imf2 = residual_1 - mean(residual_2);
+    
+    extracted_watermark = imf2(1:watermark_length, 1) - imf1(1:watermark_length, 1);
+end

ssProject/adio.wav

@@ -0,0 +1,15 @@
+function sig = audio_watermarking(host_audio_file, output_file, watermark_length)
+    
+    [host_audio, sample_rate] = audioread(host_audio_file);
+    
+    watermark = generate_watermark(watermark_length);
+    
+    watermarked_audio = host_audio;
+    
+    watermark = reshape(watermark, size(watermarked_audio(1:watermark_length)));
+
+    watermarked_audio(1:watermark_length) = watermarked_audio(1:watermark_length) + watermark;
+
+    audiowrite(output_file, watermarked_audio, sample_rate);
+    sig = host_audio;
+end

ssProject/adio1.wav

@@ -0,0 +1,29 @@
+function processed_audio = audio_watermarking_emd(host_audio, watermark, output_file)
+
+    [audio_sig, sample_rate] = audioread(host_audio);
+    fs = length(audio_sig)/sample_rate;
+    audio_sig = audio_sig(:);
+    imfs = []; imf = [];
+    residual = audio_sig;
+    while any(residual)
+        imf = residual - mean(residual);
+        padding = zeros(length(imf) - length(watermark));
+        watermark_padded = [watermark padding];
+        imf = imf + watermark_padded;
+        [~, h] = freqz(imf, 1, fs);
+        if abs(h(1)) > abs(h(end))
+            imf = -imf;
+        end
+        imfs = [imfs imf];
+        residual = residual - imf;
+    end
+    
+    signal = zeros(length(imfs(1)), 1);
+
+    for i = 1:length(imfs)
+        signal = signal + imfs(i);
+    end
+
+    processed_audio = audiowrite(output_file, signal, sample_rate);
+    
+end

ssProject/adio2.wav

@@ -0,0 +1,115 @@
+classdef emdOptions %#codegen
+    %EMDOPTIONS Create option class for the EMD command
+    %
+    %   OPT = emdOptions returns the default options for EMD.
+    %
+    %   OPT = emdOptions('Properties1',Value1,'Properties2',Value2,...) uses
+    %   name/value pairs to override the default values for 'Properties1','Properties2',...
+    % This function is only for internal use
+    
+    %   Copyright 2017 The MathWorks, Inc.
+    
+    properties (Constant)
+        defaultInterpolation    =   'spline';
+    end
+    
+    properties (Access = public)
+        Interpolation
+        SiftStopCriterion = struct('SiftRelativeTolerance', 0.2, 'SiftMaxIterations', 100);
+        DecompStopCriterion = struct('MaxNumIMF', 10, 'MaxNumExtrema', 1, 'MaxEnergyRatio', 20);
+        Display = 0;
+    end
+    
+    methods
+        %----------------------------------------------------------------------
+        function this = emdOptions(varargin)
+            if(~isempty(varargin))
+                if coder.target('MATLAB')
+                    this = this.parseInputsMatlab(varargin{:});
+                else
+                    this = this.parseInputsCodegen(varargin{:});
+                end
+            else
+                % assign default values
+                this.Interpolation = this.defaultInterpolation;
+            end
+        end
+        %----------------------------------------------------------------------
+        function [this] = parseInputsMatlab(this, varargin)
+            p = inputParser;
+            addParameter(p,'Interpolation',         this.defaultInterpolation);
+            addParameter(p,'SiftRelativeTolerance', this.SiftStopCriterion.SiftRelativeTolerance);
+            addParameter(p,'SiftMaxIterations',     this.SiftStopCriterion.SiftMaxIterations);
+            addParameter(p,'MaxNumIMF',             this.DecompStopCriterion.MaxNumIMF);
+            addParameter(p,'MaxNumExtrema',         this.DecompStopCriterion.MaxNumExtrema);
+            addParameter(p,'MaxEnergyRatio',        this.DecompStopCriterion.MaxEnergyRatio);
+            addParameter(p,'Display',               this.Display);
+            
+            % parse input
+            if(mod(length(varargin),2)~=0)
+                error(message('shared_signalwavelet:emd:general:unpairedNameValue'));
+            else
+                parse(p,varargin{:});
+            end
+            
+            % assign value
+            fieldNames = fields(p.Results);
+            for i = 1:length(fields(p.Results))
+                if(isfield(this.DecompStopCriterion, fieldNames{i}))
+                    this.DecompStopCriterion.(fieldNames{i}) = p.Results.(fieldNames{i});
+                elseif(isfield(this.SiftStopCriterion, fieldNames{i}))
+                    this.SiftStopCriterion.(fieldNames{i}) = p.Results.(fieldNames{i});
+                else
+                    this.(fieldNames{i}) = p.Results.(fieldNames{i});
+                end
+            end
+        end
+        %----------------------------------------------------------------------
+        function [this] = parseInputsCodegen(this, varargin)
+            defaultSiftRelativeTolerance  =   this.SiftStopCriterion.SiftRelativeTolerance;
+            defaultSiftMaxIterations      =   this.SiftStopCriterion.SiftMaxIterations;
+            defaultMaxNumIMF              =   this.DecompStopCriterion.MaxNumIMF;
+            defaultMaxNumExtrema          =   this.DecompStopCriterion.MaxNumExtrema;
+            defaultMaxEnergyRatio         =   this.DecompStopCriterion.MaxEnergyRatio;
+            defaultDisplay                =   this.Display;
+            
+            parms = struct( 'Interpolation',            uint32(0), ...
+                'SiftRelativeTolerance',    uint32(0), ...
+                'SiftMaxIterations',        uint32(0), ...
+                'MaxNumIMF',                uint32(0), ...
+                'MaxNumExtrema',            uint32(0),...
+                'MaxEnergyRatio',           uint32(0),...
+                'Display',                  uint32(0));
+            
+            pstruct = eml_parse_parameter_inputs(parms,[],varargin{:});
+            this.Interpolation = eml_get_parameter_value(pstruct.Interpolation,...
+                this.defaultInterpolation,varargin{:});
+            this.SiftStopCriterion.SiftRelativeTolerance = eml_get_parameter_value( pstruct.SiftRelativeTolerance,...
+                defaultSiftRelativeTolerance,varargin{:});
+            this.SiftStopCriterion.SiftMaxIterations = eml_get_parameter_value( pstruct.SiftMaxIterations,...
+                defaultSiftMaxIterations,varargin{:});
+            this.DecompStopCriterion.MaxNumIMF = eml_get_parameter_value( pstruct.MaxNumIMF,...
+                defaultMaxNumIMF,varargin{:});
+            this.DecompStopCriterion.MaxNumExtrema = eml_get_parameter_value( pstruct.MaxNumExtrema,...
+                defaultMaxNumExtrema,varargin{:});
+            this.DecompStopCriterion.MaxEnergyRatio = eml_get_parameter_value( pstruct.MaxEnergyRatio,...
+                defaultMaxEnergyRatio,varargin{:});
+            this.Display = eml_get_parameter_value( pstruct.Display, defaultDisplay,varargin{:});
+        end
+    end
+    %----------------------------------------------------------------------
+    methods(Static)
+        function props = matlabCodegenNontunableProperties(classname)
+            props = {'Interpolation'};
+        end
+    end
+end
+
+
+
+
+
+
+
+
+

ssProject/audio_watermark_extraction.m

@@ -0,0 +1,3 @@
+function watermark = generate_watermark(length)
+    watermark = randi([0, 1], 1, length);
+end

ssProject/audio_watermark_extraction_emd.m

@@ -0,0 +1,29 @@
+function ownerInfo = identify_owner(output, mappingFile, host_file, len)
+   
+    watermark = audio_watermark_extraction(output, len, host_file);
+    
+    watermark_str = num2str(watermark);
+    watermark_str = reshape(watermark_str, 1, []);
+    load(mappingFile, 'person_watermark_map');
+    
+    person_ids = keys(person_watermark_map);
+    watermarks = values(person_watermark_map);
+    
+    for i = 1:length(watermarks)
+        
+        watermark1 = strrep(watermarks{i}, ' ', '');
+        watermark2 = strrep(watermark_str, ' ', '');
+        
+        
+        
+        if strcmp(watermark1, watermark2)
+            ownerInfo = person_ids{i};
+            disp('Owner Information:');
+            disp(ownerInfo);
+            return; 
+        end
+    end
+    
+    disp('No matching owner found for the watermarked audio.');
+    ownerInfo = '';
+end

ssProject/audio_watermarking.m

@@ -0,0 +1,18 @@
+function load_mapfile()
+    mapping_file = 'watermark_mapping.mat';
+    load(mapping_file, 'person_watermark_map');
+
+    keys = person_watermark_map.keys;
+    values = person_watermark_map.values;
+
+    disp('Person-Watermark Mapping:');
+    for i = 1:length(keys)
+        disp(['Person ID: ', keys{i}]);
+        disp('Watermark:');
+        for j = 1:numel(values{i})
+            fprintf('%s ', num2str(values{i}(j)));
+        end
+        disp(' ');
+        disp('-------------------------');
+    end
+end

ssProject/audio_watermarking_emd.m

@@ -0,0 +1,45 @@
+len = 50;
+
+mode = input('Enter the mode (search or store): ', 's');
+
+if strcmp(mode, 'search')
+
+    [host_file, input_path] = uigetfile('*.wav', 'Select the input audio file');
+  
+    [output, output_path] = uigetfile('*.wav', 'Enter the output audio file name');
+
+    sig1 = audio_watermarking(host_file, output, len);
+    sig2 = audio_watermark_extraction(output, len, host_file);
+    person_watermark_mapping(prop, sig2);
+    identify_owner(outstoreput, 'watermark_mapping.mat',host_file, len);
+    %load_mapfile();
+    
+elseif strcmp(mode, 'store')
+    prop = input('Enter Owner Name: ', 's');
+
+   
+    [host_file, input_path] = uigetfile('*.wav', 'Select the input audio file');
+    
+    
+    [output, output_path] = uigetfile('*.wav', 'Enter the output audio file name');
+
+    sig1 = audio_watermarking(host_file, output, len);
+    sig2 = audio_watermark_extraction(output, len, host_file);
+    person_watermark_mapping(prop, sig2);
+    identify_owner(output, 'watermark_mapping.mat', host_file, len);
+    load_mapfile();
+    subplot(3,1,1);
+    plot(sig1);
+    title('Host signal');
+    grid;
+    subplot(3,1,2);
+    plot(sig2);
+    title('Watermark signal');
+    axis([0, 100, -2, 2])
+    grid;
+    subplot(3,1,3);
+    plot(audioread(output));
+    title('Output File');
+    grid;
+    
+end

ssProject/emdOptions.m

@@ -0,0 +1,23 @@
+function person_watermark_mapping(person_id, watermark)
+    
+    mapping_file = 'watermark_mapping.mat';
+
+    if exist(mapping_file, 'file') == 2
+        load(mapping_file, 'person_watermark_map');
+        
+        watermark_str = num2str(watermark);
+        
+        person_watermark_map(person_id) = watermark_str;
+        
+    else
+        person_watermark_map = containers.Map('KeyType','char','ValueType','char');
+        
+        watermark_str = num2str(watermark);
+        
+        person_watermark_map(person_id) = watermark_str;
+    end
+
+    save(mapping_file, 'person_watermark_map');
+    
+    disp('Person-watermark mapping saved successfully.');
+end

ssProject/generate_watermark.m

@@ -0,0 +1,22 @@
+host_file = 'adio.wav';
+prop = 'adsad';
+output = 'adio1.wav';
+len = 50;
+sig1 = audio_watermarking(host_file, output, len);
+sig2 = audio_watermark_extraction(output, len, host_file);
+person_watermark_mapping(prop, sig2);
+identify_owner('w1.wav', 'watermark_mapping.mat', 'w.wav', len);
+load_mapfile();
+subplot(3,1,1);
+plot(sig1);
+title('Host signal');
+grid;
+subplot(3,1,2);
+plot(sig2);
+title('Watermark signal');
+axis([0, 100, -2, 2])
+grid;
+subplot(3,1,3);
+plot(audioread(output));
+title('Ouput File');
+grid;