想问问大家这个程序怎么run通??
x=xlsread('2006-2016日均数据.xlsx','Sheet1','A1:T365');
function [k,v]=findpeak(x,m,w)
if nargin<2
m='10 ';
end
nx=length(x);
if any(m=='v')
x=-x('A1:T365'); % invert x if searching for valleys
else
x=x(:); % force to be a column vector
end
dx=x(2:end)-x(1:end-1);
r=find(dx>0);
f=find(dx<0);
if ~isempty(r) && ~isempty(f) % we must have at least one rise and one fall
dr=r;
dr(2:end)=r(2:end)-r(1:end-1);
rc=ones(nx,1);
rc(r+1)=1-dr;
rc(1)=0;
rs=cumsum(rc); % = time since the last rise
df=f;
df(2:end)=f(2:end)-f(1:end-1);
fc=ones(nx,1);
fc(f+1)=1-df;
fc(1)=0;
fs=cumsum(fc); % = time since the last fall
rp=repmat(-1,nx,1);
rp([1; r+1])=[dr-1; nx-r(end)-1];
rq=cumsum(rp); % = time to the next rise
fp=repmat(-1,nx,1);
fp([1; f+1])=[df-1; nx-f(end)-1];
fq=cumsum(fp); % = time to the next fall
k=find((rs
v=x(k);
if any(m=='q') % do quadratic interpolation
b=0.5*(x(k+1)-x(k-1));
a=x(k)-b-x(k-1);
j=(a>0); % j=0 on a plateau
v(j)=x(k(j))+0.25*b(j).^2./a(j);
k(j)=k(j)+0.5*b(j)./a(j);
k(~j)=k(~j)+(fq(k(~j))-rs(k(~j)))/2; % add 0.5 to k if plateau has an even width
end
% now purge nearby peaks
if nargin>2
j=find(k(2:end)-k(1:end-1)<=w);
while any(j)
j=j+(v(j)>=v(j+1));
k(j)=[];
v(j)=[];
j=find(k(2:end)-k(1:end-1)<=w);
end
end
else
k=[];
v=[];
end
if any(m=='v')
v=-v; % invert peaks if searching for valleys
end
if ~nargout
if any(m=='v')
x=-x; % re-invert x if searching for valleys
ch='v';
else
ch='^';
end
plot(1:nx,x,'-',k,v,ch);
end