爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 6221|回复: 8

[讨论] 如何找出点状图中 三条曲线的交点

[复制链接]

新浪微博达人勋

发表于 2013-2-24 17:16:09 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
data=load('g:\18002\hunan\pdf\cpd.txt');
data1=data(:,1);
data2=data(:,2);
data3=data(:,3);
data4=data(:,4);
data5=data(:,5);
data6=data(:,6);
data7=data(:,7);
plot(data2,data3,'--bd')
hold on
plot(data4,data5,'-.r*')
hold on
plot(data6,data7,':kh')
222222222222222222222222222222.png

没有进行拟合
想找出它们的交点(多个)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-2-25 11:59:16 | 显示全部楼层
你这个data2,data4,data6是一样的吗?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-2-25 12:26:54 | 显示全部楼层
如果一样的话,可以这样:
y=find(abs(data3-data5)<1e-5);  %找到data3和data5数据接近的点
y=find(abs(data7(y)-data5(y))<1e-5); %找到和data7数据接近的点
plot(data2(y),data7(y),' ') %假设data2 data4 data6数值一样

如果不一样的话,就比较麻烦了吧,我想可以通过分段线性插值把点插密,然后再用上述方法找相近的点。
函数interp(x,y,x1,'linear')
自己写的一段程序,没运行过,你试试:
data3_new=interp(data2,data3,min(data2):0.01:max(data2),'linear'); %步长你可以自己设置
data5_new=interp(data4,data5,min(data4):0.01:max(data4),'linear');
data7_new=interp(data6,data7,min(data6):0.01:max(data6),'linear');
y=find(abs(data3_new-data5_new)<1e-5);
y=find(abs(data7_new(y)-data5_new(y))<1e-5);
plot(data2(y),data7(y),'  ');
当然数据大小要一致才行,不一致就找三条曲线x公共的就行了。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-2-25 23:41:14 | 显示全部楼层
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-2-26 08:20:47 | 显示全部楼层
本帖最后由 爱的侍者 于 2013-2-26 08:47 编辑
Aires 发表于 2013-2-25 11:59
你这个data2,data4,data6是一样的吗?


不是一样的,但是相差不大,也有差别的~横坐标
数据大小一致是指??
横坐标数据的个数不一定相同

??? Operands to the || and && operators must be convertible to logical scalar values.
Error in ==> interp at 36
if l < 1 || r < 1 || cutoff <= 0 || cutoff > 1

报出这个
data1=data(:,1);
data2=data(:,2);
data3=data(:,3);
data4=data(:,4);
data5=data(:,5);
data6=data(:,6);
data7=data(:,7);
data3_new=interp(data2,data3,min(data2):0.01:max(data2),'linear'); %步长你可以自己设置
data5_new=interp(data4,data5,min(data4):0.01:max(data4),'linear');
data7_new=interp(data6,data7,min(data6):0.01:max(data6),'linear');
y=find(abs(data3_new-data5_new)<1e-5);
y=find(abs(data7_new(y)-data5_new(y))<1e-5);
figure
plot(data2(y),data7(y),'  ');

data1 没用 data2 data4 data6 为横坐标
我的数据中有NaN的缘故??

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-2-26 16:37:43 | 显示全部楼层
kongfeng0824 发表于 2013-2-25 23:41
版主你指的一样是数据个数一样吗?

不是,我指的是值是不是一样的,也就是x坐标是不是一样的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-2-26 16:47:11 | 显示全部楼层
爱的侍者 发表于 2013-2-26 08:20
不是一样的,但是相差不大,也有差别的~横坐标
数据大小一致是指??
横坐标数据的个数 ...

我这个data3_new data5_new data7_new 数组中元素个数是一样的,你这个我觉得真的听麻烦的。不行你插值完之后找到x的公共部分然后再找三条线的交点,我觉得会比较麻烦
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-2-26 22:51:43 | 显示全部楼层
Aires 发表于 2013-2-26 16:47
我这个data3_new data5_new data7_new 数组中元素个数是一样的,你这个我觉得真的听麻烦的。不行你插值完 ...

到底咋求嘛~!~
你的data3_new data5_new data7_new  元素个数为啥一样,对应 data2 data 4 data6 0.01 的个数啊~
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-3-1 11:45:28 | 显示全部楼层
爱的侍者 发表于 2013-2-26 22:51
到底咋求嘛~!~
你的data3_new data5_new data7_new  元素个数为啥一样,对应 data2 data 4 data6 0.01  ...

我是说data3_new data5_new data7_new如果个数和数值都一样,会比较好求。你这个我觉得会麻烦很多,你自己想想试试吧,图形上看三条曲线很相似啊,不行就手动找吧
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表