-
Notifications
You must be signed in to change notification settings - Fork 0
/
MIxnyn.m
64 lines (47 loc) · 1.64 KB
/
MIxnyn.m
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
%This file is part of the MutualInformationICA package
%Copyright 2005, 2008, 2019 Alexander Kraskov, Harald Stoegbauer, Sergey Astakhov, Peter Grassberger
%This program is free software: you can redistribute it and/or modify
%it under the terms of the GNU General Public License as published by
%the Free Software Foundation, either version 3 of the License, or
%(at your option) any later version.
%This program is distributed in the hope that it will be useful,
%but WITHOUT ANY WARRANTY; without even the implied warranty of
%MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
%GNU General Public License for more details.
%You should have received a copy of the GNU General Public License
%along with this program. If not, see <https://www.gnu.org/licenses/>.
function miout=MIxnyn(x,y,kneig);
% Calculate MI value between 2 vector of any dimension (rectangular
% version)
% x....input data mxn m...channelnummer n...sampling points m<<n
% kneig... k nearest neigbor for MI algorithm
%default-values
if ~exist('kneig'), kneig=6; end
% check input data if format is correct
[Ndx,Nx]=size(x);
if Ndx>Nx
x=x';
[Ndx,Nx]=size(x);
end
[Ndy,Ny]=size(y);
if Ndy>Ny
y=y';
[Ndy,Ny]=size(y);
end
if Nx~=Ny
if Nx>Ny
N=Ny;
else
N=Nx;
end
fprintf('The two input vectors must have the same length !!!!');
fprintf('Caluculation using the %d datapoints',N);
else
N=Nx;
end
% save data for C-Programm
zwsp=[x;y]';
save zwspMIxnyn.txt zwsp -ASCII
% execute C Programm
[a unout]=unix(['MIxnyn zwspMIxnyn.txt ',num2str(Ndx),' ',num2str(Ndy),' ',num2str(N),' ',num2str(kneig)]);
miout=str2num(unout);