外挿値の取得について

x軸にc, y軸にDをプロットし（実線）、y軸のデータがない部分を外挿でプロットしました（点線）。

外挿のコードは以下の通りです。

その際、以下のプロット図の点線のxy数値を取得したいのですが、どうすればよいでしょうか？

N = (linspace(1450,1460,100))';
c_ext = interp1(c, D, N, 'spline', 'extrap')

Answer 1

[Hernia Baby]

論理演算子を活用するのはいかがでしょうか？

以下、例を示します

x = linspace(0,10,100)';
y = sin(x) + rand(length(x),1);

xが5より小さいものだけを抜き出します

x1 = x(x<5);
y1 = y(x<5); 

プロットします

plot(x,y,':',x1,y1,'--')
legend({'生データ','切り取り'});

[Hernia Baby]

外装前のデータを利用して

c_min = min(c);
N_ext_only = N(N<c_min);
c_ext_only = c_ext(N<c_min);

等で解決すると思います。

[Yu]

早速コメントくださりありがとうございます。

申し訳ありません、こちらの説明が足りておりませんでした。

実際には、2セットのデータ（（c,D）および（c_ext,N））からそれぞれプロットを作成し、プロットを重ね合わせています。

N、c_extの中身を見ても、c やDとは全くことなるデータ群であり、プロットで見られるようなc, D値を取得したいと思っています。

参考に「samples.mat」にデータを添付いたしました。

ご存知でしたらご教示いただけますと幸いです。

   %c, D = 生データ（式は省略）
    
    N = (linspace(1450,1460,100))';
    c_ext = interp1(c, D, N, 'spline', 'extrap')

    plot(c,D)
    
    hold on
    plot(N,c_ext,'-.b')
    hold off
    

[Hernia Baby]

つまり繋げたいってことですかね

clc,clear,close all;
load('samples.mat');

ここで繋げてます

x_new = [N; c];
y_new = [c_ext; D];
[x_new_sort, idx] = sort(x_new);
y_new_sort = y_new(idx);

プロットして軸反転など行ってます

plot(x_new_sort,y_new_sort)
axis ij
xlim([1450 1540])
ylim([0 500])

[Yu]

sortとidxの組み合わせ方はとても勉強になりました。

お陰様で解決しました。

丁寧に説明くださりありがとうございました！