Why do I obtain identical graphs while for loop is accessing different sets of data stored in different column?

Question

Soeun Lee 2021-8-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/893447-why-do-i-obtain-identical-graphs-while-for-loop-is-accessing-different-sets-of-data-stored-in-differ

评论： Peter Perkins 2021-8-6

Hi all,

I have an array CHANNEL [10 58 27] and neural.spikes contain three columns each with different sets of data. Using the second for loop, I'm trying to plot the three sets of data separately, which should look different. But I'm obtaining three identical graphs while idx iterates from 1 to 3, which should allow me to plot the three columns of data stored in neural.spikes. Can anyone please tell me what is wrong with my code?

Thank you in advance for your help!

for di=CHANNEL(1,:)
    %load (append('210603_bare-210727-182826','ch',string(di)))
    neural = timetable(RawData,'SampleRate',fs);
    neural.Properties.VariableNames{1} = 'raw';
    neural.Properties.VariableUnits{1} = 'Volts';
    HighPass = 300;
    LowPass = 5000;
    [Z, P, K] = butter(5, [HighPass LowPass]/(fs/2), "bandpass");
    sos = zp2sos(Z, P, K);
    neural.raw = double(neural.raw);
    neural.spikes = sosfilt(sos,RawData);
    
    for idx=1:length(CHANNEL)
        figure;
        nexttile;
        plot(neural.Time,neural.spikes(:,idx))
    end
    
    xlim([duration(0,init_mn,init_s,init_ms) duration(0,end_mn,end_s,end_ms)])
    ylim([-1e-4 1e-4])
    
end

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

请先登录，再进行评论。

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

Answer 1

M.B 2021-8-6

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/893447-why-do-i-obtain-identical-graphs-while-for-loop-is-accessing-different-sets-of-data-stored-in-differ#answer_762162

在 MATLAB Online 中打开

Check the following:

- length(CHANNEL) is not equal to 1

- data in neural.spikes(:,1) are different from neural.spikes(:,2)?

Set a break point before "for idx=1:length(CHANNEL)" and run this code:

figure(123);
plot(neural.Time,neural.spikes(:,1)); hold on;
plot(neural.Time,neural.spikes(:,2));
plot(neural.Time,neural.spikes(:,3));

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Peter Perkins 2021-8-6

Yes, my guess is that neural.spikes = sosfilt(sos,RawData) isn't doing what you think it is. Set a break point at that line and check in the inputs and outputs.

请先登录，再进行评论。