%Demonstracja wyznaczania prostych aproksymujacych
%w zadaniu klasyfikacji
%Metoda analityczna i numeryczna - gradientowa
osie=[-10 10 -20 20];
u1=[-4 0];
u2=[5 2];
u3=[2 -5];
u1=[1 u1];
u2=[1 u2];
u3=[1 u3];
U=[u1;u2;u3];
z=[1;1;-1];
disp('Rozwiązanie analityczne')
%det(U'*U)
w=inv(U'*U)*U'*z
wp=w;
r=U*w-z;
disp('Błąd aproksymacji')
r=r'*r
figure(1);clf;
subplot(2,2,1);
newplot;axis(osie);
hold on;
N=length(z);
for licz=1:N
if z(licz)==1
plot(U(licz,2),U(licz,3),'bo');
elseif z(licz)==-1
plot(U(licz,2),U(licz,3),'gs');
else
plot(U(licz,2),U(licz,3),'kx');
end;
end;
pause
kolor=['b-'];
wrys=[w(1)-1;w(2:3)];
rprosta(wrys,osie,kolor);
kolor=['g-'];
wrys=[w(1)+1;w(2:3)];
rprosta(wrys,osie,kolor);
kolor=['m-'];
wrys=w;
rprosta(wrys,osie,kolor);
title('Rozwiazanie analityczne');
hold off;
pause
subplot(2,2,2);
newplot;
hold on;
plot(U(1:2,2),U(1:2,3),'bo');
plot(U(3,2),U(3,3),'gs');
wrys=w;
rprosta(wrys,osie,'m-');
rr=[];
disp('Rozwiązanie metodą gradientową')
disp('Obserwuj proces obliczeń na wykresach, liczba kroków')
NN=20
%wartość początkowa wektora wag
w=[1;-1;-1];
w=10*[rand-0.5;rand-0.5;rand-0.5]
weta=0.02;
osi=[1 NN -100 20];
for k=1:NN;
k
deltaJ=U'*(U*w-z);
w=w-weta*deltaJ
r=U*w-z;
r=r'*r;
rr=[rr 20*log10(r)];
subplot(2,2,2);
kolor=['r-'];
wrys=w;
rprosta(wrys,osie,'r-');
title('Historia');
hold on;
subplot(2,2,3);
newplot;
hold on;
plot(U(1:2,2),U(1:2,3),'bo');
plot(U(3,2),U(3,3),'gs');
wrys=w;
kolor=['r-'];
wrys=w;
rprosta(wrys,osie,'r-');
title('Rozwiazanie aktualne');
hold off;
grid on;
subplot(2,2,4);
plot(1:k,rr);
title('Błąd w [dB]');
axis(osi);
grid on;
disp('Pauza - kliknij')
pause
end;
subplot(2,2,2);
hold off;
subplot(2,2,4);
hold off;
disp('Porównanie wyników: analitycznego i numerycznego')
[wp w]
Wyszukiwarka
Podobne podstrony:
asystent osoby niepelnosprawnej46[02] z1 02 u311[15] Z1 01 Wykonywanie pomiarów warsztatowychfotograf13[05] z1 10 uwzory protokołów pomiarowych zap1102012 z1fotograf13[05] z1 01 nlakiernikq4[03] z1 02 nmonter systemow rurociagowychq3[04] z1 01 nfototechnik13[01] z1 01 nTechnik?zpieczenstwa i higieny pracy15[01] Z1 01 nmonter instrumentow muzycznychs1[02] z1 01 uasystent osoby niepelnosprawnej46[02] z1 01 u (2)monter sieci komunalnychq3[03] z1 2 3 04 uJak działa mózg neuronu –z1 podejmowanie?cyzji w warunkach ryzykaneuronalne korelaty obecnosciTech tech chem11[31] Z1 01 uwięcej podobnych podstron