%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]