我编写了一段数据处理程序，如下：

clc;

clear all;

close all;

load F:\科研立项\Data\kuai2_e14.mat

s=data1(1*10^5:4*10^5)-mean(data1(1*10^5:4*10^5));

clear data1

n=length(s);                      %样本长度

fs=25000;                         %采样频率

s=s/max(abs(s));                  %归一化

pp=downsample(s,5);               %降采样

nn=length(s);

fss=fs/5;

nn=fss;                           %fix(length(pp)/5);

s=pp(nn+1:2*nn)';

max_tau=500;

%[ a,R ] = ComputeRalpha(s',20);

%[ a,R ] = Blanchit(alpha,R);

%I = floor(0.50*length(R)):length(R);

T=1000;

alpha=1/T;%1./(1:T);%Ta(floor(0.50*length(R))+1:length(R));

%clear a;

%nm=length(alpha);

%rb=zeros(nm,max_tau*2+1,T+1);

%for k=1:nm

figure;

R=frac_delay_cyclic_ac(s,T+1,alpha,max_tau);%R=cyclic_cross_covariance(s,s,alpha,max_tau);

plot(abs(R));

在其中我分析了一个特定的数据，由于数据太长，我进行了截取(s=data1(1*10^5:4*10^5)-mean(data1(1*10^5:4*10^5)); )  ，之后运行正常，绘制出了正确的谱图。

然后我想做一个GUI，有导入不同数据的功能。然后同样的进行数据截取，再进行分析、绘图。

我用菜单项LOAD来完成数据的导入，Load的callback如下：

[ fpath, fname ] =uigetfile( { '*.mat', 'Mat File' }, '请选择数据文件' );

if( ~ischar(fname) )

return;                             % 没有选择，直接返回

end

file_name = [ fname, fpath ];

handles.mydata = load( file_name );

guidata(hObject, handles);

这里应该没有错误，导入的数据在handles.mydata中。

然后我同样复制了以上数据处理部分的程序，作为一个pushbotton的callback，如下：

ha1 = findobj('tag','axes1');     %show figure in axes1

axes(ha1);

data1 = handles.mydata;

s=data1(1*10^5:4*10^5)-mean(data1(1*10^5:4*10^5));

clear data1

n=length(s);                      %样本长度

fs=25000;                         %采样频率

s=s/max(abs(s));                  %归一化

pp=downsample(s,5);               %降采样

nn=length(s);

fss=fs/5;

nn=fss;%fix(length(pp)/5);

s=pp(nn+1:2*nn)';

max_tau=500;

T=1000;

alpha=1/T;

R=frac_delay_cyclic_ac(s,T+1,alpha,max_tau);%R=cyclic_cross_covariance(s,s,alpha,max_tau);

plot(abs(R));

由于导入的数据可能长短不一，因此我在测试的时候，使用Load菜单项，依然导入的是我之前分析过的F:\科研立项\Data\kuai2_e14.mat文件，这是为了避免在截取时(s=data1(1*10^5:4*10^5)-mean(data1(1*10^5:4*10^5));)数据长度不够4*10^5个，因为已经利用m文件作过分析，因此应该画出同样的图来。但是我点击pushbotton进行数据分析时报错：

??? Index exceeds matrix dimensions.

Error in ==> Test>R_Callback at 119

s=data1(1*10^5:4*10^5)-mean(data1(1*10^5:4*10^5));

Error in ==> gui_mainfcn at 96

feval(varargin{:});

Error in ==> Test at 42

gui_mainfcn(gui_State, varargin{:});

Error in ==>

guidemfile>@(hObject,eventdata)Test('R_Callback',hObject,eventdata,guidata(hObject))

??? Error while evaluating uicontrol Callback

请问是怎么回事呢？我已经调了好久了，比较着急啊！