0001 function [PSD,operador,coef_ar] = ajusta_bloque_ar4(observado) 0002 0003 M = size(observado,1); 0004 0005 Y=[]; 0006 rhs=[]; 0007 0008 bloque = observado; 0009 0010 for jj=2:M-1 0011 for kk=2:M-1 0012 Y=[Y; bloque(jj+1,kk) bloque(jj,kk+1) bloque(jj-1,kk) bloque(jj,kk-1)]; 0013 rhs =[rhs;bloque(jj,kk)]; 0014 end 0015 end 0016 0017 coef_ar = Y\rhs; 0018 0019 0020 ptos=M; 0021 0022 w=linspace(-pi,pi,ptos); 0023 0024 z1=exp(sqrt(-1)*w); 0025 z2=exp(sqrt(-1)*w); 0026 0027 0028 F = zeros(ptos); 0029 espectro = F; 0030 0031 for jj=1:ptos 0032 for kk=1:ptos 0033 v1 = coef_ar(1)*z1(jj)^1; 0034 v2 = coef_ar(2)*z2(kk)^1; 0035 v3 = coef_ar(3)*z1(jj)^(-1); 0036 v4 = coef_ar(4)*z2(kk)^(-1); 0037 0038 val = 1 - (v1 + v2 + v3 + v4); 0039 F(jj,kk) = val; 0040 PSD(jj,kk) = 1 / abs(val)^2; 0041 end 0042 end 0043 0044 operador = F;