1  % 'signal_band_filter': example of passband filter (using the signal toolbox fct fdesign)


2  % frequency and bandwidth need to be modified in the function


3  %


4  % OUTPUT:


5  % DataOut: Matlab structure representing the output (filtered) field


6  %


7  %INPUT:


8  % DataIn: Matlab structure representing the output field


9 


10  %=======================================================================


11  % Copyright 20082016, LEGI UMR 5519 / CNRS UGA GINP, Grenoble, France


12  % http://www.legi.grenobleinp.fr


13  % Joel.Sommeria  Joel.Sommeria (A) legi.cnrs.fr


14  %


15  % This file is part of the toolbox UVMAT.


16  %


17  % UVMAT is free software; you can redistribute it and/or modify


18  % it under the terms of the GNU General Public License as published


19  % by the Free Software Foundation; either version 2 of the license,


20  % or (at your option) any later version.


21  %


22  % UVMAT is distributed in the hope that it will be useful,


23  % but WITHOUT ANY WARRANTY; without even the implied warranty of


24  % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the


25  % GNU General Public License (see LICENSE.txt) for more details.


26  %=======================================================================


27 


28  function DataOut=signal_bandpass_filter(DataIn,Param)


29 


30  %% request input parameters


31  if isfield(DataIn,'Action') && isfield(DataIn.Action,'RUN') && isequal(DataIn.Action.RUN,0)


32  VarNbDim=cellfun('length',DataIn.VarDimName);


33  [tild,rank]=sort(VarNbDim,2,'descend');% sort the list of input variables, putting the ones with higher dimensionality first


34  ListVarName=DataIn.ListVarName(rank);


35  VarDimName=DataIn.VarDimName(rank);


36  InitialValue=1;%default choice


37  if isfield(Param,'TransformInput') && isfield(Param.TransformInput,'VariableName')


38  val=find(strcmp(Param.TransformInput.VariableName,ListVarName));


39  if ~isempty(val);


40  InitialValue=val;


41  end


42  end


43  [s,OK] = listdlg('PromptString','Select the variable to process:',...


44  'SelectionMode','single','InitialValue',InitialValue,...


45  'ListString',ListVarName);


46  if OK==1


47  VarName=ListVarName{s};


48  DataOut.TransformInput.VariableName=VarName;


49  dlg_title = [mfilename ' filter signal along first dim ' VarDimName{s}{1}];% title of the input dialog fig


50  prompt = {'central filtering frequency';'bandwith (relative to the central frequency) '};% titles of the edit boxes


51  %default input:


52  def={'';'0.1'};% filtering frequency and relative bandwidth


53 


54  if isfield(Param,'TransformInput')% if parameters have been memorised


55  if isfield(Param.TransformInput,'CentralFrequency')


56  def{1}=num2str(Param.TransformInput.CentralFrequency);


57  end


58  if isfield(Param.TransformInput,'BandWidth')


59  def{2}=num2str(Param.TransformInput.BandWidth);


60  end


61  end


62  num_lines= 1;%numel(prompt);


63  % open the dialog fig


64  answer = inputdlg(prompt,dlg_title,num_lines,def);


65  DataOut.TransformInput.CentralFrequency=str2num(answer{1});


66  DataOut.TransformInput.BandWidth=str2num(answer{2});


67  end


68  return


69  end


70 


71 


72 


73  %% set GUI config


74  DataOut=[];


75  if strcmp(DataIn,'*')


76  DataOut.InputFieldType='1D';


77  return


78  end


79  %%%%%%%%%%%%%%%%%%%%%%%%%


80  frequency=Param.TransformInput.CentralFrequency;% frequency at which the signal needs to be filetered


81  Bw=Param.TransformInput.BandWidth;


82  DataOut=DataIn;


83  %d=fdesign.bandpass(0.8*frequency, 0.95*frequency, 1.05*frequency, 1.2*frequency, 60, 1, 60);


84  d=fdesign.bandpass((12*Bw)*frequency, (1Bw)*frequency, (1+Bw)*frequency, (1+2*Bw)*frequency, 60, 1, 60);


85  Hd=design(d);% command fvtool(Hd) to visualize the filter frequency response


86 


87  if isfield(DataIn,'U')


88  DataOut.U = filter(Hd,DataIn.U);


89  end


90  if isfield(DataIn,'V')


91  DataOut.V = filter(Hd,DataIn.V);


92  end


93 

