z2 learnWH cykl.M, WAT, SEMESTR VIII, Metody i zadania rozpoznawania wzorców
[ Pobierz całość w formacie PDF ]
%Aproksymacja liniowa%Demonstracja wyznaczania prostych aproksymujacych%w zadaniu klasyfikacji%Algorytm Widrowa-Hoffaosie=[-2 2 -2 2];L=4;%L-liczba punktówu1=[0 0];u2=[0 1];u3=[1 0];u4=[1 1];u1=[1 u1];u2=[1 u2];u3=[1 u3];u4=[1 u4];U=[u1;u2;u3;u4];det(U'*U);z=[1;1;-1;-1];w=inv(U'*U)*U'*z;wp=w;r=U*w-z;r=r'*r;figure(1);clf;subplot(2,2,1);newplot;hold on;N=length(z);for licz=1:Nif z(licz)==1plot(U(licz,2),U(licz,3),'bo');elseif z(licz)==-1plot(U(licz,2),U(licz,3),'gs');elseplot(U(licz,2),U(licz,3),'kx');end;end;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-'];rprosta(wp,osie,kolor);hold off;rr=[];NN=50;w=[0;0;1];w=10*[rand-0.5;rand-0.5;rand-0.5]weta=0.5;osi=[1 NN -150 20];subplot(2,2,2);newplot;rprosta(wp,osie,kolor);hold on;N=length(z);for licz=1:Nif z(licz)==1plot(U(licz,2),U(licz,3),'bo');elseif z(licz)==-1plot(U(licz,2),U(licz,3),'gs');elseplot(U(licz,2),U(licz,3),'kx');end;end;rprosta(wp,osie,kolor);for k=1:NN;k%cykliczny wybór punktu ze zbioru uczącego');nr=mod(k-1,L)+1;x=U(nr,:);x=x';klasa=z(nr);delta=w'*x-z(nr);w=w-weta*delta*x;r=delta'*delta';bw=wp-w;bw=bw'*bw;wiersz=[20*log10(r) 20*log10(bw)];rr=[rr;wiersz];subplot(2,2,2);hold on;rprosta(w,osie,'r-');hold on;subplot(2,2,3);newplot;hold on;N=length(z);for licz=1:Nif z(licz)==1plot(U(licz,2),U(licz,3),'bo');elseif z(licz)==-1plot(U(licz,2),U(licz,3),'gs');elseplot(U(licz,2),U(licz,3),'kx');end;end;if klasa==1plot(x(2),x(3),'ro');elseplot(x(2),x(3),'rs');endrprosta(w,osie,'r-');wrys=[w(1)-1;w(2:3)];rprosta(wrys,osie,'b-');wrys=[w(1)+1;w(2:3)];rprosta(wrys,osie,'g-');hold off;grid on;subplot(2,2,4);plot(1:k,rr(:,1),'b-',1:k,rr(:,2),'g-');axis(osi);grid on;pauseend;subplot(2,2,2);hold off;subplot(2,2,4);hold off;disp('Porównanie wyniku analitycznego i numerycznego');[wp w]
[ Pobierz całość w formacie PDF ]