How obtain the Y value that corresponds to the largest X value for a number of plots?

Question

Ismail Qeshta 2019-2-18

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/445492-how-obtain-the-y-value-that-corresponds-to-the-largest-x-value-for-a-number-of-plots

评论： Yasasvi Harish Kumar 2019-2-18

在 MATLAB Online 中打开

Hi,

I am trying to write a code that gives me the Y value that corresponds to the largest X value for a number of plots.

I have a code that interpolates X values from Y. I just need the oposite and for the largest X values instead.

clear; clc;
Folder = cd;
N=170;
x2    = zeros(N, 10);
for k = 1:N;
  Driftt  = sprintf('Disp%d.out', k);
  Reactt  = sprintf('React%d.out', k);
  matDrift     = importdata(fullfile(Folder, Driftt));
  matReact = importdata(fullfile(Folder, Reactt));
  x1= matDrift(:,1);
  y1= matDrift(:,2);
[x3, ix] = unique(x1);
y3       = y1(ix);
A=dlmread('Time.txt');
for i=1:size(A,2)
    x2 (k,:) = interp1(y3, x3, A(:,i), 'linear');
    temp=x2(k,:);
temp(isnan(temp))=0.05;
x2(k,:)=temp;
    fid=fopen(['result_' num2str(i) '.txt'],'a');
    fprintf(fid,'%f\n',x2(k,:)); 
    fclose(fid);
end 

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Yasasvi Harish Kumar 2019-2-18

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/445492-how-obtain-the-y-value-that-corresponds-to-the-largest-x-value-for-a-number-of-plots#answer_361509

编辑：madhan ravi 2019-2-18

在 MATLAB Online 中打开

Hey,

I think the following should help. The variable big is an array of the largest x value in each row and the variable bigy is contains corresponding y values.

clear; clc;
Folder = cd;
N=170;
big = zeros(N);
bigy = zeros(N);
x2    = zeros(N, 10);
for k = 1:N;
  Driftt  = sprintf('Disp%d.out', k);
  Reactt  = sprintf('React%d.out', k);
  matDrift     = importdata(fullfile(Folder, Driftt));
  matReact = importdata(fullfile(Folder, Reactt));
  x1= matDrift(:,1);
  y1= matDrift(:,2);
[x3, ix] = unique(x1);
y3       = y1(ix);
A=dlmread('Time.txt');
for i=1:size(A,2)
    x2 (k,:) = interp1(y3, x3, A(:,i), 'linear');
    temp=x2(k,:);
temp(isnan(temp))=0.05;
x2(k,:)=temp;
for z = 1:10
    if x2(k,z)>big(z)
        big(z) = x2(k,z);
        bigy(z) = A(z,i);
    end
end
    fid=fopen(['result_' num2str(i) '.txt'],'a');
    fprintf(fid,'%f\n',x2(k,:)); 
    fclose(fid);
end 

I hope I was of some help.

Regards

2 个评论
显示无隐藏无

Ismail Qeshta 2019-2-18

编辑：Ismail Qeshta 2019-2-18

在 MATLAB Online 中打开

Hi Yasasvi,

Many thanks for your answer and suggested code. I actually keep getting the following error message:

Error using griddedInterpolant
The grid vectors must contain unique points.
Error in interp1 (line 151)
        F = griddedInterpolant(X,V,method);
Error in Interpolation (line 44)
    x2 (k,:) = interp1(y3, x3, A(:,i), 'linear');

Yasasvi Harish Kumar 2019-2-18

To interpolate you will need unique y3 and x3 values. Is there any repeating value in x3 and y3?

请先登录，再进行评论。

Answer 2

KSSV 2019-2-18

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/445492-how-obtain-the-y-value-that-corresponds-to-the-largest-x-value-for-a-number-of-plots#answer_361494

在 MATLAB Online 中打开

Read about the function max

Let X,Y be your data.

[max_x,idx] = max(X) ;
max_y = Y(idx) ;

7 个评论
显示 5更早的评论隐藏 5更早的评论

KSSV 2019-2-18

YOu need not to use interp1..you got your maximum vlaues.

Ismail Qeshta 2019-2-18

编辑：Ismail Qeshta 2019-2-18

在 MATLAB Online 中打开

Thanks. This is my latest code, but it still does not work. I get all zeroes in my output "Disp.txt" file.

clear; clc;
Folder = cd;
N=170;
x2    = zeros(N, 10);
for k = 1:N;
  X = sprintf('SaadC2DynDisp%d.out', k);
  Y  = sprintf('SaadC2DynReact%d.out', k);
  matDrift     = importdata(fullfile(Folder, X));
  matReact = importdata(fullfile(Folder, Y));
  x1= matDrift(:,1);
  y1= matDrift(:,2);
[max_x,idx] = max(X) ;
max_y = Y(idx) ;
[x3, ix] = unique(x1);
y3       = y1(ix);
A=dlmread('Time.txt');
    fid=fopen(['Disp_' num2str(i) '.txt'],'a');
    fprintf(fid,'%f\n',x2(k,:)); 
    fclose(fid);
end