Info

This question is closed. Reopen it to edit or answer.

Indexing matrix to get corresponding values for condition,store in a new matrix each time?

1 view (last 30 days)
Sorry for the perhaps confusing title...
Basically I have a 3x3 matrix containing elevation angle, azimuth angle and range. I want to generate new matrices each time elevation >5 deg. There are usually about 5 segments that have this data and I want to separate each one into a new matrix.
I know how to index but not sure how to put this condition in...
Thanks
sat_tcs=llh2tcsT(sat_llh,station_llh);
sat_elev=atan2(sat_tcs(3,:),sqrt(sat_tcs(1,:).^2+sat_tcs(2,:).^2)); sat_azim=atan2(-sat_tcs(2,:),sat_tcs(1,:));
range=sqrt(sat_tcs(1,:).^2+sat_tcs(2,:).^2+sat_tcs(3,:).^2);` sat_elev(sat_elev < 5*deg2rad) = NaN;
sat_look_tcs=[sat_elev;sat_azim;range];

Answers (1)

David Sanchez
David Sanchez on 22 Aug 2014
It is not very clear what you want, but I think it has to be something like this:
sat_look_tcs=[sat_elev;sat_azim;range];
elev_5 = find(sat_elev >5);
sat_look_tcs_5 = zeros(length(elev_5, 3, 3)); % initialize 3D matrix to hold your submatrices
for k=1:(length(elev_5)-1)
sat_look_tcs_5(k+1,:,:) = sat_look_tcs((elev_5(k)+1):elev_5(k+1),:,:);
end
And you access your matrix with
sat_look_tcs_5(1,:,:)
sat_look_tcs_5(2,:,:)
etc...
  3 Comments
Michael Haderlein
Michael Haderlein on 22 Aug 2014
There's a typo. It should read
sat_look_tcs_5 = zeros(length(elev_5),3,3);
I didn't check the rest of the code as I also didn't really understand the question. But maybe that's already what you need.
Buzz
Buzz on 8 Sep 2014
Hi, yeah thanks. The problem I have is that I am getting a 3x3 with the elevation angle in each column, rather than the corresponding azimuth and range for each value.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!