townini88888

PIERWSZE:

close all
clear all
%podpunkt a
X=magic(4);
%dekompozycja QR - start
[Q,R]=qr(X);
for i=1:10
    [Q,R]=qr(R*Q);
end
lambda=diag(R*Q)
%dekompozycja QR - koniec
%funkcja matlabowska eig
eig(X)

%podpunkt b
A=X;
x=[1; 1; 1; 1];
for i=1:10
    y=A*x;
    [ym,i]=max(abs(y));
    x=y/y(i);
    l=y(i);
end
eig(A)

..--..--..--..--..--..-..--..--..--..--..---.-.-.-.--.-.-.--.-.--..-.--.....................................

DRUGIE:

%Wczytywanie plików
mix1 = wavread('mix1.wav');
mix2 = wavread('mix2.wav');
%Budowa macierzy kowariancji
kow = cov(mix1, mix2);

%Obliczanie wektorów własnych macierzy
[V, D]=eig(kow);
aaa=[mix1 mix2]*V;
plot(aaa)

----------------------------------------------------------------------------------------------

TRZECIE:

xx=double(X);
[U,S,V]= svd(xx);
k=50;
s=zeros(k,k);
for n=1:1:k
    s(n,n)=S(n,n);
end
u=zeros(512,k);
v=zeros(512,k);
for n=1:1:k
    for i=1:1:512
        u(i,n)=U(i,n);
        v(i,n)=V(i,n);
    end
end
x=u*s*v';
mse=0;
for i=1:512
    for j=1:512
        mse=mse+(xx(i,j)-x(i,j))^2;
    end
end
mse=mse/(512^2);
psnr=10*log10((255^2)/mse);