audio_toolbox/filterdBA.m at master · amilo/audio_toolbox · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100

inputdir = '/Users/amilo/Desktop/karl/P01/yourfolder/filt/';
S = dir(fullfile(inputdir,'*.wav'));

Fs=44100;
[HdAw, HdBw, HdCw] = ABCdsgn(Fs);

for k = 1:numel(S)
    fnm = fullfile(inputdir,S(k).name);
    disp(fnm);
    [x,Fs] = audioread(fnm);
    y = filter(HdAw, x);
    filename = strcat(fnm(1:end-4), '_A.wav');
    disp(strcat('filtered file: ', filename));
    audiowrite(filename,y,Fs);
    ...
end


%
% for k = 1:numel(S)
%     fnm = fullfile(inputdir,S(k).name);
% %     txt = load(fnm);
%     filter(fnm)
%     disp(fnm);
%     ...
% end
% file='/Users/amilo/Desktop/karl/participants/yourfolder/251119_P1-left.wav';
%  filter(file);
%
%  function filter = filterA(file)


% function filter(file)


% HdAw

% B = coeffs(HdAw)

% B.Stage1.Stage1
% B.Stage1.Stage2
% B.Stage2


% Create a time record
% buf=randn(1, 44100);
% t=1/Fs*(0:(Fs-1));
% buf2=buf;


% t=0:length(y);
% length(y)

% Apply the A-weighting filters
% buf = filter(HdAw, buf );


% Plot the results!
% figure(1); plot(t(1:15), buf2(1:15), 'k'); hold on; plot(t(1:15), buf(1:15), 'g'); legend('Linear-weighted Time Record', 'A-weighted Time Record');
% figure(1); plot(t(1:512), x(1:512), 'k'); hold on; plot(t(1:512), y(1:512), 'g'); legend('Linear-weighted Time Record', 'A-weighted Time Record');

% end


%
% groupArray = [group1,group2,group3];
%
% index=0;
%
% countMatrix=[14500,3];
%
% for SUF = groupArray %[subgroup] %SUF
%     SUF = num2str(SUF);
%     index=index+1;
%
% %SUF='n';
% run(['sensorLog18Rev_' SUF '.m'])
%
% filename = ['audio_' SUF '.bin'];
%
%  end
%
%
%
% average=std(countMatrix ,0,2);
% %average=mean(countMatrix ,2);
%
%
%
% size(average)
% size(Tred)
%