%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% AGH - EAIiE - M.I.S. III rok %
% cw_1 %
% autor: Mariusz Dyba�a 9.04.2010 %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all;
disp('Czy chcesz kreslic charakterystyki b0,b1 od K czy od N?');
disp('Wpisz 1 jesli N, lub 0 jesli K');
wybor=input('Wybor: ');
switch wybor
case 1
beta0 = input('podaj parametr beta 0: ');
beta1 = input('podaj parametr beta 1: ');
disp(' zakres argumentow x');
x1 =input('dolna granica: ');
x2 = input('gorna granica: ');
x3 = input('krok: ');
x = x1:x3:x2;
K1=5;
x4= x2-x1;
y = 2+3*x+K1*randn(1,length(x));
M = length(x);
N=1;
sumaxy = x*(y');
sumax = sum(x);
sumay = sum(y);
sumaxkwadrat = x*(x');
N = x1;
figure(1)
for k=1:x4
b1 = ((N*sumaxy) - (sumax*sumay))/((N*sumaxkwadrat) - (sumax^2));
b0 = (1/N)*(sumay -( b1*sumax));
blad1 = (b1 - beta1)/beta1;
blad0 = (b0 - beta0)/beta0;
y1 = (b1*x) + b0;
plot(N,b0,'r*',N,b1,'b*',N,blad1,'ro',N,blad0,'go')
legend('b0', 'b1', 'bladb1','bladb2')
title('zmiana wspolczynnikow b0 i b1 i ich bledow bezwzglednych w zaleznosci od N, przy K=5')
hold on
N=N+x3;
end
figure(3)
plot(x,y,x,y1)
title('sprawdzenie poprawnosci algorytmu');
case 0
beta0 = input('podaj parametr beta 0: ');
beta1 = input('podaj parametr beta 1: ');
disp (' podaj zakres wspolczynnika K: ');
k1 =input('dolna granica: ');
k2 = input('gorna granica: ');
k3 = input('krok: ');
k= k1:k3:k2;
x=1:1:50;
M = length(x);
sumax = sum(x);
sumaxkwadrat = x*(x');
for l=k1:k3:k2;
y = 2+3*x+l*randn(1,length(x));
sumaxy = x*(y');
sumay = sum(y);
b1(l) = ((M*sumaxy) - (sumax*sumay))/((M*sumaxkwadrat) - (sumax^2));
b0(l) = (1/M)*(sumay -( b1(l)*sumax));
%blad1 = (b1- beta1)/beta1;
%blad0 = (b0- beta0)/beta0;
end
figure (3)
plot(b1,'rx')
title('zmiana wspolczynnika b1 w zaleznosci od K')
legend('b1')
figure(4)
plot(b0,'bx')
title('zmiana wspolczynniks b0 w zaleznosci od K')
legend('b0')
otherwise
disp('zly wybor');
end