tfpm_simu_defpre_tfma.m

来自「用于模拟时变非平稳的ARMA过程」· M 代码 · 共 11 行

SPRE= [
sprintf('                                                                                                                                                                                                                                                                                                 '); 
sprintf('[AAlong, BBlong]= tfar_est_tfywu(Ayy(N/2-1*3*LMA+1:N/2+1*3*LMA+1, N/2-2*MMA+1:N/2+2*MMA+1).*Psilong(N/2-1*3*LMA+1:N/2+1*3*LMA+1, N/2-2*MMA+1:N/2+2*MMA+1), N);e= tfarma_inv(y, AAlong(:, :, end), BBlong(:, :, end));Rey=corr_est(e, y, MMA, 1/2);Aey=fft(Rey);Aey=[Aey(N/2+1:N,:);Aey(1:N/2,:)];'); 
sprintf('                                                                                                                                                                                                                                                                                                 '); 
sprintf('                                                                                                                                                                                                                                                                                                 '); 
sprintf('[AAlong, BBlong]= tfar_est_tfywu(Ayy(N/2-2*3*LMA+1:N/2+2*3*LMA+1, N/2-2*MMA+1:N/2+2*MMA+1).*Psilong(N/2-2*3*LMA+1:N/2+2*3*LMA+1, N/2-2*MMA+1:N/2+2*MMA+1), N);e= tfarma_inv(y, AAlong(:, :, end), BBlong(:, :, end));Rey=corr_est(e, y, MMA, 1/2);Aey=fft(Rey);Aey=[Aey(N/2+1:N,:);Aey(1:N/2,:)];'); 
sprintf('                                                                                                                                                                                                                                                                                                 '); 
      
];