Notes

function La13Nastula_urn
t0=clock; % poczatek obliczen
rok=int2str(t0(1)); miesiac=int2str(t0(2)); dzien=int2str(t0(3));
godz=int2str(t0(4));min=int2str(t0(5)); sek=int2str(t0(6));
% tworzenie nazwy pliku ze znacznikiem czasu
nazwa=['La13Nastula_urnOUT',dzien,'-',miesiac,'-',rok,'_',godz,'h',min,'min.m', sek, 'sek'];
[nazwa,sciezka] =...
uiputfile(nazwa,' UWAGA! Zakceptuj nazwe ... albo wpisz nowa');
fnames=strcat(sciezka,nazwa);
[fd,msg]=fopen(fnames,'wt');
fprintf(fd,'n Otwarto plik tekstowy o nazwie ze znacznikiem czasowym n %s',...
fnames);
fprintf(fd,'n Data: %5d-%2d-%2d godz. %2d, %2dmin, %2.0fs',clock);
fprintf( 'n Otwarto plik tekstowy o nazwie ze znacznikiem czasowym n %s',...
fnames);
fprintf( 'n Data: %5d-%2d-%2d godz. %2d, %2dmin, %2.0fs',clock);
fprintf(fd,'n Uklad rownan zdefiniowny w funkcji f=urn7par(x,p1,p2,p3,p4,p5,p6,p7)');
fprintf(fd,'n f(1)=a* x(1)^2 + b*x(2)^2 +c*x(3)^2-4;');
fprintf(fd,'n f(2)=d*x(1)*x(2)-1;');
fprintf(fd,'n f(3)=e*x(1)*x(2)+f*x(1)*x(3)+g*x(2)*x(3)-2;');
fprintf(fd,'n zostanie rozwiazany za pomoca fsolve z uzyciem wskaznika do funkcji');
fprintf(fd,'n [x,fval,exitflag,output]=fsolve(@(x) urn7par(x,a,b,c,d,e,f,g),x0,options);');
%
wdold=cd;
[fname,sciezka]=uigetfile('La13Nastula_urnDAT.m','Wybierz z m-pliku ...n');
eval(['cd(''',sciezka,''')']);
datafile=strtok(fname,'.');
fprintf('n ... CZYTANIE danych tekstowych ...');
fprintf(fd, 'n ... wybrano sciezke dostepu: %s',sciezka);
fprintf( 'n ... wybrano sciezke dostepu: %s',sciezka);
fprintf(fd, 'n ... nazwa wybranego m-pliku: %s.m',datafile);
fprintf( 'n ... nazwa wybranego m-pliku: %s.m',datafile);
[imie,nazwisko,album,a]=feval(datafile); % czytanie danych z wybranego m-pliku
function f = La13Nastula_urnf(x,par) %uklad rownan nieliniowych
% uklad rownan
a=par(1); b=par(2); c=par(3); d=par(4);
f(1)=a*x(1).^2+ a*x(2).^2+ a*x(3).^2-b;
f(2)=x(1)*x(2)-c;
f(3)=x(1)*x(2)+x(1)*x(3)+x(2)*x(3)-d;
end
%wyznaczenie wartości a
k=length(album);
a=length(imie)+length(nazwisko);
while(k>0)
a=a+str2double(album(k));
k=k-1;
end
b=4*a;
c=a/b;
d=c/4;
e=1;
f=1;
g=1;
%par=[a b c d e f g];
options=optimset('Display','on','LargeScale','off','MaxIter',100,'TolCon',1e-4,'TolFun',1e-4,'TolX',1e-4);
%
fprintf(fd,'nn options=optimset(''Display'',''off'',''LargeScale'',''off'',''MaxIter'',100,''TolCon'',1e-4,''TolFun'',1e-4,''TolX'',1e-4);');

% 1-szy p. startowy
x0=[1 22 333]; x=x0; % punkt startowy obliczen
fprintf('nn 1-szy punkt startowy obliczen: x0 = [%.2f %.2f %.2f]',x0(1),x0(2),x0(3));
[x,fval,exitflag,output]=fsolve(@(x) urn7par(x,a,b,c,d,e,f,g),x0,options);
opisrozw(fd,x0,x,fval,exitflag,output);
%
%2-gi punkt startowy
x0=[-1 -22 -333]; x=x0; % punkt startowy obliczen
fprintf('nnn 2-gi punkt startowy obliczen: x0 = [%.2f %.2f %.2f]',x0(1),x0(2),x0(3));
[x,fval,exitflag,output]=fsolve(@(x) urn7par(x,a,b,c,d,e,f,g),x0,options);
opisrozw(fd,x0,x,fval,exitflag,output);

%3-ci punkt startowy
x0=[1 2 -3]; x=x0; % punkt startowy obliczen
fprintf('nnn 2-gi punkt startowy obliczen: x0 = [%.2f %.2f %.2f]',x0(1),x0(2),x0(3));
[x,fval,exitflag,output]=fsolve(@(x) urn7par(x,a,b,c,d,e,f,g),x0,options);
opisrozw(fd,x0,x,fval,exitflag,output);

%4-ty punkt startowy
x0=[1.1 0.9 1.0]; x=x0; % punkt startowy obliczen
fprintf('nnn 2-gi punkt startowy obliczen: x0 = [%.2f %.2f %.2f]',x0(1),x0(2),x0(3));
[x,fval,exitflag,output]=fsolve(@(x) urn7par(x,a,b,c,d,e,f,g),x0,options);
opisrozw(fd,x0,x,fval,exitflag,output);
%
et=etime(clock,t0); % czas obliczen
fprintf(fd,'n Czas obliczen programem urnsolve tobl = %.2f sekund !', et);
fprintf( 'n Czas obliczen programem urnsolve tobl = %.2f sekund !', et);
fclose('all');
fprintf('nn ***Wyniki w pliku La13Nastula_urnOUT.m *** nn');



function f = urn7par(x,p1,p2,p3,p4,p5,p6,p7)
% uklad rownan
%p1=par(1); p2=par(2); p3=par(3); p4=par(4);
%p5=par(5); p6=par(6); p7=par(7);
f(1)= p1*x(1)^2 + p2*x(2)^2 +p3*x(3)^2-4;
f(2)=p4*x(1)*x(2)-1;
f(3)=p5*x(1)*x(2)+p6*x(1)*x(3)+p7*x(2)*x(3)-2;
end




function opisrozw(fd,x0,x,fval,exitflag,output);
% opis rozwiazania
fprintf(fd,'Michal Tyranski, album=249918, a=47');
fprintf(fd,'nn O P I S R O Z W I A Z A N I A n');
fprintf( 'nn O P I S R O Z W I A Z A N I A n');
fprintf(fd,'n zmienna PunktStartowy Rozwiazanie');
fprintf( 'n zmienna PunktStartowy Rozwiazanie');
n=length(x);
for i=1:n
fprintf(fd,'n %2d %8.4f %8.4f',i,x0(i),x(i));
fprintf( 'n %2d %8.4f %8.4f',i,x0(i),x(i));
end
fprintf(fd,'n fval = %8.4e',fval);
fprintf( 'n fval = %8.4e',fval);
fprintf(fd,'n exitflag = %d',exitflag);
fprintf( 'n exitflag = %d',exitflag);
switch exitflag
case -4
fprintf(fd,'n -4 - nie mozna zmnieszyc kroku w kierunku poszukiwan');
fprintf( 'n -4 - nie mozna zmnieszyc kroku w kierunku poszukiwan');
case -3
fprintf(fd,'n -3 - promien poszukiwan zbyt maly');
fprintf( 'n -3 - promien poszukiwan zbyt maly');
case -2
fprintf(fd,'n -2 - algorytm zakonczyl obliczenia, ale x nie jest rozwiazaniem');
fprintf( 'n -2 - algorytm zakonczyl obliczenia, ale x nie jest rozwiazaniem');
case -1
fprintf(fd,'n -1 - obliczenia zakonczone przez funkcje wyjscia');
fprintf( 'n -1 - obliczenia zakonczone przez funkcje wyjscia');
case 0
fprintf(fd,'n 0 - przekroczona maksymalna liczba iteracji');
fprintf( 'n 0 - przekroczona maksymalna liczba iteracji');
case 1
fprintf(fd,'n 1 - x rozwiazanie zbiezne');
fprintf( 'n 1 - x rozwiazanie zbiezne');
case 2
fprintf(fd,'n 2 - proces iteracyjny zbiezny do punktu nie bedascego rozwiazaniem');
fprintf( 'n 2 - proces iteracyjny zbiezny do punktu nie bedascego rozwiazaniem');
case 3
fprintf(fd,'n 3 - promien obszary zbyt maly');
fprintf( 'n 3 - promien obszary zbyt maly');
case 4
fprintf(fd,'n 4 - krok w kierunku szukania mniejszy od dokladnosci');
fprintf( 'n 4 - krok w kierunku szukania mniejszy od dokladnosci');
otherwise
fprintf(fd,'n inny przypadek');
fprintf( 'n inny przypadek');
end
fprintf(fd,'n output.iterations - liczba iteracji = %4d',output.iterations);
fprintf(fd,'n output.funcCount - liczba obliczania funkcji = %4d',output.funcCount);
fprintf(fd,'n output.algorithm - zastosowany algorytm = %s',output.algorithm);
%fprintf(fd,'n output.cgiterations - liczba iteracji CG = %d',output.cgiterations); % tylko dla 'LargeScale'='on'
fprintf(fd,'n output.firstorderopt - koncowe podstawienie w algorytmie Levenberg-Marquardt = %8.5e',output.firstorderopt);
fprintf(fd,'n output.message - komunikat koncowy = %s',output.message);
fprintf(fd,'n K O N I E C O P I S U n');

fprintf( 'n output.iterations - liczba iteracji = %4d',output.iterations);
fprintf( 'n output.funcCount - liczba obliczania funkcji = %4d',output.funcCount);
fprintf( 'n output.algorithm - zastosowany algorytm = %s',output.algorithm);
%fprintf( 'n output.cgiterations - liczba iteracji CG = %d',output.cgiterations); % tylko dla 'LargeScale'='on'
fprintf( 'n output.firstorderopt - koncowe podstawienie w algorytmie Levenberg-Marquardt = %8.5e',output.firstorderopt);
fprintf( 'n output.message - komunikat koncowy = %s',output.message);
fprintf( 'n K O N I E C O P I S U n');

% okrag x^2+y^2=r^2 w ukladzie wsp. x=-4:0.1:4, y=[-4 ]
r=2; % promien
x=-4:0.1:4; % zakres zmiennosci x
y1p=sqrt(-x.^2+r^2); % 1-szy polokrag
x1okrag=x;
x1okrag(imag(y1p)~=0) = []; % dla braku punktow okregu
y1p(imag(y1p)~=0) = []; % bez czesci urojonych
y1m=-y1p; % 2gi polokrag
plot(x1okrag,y1p, '-b',x1okrag,y1m, '-r');
axis([-4 4 -4 4]); hold on;
% hiperbola y=a/x, y=b/x
b=3;
y2=a*x.^(-1); y3=b*x.^(-1);
plot(x, y2, '--g', x, y3, ':k');
axis([-4 4 -4 4]);
grid on;
grupa='grupa C godz. 11-12.30';
nazw={[imie,' ' ,nazwisko,' ' , album, ' ' ,'a=' num2str(a),' ', grupa]};
title(nazw);
legend(' Polokrag minus', ' Polokrag plus', ' Hiperbola 1', ' Hiperbola 2');
xlabel('x'); ylabel('y');
fprintf('nn ... czekaj na zamkniecie okna graficznego ...n');
pause(10);
close;

end
end