[未答复] 这段Gabor滤波器代码什么意思

[复制链接]
不多sali 发表于 6 天前
本帖最后由 不多sali 于 2021-1-14 17:32 编辑
  1. function out = auditory_RF(sideband,postinh)

  2. % function out = auditory_RF(LAT,BW,sideband,postinh)
  3. %
  4. % Computes and auditory RF model consisting of a Gabor, with or
  5. % without sidebands in frequency direction and a possible post-inhibition.
  6. % the parameters are:

  7. % LAT = [Latency Gabor, Latency post-inh]潜伏
  8. % BW  = [Bandwidth Gabor, Bandwidth post-inh] 带宽
  9. % sideband if 1 sidebands appear
  10. % postinh if 1 a post-inhibition is added
  11. %
  12. % time in 1 ms steps and frequency in quarter tones!

  13. % the parameters
  14. LAT = [24,50];
  15. BW = [0.035,0.04];
  16. Freq = 1;
  17. DUR = 5;

  18. if sideband
  19.   BW(1) = 0.08;
  20. else
  21.   BW(2) = 0.035;
  22. end

  23. % this is the frequency vector
  24. OCTAVES = round(1.3458*2.^[5:1/8:13]);

  25. % time in 1 ms steps and frequency in quarter tones!
  26. Fax = [0.3:1/32:1.7-0.000001];
  27. Tax = [1:80];
  28. [T,F] = meshgrid(Tax,Fax);%变成T,F两个坐标矩阵
  29. LAT = LAT;
  30. BW = BW*4;

  31. % the gabor
  32. env = exp( -( ((T-LAT(1)).^2)/(2*DUR^2) + ((F-1).^2)/(2*BW(1)^2)));
  33. osc = cos(2*pi*F*Freq);
  34. out = osc.*env;

  35. if postinh
  36.   DUR2 = DUR*1.3;
  37.   env = exp( -( ((T-LAT(2)).^2)/(2*DUR2^2) + ((F-1).^2)/(2*BW(2)^2)));

  38. out = out - osc.*env/2;  
  39. end
  40. out = out/abs(sum(out(:)));

  41. return;
复制代码

您需要登录后才可以回帖 登录 | 注册

本版积分规则

热门教程
站长推荐
快速回复 返回顶部 返回列表