pause;
subplot(2,1,1); plot(wx/pi,Wx/max(Wx)); grid; title('PSD of the input time-series');xlabel('unormowana czestotliwosc');ylabel('wartosc');
subplot(2,1,2); plot(wrho/pi,Wrho4/max(Wrho4)); grid; title('Amplitude spectrum of the time-varying Schur coefficient');xlabel('unormowana czestotliwosc');ylabel('wartosc');


hold on;
for m = 2:nrwsp
    [Wrho4,wrho] = wgmx(rho(:,m));
    plot(wrho/pi,Wrho4/max(Wrho4)); grid on; title('Amplitude spectrum of the time-varying Schur coefficient');xlabel('unormowana czestotliwosc');ylabel('wartosc');
end

Legend = cell(nrwsp,1);
for i=1:nrwsp
    Legend{i} = strcat('nrwsp = ',num2str(i));
end
legend(Legend)

hold off;
pause;