Home > PR > ajusta_bloque_ar8.m

ajusta_bloque_ar8

PURPOSE ^

SYNOPSIS ^

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

DESCRIPTION ^

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

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

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