Jump to content


Photo
- - - - -

Pomoc oko zadataka iz hemije, fizike, matematike ....


This topic has been archived. This means that you cannot reply to this topic.
68 replies to this topic

#61 tomas.hokenberi

tomas.hokenberi
  • Members
  • 4,349 posts

Posted 27 March 2014 - 18:48

Mozes da postavis ogranicenja tipa nejednakosti koliko god zelis - bitno je samo da nisu nesaglasne;

 

Daklem  ogranicenja su kao sto si ih postavio

 

x1+x2+x3+x4+x5+x6=80

p1*x1+p2*x2+...+p6*x6>=1600

p1*x1+p2*x2+...+p6*x6<=2000 (gde je pi prinos tona/ha kulture i)

p1*x1<=300; p2*x2>=100, itd



#62 AgroUrban

AgroUrban
  • Members
  • 4,027 posts

Posted 28 March 2014 - 21:58

Ok, hvala Tomase. :)

 

Tek sam poceo da radim u solveru, nisam siguran da li moze da prepozna greske u odnosu kao kad se radi rucno na papiru, pa tacno vidis da si zajebao nesto recimo ako si dobio optimalno resenje a vestacke promenljive ostale u bazi. :D



#63 tomas.hokenberi

tomas.hokenberi
  • Members
  • 4,349 posts

Posted 28 March 2014 - 23:45

Nego nije mi jasno zasto ne probas malo mocniji alat od Excela, npr imas open-source Octave:

 

http://wiki.octave.o...ctaveforge_pkgs

http://lpsolve.sourc.../5.1/Octave.htm

http://nadiah.org/blog/?p=74

 

## agrourban
## Author: ...
## Created: 2014-03-28   %Maksimiziras profit
P=[5000 3000 4000 1500 1800 2500];
%Pojedinacna ogranicenja su
O1=[1, 1, 1, 1, 1, 1]; O1g=80; O1t="S"; %X1+...+X6 = 80
O2=[20, 50,50,60, 25,20]; O2g=1600; O2t="L"; %a1X1+...+a6X6 >= 1600
O3=[20, 50,50,60, 25,20]; O3g=2000; O3t="U"; %a1X1+...+a6X6 <= 2000
O4=[1.9,4.8,5.1,5.8,2.3,2.1]; O4g=200; O4t="U"; %r1X1+...+r6X6 <= 200
O5=[20,0,0,0,0,0]; O5g=300; O5t="U"; %paradajz
O6=[0,50,0,0,0,0]; O6g=200; O6t="L"; %paprika
O7=[0,0,50,0,0,0]; O7g=100; O7t="L"; %krastavac
O8=[0,0,0,60,0,0]; O8g=30; O8t="L"; %kupus
O9=[0,0,0,0,25,0]; O9g=50; O9t="L"; %luk
O10=[0,0,0,0,0,20]; O10g=420; O10t="U"; %sargrepa
%Sada spakujes ogranicenja u matricu
OGR=[O1;O2;O3;O4;O5;O6;O7;O8;O9;O10];
BDS=[O1g;O2g;O3g;O4g;O5g;O6g;O7g;O8g;O9g;O10g];
TYP=[O1t;O2t;O3t;O4t;O5t;O6t;O7t;O8t;O9t;O10t];
VTP=["C";"C";"C";"C";"C";"C"]; % Setujes da je X kontinualno
LBX=[0;0;0;0;0;0]; %Donja granica po X - NEMA negativne proizvodnje
UBX=[80;80;80;80;80;80]; %Gornja granica po X - 80 ha
sense=-1; %Maksimizacija
[xopt,zmx]=glpk(P,OGR,BDS,LBX,UBX,TYP,VTP,sense)

 

i kao rezultat dobijes:

 

octave-3.6.4.exe:16> agrourban
xopt =
   15.00000
    4.00000
    2.50000
    0.50000
   37.00000
   21.00000
zmx =  216850

 

 

Edit: Kako su tvoja ogranicenja "po pojedinacnom povrcu" zapavo ogranicenja po Xi  mozes i jednostavnije:

 

%Maksimiziras profit
P=[5000 3000 4000 1500 1800 2500];
%Pojedinacna ogranicenja su
O1=[1, 1, 1, 1, 1, 1]; O1g=80; O1t="S"; %X1+...+X6 = 80
O2=[20, 50,50,60, 25,20]; O2g=1600; O2t="L"; %a1X1+...+a6X6 >= 1600
O3=[20, 50,50,60, 25,20]; O3g=2000; O3t="U"; %a1X1+...+a6X6 <= 2000
O4=[1.9,4.8,5.1,5.8,2.3,2.1]; O4g=200; O4t="U"; %r1X1+...+r6X6 <= 200
%Sada spakujes ogranicenja u matricu
OGR=[O1;O2;O3;O4];
BDS=[O1g;O2g;O3g;O4g];
TYP=[O1t;O2t;O3t;O4t];
VTP=["C";"C";"C";"C";"C";"C"]; % Setujes da je X kontinualno
%Donja granica po X(granice za luk,papriku itd.  0 je "prirodna" granica tamo gde nije zadato drugacije
% 4=200/50, 2=100/50, 0.5=30/60, 2=50/25)
LBX=[0;4;2;0.5;2;0];
%Gornja granica  po X(granice za luk,papriku itd.  80 je maksimum u svakom slucaju gde drugacije nije napisano. 15=300/20, 21=420/20)
UBX=[15;80;80;80;80;21];
sense=-1; %Maksimizacija
[xopt,zmx]=glpk(P,OGR,BDS,LBX,UBX,TYP,VTP,sense)


Edited by tomas.hokenberi, 29 March 2014 - 10:35.


#64 AgroUrban

AgroUrban
  • Members
  • 4,027 posts

Posted 29 March 2014 - 23:38

Nego nije mi jasno zasto ne probas malo mocniji alat od Excela, npr imas open-source Octave:

 

http://wiki.octave.o...ctaveforge_pkgs

http://lpsolve.sourc.../5.1/Octave.htm

http://nadiah.org/blog/?p=74

 

 

## agrourban
## Author: ...
## Created: 2014-03-28   %Maksimiziras profit
P=[5000 3000 4000 1500 1800 2500];
%Pojedinacna ogranicenja su
O1=[1, 1, 1, 1, 1, 1]; O1g=80; O1t="S"; %X1+...+X6 = 80
O2=[20, 50,50,60, 25,20]; O2g=1600; O2t="L"; %a1X1+...+a6X6 >= 1600
O3=[20, 50,50,60, 25,20]; O3g=2000; O3t="U"; %a1X1+...+a6X6 <= 2000
O4=[1.9,4.8,5.1,5.8,2.3,2.1]; O4g=200; O4t="U"; %r1X1+...+r6X6 <= 200
O5=[20,0,0,0,0,0]; O5g=300; O5t="U"; %paradajz
O6=[0,50,0,0,0,0]; O6g=200; O6t="L"; %paprika
O7=[0,0,50,0,0,0]; O7g=100; O7t="L"; %krastavac
O8=[0,0,0,60,0,0]; O8g=30; O8t="L"; %kupus
O9=[0,0,0,0,25,0]; O9g=50; O9t="L"; %luk
O10=[0,0,0,0,0,20]; O10g=420; O10t="U"; %sargrepa
%Sada spakujes ogranicenja u matricu
OGR=[O1;O2;O3;O4;O5;O6;O7;O8;O9;O10];
BDS=[O1g;O2g;O3g;O4g;O5g;O6g;O7g;O8g;O9g;O10g];
TYP=[O1t;O2t;O3t;O4t;O5t;O6t;O7t;O8t;O9t;O10t];
VTP=["C";"C";"C";"C";"C";"C"]; % Setujes da je X kontinualno
LBX=[0;0;0;0;0;0]; %Donja granica po X - NEMA negativne proizvodnje
UBX=[80;80;80;80;80;80]; %Gornja granica po X - 80 ha
sense=-1; %Maksimizacija
[xopt,zmx]=glpk(P,OGR,BDS,LBX,UBX,TYP,VTP,sense)

 

 

i kao rezultat dobijes:

 

 

octave-3.6.4.exe:16> agrourban
xopt =
   15.00000
    4.00000
    2.50000
    0.50000
   37.00000
   21.00000
zmx =  216850

 

 

 

Edit: Kako su tvoja ogranicenja "po pojedinacnom povrcu" zapavo ogranicenja po Xi  mozes i jednostavnije:

 

 

%Maksimiziras profit
P=[5000 3000 4000 1500 1800 2500];
%Pojedinacna ogranicenja su
O1=[1, 1, 1, 1, 1, 1]; O1g=80; O1t="S"; %X1+...+X6 = 80
O2=[20, 50,50,60, 25,20]; O2g=1600; O2t="L"; %a1X1+...+a6X6 >= 1600
O3=[20, 50,50,60, 25,20]; O3g=2000; O3t="U"; %a1X1+...+a6X6 <= 2000
O4=[1.9,4.8,5.1,5.8,2.3,2.1]; O4g=200; O4t="U"; %r1X1+...+r6X6 <= 200
%Sada spakujes ogranicenja u matricu
OGR=[O1;O2;O3;O4];
BDS=[O1g;O2g;O3g;O4g];
TYP=[O1t;O2t;O3t;O4t];
VTP=["C";"C";"C";"C";"C";"C"]; % Setujes da je X kontinualno
%Donja granica po X(granice za luk,papriku itd.  0 je "prirodna" granica tamo gde nije zadato drugacije
% 4=200/50, 2=100/50, 0.5=30/60, 2=50/25)
LBX=[0;4;2;0.5;2;0];
%Gornja granica  po X(granice za luk,papriku itd.  80 je maksimum u svakom slucaju gde drugacije nije napisano. 15=300/20, 21=420/20)
UBX=[15;80;80;80;80;21];
sense=-1; %Maksimizacija
[xopt,zmx]=glpk(P,OGR,BDS,LBX,UBX,TYP,VTP,sense)

 

Uh, kakve su ovo maDJije? :D Moracu da proucim lepo pa da impresioniram profesorku. :D

Zasto radimo u solveru? Paaa, tako na faxu radimo. Jedne nedelje radimo primer neke proizvodnje a sledece radimo test na racunaru u solveru. :) Tako da svakako na njemu moram da vezbam. Ovo je iz predmeta farm menadzment (organizacija malih poljoprivrednih gazdinstava).

 

,,%Maksimiziras profit
P=[5000 3000 4000 1500 1800 2500],,

Ovde bi trebalo da idemo umesto 5000, 5000x20 jer je 5000 dinara prihod po toni, a imamo 20 tona po hektaru.  

 

Moracu da proucim ove programe bolje.

 

Prosle godine sam slusao operaciona istrazivanja a sada imam par predmeta gde imamo njihovu primenu.

Jako zanimljiv predmet i izgledamo da je jako koristan. 

 

Hvala T.H. :)



#65 tomas.hokenberi

tomas.hokenberi
  • Members
  • 4,349 posts

Posted 30 March 2014 - 13:05

Uh, kakve su ovo maDJije? :D Moracu da proucim lepo pa da impresioniram profesorku. :D

 

Magija je open source alat za Octave ciji je "jezik" kompatibilan sa Matlabom

 

GNU Octave is a high-level programming language, primarily intended for numerical computations. It provides a command-line interface for solving linear and nonlinear problems numerically, and for performing other numerical experiments using a language that is mostly compatible with MATLAB. It may also be used as a batch-oriented language. As part of the GNU Project, it is free software under the terms of the GNU General Public License.

 

Sto bi rekli moje kolege inzenjeri "Excel je za sekretarice" ;).

 

Sto se tice profita, dobija se identican rezultat ako stavis 5000x20, ... itd, osim naravno sto je suma veca:

 

%Maksimiziras profit
P=[5000*20 3000*50 4000*50 1500*60 1800*25 2500*20];
%Pojedinacna ogranicenja su
O1=[1, 1, 1, 1, 1, 1]; O1g=80; O1t="S"; %X1+...+X6 = 80
O2=[20, 50,50,60, 25,20]; O2g=1600; O2t="L"; %a1X1+...+a6X6 >= 1600
O3=[20, 50,50,60, 25,20]; O3g=2000; O3t="U"; %a1X1+...+a6X6 <= 2000
O4=[1.9,4.8,5.1,5.8,2.3,2.1]; O4g=200; O4t="U"; %r1X1+...+r6X6 <= 200
%Sada spakujes ogranicenja u matricu
OGR=[O1;O2;O3;O4];
BDS=[O1g;O2g;O3g;O4g];
TYP=[O1t;O2t;O3t;O4t];
VTP=["C";"C";"C";"C";"C";"C"]; % Setujes da je X kontinualno
%Donja granica po X(granice za luk,papriku itd.  0 je "prirodna" granica tamo gde nije zadato drugacije
% 4=200/50, 2=100/50, 0.5=30/60, 2=50/25)
LBX=[0;4;2;0.5;2;0];
%Gornja granica  po X(granice za luk,papriku itd.  80 je maksimum u svakom slucaju gde drugacije nije napisano. 15=300/20, 21=420/20)
UBX=[15;80;80;80;80;21];
sense=-1; %Maksimizacija
[xopt,zmx]=glpk(P,OGR,BDS,LBX,UBX,TYP,VTP,sense)

 

Rezultat:

octave-3.6.4.exe:2> agrourban2
xopt =
   15.00000
    4.00000
    2.50000
    0.50000
   37.00000
   21.00000
zmx =  5360000

 

Jasno ti je naravno da xopt vektor zasejanih povrsina po "kulturama", a zmx maksimizirani profit. E sad ja se nisam bas trudio da proveravam da li sam dobro uneo tvoje podatke, to ostavljam tebi da uradis kada i ako instaliras i isprobas Octave.

 

Inace postoji cak i verzija za Android platformu (ja ga imam instaliranog na Android tabletu).


Edited by tomas.hokenberi, 30 March 2014 - 13:17.


#66 AgroUrban

AgroUrban
  • Members
  • 4,027 posts

Posted 31 March 2014 - 10:09

Magija je open source alat za Octave ciji je "jezik" kompatibilan sa Matlabom

 

 

 

Sto bi rekli moje kolege inzenjeri "Excel je za sekretarice" ;).

 

Sto se tice profita, dobija se identican rezultat ako stavis 5000x20, ... itd, osim naravno sto je suma veca:

 

 

%Maksimiziras profit
P=[5000*20 3000*50 4000*50 1500*60 1800*25 2500*20];
%Pojedinacna ogranicenja su
O1=[1, 1, 1, 1, 1, 1]; O1g=80; O1t="S"; %X1+...+X6 = 80
O2=[20, 50,50,60, 25,20]; O2g=1600; O2t="L"; %a1X1+...+a6X6 >= 1600
O3=[20, 50,50,60, 25,20]; O3g=2000; O3t="U"; %a1X1+...+a6X6 <= 2000
O4=[1.9,4.8,5.1,5.8,2.3,2.1]; O4g=200; O4t="U"; %r1X1+...+r6X6 <= 200
%Sada spakujes ogranicenja u matricu
OGR=[O1;O2;O3;O4];
BDS=[O1g;O2g;O3g;O4g];
TYP=[O1t;O2t;O3t;O4t];
VTP=["C";"C";"C";"C";"C";"C"]; % Setujes da je X kontinualno
%Donja granica po X(granice za luk,papriku itd.  0 je "prirodna" granica tamo gde nije zadato drugacije
% 4=200/50, 2=100/50, 0.5=30/60, 2=50/25)
LBX=[0;4;2;0.5;2;0];
%Gornja granica  po X(granice za luk,papriku itd.  80 je maksimum u svakom slucaju gde drugacije nije napisano. 15=300/20, 21=420/20)
UBX=[15;80;80;80;80;21];
sense=-1; %Maksimizacija
[xopt,zmx]=glpk(P,OGR,BDS,LBX,UBX,TYP,VTP,sense)

 

 

Rezultat:

 

octave-3.6.4.exe:2> agrourban2
xopt =
   15.00000
    4.00000
    2.50000
    0.50000
   37.00000
   21.00000
zmx =  5360000

 

 

Jasno ti je naravno da xopt vektor zasejanih povrsina po "kulturama", a zmx maksimizirani profit. E sad ja se nisam bas trudio da proveravam da li sam dobro uneo tvoje podatke, to ostavljam tebi da uradis kada i ako instaliras i isprobas Octave.

 

Inace postoji cak i verzija za Android platformu (ja ga imam instaliranog na Android tabletu).

 

Ajooj, jeste meni krajnja titula inzenjer agroekonomije, ali daleko sam ja od inzenjera, ja sam samo prost sumadijski paor. :)

 

Kad mi prodje ispit u petak videcu da proucim kako se instalira ovo cudo, mada za mene je i instaliranje solvera bio veliki uspeh jer sam duduk za te stvari. :( Sto se same tvoje postavke tice, to sam ispratio, to mi je ok.

Ali pre toga, zanima me da li on moze da odradi i analizu senzitivnosti, znaci treba da mi pokaze podatke o konkurentnosti proizvoda, cene u senci, limits, itd...

 

A sto se tice Androida, jbg, koristim Lumiu 1520, (Zarazio me miodrage sa WP8. :P )



#67 tomas.hokenberi

tomas.hokenberi
  • Members
  • 4,349 posts

Posted 31 March 2014 - 19:18

Ajooj, jeste meni krajnja titula inzenjer agroekonomije, ali daleko sam ja od inzenjera, ja sam samo prost sumadijski paor. :)

 

Kad mi prodje ispit u petak videcu da proucim kako se instalira ovo cudo, mada za mene je i instaliranje solvera bio veliki uspeh jer sam duduk za te stvari. :( Sto se same tvoje postavke tice, to sam ispratio, to mi je ok.

Ali pre toga, zanima me da li on moze da odradi i analizu senzitivnosti, znaci treba da mi pokaze podatke o konkurentnosti proizvoda, cene u senci, limits, itd...

 

A sto se tice Androida, jbg, koristim Lumiu 1520, (Zarazio me miodrage sa WP8. :P )

 

Na PC-u se instalira lako. Skines instalacionu datoteku, pokrenes setup i "next, next, ...".

 

Sto se tice analize osetljivosti naravno da moze. Sve sto je potrebno jeste da ti razumes matematicki model problema i da umes da ga jasno postavis.

 

Za ono sto tebi treba ovaj tutorijal ispod je prava stvar:

 

http://lipas.uwasa.f...with_octave.pdf

 

P.S. Miodrage sa mobilnog sveta?



#68 AgroUrban

AgroUrban
  • Members
  • 4,027 posts

Posted 31 March 2014 - 19:31

Na PC-u se instalira lako. Skines instalacionu datoteku, pokrenes setup i "next, next, ...".

 

Sto se tice analize osetljivosti naravno da moze. Sve sto je potrebno jeste da ti razumes matematicki model problema i da umes da ga jasno postavis.

 

Za ono sto tebi treba ovaj tutorijal ispod je prava stvar:

 

http://lipas.uwasa.f...with_octave.pdf

 

P.S. Miodrage sa mobilnog sveta?

Ono sto sam ja citao iz onih tvojih linkova, fazon, jedna recenica objasnjenja, pa milion nepovezanih slova. :D Ako je samo instalaciona datoteka onda ekstra. :D

 

Zanimljivo je da su na jednom linku objasnjavali primenu bas na setvenoj strukturi psenice i jecma, i cak ima graficki prikaz. Doduse on nesto znaci samo u slucaju kad imamo samo dve promenljive.

Razumem ja matematicki model, prosle godine sam imao klot operaciona istrazivanja sa nekim random primerima pa se razmisljaj, ovako kada radimo konkretne primere pa kad ti kaze: Imas toliko i toliko radnih casova na raspolaganju, normalno da je to ogranicenje

 

. Za analizu mi treba jos pored strukture setve i maksimalnog prihoda konkurentnost proizvoda, dokaz da je model pouzdan, cene u senci, i ,,limits'', ne znam kako se to zove kod nas, to je ono (na primeru) ako imas 3 parcele i 5 kultura, i dobijemo da se psenica gaji na prvoj parceli, kolika bi bila funcija cilja da se psenica ne gaji na toj parceli. Ne znam u teoriji sta sta to ispituje. :-/

 

 

EDIT: Da da, MioVrag stari. :)


Edited by AgroUrban, 31 March 2014 - 19:33.


#69 AgroUrban

AgroUrban
  • Members
  • 4,027 posts

Posted 22 April 2014 - 12:47

Da li neko ima neku literaturu iz koje bi mogao da mi napise definiciju operacione krive koja se tice statisticke kontrole kvaliteta. Kupio sam knjigu ali sam izgleda slep ili nema te definicije a potrebna mi je.