Home > PR > ajusta_bloque_ar12.m

ajusta_bloque_ar12

PURPOSE ^

SYNOPSIS ^

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

DESCRIPTION ^

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

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;

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