Explotoooooooooooooooooo
What I did this afternoon
PS: If adiviniais that is the image of the explosion I shall give a prize
PD2: I put the post datas pq I doubt anyone here read the trash below XD
x = 100000;
num_muestra = input ('with that patient wants to work? 26 32 34 42 46 47 53 56''\\ n');
num_muestra = num_muestra * 1;
matrix_muestra
% = [1 26, February 1932, March 1934, April 1942, May 1946, June 1947, July 1953, August 1956];
% num_muestra = matrix_muestra (n_muestra, 2);
if num_muestra
== 26 & amp; nbsp; fid = fopen ('26 xls ',' r '.)
[vector, txt, raw] = xlsread ('26 xls.');
fid = fopen ('P26Sig.dat', 'r');
elseif num_muestra == 32
fid = fopen ('32 xls ','. r ');
[vector, txt, raw] = xlsread ('32 xls.'), fid = fopen
('P32Sig.dat', ' r ');
elseif num_muestra == 34
fid = fopen ('34 xls.', 'r');
[vector, txt, raw] = xlsread ('34 xls. ');
fid = fopen (' P34Sig.dat ','r ');
elseif num_muestra == 42
fid = fopen ('42 xls.', 'r');
[vector, txt, raw] = xlsread ('42 xls. ');
fid = fopen (' P42Sig.dat ',' r ');
elseif num_muestra == 46
fid = fopen ('46 xls', 'r'.)
[vector, txt, raw] = xlsread ('46 xls. '), fid = fopen
(' P46Sig.dat ',' r ');
elseif num_muestra == 47
fid = fopen ('47.xls ',' r ');
[vector, txt, raw] = xlsread ('47 xls.');
; fid = fopen ('P47Sig.dat', 'r');
elseif num_muestra == fid = 53
fopen ('53 xls ',' r '.)
[vector, txt, raw] = xlsread ('53 xls.');
fid = fopen ('P53Sig.dat', 'r');
elseif num_muestra == 56
; fid =fopen ('56. xls ',' r ');
[vector, txt, raw] = xlsread ('56. xls');
fid = fopen ('P56Sig.dat', 'r');
end
num_instante vector = input ('enter the moment''\\ n ');
inicio_segmento = vector (num_instante, 2);
final_segmento = vector (num_instante, 3);
% n = num_instante inicio_segmento
% n1 = n2 = final_segmento
%
r1 = rem (inicio_segmento, 100000);
r2 = rem (final_segmento, 100000);
counter = fix ((inicio_segmento) / 100000);
% counter = fix ((n1-r1) / 100000), -> Because it is not necessary to fix the
% remains of the division
% to avoid collapse of memory we calculate the vector
% 100,000 100,000 in figures. counter here would be the integer part of n1/100.000
for i = 1: counter
a = fread (fid, x, 'int16');
end mean = mean (a);
a2 = a-half;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%% PLOT TO SEGMENT WE WANT %%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %
bi = 1;
b = zeros (r2-r1, 1);
for i = r1, r2
b (bi, 1) = a2 (i, 1);
bi = bi + 1;
end
%% b is the vector of the segment that we chose
% plot (b, 'k');
% axis ([xmin xmax 32 768 -32 768]);
% xlabel ('time (s)');
% desviacion_tipica = std (b);
figure (1)
if vector (num_instante, 4) ~ = 999
subplot (2,1,1) hold on
%
vector_ejex_absoluto = (counter * 100 000 +1:1: (counter +1) * 100000) ';
% plot (vector_ejex_absoluto, a2,' k ');
plot (a2,' k ');
title ([, ['segment', int2str (counter)], 'test 00000']);
% axis ([r1 r2 32 768 -32 768]);
xlabel ('time (s)') ;
plot ([r1, r1], [min (b) max (b)], 'r');
plot ([r2, r2], [min (b) max (b) ], 'r');
plot ([r1, r2], [min (b), min (b)], 'r');
plot ([r1, r2], [max (b) max (b)], 'r ');
hold off subplot (2,1,2)
plot (b,' k ') title
([,[' Segment from', int2str (inicio_segmento )], 'a', ['' int2str (final_segmento)], 'test']);
% axis ([n1 n2 32 768 -32 768]);
xlabel ('time (s) ');
else
plot (b,' k ');
([,[' Segment title', int2str (counter)], 'test 00000 32' ]);
axis ([xmin xmax 32 768 -32 768]);
xlabel ('time (s)');
end
%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% CHTML %%%%%%%%%%a = PX ';
else data = PX;
end;
tamano_PX = size (PX) M = tamano_PX
(1.1);
N = tamano_PX (1.2);
% New PSD parameters Probabilistic
%
mean index = 1: length (data);
pmean = sum (index. * data) / sum (data);
Probabilistic
%
index2 = standard deviation (index-pmean). ^ 2;
var = sum (index2. * data) / sum (data);
var = sqrt (var);
% Symmetry parameter (3rd moment)
index3 = (index-pmean). ^ 3;
Symmetry = sum (index3. * data) / sum (data);
% Median (f50%) and frequency spread (f95%)
datacum = cumSum (data);
f50 = index (datacum / datacum (length (datacum) )> 0.5);
f50 = f50 (1);
f95 = index (datacum / datacum(length(datacum)) > 0.95);
f95 = f95(1);
%Flatness measure (in f95%)
flatness = datacum(f95-1)/(max(data)*(f95-1));
%% NEW PARAMETERS
%Ratio of power above/below 800Hz
k800 = floor(800/deltaf);
pab800 = sum(data(k800:length(data)));
pbe800 = sum(data(1:k800-1));
ratio800 = pab800 / pbe800;
%Fraction of power in the 100-500Hz band
k100 = floor(100/deltaf);
k500 = floor(500/deltaf);
p100_500 = sum(data(k100:k500-1)) / sum(data);
%Same, normalized to the power out of 100-500 band
ratio100_500 = sum(data(k100:k500-1)) / ( sum(data(1:k100-1)) + sum(data (k500: length (data))));
% Below 500Hz Fraction of power
k500 = floor (500/deltaf)
pbe500 = sum (data (1: k500-1)) / sum (data); Same
%, normalized to the power Above 500Hz
ratio500 = sum (data (1: k500-1)) / sum (data (k500: length (data)));
% NEW! FPEAK
index = 1: length (data);
[m, i] = max (data);
FPEAK = i-1;
% return pars Results (1) = pmean * deltaf;
pars (2) = var * deltaf;
pars (3) = (abs (Symmetry)) ^ (1 / 3) * deltaf;
pars (4) = flatness * 1000;
pars (5) = f50 * deltaf;
pars (6) = f95 * deltaf;
pars (7) = ratio800 * 1000;
pars (8) = p100_500;
pars (9) = ratio100_500;
pars (10) = pbe500;
pars (11) = ratio500;
pars (12) = FPEAK * deltaf;
% PLOTS
% Calculate
frequency vector f = 0: (Fs / 2) / length (data) : (Fs / 2) / length (data) * (length (data) -1);
% Plot frequency waveform
figure (FIG);
subplot (M, N, SUBFIG)
plot (f, data, 'r');
% plot (f, 10 * log10 (data));
grid on a = axis;
axis ([f (1) f (length (f)) to (3) 1.2 * max (data)]);
xlabel ('frequency (Hz)');
title ('title');
drawing
% Parameters (MATLAB 4.2)
a = axis;
p = num2str (pars (1), 5);
p = ['fMean =' p 'Hz'];
text (0.5 * (a (2)-a (1)), a (4) -0.08 * (a (4)-a (3 )), p);
p = num2str (paras (2), 5);
p = ['Stdev = 'p' Hz '];
text (0.5 * (a (2)-a (1)), a (4) -0.16 * (a (4)-a (3)), p);
p = num2str (pars (3), 5);
p = ['Symmetry =' p];
text (0.5 * (a (2)-a (1)), a (4) -0.24 * (a ( 4)-a (3)), p);
p = num2str (pars (4), 5);
p = ['Flatness =' p];
text (0.5 * (a (2)-a ( 1)), a (4) -0.32 * (a (4)-a (3)), p);
p = num2str (pars (5), 5);
p = ['f (50%) = 'p' Hz '];
text (0.5 * (a (2)-a (1)), a (4) -0.40 * (a (4)-a (3)), p);
p = num2str (pars (6), 5);
p = ['f (95%) =' p 'Hz'];
text (0.5 * (a (2)-a (1)), a (4) -0.48 * (a (4)-a (3)), p);
p = num2str (pars (7), 5);
p = ['ratio800 =' p];
text (0.5 * (a (2)-a (1)), a (4) -0.56 * (a (4)-a (3)), p);
p = num2str (pars (8), 5);
p = [ 'p100-500 =' p];
text (0
num_muestra = input ('with that patient wants to work? 26 32 34 42 46 47 53 56''\\ n');
num_muestra = num_muestra * 1;
matrix_muestra
% = [1 26, February 1932, March 1934, April 1942, May 1946, June 1947, July 1953, August 1956];
% num_muestra = matrix_muestra (n_muestra, 2);
if num_muestra
== 26 & amp; nbsp; fid = fopen ('26 xls ',' r '.)
[vector, txt, raw] = xlsread ('26 xls.');
fid = fopen ('P26Sig.dat', 'r');
elseif num_muestra == 32
fid = fopen ('32 xls ','. r ');
[vector, txt, raw] = xlsread ('32 xls.'), fid = fopen
('P32Sig.dat', ' r ');
elseif num_muestra == 34
fid = fopen ('34 xls.', 'r');
[vector, txt, raw] = xlsread ('34 xls. ');
fid = fopen (' P34Sig.dat ','r ');
elseif num_muestra == 42
fid = fopen ('42 xls.', 'r');
[vector, txt, raw] = xlsread ('42 xls. ');
fid = fopen (' P42Sig.dat ',' r ');
elseif num_muestra == 46
fid = fopen ('46 xls', 'r'.)
[vector, txt, raw] = xlsread ('46 xls. '), fid = fopen
(' P46Sig.dat ',' r ');
elseif num_muestra == 47
fid = fopen ('47.xls ',' r ');
[vector, txt, raw] = xlsread ('47 xls.');
; fid = fopen ('P47Sig.dat', 'r');
elseif num_muestra == fid = 53
fopen ('53 xls ',' r '.)
[vector, txt, raw] = xlsread ('53 xls.');
fid = fopen ('P53Sig.dat', 'r');
elseif num_muestra == 56
; fid =fopen ('56. xls ',' r ');
[vector, txt, raw] = xlsread ('56. xls');
fid = fopen ('P56Sig.dat', 'r');
end
num_instante vector = input ('enter the moment''\\ n ');
inicio_segmento = vector (num_instante, 2);
final_segmento = vector (num_instante, 3);
% n = num_instante inicio_segmento
% n1 = n2 = final_segmento
%
r1 = rem (inicio_segmento, 100000);
r2 = rem (final_segmento, 100000);
counter = fix ((inicio_segmento) / 100000);
% counter = fix ((n1-r1) / 100000), -> Because it is not necessary to fix the
% remains of the division
% to avoid collapse of memory we calculate the vector
% 100,000 100,000 in figures. counter here would be the integer part of n1/100.000
for i = 1: counter
a = fread (fid, x, 'int16');
end mean = mean (a);
a2 = a-half;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%% PLOT TO SEGMENT WE WANT %%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %
bi = 1;
b = zeros (r2-r1, 1);
for i = r1, r2
b (bi, 1) = a2 (i, 1);
bi = bi + 1;
end
%% b is the vector of the segment that we chose
% plot (b, 'k');
% axis ([xmin xmax 32 768 -32 768]);
% xlabel ('time (s)');
% desviacion_tipica = std (b);
figure (1)
if vector (num_instante, 4) ~ = 999
subplot (2,1,1) hold on
%
vector_ejex_absoluto = (counter * 100 000 +1:1: (counter +1) * 100000) ';
% plot (vector_ejex_absoluto, a2,' k ');
plot (a2,' k ');
title ([, ['segment', int2str (counter)], 'test 00000']);
% axis ([r1 r2 32 768 -32 768]);
xlabel ('time (s)') ;
plot ([r1, r1], [min (b) max (b)], 'r');
plot ([r2, r2], [min (b) max (b) ], 'r');
plot ([r1, r2], [min (b), min (b)], 'r');
plot ([r1, r2], [max (b) max (b)], 'r ');
hold off subplot (2,1,2)
plot (b,' k ') title
([,[' Segment from', int2str (inicio_segmento )], 'a', ['' int2str (final_segmento)], 'test']);
% axis ([n1 n2 32 768 -32 768]);
xlabel ('time (s) ');
else
plot (b,' k ');
([,[' Segment title', int2str (counter)], 'test 00000 32' ]);
axis ([xmin xmax 32 768 -32 768]);
xlabel ('time (s)');
end
%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% CHTML %%%%%%%%%%a = PX ';
else data = PX;
end;
tamano_PX = size (PX) M = tamano_PX
(1.1);
N = tamano_PX (1.2);
% New PSD parameters Probabilistic
%
mean index = 1: length (data);
pmean = sum (index. * data) / sum (data);
Probabilistic
%
index2 = standard deviation (index-pmean). ^ 2;
var = sum (index2. * data) / sum (data);
var = sqrt (var);
% Symmetry parameter (3rd moment)
index3 = (index-pmean). ^ 3;
Symmetry = sum (index3. * data) / sum (data);
% Median (f50%) and frequency spread (f95%)
datacum = cumSum (data);
f50 = index (datacum / datacum (length (datacum) )> 0.5);
f50 = f50 (1);
f95 = index (datacum / datacum(length(datacum)) > 0.95);
f95 = f95(1);
%Flatness measure (in f95%)
flatness = datacum(f95-1)/(max(data)*(f95-1));
%% NEW PARAMETERS
%Ratio of power above/below 800Hz
k800 = floor(800/deltaf);
pab800 = sum(data(k800:length(data)));
pbe800 = sum(data(1:k800-1));
ratio800 = pab800 / pbe800;
%Fraction of power in the 100-500Hz band
k100 = floor(100/deltaf);
k500 = floor(500/deltaf);
p100_500 = sum(data(k100:k500-1)) / sum(data);
%Same, normalized to the power out of 100-500 band
ratio100_500 = sum(data(k100:k500-1)) / ( sum(data(1:k100-1)) + sum(data (k500: length (data))));
% Below 500Hz Fraction of power
k500 = floor (500/deltaf)
pbe500 = sum (data (1: k500-1)) / sum (data); Same
%, normalized to the power Above 500Hz
ratio500 = sum (data (1: k500-1)) / sum (data (k500: length (data)));
% NEW! FPEAK
index = 1: length (data);
[m, i] = max (data);
FPEAK = i-1;
% return pars Results (1) = pmean * deltaf;
pars (2) = var * deltaf;
pars (3) = (abs (Symmetry)) ^ (1 / 3) * deltaf;
pars (4) = flatness * 1000;
pars (5) = f50 * deltaf;
pars (6) = f95 * deltaf;
pars (7) = ratio800 * 1000;
pars (8) = p100_500;
pars (9) = ratio100_500;
pars (10) = pbe500;
pars (11) = ratio500;
pars (12) = FPEAK * deltaf;
% PLOTS
% Calculate
frequency vector f = 0: (Fs / 2) / length (data) : (Fs / 2) / length (data) * (length (data) -1);
% Plot frequency waveform
figure (FIG);
subplot (M, N, SUBFIG)
plot (f, data, 'r');
% plot (f, 10 * log10 (data));
grid on a = axis;
axis ([f (1) f (length (f)) to (3) 1.2 * max (data)]);
xlabel ('frequency (Hz)');
title ('title');
drawing
% Parameters (MATLAB 4.2)
a = axis;
p = num2str (pars (1), 5);
p = ['fMean =' p 'Hz'];
text (0.5 * (a (2)-a (1)), a (4) -0.08 * (a (4)-a (3 )), p);
p = num2str (paras (2), 5);
p = ['Stdev = 'p' Hz '];
text (0.5 * (a (2)-a (1)), a (4) -0.16 * (a (4)-a (3)), p);
p = num2str (pars (3), 5);
p = ['Symmetry =' p];
text (0.5 * (a (2)-a (1)), a (4) -0.24 * (a ( 4)-a (3)), p);
p = num2str (pars (4), 5);
p = ['Flatness =' p];
text (0.5 * (a (2)-a ( 1)), a (4) -0.32 * (a (4)-a (3)), p);
p = num2str (pars (5), 5);
p = ['f (50%) = 'p' Hz '];
text (0.5 * (a (2)-a (1)), a (4) -0.40 * (a (4)-a (3)), p);
p = num2str (pars (6), 5);
p = ['f (95%) =' p 'Hz'];
text (0.5 * (a (2)-a (1)), a (4) -0.48 * (a (4)-a (3)), p);
p = num2str (pars (7), 5);
p = ['ratio800 =' p];
text (0.5 * (a (2)-a (1)), a (4) -0.56 * (a (4)-a (3)), p);
p = num2str (pars (8), 5);
p = [ 'p100-500 =' p];
text (0
0 comments:
Post a Comment