[[email protected]]

x=-2:0.001:2;
y1 = 0.2.*exp(-0.5).*x.*cos(4.*pi.*x);
y2 = 2.*exp(-0.5).*x.*cos(pi.*x);
plot(x,y1,'-r',x,y2,'g:')
legend('y1','y2');
hold on
a=find(abs(y1-y2)<=0.005);
x3=a.*0.001-2;
y3= 0.2.*exp(-0.5).*x3.*cos(4.*pi.*x3);
plot(x3,y3,'+')
hold off

这是网上找的一种方法，但是相近点会重复标记。哪位高手，有没有更好的方法？