如何找到从阵列创建线条时的交点

下面推荐的十字路口功能对于数组达到8000个值非常适用，但如果我有一个100000个或更多数值的数组，我会耗尽内存（和I有16g内存），这很可能是由于具有交集函数的repmat命令引起的。如何找到从阵列创建线条时的交点

我试图找到从数组创建的线的交点。但不断收到错误“fzero：不是一个有效的初始包围”我正在使用八度3.8.1（这是一个开放源代码版本的matlab）下面的图片是我想用黑色在交点处圈。我是否需要在for循环中使用fzero来遍历x和y值的数组？

clear all,clf, clc,tic 

%freq array here 
x1=[20,30,40,50,60,70,80]'; 
x2=[20,30,40,50,60,70,80]'; 
y1=[2,4,3,7,1,8,4]'; 
y2=abs(y1-max(y1)); %used to switch high and low amplitude of freq 

%fit linear polynomial 
p1 = polyfit(x1,y1,1); 
p2 = polyfit(x2,y2,1); 

%calculate intersection 
x_intersect = fzero(@(x) polyval(p1-p2,x),3); 
y_intersect = polyval(p1,x_intersect); 

line(x1,y1); 
hold on; 
line(x2,y2); 
plot(x_intersect,y_intersect,'r*')

交点运作建议采取以下为数组高达8000个值伟大的工作，但如果我有10万个值的数组或更多我会耗尽内存，（我有RAM的16gig）这个最有可能是由于具有交叉函数的repmat命令引起的。

So now I'm trying to: 
1) cycle though each row in the array which represents a line 

linea1-6 xvalues = 20 to 30, 30 to 40, 40 to 50, 50 to 60, 60 to 70, 70 to 80 
linea1-6 yvalues =2 to 4, 4 to 3, 3 to 7, 7 to 1, 1 to 8, 8 to 4 
lineb1-6 xvalues = 20 to 30, 30 to 40, 40 to 50, 50 to 60, 60 to 70, 70 to 80 
lineb1-6 yvalues =6 to 4, 4 to 5, 5 to 1, 1 to 7, 7 to 0, 0 to 4 

**I'm having problems coding the for loop to work with polyfit and fzero** 


2) store the intersection values found for each line into an array. 
This should solve running out of memory issues when using large arrays

enter image description here

来源

2014-09-26 Rick T

我不熟悉你使用的语言;不过，在声明x1，x2和y1的行末尾有单引号。也许这与你的问题有关？ – Bobort 2014-09-26 14:47:06

@Ander Biguri不，它不是重复的，请阅读更新的问题原始答案不适用于大型数组。 – 2014-09-26 15:29:46

我认为这将有助于如果你能更好地描述你的输入数据。 X2总是等于x1，并且它们是等间隔和单调的吗？仅在您的示例中，y2 = abs（y1-max（y1））。如果我们能够做出这些假设，这将大大简化计算。在这种情况下不需要fzero或polyval。 – Andy 2014-09-28 14:51:22

我不知道你为什么不使用提供给您刚才的问题Finding where plots may cross with octave/matlab

但这里是这里发生了什么（从DOC）的解决方案：

如果X0是单个标量，那么探测多个附近和远处的值是，试图获得有效包围。如果这不是成功，则该功能失败。

所以发生的是你最初的猜测3离解决方案太远了。尝试改为：

>> x_intersect = fzero(@(x) polyval(p1-p2,x),30) 
x_intersect = 46.667

但是，我不知道你正在尝试做的拟合第一多项式对您的数据，它没有任何意义，我...

来源

2014-09-26 14:48:49 am304

下面推荐的十字路口功能对于数组高达8000个值是非常好的，但是如果我有一个数值为100000或更多的数组，我会耗尽内存（并且我有16g内存），这很可能是由于repmat命令与交叉口功能。所以现在我试图： 1）虽然循环表示行数组中的每一行2）将交叉点值存储在一个数组中。这应该解决使用大型数组的内存不足问题 – 2014-09-26 15:05:59

如何找到从阵列创建线条时的交点

回答

相关问题