Home > PR > ajusta_bloque_ar4.m

ajusta_bloque_ar4

PURPOSE ^

SYNOPSIS ^

function [PSD,operador,coef_ar] = ajusta_bloque_ar4(observado)

DESCRIPTION ^

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

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;

Generated on Wed 29-Nov-2006 16:19:19 by m2html © 2003