Matlab內插一維散亂數據

Question

我有一些簡單地組織爲二維（x，y）座標的數據。我有大量的這些數據，雖然X軸的數據對於所有點具有相同的範圍，但它不使用數據集之間相同的精確X點。我想插入每組數據，然後爲每個數據集獲取相同的X點。每當我使用Matlab嘗試插入數據時遇到問題。

所以，我的數據是這樣的：

x = [0 1 2 3 4 5 6 7] 
y = [2.2 3.7 3.9 4.1 4.2 8.9 9.1 9.3] 
xq = [0.5 0.75 2 2.25] 

其中x和y是我記錄的數據值和XQ是我想從內插數據，以獲得新的x點。

通常我使用interp1命令這樣：

f = interp1(x,y,xq); 

不幸該命令產生錯誤：「網格向量不是嚴格單調遞增」如果我使用griddata或相關命令，會發生這種情況。我知道這是因爲我的數據不是來自函數，因此我需要一種處理分散數據的技術。所以我試圖使用scatterInterpolant，但似乎這個函數似乎不適合這種類型的數據，因爲它需要x，y和一個v（值）矩陣，它比我的尺寸更大。

我不知道如何繼續，建議和建議將不勝感激。

Answer 1

雖然您確實需要提供一個真實的代碼示例和引發錯誤的輸入數據的實際版本，但我會對此做一個嘗試。代碼和數據，你投遞工作就好了：

f = 

     2.95  3.325   3.9   3.95 

你得到的錯誤意味着實際x向量不排序，因爲它是在你的例子這裏，否則它包含重複的值（即x = [0 1 1 2 3];） 。您可以使用issorted(x)命令測試第一種情況，使用any(diff(x) == 0)測試第二種情況。

第一種情況很容易解決：

[x,ix] = sort(x); 
y = y(ix); 
xq = sort(xq); 
yq = interp1(x,y,xq); 

有幾個方法來處理第二種情況，根據您的應用程序。您可以搜索重複項並將它們移動±eps，將它們平均或丟棄。

雖然沒有重複您的錯誤的數據，但我們都只是猜測。

更新：

看您發佈的數據（和對子孫後代的緣故，你應該直接編輯成你的問題），你有兩個問題。您的數據沒有排序，並有重複。正如我所說，排序很簡單。 排序如上之後，這是我在過去所做的那樣採取重複的條目的護理：

xu = unique(x); 
yu = y; 
duplicated = xu(histc(x,xu) > 1); % find the duplicated entries 
discardIndex = []; 
for k = 1:length(duplicated); 
    dupIndex = sort(find(x == duplicated(k))); % look for duplicated entries 
    keepIndex = dupIndex(1); % keep only the first one 
    discardIndex = [discardIndex dupIndex(2:end)]; % add the rest to a list 
    yu(keepIndex) = mean(y(:,dupIndex)); % take the mean of the y values at the duplicated x values 
end 
yu(discardIndex) = []; % after all is said and done, delete the duplicated entries. 

我寫了很長一段時間以前，幾乎可以肯定是不這樣做的最有效方法，但它會起作用。