Проблема с построением графиков

[a,path]=uigetfile('93_630_6_P_1.ob','Открытие файла данных')
fid=fopen(a,'r');
[a,count]=fscanf(fid,'%d',inf);
fclose(fid);
choice=menu('Математические основы системы диагностики','Пункт_2.1 Математическое ожидание','Дисперсия','Пункт_2.2 График нормально распределённых данныx','Гистограмма','Экспериментально полученная функция плотности распределения вероятности','Теоретическая функция распределения вероятности','Функция распределения вероятности','Пункт 2.3 Автокорреляционная функция','Операция нормирования','Вычисление спектральной плотности','Вычисление автоспектра от временной реализации cm','Совмещение графиков','Пункт_2.4 Спектральный анализ','Пункт_2.5 Синхронное и синфазное усреднение');
%2.1
%Вычисление математического ожидания
M=mean(a);
%Вычисление дисперсии
D=cov(a);
%2.2
N=count; %кол-во элементов в массиве
% формируем n нормально распределенных данных
X=normrnd(215,18,N,1);
%построим график этих данных
plot(X); 
%построим гистограмму плотности распределения вероятности
n=19 %кол-во карманов
gistogramm=hist(X,n);
%экспериментально полученная функция плотности распределения вероятности
P=hist(X,n);
z=bar(P/N);
%теоретическая функция распределения вероятности
x=min(X):1:max(X); %Теоритическая
y=1/(std(X)*sqrt(2*pi))*exp(-(x-mean(X)).*(x-mean(X))/(2*std(X)*std(X)));
GP=hist(X,length(x)); %Экспериментальная    
plot(x,y,'k',x,GP/N,'r');
%Функция распределения вероятности 
P=hist(X,n);
F(1)=0;
for L=1:(n-1)
       F(L+1)=F(L)+P(L)/N;
end
bar(F);
%2.3
for K=0:N/4  %N – количество данных
    AK(K+1)=0;
       for L=1:(N-N/4) %шаг по массиву 
              AK(K+1)=AK(K+1)+a(L)*a(K+L);  %центрированный массив (здесь у меня может быть ошибка а-должен быть центрированный массив)
       end
          AK(K+1)=AK(K+1)/(N-N/4);
end
plot(AK(1:N/4));
%Операция нормирования
AKN=AK/cov(a); (здесь тоже возможна ошибка, т.к. а-должен быть центрированный массив.)
plot(AKN(1:N/4));
%Вычисление спектральной плотности
Sxx=abs(fft(AK));
NS=length(Sxx)/2; %пол длины Sxx
dt=360/1250;%Временной шаг между данными, сек.    
F=1/2/dt;  %частота Найквиста, Гц
dF=F/NS;
Sxx=Sxx.*cov(a)/sum(Sxx)./dF; %Коррекция Sxx по теореме Парсеваля
D0=dF*sum(Sxx); %Проверка Sxx по теореме Парсеваля
disp(D) %Дисперсия временной реализации
disp(D0) %Дисперсия как интеграл от спектра
KK=0:NS;
plot(dF*KK,Sxx(KK+1),'k') %График спектральной плотности в частотном диапазане [0..F]
%Вычисление автоспектра от временной реализации cm
Gxx=abs(fft(a))./(NS*dt);
ns=length(Gxx)/2;   %Половина элементов в Gxx
df=F/ns;            %Шаг по частоте, Гц
Gxx(2*ns+1)=Gxx(1); %Добавим один элемент для симметрии
Gxx=Gxx.*cov(a)/sum(Gxx)./df;  %Коррекция Gxx по теореме Парсиваля
D0=df*sum(Gxx);
disp(D)            %Дисперсия временной реализации
disp(D0)           %Дисперсия как интеграл от спектра
LL=0:ns;
plot(df*LL,Gxx(LL+1),'-m')  %График спектральной плотности в частотном диапозоне [0..F]
%Совмещение графиков
plot(dF*KK,Sxx(KK+1),'-k',df*LL,Gxx(LL+1),'-m');  
%2.4 Спектральный анализ
fn=1/2/360*1250; %частота Найквиста  
chM=abs(fft(a(3001:5048)));  %a – исходная центрированная матрица
xsh=0:fn/1024:(fn-fn/1024);
plot(xsh,chM(1:1024),'r');
%2.5 Синхронное и синфазное усреднение
%ob=fix(N/1250);   %ob – количество оборотов, N – размерность исходного массива
%for K=1:1250
 %  um(K,1)=0;
  % for L=1:ob
   %   um(K+1)=um(K,1)+mat(K+1250*(L-1),1);
   %end
%end
%um=um/ob;   
%Получаем um – массив усредненных элементов.
if choice==1
    result=M;
elseif choice==2
    result=D;
elseif choice==3
    result=plot(X);
elseif choice==4
    result=gistogramm;
elseif choice==5
    result=z;
elseif choice==6
    result=plot(x,y,'k',x,GP/N,'r'); 
elseif choice==7
    result=bar(F);
elseif choice==8
    result=plot(AK(1:N/4));
elseif choice==9
    result=plot(AKN(1:N/4));
elseif choice==10
    result=plot(dF*KK,Sxx(KK+1),'k');
elseif choice==11
    result=plot(df*LL,Gxx(LL+1),'-m');
elseif choice==12
    result=plot(dF*KK,Sxx(KK+1),'-k',df*LL,Gxx(LL+1),'-m');
elseif choice==13
    result=plot(xsh,chM(1:1024),'r');
%else choice==14
 %   result=um;
end
H=msgbox(['Результат= ' num2str(result)]);