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