jizhibaoluo.m

自动读取记录在文件中的数据

%%%%%%%%%求极大值，极小值，并取包络
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%
clear all;
close all
[FileName,PathName]=uigetfile('*.*');
lujing=[PathName,FileName];%输入数据文件名路径：'路径'
shuju=dlmread(lujing);%读取数据

a=shuju(:,2);

%求极大值
b=[1:1]*a(1);
num=[1:1]*1;
j=1+1;
mark=0;

for i=2:length(a)-1
    if (a(i)>a(i-1)&a(i)>a(i+1))
        b(j)=a(i);
        num(j)=i;
        mark=1;
    end
    if mark==1
        j=j+1;
    end
    mark=0;
end


%求极小值
b2=[1:1]*a(1);
num2=[1:1]*1;
j=1+1;
mark=0;

for i=2:length(a)-1
    if (a(i)<a(i-1)&a(i)<a(i+1))
        b2(j)=a(i);
        num2(j)=i;
        mark=1;
    end
    if mark==1
        j=j+1;
    end
    mark=0;
end
% %原始数据
% plot(shuju);hold on;
% %极大值数据
% plot(num,b);hold on;
% %极小值数据
% plot(num2,b2);hold on
% 
% figure

cc=fit(num',b','spline');%拟合
plot(cc);hold on
cc2=fit(num2',b2','spline');
aaa=plot(cc2);hold on
plot(shuju)