Identificarea Semnalelor

%criteriul aplatizarii

function [Vmin]=Aplatizare(p,y)

N=length(y);

M=N;

%for p=1:N

%Se calculeaza r

r=calc_r(p,y);

% Se calculeaza R

R=calculez_R(p,y);

% Se calculeaza B

B=zeros(p+1,p+1);

for i=1:1:p+1

B(i,i)=1/(sqrt(M)*M^(i-1));

end

theta_T=B*inv(B*R*B)*B*r';

% Se calculeaza V

%V(p)=0;

for n=1:N

yT(n)=0;

V(n)=0;

for i=1:p+1

yT(n)=yT(n)+theta_T(i)*(n^(i-1));

end

Ystat(n)=y(n)-yT(n);

V(n)=V(n)+Ystat(n)^2;

V(n)=V(n)*N/(N-p);

end

%V(n)=V(n)*N/(N-p);

%end

figure(3);

plot(1:N,V);

title('Criteriul aplatizarii');

axis([0 26 0 20]);

Vmin=min(V);

for p=1:N

if (Vmin==V(p)) break;

end

end

4 teme la identificarea semnalelor IS