oard end fprintf ( Ітерація:% g з% g t кращий боковий пік АКФ:% g t ВКФ: % g n laquo ;, iter, iterations, best_lm_AKF, best_lm_VKF) keyboard end% fprintf ( Ітерація:% g з% g t кращий боковий пік АКФ:% g t ВКФ:% g n laquo ;, iter, iterations , best_lm_AKF, best_lm_VKF)% keyboard endІщем зі знайдених кращих ВКФ кращий перший бічний пік АКФbest_lm_AKF_first=1; iter_int=1; for h=1: length (bests_kk_VKF) z=[q (1,1). * exp (i * phi (kk1,1)) q (2,1). * exp (i * phi (kk2,1)) q (3,1). * exp (i * phi (kk3,1)) q (4,1). * exp (i * phi (kk4,1)) q (5,1). * exp (i * phi (kk5,1)) q (6,1). * exp (i * phi (kk6,1)) q (7,1). * exp (i * phi (kk7,1)) q (8,1 ). * exp (i * phi (kk8,1)) q (9,1). * exp (i * phi (kk9,1)) q (10,1). * exp (i * phi (kk10,1 )) q (11,1). * exp (i * phi (kk11,1)) q (12,1). * exp (i * phi (kk12,1)) q (13,1). * exp ( i * phi (kk13,1)) q (14,1). * exp (i * phi (kk14,1)) q (15,1). * exp (i * phi (kk15,1)) q (16 , 1). * exp (i * phi (kk16,1)) q (17,1). * exp (i * phi (kk17,1)) q (18,1). * exp (i * phi (kk18 , 1)) q (19,1). * exp (i * phi (kk19,1)) q (20,1). * exp (i * phi (kk20,1)) q (21,1). * exp (i * phi (kk21,1)) q (22,1). * exp (i * phi (kk22,1)) q (23,1). * exp (i * phi (kk23,1)) q (24,1). * exp (i * phi (kk24,1)) q (25,1). * exp (i * phi (kk25,1)) q (26,1). * exp (i * phi (kk26,1)) q (27,1). * exp (i * phi (kk27,1)) q (28,1). * exp (i * phi (kk28,1)) q (29,1). * exp (i * phi (kk29,1)) q (30,1). * exp (i * phi (kk30,1)) q (31,1). * exp (i * phi (kk31,1) )] raquo ;; ERMIT=diag (z, 1) + diag (conj (z), - 1); [V, D]=eig (ERMIT); for k=1: n [AKF (1: 2 * n - 1, k), lagsAKF]=xcorr (V (1: n, k), coeff ); еnd% AKF (n, 1: n)=0; % для пошуку бічного максимуму обнулив n-й рядок max_AKF=max (abs (AKF (n - 1,1: n))); maximums_AKF_first (h, 1)=max_AKF; if max_AKF lt; best_lm_AKF_first best_lm_AKF_first=max_AKF; best_kk_AKF_first=[kk1; kk2; kk3; kk4; kk5; kk6; kk7; kk8; kk9; kk10; kk11; kk12; kk13; kk14; kk15; kk16; kk17; kk18; kk19; kk20; kk21; kk22; kk23; kk24; kk25; kk26; kk27; kk28; kk29; kk30; kk31]; end if max_AKF lt; =. 2 501 bests_lm_AKF_first (iter_int, 1)=max_AKF; bests_kk_AKF_first (iter_int, 1)=kk1; bests_kk_AKF_first (iter_int, 2)=kk2; bests_kk_AKF_first (iter_int, 3)=kk3; bests_kk_AKF_first (iter_int, 4)=kk4; bests_kk_AKF_first (iter_int, 5)=kk5; bests_kk_AKF_first (iter_int, 6)=kk6; bests_kk_AKF_first (iter_int, 7)=kk7; bests_kk_AKF_first (iter_int, 8)=kk8; bests_kk_AKF_first (iter_int, 9)=kk9; bests_kk_AKF_first (iter_int, 10)=kk10; bests_kk_AKF_first (iter_int, 11)=kk11; bests_kk_AKF_first (iter_int, 12)=kk12; bests_kk_AKF_first (iter_int, 13)=kk13; bests_kk_AKF_first (iter_int, 14)=kk14; bests_kk_AKF_first (iter_int, 15)=kk15; bests_kk_AKF_first (iter_int, 16)=kk16; bests_kk_AKF_first (iter_int, 17)=kk17; bests_kk_AKF_first (iter_int, 18)=kk18; bests_kk_AKF_first (iter_int, 19)=kk19; bests_kk_AKF_first (iter_int, 20)=kk20; bests_kk_AKF_first (iter_int, 21)=kk21; bests_kk_AKF_first (iter_int, 22)=kk22; bests_kk_AKF_first (iter_int, 23)=kk23; bests_kk_AKF_first (iter_int, 24)=kk24; bests_kk_AKF_first (iter_int, 25)=kk25; bests_kk_AKF_first (iter_int, 26)=kk26; bests_kk_AKF_first (iter_int, 27)=kk27; bests_kk_AKF_first (iter_int, 28)=kk28; bests_kk_AKF_first (iter_int, 29)=kk29; bests_kk_AKF_first (iter_int, 30)=kk30; bests_kk_AKF_first (iter_int, 31)=kk31; iter_int=iter_int + 1; end end for h=1: length (bests_kk_VKF) bests_phi_radian (h, 1: n - 1)=[phi (bests_kk_VKF (h, 1), 1) phi (bests_kk_VKF (h, 2), 1) phi (bests_kk_VKF (h , 3), 1)]; end bests_phi_degree=bests_phi_radian. * 180/pi;
