MRESPA Get measurement response from an averaging kernel matrix FORMAT [mresp,As] = mrespA( A, ji ) OUT mresp Measurement response. As Array of matrices. One matrix for each retrieval quantity. As returned by *splitA*. IN A Complete averaging kernel matrix. ji Matrix of two columns, where the columns hold start and end index for each retrieval quantity.
0001 % MRESPA Get measurement response from an averaging kernel matrix 0002 % 0003 % FORMAT [mresp,As] = mrespA( A, ji ) 0004 % 0005 % OUT mresp Measurement response. 0006 % As Array of matrices. One matrix for each retrieval quantity. 0007 % As returned by *splitA*. 0008 % IN A Complete averaging kernel matrix. 0009 % ji Matrix of two columns, where the columns hold start and end 0010 % index for each retrieval quantity. 0011 0012 % 2006-09-29 Created by Patrick Eriksson. 0013 0014 0015 function [mresp,As] = mrespA(A,ji) 0016 0017 % Checks are done in *splitA* 0018 0019 As = splitA( A, ji ); 0020 0021 mresp = zeros( size(A,1), 1 ); 0022 0023 i0 = 0; 0024 0025 for i = 1:length(As) 0026 ind = i0+1:i0+size(As{i},1); 0027 mresp(ind) = sum( A(ind,ind)' ); 0028 i0 = ind(end); 0029 end 0030 0031 0032