i want to find the values in the x axis for the respective y axis
10 次查看(过去 30 天)
显示 更早的评论
i have written a code for autocorrelation and got the plot for .nc file . Now i want to find the value of x axis for the respective y axis(at 0.5 and 0 )i.e.value for the first intersection of yaxis at 0.5 and 0 with the x axis.
can anyone please help me in getting this value using code...
thank you in advance.
Y = xlsread('mat');
% m - columns
% n - rows
[n,m] = size(Y);
Ym = zeros(n,m);
N = sum(Y==Y);
lag = min(N) - 1; % lag is 20
acf = zeros(lag+1,m); % auto correlation results
for i = 1:m
n = N(i);
ymean = mean(Y(1:n,i));
Ym(1:n,i) = Y(1:n,i)-ymean;
for k = 0:lag
tot = 0;
idivisor = 0;
y = Y(:,i);
for j = 1:(n-k)
tot = tot + Ym(j,i)*Ym(j+k,i);
end
for j = 1:n
idivisor = idivisor + (Ym(j,i)*Ym(j,i));
end
acf(k+1,i) = tot/idivisor;
end
end
for i = 1:m
plot(acf(:,i)),xlabel('Lag'),ylabel('SWH'),title('Autocorrelation');
grid minor;
set(gca,'xminortick','on');
set(gca,'xtick',0:50:1200);
filename = sprintf('%d',i);
print(filename,'-dpng');
end
0 个评论
回答(1 个)
Jakob B. Nielsen
2020-2-4
You seem to do your autocorrelation in steps of 1. What is the precision you are after? Is it simply "the first time SWH drops below 0.5 or 0", respectively? Because then, a simple way is to just use the find function.
Example:
xxx=1:1:419; %your x values
yyy=acf(:,1); %For the sake of example I just take the SWH values from your first plot.
index1=find(yyy < 0.5); %This will give a vector of all the indexes of yyy, where the value is smaller than 0.5
%In this example, index1(1) is 20, which means yyy(20) is the first value of yyy below 0.5. The corresponding x value is simply xxx(20).
%You can also just type yyy(index1(1)) and xxx(index1(1)) to get the values out.
If you need "more precision", ie. you need to know the actual intersect point and not just the step, then combine the point you find above with the point immediately before that, make a simple linear fit with those two x,y pairs, and find the values that satisfy your condition, but that takes a few lines of coding extra.
0 个评论
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!