% 拟合曲线找任一点
clear all
open('干扰ul100%.fig');
lh = findall(gca, 'type', 'line');% 如果图中有多条曲线,lh为一个数组
xc = get(lh, 'xdata'); % 取出x轴数据,xc是一个元胞数组
yc = get(lh, 'ydata'); % 取出y轴数据,yc是一个元胞数组
%如果想取得第2条曲线的x,y坐标
x2=xc{2};
y2=yc{2};
x1=xc{1};
y1=yc{1};
a = find(y1<5,1);
a1 = [y1(a-1),y1(a)];
b1 = [x1(a-1),x1(a)];
c = find(y2<5,1);
a2 = [y2(c-1),y2(c)];
b2 = [x2(c-1),x2(c)];
plot(a1,b1)
p(1)=spline(a1,b1);
T=[5];%数值可以输入自己需要的任意值
q1=ppval(p(1),T);
hold on
plot(a2,b2)
p(1)=spline(a2,b2);
T=[5];%数值可以输入自己需要的任意值(横坐标点)
q2=ppval(p(1),T);%纵坐标点