matlab共振峰的提取

利用matlab提取语音信号中的F0,F1,F2即前三个共振峰的属性

% ANISODIFF - Anisotropic diffusion.
%
% Usage:
%  diff = anisodiff（im niterations kappa lambda option）
%
%         im  - input image
%         nititerations 
%         kappa - conduction coefficient 20-100 ?
%         lambda - max value of .25 for stability
%         option - 1 Perona Malik diffusion equation No 1
%                  2 Perona Malik diffusion equation No 2
%
% Reference: 
% P. Perona and J. Malik. 
% Scale-space and edge detection using ansotropic diffusion.
% IEEE Transactions on Pattern Analysis and Machine Intelligence 
% 12（7）:629-639 July 1990.
%
% Author: Peter Kovesi   pk@cs.uwa.edu.au
% Department of Computer Science & Software Engineering
% The University of Western Australia
%
% June 2000         

function diff = anisodiff2（im niterations kappa lambda option）

im = double（im）;
[rowscols] = size（im）;
diff = im;
  
for i = 1:niterations
  fprintf（‘\rIteration %d‘i）;

  % Construct diffl which is the same as diff but
  % has an extra padding of zeros around it.
  diffl = zeros（rows+2 cols+2）;
  diffl（2:rows+1 2:cols+1） = diff;

  % North South East and West differences
  deltaN = diffl（1:rows2:cols+1） - diff;
  deltaS = diffl（3:rows+22:cols+1） - diff;
  deltaE = diffl（2:rows+13:cols+2） - diff;
  deltaW = diffl（2:rows+11:cols） - diff;

  % Conduction

  if option == 1

    cN = exp（-（deltaN/kappa）.^2）;
    cS = exp（-（deltaS/kappa）.^2）;
    cE = exp（-（deltaE/kappa）.^2）;
    cW = exp（-（deltaW/kappa）.^2）;

  elseif option == 2

    cN = 1./（1+exp（-（deltaN/kappa）.^2））;
    cS = 1./（1+exp（-（deltaS/kappa）.^2））;
    cE = 1./（1+exp（-（deltaE/kappa）.^2））;
    cW = 1./（1+exp（-（deltaW/kappa）.^2））;
  end


  diff = diff + lambda*（cN.*deltaN + cS.*deltaS + cE.*deltaE + cW.*deltaW）;

%  Uncomment the following to see a progression of images
%  subplot（ceil（sqrt（niterations））ceil（sqrt（niterations）） i）
%  imagesc（diff） colormap（gray） axis image

end
fprintf（‘\n‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       4155  2002-07-15 16:28  F0F1F2F3\ftrack.m

     文件        290  2007-05-24 20:27  F0F1F2F3\demo01.m

     文件        414  2007-05-19 11:10  F0F1F2F3\nanmax.m

     文件        741  2007-05-19 11:06  F0F1F2F3\nanmean.m

     文件        414  2007-05-19 11:09  F0F1F2F3\nanmin.m

     文件     812656  2002-03-11 14:27  F0F1F2F3\wioioi01.mat

     文件       1718  2002-06-18 01:47  F0F1F2F3\lpcsgram.m

     文件       1069  2002-03-11 15:20  F0F1F2F3\frmnts1.m

     文件       1928  2001-11-20 16:43  F0F1F2F3\anisodiff.m

     目录          0  2007-05-26 14:30  F0F1F2F3

----------- ---------  ---------- -----  ----

               823603                    11