0001 function [PSD,operador,coef_ar] = ajusta_bloque_ar12(observado) 0002 0003 M = size(observado,1); 0004 Y=[]; 0005 rhs=[]; 0006 0007 bloque = observado; 0008 0009 for jj=3:M-2 0010 for kk=3:M-2 0011 Y=[Y; bloque(jj,kk-1) bloque(jj,kk-2) bloque(jj-1,kk+2:-1:kk-2) bloque(jj-2,kk+2:-1:kk-2)]; 0012 rhs =[rhs;bloque(jj,kk)]; 0013 end 0014 end 0015 0016 coef_ar = Y\rhs; 0017 0018 0019 ptos=M; 0020 0021 w=linspace(-pi,pi,ptos); 0022 0023 z1=exp(sqrt(-1)*w); 0024 z2=exp(sqrt(-1)*w); 0025 0026 0027 F = zeros(ptos); 0028 0029 0030 for jj=1:ptos 0031 for kk=1:ptos 0032 v1 = coef_ar(1)*z1(jj)^0*z2(kk)^(-1); 0033 v2 = coef_ar(2)*z1(jj)^0*z2(kk)^(-2); 0034 v3 = coef_ar(3)*z1(jj)^(-1)*z2(kk)^2; 0035 v4 = coef_ar(4)*z1(jj)^(-1)*z2(kk)^1; 0036 v5 = coef_ar(5)*z1(jj)^(-1)*z2(kk)^(0); 0037 v6 = coef_ar(6)*z1(jj)^(-1)*z2(kk)^(-1); 0038 v7 = coef_ar(7)*z1(jj)^(-1)*z2(kk)^(-2); 0039 v8 = coef_ar(8)*z1(jj)^(-2)*z2(kk)^2; 0040 v9 = coef_ar(9)*z1(jj)^(-2)*z2(kk)^1; 0041 v10 = coef_ar(10)*z1(jj)^(-2)*z2(kk)^(0); 0042 v11 = coef_ar(11)*z1(jj)^(-2)*z2(kk)^(-1); 0043 v12 = coef_ar(12)*z1(jj)^(-2)*z2(kk)^(-2); 0044 0045 val = 1 - (v1 + v2 + v3 + v4 + v5 + v6 + v7 + v8 +v9 + v10 + v11 + v12); 0046 F(jj,kk) = val; 0047 PSD(jj,kk) = 1 / abs(val)^2; 0048 end 0049 end 0050 0051 operador = F;