t=0:1/500:5;
s_sin=sin(2*pi*1.5*t);
s_940nm=s_sin*0.2+0.75;
s_660nm=s_sin*0.1+0.8;
plot(t,s_sin);axis([0 5 -1 1]);
figure;plot(t,s_940nm);axis([0 5 0 1]);
figure;plot(t,s_660nm);axis([0 5 0 1]);
max_total=0;
max_num=0;
min_total=0;
min_num=0;
b_value=0;
b_f=0;
n_value=0;
n_f=0;
for i=1:2500
b_value=n_value;
b_f=n_f;
n_value=s_940nm(i);
if n_value>b_value
n_f=1;
elseif n_value<b_value
n_f=-1;
end
if (n_f==-1)&&(b_f==1)
max_total=max_total+s_940nm(i);
max_num=max_num+1;
end
if (n_f==1)&&(b_f==-1)
min_total=min_total+s_940nm(i);
min_num=min_num+1;
end
end
max_940nm=max_total/max_num;
min_940nm=min_total/min_num;
max_total=0;
max_num=0;
评论0