114

1 14 Rozdział 9

fprintf(fd,'\n -10 - xl*x2 <= 0 ' ) ;

fprintf(fd,'\n - nieliniowych równościowych');

fprintf(fd,'\n -46 + xl~2 + x2~2 + x3~2 = 0 ' ) ;

fprintf(fd,'\n oraz ograniczeniach nakładanych na sterowania');

fprintf(fd,'\n    -10    <=    xl    <=    1');

fprintf(fd,'\n    -5    <=    x2    <=    5 ' ) ;

fprintf(fd,'\n    0    <=    x3    <=    inf');

fprintf(fd,'\n\n\n*** OPTYMALIZACJA ***');

fprintf(fd,'\n punkt startowy    x0 = %f',x0);

[f0,fgrad0] = optnlcfg (x0);

fprintf(fd, '\n Funkcja celu przed optymalizacja fO = %f' ,f0) ;

fprintf(fd,'\n

===============================================');

fprintf(fd,'\n rozwiązanie optymalne x = %f',x); fprintf(fd,'\n Funkcja celu po optymalizacji f = %f',fval); fprintf(fd,'\n exitflag = %d',exitflag);

%UWAGA !

% zmienna output oraz lambda to tablice strukturalne

%

fprintf(fd,'\n output.iterations = %d',output.iterations); fprintf(fd,'\n output.funcCount = %d',output.funcCount); fprintf(fd, '\n output.algorithm = %s', output.algorithm) ; fprintf(fd,'\n lambda.eąnonlin = %f', lambda.eąnonlin); laml = lambda . eąnonlin (1) ,-

fprintf(fd,'\n lambda.eqlin = %f', lambda.eąlin); lam.2 = 1 ambda . eqlin (1) ;

fprintf(fd,'\n lambda.ineąnonlin = %f', lambda. ineąnonlin) ,-1am3=lambda.ineąnonlin(1); lam4=lambda.ineąnonlin(2); fprintf(fd,¹\n lambda.ineąlin = %£', lambda.ineąlin);

lam5=lambda.ineąlin(1);

fprintf(fd, '\n lambda.lower = %f' , lambda.lower); lam6-lambda.]ower(1);lam7=lambda.lower(2);

1am8=1ambda.1ower(3);

fprintf(fd, '\n lambda.upper = %f' , lambda.upper);

1am9=lambda.upper(1);    1aml0=lambda.upper(2);

fprintf(fd,'\n ==============================\n\n');

fprintf(fd,'\n *** Współczynniki Lagrangea ***');

!prir.t.f(fd, ' \n - równościowe nieliniowe'); i f 1ami

f printf ( id, ¹ \n poniższe ograniczenie jest aktywne: ' ) ,- end f pr i nt. f ( f d , ' \ n x IG 2 +x2 ~ 2 +x3 ^ 2 - 4 6 = 0 , laml = % f ' , laml) ; fprintf. (td, ' \n - równościowe liniowe'); if lam2

fprintf(fd,'\n poniższe ograniczenie jest aktywne:'); end fprintf(fd,'\n xl+x2+x3-10=0,    Iam2=%f',lam2);

f pru nt f ( f d , ' \n - n i erownosc .iowych nieliniowe'); if .laml

; end

fprintf(fd,'\n poniższe ograniczenie jest aktywne:')