基于MATLAB的edf脑电信号阅读器，仅支持edf格式信号；可以生成时域图，时频图及频谱图；可以进行简单的信号预处理，

classdef ECoGAnalysisApp < matlab.apps.AppBase % Properties that correspond to app components properties (Access = public) UIFigure matlab.ui.Figure GridLayout matlab.ui.container.GridLayout % 左侧控制面板元素 ControlPanel matlab.ui.container.Panel LoadECoGButton matlab.ui.control.Button ECoGFilePathLabel matlab.ui.control.Label LoadEMGButton matlab.ui.control.Button EMGFilePathLabel matlab.ui.control.Label ChannelListBoxLabel matlab.ui.control.Label ChannelListBox matlab.ui.control.ListBox % 将支持多选 AnalyzeChannelButton matlab.ui.control.Button % 新增：处理操作面板 (未来用于滤波和参考) ProcessingPanel matlab.ui.container.Panel FilterTypeDropDownLabel matlab.ui.control.Label FilterTypeDropDown matlab.ui.control.DropDown CutoffFreq1Label matlab.ui.control.Label CutoffFreq1EditField matlab.ui.control.NumericEditField CutoffFreq2Label matlab.ui.control.Label CutoffFreq2EditField matlab.ui.control.NumericEditField FilterOrderLabel matlab.ui.control.Label FilterOrderEditField matlab.ui.control.NumericEditField ApplyFilterButton matlab.ui.control.Button ApplyCARButton matlab.ui.control.Button ResetProcessingButton matlab.ui.control.Button ProcessingStatusLabel matlab.ui.control.Label % 右侧显示面板元素 (使用 TabGroup) DisplayTabGroup matlab.ui.container.TabGroup TimeDomainTab matlab.ui.container.Tab TimeDomainGridLayout matlab.ui.container.GridLayout ECoGAxes matlab.ui.control.UIAxes EMGAxes matlab.ui.control.UIAxes TimeFrequencyTab matlab.ui.container.Tab SpectrogramGridContainer matlab.ui.container.GridLayout % 新属性：可滚动的父grid PowerSpectrumTab matlab.ui.container.Tab PSDGridContainer matlab.ui.container.GridLayout % 新属性：可滚动的父grid end % Properties to store data properties (Access = private) ecog_data_raw % 存储原始加载的ECoG数据 ecog_data_filtered % 存储滤波后的ECoG数据 ecog_data % 当前用于显示和分析的ECoG数据 (可能是原始、滤波后或参考后) ecog_srate ecog_chanlocs % Cell array of channel labels emg_data emg_srate emg_chanlocs % Cell array of channel labels isFilterApplied matlab.lang.OnOffSwitchState = 'off' isCARApplied matlab.lang.OnOffSwitchState = 'off' lastFilterParams struct % 存储上次应用的滤波器参数 end % Callbacks that handle component events methods (Access = private) % Button pushed function: LoadECoGButton function LoadECoGButtonPushed(app, event) [filename, filepath] = uigetfile({'*.edf;*.EDF', 'EDF Files (*.edf, *.EDF)'; '*.*', 'All Files (*.*)'}, ... 'Select ECoG EDF File'); if isequal(filename, 0) || isequal(filepath, 0) uialert(app.UIFigure, 'User cancelled ECoG file selection.', 'Info'); return; end fullpathname = fullfile(filepath, filename); app.ECoGFilePathLabel.Text = ['ECoG: ', filename]; try [data_temp, Hdr] = sload(fullpathname); if isempty(data_temp) || isempty(Hdr) error('sload did not return data or header.'); end if size(data_temp,1) == Hdr.NS raw_data = data_temp; elseif size(data_temp,2) == Hdr.NS raw_data = data_temp'; else error('Mismatch between data dimensions and number of channels in header from sload.'); end app.ecog_data_raw = raw_data; % Store raw data app.ecog_data_filtered = app.ecog_data_raw; % Initially filtered is same as raw app.ecog_data = app.ecog_data_raw; % Active data is initially raw if isfield(Hdr, 'SampleRate') && ~isempty(Hdr.SampleRate) if length(Hdr.SampleRate) == Hdr.NS app.ecog_srate = Hdr.SampleRate(1); if ~all(Hdr.SampleRate == app.ecog_srate) uialert(app.UIFigure, 'ECoG channels have different sampling rates. Using the first one. Processing might be incorrect.', 'Warning', 'Interpreter', 'none'); end elseif length(Hdr.SampleRate) == 1 app.ecog_srate = Hdr.SampleRate; else error('Cannot determine ECoG sampling rate from sload header (Hdr.SampleRate format unknown).'); end elseif isfield(Hdr, 'SPR') if isfield(Hdr, 'Dur') && Hdr.Dur(1) ~= 0 app.ecog_srate = Hdr.SPR(1) / Hdr.Dur(1); if ~all( (Hdr.SPR ./ Hdr.Dur) == app.ecog_srate) uialert(app.UIFigure, 'ECoG channels have different sampling rates calculated from SPR/Dur. Using the first one.', 'Warning', 'Interpreter', 'none'); end else app.ecog_srate = Hdr.SPR(1); uialert(app.UIFigure, 'ECoG sampling rate determined from SPR without Duration. May be inaccurate if record duration is not 1s.', 'Warning', 'Interpreter', 'none'); end else error('Cannot determine ECoG sampling rate from sload header (SampleRate or SPR not found).'); end if isfield(Hdr, 'Label') && ~isempty(Hdr.Label) app.ecog_chanlocs = cellstr(Hdr.Label); else app.ecog_chanlocs = arrayfun(@(x) sprintf('Ch %d', x), 1:size(app.ecog_data,1), 'UniformOutput', false); end app.ecog_chanlocs = strtrim(app.ecog_chanlocs); emptyLabels = cellfun('isempty', app.ecog_chanlocs); for k_idx = find(emptyLabels)' app.ecog_chanlocs{k_idx} = sprintf('Ch %d', k_idx); end app.isFilterApplied = 'off'; % Reset states app.isCARApplied = 'off'; app.updateProcessingStatusLabel(); app.updateChannelListBox(); app.plotTimeDomainECoG(); uialert(app.UIFigure, 'ECoG data loaded successfully using sload!', 'Success', 'Interpreter', 'none'); catch ME uialert(app.UIFigure, ['Error loading ECoG with sload: ', ME.message], 'Error', 'Interpreter', 'none'); fprintf('Error loading ECoG with sload. Details:\n%s\n', ME.getReport('extended', 'hyperlinks','off')); app.ECoGFilePathLabel.Text = 'ECoG: Failed to load'; app.ecog_data_raw = []; app.ecog_data_filtered = []; app.ecog_data = []; app.ecog_chanlocs = {}; app.updateChannelListBox(); cla(app.ECoGAxes); title(app.ECoGAxes, 'ECoG Signals'); end end % Button pushed function: LoadEMGButton (similar logic, not detailed for brevity now) function LoadEMGButtonPushed(app, event) [filename, filepath] = uigetfile({'*.edf;*.EDF', 'EDF Files (*.edf, *.EDF)'; '*.*', 'All Files (*.*)'}, ... 'Select EMG EDF File'); if isequal(filename, 0) || isequal(filepath, 0) uialert(app.UIFigure, 'User cancelled EMG file selection.', 'I