dan@0: function [dddlss ddlss dlss A_der2_plus A_der_plus A_plus] = ... dan@0: ddm_der2_lin_sys_sol(A_sys,b_sys,A_der_sys,b_der_sys,A_der2_sys,b_der2_sys,N_fft,R,Q) dan@0: dlss = zeros(Q+1, N_fft-R+1); dan@0: ddlss = zeros(Q+1, N_fft-R+1); dan@0: dddlss = zeros(Q+1, N_fft-R+1); dan@0: dan@0: for l=1:N_fft-R+1 dan@0: A = A_sys(:,:,l); dan@0: b = b_sys(:,:,l); dan@0: dan@0: AA_inv = inv(A'*A); dan@0: A_plus = AA_inv*A'; dan@0: ddm_est_ = fliplr((A_plus * b).'); dan@0: dlss(:,l) = [ddm_est_ 0]; dan@0: dan@0: A_der = A_der_sys(:,:,l); dan@0: b_der = b_der_sys(:,:,l); dan@0: AA_mix = ( A_der' * A + A' * A_der ); dan@0: A_der_plus = (-AA_inv * AA_mix * AA_inv * A' + AA_inv * A_der'); dan@0: ddm_est_der_ = fliplr((A_der_plus *b + A_plus * b_der).'); dan@0: ddlss(:,l) = [ddm_est_der_ 0]; dan@0: dan@0: A_der2 = A_der2_sys(:,:,l); dan@0: b_der2 = b_der2_sys(:,:,l); dan@0: dan@0: AA_inv_der = - AA_inv * AA_mix * AA_inv; dan@0: AA_mix_der = A_der2'*A + 2*A_der'*A_der + A'*A_der2; dan@0: A_der2_plus = ... dan@0: - AA_inv_der * AA_mix * AA_inv * A'... dan@0: - AA_inv * AA_mix_der * AA_inv * A'... dan@0: - AA_inv * AA_mix * AA_inv_der * A'... dan@0: - AA_inv * AA_mix * AA_inv * A_der'... dan@0: + AA_inv_der * A_der' + AA_inv * A_der2'; dan@0: ddm_est_der2_ = fliplr((A_der2_plus *b +2 * A_der_plus * b_der + A_plus*b_der2).'); dan@0: dddlss(:,l) = [ddm_est_der2_ 0]; dan@0: dan@0: end dan@0: end