Loop calculation and output

Question

siri meka 2022-2-11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1647725-loop-calculation-and-output

评论： Yongjian Feng 2022-2-11

i have 170 files of eeg recording. i want output of each file to be stored in a seperate excel file. i am not sure how to do this. pls help

My code to initiate is as below:

txtFiles = dir('*.txt') ;
l=length(txtFiles);
N = l*10 ;
Bandsheet = zeros(N,8);
PSDSheet1 =zeros(N,16);
PSDSheet2 =zeros(N,16);
PSDSheet3 =zeros(N,16);
for n=1:170
    filename = ['00000460_', int2str(n), '.txt'];
    x=load(filename);
    figure;

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Yongjian Feng 2022-2-11

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1647725-loop-calculation-and-output#answer_893790

在 MATLAB Online 中打开

You meant this:

txtFiles = dir('*.txt') ;
l=length(txtFiles);
N = l*10 ;
Bandsheet = zeros(N,8);
PSDSheet1 =zeros(N,16);
PSDSheet2 =zeros(N,16);
PSDSheet3 =zeros(N,16);
for n=1:170
    filename = ['00000460_', int2str(n), '.txt'];
    x=load(filename);
    figure;
    % then you are going to update Bandsheet, PSDSheet1, etc here?
    % and you want to write them into different excel files?
    fname = ['00000460_', int2str(n), '.xlsx'];
    writetable(Bandsheet, fname, 'Sheet', 'Bandsheet');
    writetable(PSDSheet1, fname, 'Sheet', 'PSDSheet1');
    writetable(PSDSheet2, fname, 'Sheet', 'PSDSheet2');
    writetable(PSDSheet3, fname, 'Sheet', 'PSDSheet3');
end

6 个评论
显示 4更早的评论隐藏 4更早的评论

siri meka 2022-2-11

% this will give you all test files names

txtFiles = dir('*.txt') ;

l=length(txtFiles);

N = l*10 ;

Bandsheet = zeros(N,8);

PSDSheet1 =zeros(N,16);

PSDSheet2 =zeros(N,16);

PSDSheet3 =zeros(N,16);

for n=1:2

filename = ['00000460_', int2str(n), '.txt'];

x=load(filename);

figure;

x=load(filename);

ab= x2(0001:24000);

%%band power%%

delta = bandpower(ab,24000, [1 4]);

theta = bandpower(ab,24000,[4 8]);

alpha = bandpower(ab,24000,[8 12]);

low_beta = bandpower(ab,24000,[12 20]);

high_beta = bandpower(ab,1000,[20 35]);

beta = bandpower(ab,24000,[12 35]);

gamma = bandpower(ab,24000,[35 80]);

HFO= bandpower (ab,24000,[200 300]);

bandcomplete1 = [delta theta alpha low_beta high_beta beta gamma HFO]

Fs=24000;

N=ceil(log2(length(ab)));

L=length(ab);

%%Delta Filtering and Spectrum%%

delta1 = bandpass(ab,[1 4],1000);

figure

delta2= fft(delta1,2^N)/(L/2);

if isreal(ab)==0

fy=fy/2;

end

Power=delta2.*conj(delta2);

phy=angle(delta2);

% Frequency axis

f=(Fs/2^N)*(0:2^(N-1)-1);

%if nargin==4

deltapower= Power(1:2^(N-1));

% Find local maxima

maxIndices = islocalmax(deltapower,'MaxNumExtrema',3,'SamplePoints',f);

% Display results

clf

plot(f,deltapower,'Color',[109 185 226]/255,'DisplayName','Input data')

hold on

% Plot local maxima

plot(f(maxIndices),deltapower(maxIndices),'^','Color',[217 83 25]/255,...

'MarkerFaceColor', [217 83 25]/255, 'DisplayName', 'Local maxima')

hold off

deltavalue = plot(f(maxIndices),deltapower(maxIndices), '*','Color',[20 83 20]/255,...

'MarkerFaceColor', [55 83 25]/255, 'DisplayName', '.')

legend

x=get(deltavalue,'XData') ;

y=get(deltavalue,'YData') ;

disp(x)

disp(y)

t= [x;y]

sorted_delta = sortrows(t',2,'descend')

AllPSD2= [sorted_delta sorted_theta sorted_alpha sorted_lowbeta sorted_highbeta sorted_beta sorted_gamma sorted_hfo ]

psd11 = AllPSD2 (1,:)

psd12 = AllPSD2 (2,:)

psd13 = AllPSD2(3,:)

Bandsheet (n,:) = bandcomplete2;

PSDSheet1 (n,:) =psd21;

PSDSheet2 (n,:) =psd22;

PSDSheet3 (n,:) =psd23;

ab= x2(24001:48000);

%%band power%%

delta = bandpower(ab,24000, [1 4]);

theta = bandpower(ab,24000,[4 8]);

alpha = bandpower(ab,24000,[8 12]);

low_beta = bandpower(ab,24000,[12 20]);

high_beta = bandpower(ab,1000,[20 35]);

beta = bandpower(ab,24000,[12 35]);

gamma = bandpower(ab,24000,[35 80]);

HFO= bandpower (ab,24000,[200 300]);

disp (delta);

disp (theta);

disp (alpha);

disp (low_beta);

disp (high_beta);

disp (beta);

disp (gamma);

disp (HFO);

bandcomplete2 = [delta theta alpha low_beta high_beta beta gamma HFO]

%%Delta Filtering and Spectrum%%

delta1 = bandpass(ab,[1 4],1000);

figure

delta2= fft(delta1,2^N)/(L/2);

if isreal(ab)==0

fy=fy/2;

end

Power=delta2.*conj(delta2);

phy=angle(delta2);

% Frequency axis

f=(Fs/2^N)*(0:2^(N-1)-1);

%if nargin==4

deltapower= Power(1:2^(N-1));

% Find local maxima

maxIndices = islocalmax(deltapower,'MaxNumExtrema',3,'SamplePoints',f);

% Display results

clf

plot(f,deltapower,'Color',[109 185 226]/255,'DisplayName','Input data')

hold on

% Plot local maxima

plot(f(maxIndices),deltapower(maxIndices),'^','Color',[217 83 25]/255,...

'MarkerFaceColor', [217 83 25]/255, 'DisplayName', 'Local maxima')

hold off

deltavalue = plot(f(maxIndices),deltapower(maxIndices), '*','Color',[20 83 20]/255,...

'MarkerFaceColor', [55 83 25]/255, 'DisplayName', '.')

legend

x=get(deltavalue,'XData') ;

y=get(deltavalue,'YData') ;

disp(x)

disp(y)

t= [x;y]

sorted_delta = sortrows(t',2,'descend')

AllPSD2= [sorted_delta sorted_theta sorted_alpha sorted_lowbeta sorted_highbeta sorted_beta sorted_gamma sorted_hfo ]

psd21 = AllPSD2 (1,:)

psd22 = AllPSD2 (2,:)

psd23 = AllPSD2(3,:)

Bandsheet (n+1,:) = bandcomplete2;

PSDSheet1 (n+1,:) =psd21;

PSDSheet2 (n+1,:) =psd22;

PSDSheet3 (n+1,:) =psd23;

siri meka 2022-2-11

thanks a ton worked

Yongjian Feng 2022-2-11

Great.

Accept the answer please then.

请先登录，再进行评论。

Loop calculation and output

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

6 个评论
显示 4更早的评论隐藏 4更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

Loop calculation and output

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

6 个评论 显示 4更早的评论隐藏 4更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

6 个评论
显示 4更早的评论隐藏 4更早的评论