clc
clear all
close all
w=120*pi;
fi=pi/6;
ft=1440;
T=1/ft;
a=w*T;
t=0:T:0.02;
f1=0.25*64*(0.5-(rand(size(t)))); %zakłócenie sygnału
f=64*sin(w*t+fi)+f1; %sygnał z zakłóceniem
%Metoda aproksymacji z użyciem modelu trygonometrycznego
for k=1:length(t)
y(:,k)=[sin(a*k);cos(a*k)];
end;
H=inv(y*y')*y*f';
aprox=H'*y;
figure(1);
plot(t,f,'--k');
xlabel('Czas'); ylabel('Amplituda');
hold on
grid on
plot(t,aprox,'k');
legend('sygnal z zakloceniem','aproksymacja ');
%Metoda aproksymacji wspomaganej SVD
[U,W,V]=svd(y');
x=inv(V*W'*W*V')*(V*W'*U')*f';
aproxSVD=x'*y;
figure(2);
plot(t,aproxSVD,'k');
xlabel('Czas'); ylabel('Amplituda');
hold on
grid on
plot(t,f,'--k');
legend('aproksymacja z SVD','sygnal z zakloceniem');
err=max(aprox)-max(f)%błąd aproksymacji z użyciem modelu trygonometrycznrego
err2=max(aproxSVD)-max(f)%błąd aproksymacji wspomaganej SVD