clear
close all
x=[0:0.2:8.8];
y=[-0.02 -0.022 -0.016 -0.016 -0.013 -0.012 -0.018 -0.023 -0.02 -0.01 -0.005 -0.009 -0.02 -0.012 -0.006 0.03 0.034 0.004 -0.019 0.017 0.185 0.536 0.924 1.158 1.12 0.845 0.312 0.179 0.019 -0.022 0.005 0.025 0.013 -0.011 -0.023 -0.01 -0.009 -0.008 -0.014 -0.021 -0.022 -0.018 -0.015 -0.019 -0.021];
xx=[0:0.05:8.8];
yy=interp1(x,y,xx,'cubic');
yym=max(yy);
y2=y./yym;
yy2=yy./yym;
plot(xx,yy2,'-');
hold on;
xlabel('x(mm)');
ylabel('I/I0')
j=1;
xmin=[];
ymin=[];
for i=2:(length(xx)-1)
if ( yy2(i)<yy2(i+1) && yy2(i)<yy2(i-1) )
xmin=[xmin xx(i)];
c(j)=i;
j=j+1;
end
end
disp(j-1);
disp (xmin)
<script src="https://mp.csdn.net/postedit"></script>
aa
最新推荐文章于 2023-08-23 21:36:22 发布