From fdcf4b4a3d33eeb9ced8c16c489119d6ce6fbadd Mon Sep 17 00:00:00 2001 From: szymon Date: Sat, 20 Jun 2026 14:17:14 +0200 Subject: [PATCH] blabla --- media/Pasted image 20260616200137.png | Bin 0 -> 1896 bytes media/Pasted image 20260617124924.png | Bin 0 -> 8676 bytes .../2_semestr/metody_numeryczne/egzamin.md | 259 ++++++++++++++++++ 3 files changed, 259 insertions(+) create mode 100644 media/Pasted image 20260616200137.png create mode 100644 media/Pasted image 20260617124924.png create mode 100644 uczelnia/2_semestr/metody_numeryczne/egzamin.md diff --git a/media/Pasted image 20260616200137.png b/media/Pasted image 20260616200137.png new file mode 100644 index 0000000000000000000000000000000000000000..655215dbee4ad64720abeb0ab010dfb78b2d014c GIT binary patch literal 1896 zcmZ`)c~sKr8vZGksiC2yP|!4)=}kt>n=804$kat~8A(k=k7eb);8yDBlwpE5!!o0A z%jLS1ljw|wnj@~Lxk73+N{MK=q?nfGP&0GRy?2^FzVA8Td!9eO^S;mfJQNo^7Oc2W z5dZ)%&dJ_YisjOCcGphnO_wN)N9=eNbYf6e zrJ&QfsJR*P^$!In@Qxm%7})OOQQBTLD)p>O+)pybcaX&rk&aFmhhg^GU?2EAc3xxI z&!xpGRmF00H*+9+Z||+C=W)KB7DqJe!=9(-pMKe480eFmo9l1uYwYDUek^VQ)pzqz zg$Pw2QqAGs9$p$8&m}gl;uz-7vpgz!MNJbGWYik-w@KE)cDc0WJ zUspD2!?z(U&9!Sk6=*x3`ZtrAUQ&YId73GSfWJ?(Id0P6Xuy{*eYV?LliGBtqcfjFYp?Ob>*rR7616u8+?Wm5Ft6A-Gd9TBS76POjN2{ zwrJq}acu(8_rFZphK3@2^h}QG>Wv zBJyft(jg0hiiBETVFqNE$vYOyJZOH4I$2s9=>$XCC6Jy4oa5AoCi!naQUY9jFpvRH zeW)ao8a{F%YplQff%%uU_@_W=8P|h1|Kk=|N5)%D|2PCLa8(_VMUR##qFe^oVxch~nWATzwD}y?JknT;-C;m?Zg&KaE(4qcHCHU<)JFD}Duc)XnTd$`HNTBtL$wXAO4vM=!#zHH| zw5R9Wr_T~jjU1m6yFz#D`PiNHnZv`fa()>w3-dVww2B_@=$VC2hk9<2KA4ZqkgxXyoh5JE+BaNVP< zBUsG1*nZ=$2K5J3eL0@!?JZAEiQU6yvt8(<_^KaeIh-A}+|qn|m5k*pq~*Xl|QRka&28dI4n&4%s}$4Z@v-h9PJqnF;fe|RQ;`oWahtv4{5%j zpOXW!jA#@bOhJ`|;m%7)!!yr`1}Vm8PBODbM$iOeUIyyv56f5D24}SR(!3>^e&Q*>FTQK z+1N#(W-|k zW?|*L;?R(z|MGMvVa*@M2zX&h%eefA=v5m!ae%-5Mizi-IJd8iDmjF*-%=B)gX3l! z@GpeIpNl9B3_(1MQW}Gvl*SAPi}}B(`83HCP;w7aPziI;lSUhrI4x3Rm5`|Qy*0Hle80y= zHqNF~bva7w9j&A}U?BC4d3yUMMBAj>42vK^?ZylY^85M!<3oB=j|Vm0GQ)gvJ93B_ z)c))KHo4mAq5)g@pM2_ve~-Y&Evzw`F*R_#wa8;YWB#jt|8Ies24VvfdoNa+`L-%2;6XwdsKg7`}!_&(+u z$dxz{1|{oEjj^w+JSqx4>mtQ)o?O&Vlv!V0Y?z|wS)8767t=22t>ZDTlZ>)%Jw8X? zsEl0}-f&USO24mt@>o%?X9GM%b-J4!xGS-r)3nY1MwBga zd~J&@*)U#BiMsOb#6 zkyhsNf{7P}2vxTf{a5A|rh8PJ;=Rq=O_6R5t;`4aODZ^TryH`9rP9Jd1s_@a>hPh&TR_yng79GKhWNUn%w~PpNsoOO|Vqkp+KUujo zNXeLSW~h?g2p==7V3tH_UgwZZqjf3+O+n=Y`}BWw9mLRi6}mT5uXbV)g!kJr`N zx%Oyk9Jww++NmY+jOx%|LamwcPA*73?;$~c701TcY0}6g%s?y8|4l3Y|8?z|p*<&T zWRSlx#GKl<$H>~>|X^QWS)rIkhPFuf>iKEDmaG2M1&X^WWcdyN5fRtmfMA~pc z)35JIxJLP!AB%URysXKDiqZB(dCZlP9cF@vdDhx2RjWpvmYI>cqFKA8b>|kPF*d zwFZF1ujuuw#rp2I-NG8abGOz@q4S%xHdtC&3)Y>PHmo%xD(&-eL$b&Vo?OzBv5t-o zadD!cl?I=_19hBlsj_KqM~9v7VLKZ8qEkAqubgl%mw-Q z*BAqn?r>%gBbDuvH>|=4;T$EUUE7H~s^u~NUi`EnmXJrEPre${h9mMmLC&8Q%jO?- z3T$pFnf2=jNS%17F(_1aWbaL*{tfI(q@zn27z3deD7@RQ5!+l}w^BU;Pw zv!|4lM(Sm<8msQ{;D(>6m}7S)whIxo`$YZcK1D0sya14yRnB2H05s{Lw;HQln=u$` zjt780k6%E{!)tjaCc^c^>b^b_%A2p=KejrVMN^ynn8)dDR-cWNx36~i{#I9Yww!g$ zdTYz-BaOt&z{PYsNqwJzfMu`o)?_*xyWFd&?+IAc;k{SeMeQ;cakioW1WVMlTiC_N zx|e>Svn?|sB?#lmQLR6Iy!4s&Kfynz{}67~AF$jfQ5S^aY0a{u3YPoSktN5wa_mLT zF)#HP81Hs?L5Vtm+weAHt%{#Y$Osen+qGA+b!5jHwQ}0>3zR#^4nMy>lEG>bn2f24 z<$hsxpDCvO#fuk9;MH-1_xMI*2RM$^PlY4ctg)<+5oBZ+q3W2Qx`&E-hmws2aawb2 zD#Escc(e7Le#y_rGUlg?oWH+4rL64u`LC@qXM5k9Bzb<`D$RIe$GNf8Wgph-)8mFV zQa@oVzI}}kbL6~^8$B-$qxmM3dC8}Chm45dNkPZ+As-i3=Z9;YsOw7xrlzJ#K!_-1 zL~D@(Q5r=VoM+4SS{*^`%l-;Kr;d+TuaedmSq*kv1)b0x7s=)9(O9~3(laJ2VM8NS z#Itta1`CsajuOH5qYoYNBO&S0mEka&Y*8=S%DqHVzZJ0MsfenA6NNC!9GhvsRq%Ab ziX3SfPb%WDI&JWMUbFq!ta9cB$u*26PdMM@8)NbK<_KfWXrK=VfVuJSNqVcxvkC5f zr5ShnpbQdBlFnt^nJVmYZD^jXnJQ z^r2urLEGR4#8Tx%Ib!*2fHSH>zgr{9D#!|ujk58%+KxjkPlfa4^LN^j0(|_^9*CXLB=4z%8u?+bA2wgbD+f{wpN1ZRj|#tWA@1KrY~RHM;?L+O2{Y_@V_ zOquQwal4T2xI2@_x#X-^%iBJGR@2quJJ;PS7CU;0lNmB`>TE^?uH&kcS^oCKn=HrL zs$*-#9wFPSV#Xfq{qN}Rs$*6$$xY7-jMw-hk?f}cgmf!HIz7jI-2BrnBy#8f9$mzL zf@Ej_ZpmBhuu3IuaM7|e!Y=-3q+rV?T##7<2F5=Kl|%F!Ph>rp=uTra5*MerK5jfV z5@pwnDnU5GB20WpO{4rH`_DsR${tYE!pgU=C=d9v{Fm|)n)>}_<3Wt1Ph8DMKtI*l z`c1!5#f>s*w77R2n_^xKaF>dr#y2N^e26~f@@LA$*huuP zN*H7V=7|mF2{YQeKR{^W$qpOq%XZ^z5So(uQbuuCiY_AJW;%fEX4p*=ZbP3R(C>T^ zP;_dJGmKV{v~nwG=N}UXC3hRjgznXLk=2W%ndRp$7AxW%!=`dR=NyeDXXLi8`JbC? zp^`Sr=RR*6?=Emv1L~dPfcjzqK04Hh;v2I_>yxX3n^%ygjgI@ugD9BBsp6P)T43Q- za0OA73(2>Ohe8Nk_89B)wiEB{sDX>(`Xm32eNixLUXx2)YRa?+chqA-Y3tx&Aws@B zO>4uS?9iR>GTxrvc_Hg+_O&x)v`Y_VJ)7t&!#g&|MUyx+q(-^2sXrG*g2Hqe-=$>B zbn-rJ+3x`S%&l?-?6eKj&)0^gYzFv@2Sj#C4hJFKIg(j6Q=fhPizlF?dG>FE>yI?L z-vJ6q>!bDqH9qAbkrQ)Y1ZgOMZ#z+b^6rX`#s9cj}1xRJKm_v;9#L1olW2U|a&+0ZVr)>PzD! zUC#Vbbct_j@=53nwRVO+$ZSx?vz;f81k}3g1zTi(^VK#u$UIQ@x1#v=Pw~-6EB88Bok(=vlFmKzlokxy zQ@8EJpuHNDk)rp+4~@kirX{Nr>59aKU%(04iiZR>-!!QVXnYA3q`~;j7Jo;a&oRE= zSpTJ`RlhLamQD$mK;d4w+c~Kh`dCIJK^^a2>N9Qm{TSj)Iuy2N(n)LWDp=>_G$#pt zHSdX+uq<+^@jjNe;BY1C{@NeV-g@ zg>&U*qV=u4Eq=i^Ca##4D~|LS&!RCYf37!3QfGS#wA*tb3`nbOi*L5z1F6?Gu=F(u(+* zUJ(FM*QxozFTJsTiV`}S=2IBA(Ap&3SQU@{3t^p#Gh9m=1)DdX5}RV9zQJx$IR2z! zBYs16=?a&UFu0NIoKp(Uc&aP6cR{(e&iVt=y{dY2GRN%h1gcH-984X_7ksLoR&YkS)|HSkAzU8T zmut1p&ZAK*q5E)FOv1I3O~HA_4@@1Y$83&92?yJYtAGz zb$DEqRlU>WGWLPs_qB@t=!&JOG^e>UqmtvVtQF?Tx67If&9bfw>{Ih{1~@Hg#eE48 zXmCl5ms;w6c;nD)_Lm(EJc1gR7Pp%x?R6@DU9fDWK5NT#k0G+xvt_S-krDV1Q2LO1 zphMv8faW^U#73V_6kw$pF)t@TTmr|qbN zQG1@{yMkQv(v;GH^6!VS%p11WLsyo)jMwtB8iuJ6_}n8;CYV<;nxyCwT8)Bvo#how zDser_FhLU=u(p&z+~e9(A7~vb5wDJK{;nq;4@yoIWLh9^!+VfliMp|%G)q(cjG8{E z%GKPot`4|1xIH*-;d{&e}9iG`>Y%q-3 zraNduO`Ea|qBbhqIDnCt<-^n_dP{JX8|$+m(QIpZ_gWvA>pqERI_IEDUbe97X8>*a zd>H4XR^u|W!usVC&%oxaVKrxAPe+%j}2o^9mV=B zeI#y?4XRE0w%Sotxfq_g*3eZZz5H>=ti{thTt@TMjI-ZnG4Ps?(KX5N&p-LsE|KuV zX%r;?8fb6i5*T3gp!a0H!^&f%pp(Njy<#c`yoO3pGa54S6)<;aVVFDuYcuvMDXh7i!17!j zPaT7?We|7$KKyO3o#6&|9oMQaZHDE2|DrZ&7qkowJdJ`|cYYjk;zI*zWLNgBe&69| zQ{W6ynbB?<^q2tS`PGqO1+-+bxcvxdeAC?7$~EExwHIIGHG<4NxE!|X_`8VQCgmPf zq)>+WB|yby^h3kE0QORWz^fTp76%w(j@{~v;S@)d$OIJ-OvRydA21 z<2>t+wm*Fbg19}?Z+S)-^1XINsU5#QQVr=x|=azWV8iN8mablF^<4t`8a}C2$u~**(XCGd(tzJpL9Hu&5dNhYpw@frU3EraqHsj3BqNF>3?_p`f%rm_FE;)dYN>} z=Y55~^FOzXYigi+FOV6L;7}EW5mvf&Ib^#qIh8^@20ypjg+9Ch19WxS%fsjPnLU@Y zH+M&h1DB(8q~JR@C<5`bJ~<>LgfmIyx$D(OzQMrp*>7x2KAuX!!&R}GbQ%;$| zmCyOa{ZE3F@IhJKHM2wts@pEmJ+Dc*0kW-d;GJMW&C}3@>35Erx%@lVAQ?S6gd^$> zp}U`BU-JR69{X=gdH$!@)V7heE)iMWGY&8bk7D^+}E9Y3@Exx;y!JYvtpN2 zz^yb6ehO~u>UymoA(aet?sDr^{VNOG=zBbTqo4*q5)Wz1qo0I{xRYWSO$NA?IYH%k zib3F`vA)n&0&Uw05gtSXt)keV_OlSH%!3R@LlZi5Jr1NZ>;-gO4RlBIfY306&b#Q#mzas_+hzx!V$DSu99?l-|9e=UCp7@h=AnCxyDa;?^wr z|CIb(pbbDSgJP@s-J{w)2bn>oT~Hsb=@wZi{#!t3%0Gr?;?>pFhObW^gxw&Ir4m$A zfQXxA00gFYS_hfa!d&l_q5zm#I*~R7xn`CZ$dXaHOD-wT%ET zc57ytr;ohQ*`Z(L-j@q7_PB#5UvmcD86@%=I@u`5prNUV_YV^kI+=}%GzUg)c*zEL zElLaT5q)>K`sUIUoy}dQ+I_zj)KM1?0htfPOd610pcXvEIDsODCBB0sU>t^)=ejW$ zs^!r5TAoeY;Lv(W8|%A2@;X;QH9XP`SbYZ(&8%M_gm=Q4{=U`EMPq3JGq^*7gZ9i} z%S2r@z_{Fg;^av2@bh#czjivqNUY^&Ie;;m^(kDd3Yr1|q))c6B}nEuHr2m4sH5OB z)issdAh$bfvojrnmY-fzUGjTveTP6iuwn-sVDNO6LkKsf zmPfNt@+w(2SC`V1CqBL^HXgp0`A0{cnyarc(%@?ZS+n=~P|gvp)!Iq7;Bk7#v5 zm2~bYgFM~B+5~2*V0oq<(SewfYC&vw5k?-V6IvXbC|GsyPwcM=1Dex{Qudk)rPZhS zVL$PqhI24#u@`0sw4~#`T8{<9^Wd&+aj06O%tljGHEXh3?UA`U_?Z2F;_OeyNEtpZ z%lGiCz!&1hPxvnFgq@^XAnm_xG->6FxOq=!uhyE=liS~?d}#}Y_WLY)hLO2>xDU@L zEg>PoNbP9P?YxU>Lcl%?3#+!1{;l_$B*ECpgk($aLuaG%5IiRU`^Y(iyN!%{esNk1w&dp_cY?~ z6^NJ3$NUGAw=)ZDAs*Ten^X_%r@IRbz|^!|zZ-hog(@-J^Nc1sJ#(H$kqPNf)m;1l z8;1hvH~A28Lk33N0CpeI*emPKN7&Ruwbo54^>`NE7fYfp4xG?hiUJR5Li>QYp{p1v zXH)6`34_%*fE7uOL}8DK07!b<@>bcw{$ITfOgm~kV@0)~^gYr29B%`ieXZBuP|Xp` zO{D!gGAh9H8sJl#oPnsO=RDdMYQjMLi;a2aK-z~cg9^zPk?V`Ac(%M(b>#D HuHFAPjqBAw literal 0 HcmV?d00001 diff --git a/uczelnia/2_semestr/metody_numeryczne/egzamin.md b/uczelnia/2_semestr/metody_numeryczne/egzamin.md new file mode 100644 index 0000000..4bce87f --- /dev/null +++ b/uczelnia/2_semestr/metody_numeryczne/egzamin.md @@ -0,0 +1,259 @@ +# Węzły +Węzły interpolacyjne to inaczej nazwane punkty przez które przechodzi interpolowana funkcja. Żeby możliwe było znalezienie rozwiązania takiego zadania w każdym węźle koordynat X musi być różny. + +# [Efekt Rungego](https://pl.wikipedia.org/wiki/Efekt_Rungego) +Efekt rungego polega na pogorszeniu się jakości interpolacji pomimo zwiększenia liczby węzłów (punktów). Początkowo im więcej węzłów tym lepsze przybliżenie, ale od pewnego momentu zaczyna się pogarszać, co widać szczególnie na krańcach przedziałów. To zachowanie dotyczy [interpolacji wielomianowej]. Aby to zjawisko wystąpiło muszą być spełnione dwie przesłanki: +- wysoki stopień wielomianu +- równa odległość pomiędzy węzłami interpolacyjnymi + +# [Macierz Vandermonda](https://pl.wikipedia.org/wiki/Macierz_Vandermonde%E2%80%99a) +Taka macierz: +$$ +A = \left(\begin{matrix} +1 & 1 & 1 & 1 \\ +1 & 2 & 4 & 8 \\ +1 & 3 & 9 & 27 \\ +1 & 4 & 16 & 64 +\end{matrix}\right) +$$ + +Albo bardziej ogólnie: +$$ +A = \left(\begin{matrix} +1 & x_1 & x^2_1 & \cdots & x^{n-1}_1 \\ +1 & x_2 & x^2_2 & \cdots & x^{n-1}_2 \\ +\vdots & \vdots & \vdots & \ddots & \vdots \\ +1 & x_n & x^2_n & \cdots & x^{n-1}_n +\end{matrix}\right) +$$ +Wyznacznik tej macierzy to specjalny **wyznacznik Vandermonde'a** i jest zapisywany w postaci: +$$ +\det A = \prod_{1 \leqslant i | $f(x_0)$ | $f[x_0, x_1]$ | $f[x_0, x_1, x_2]$ | $f[x_0, x_1, x_2, x_3]$ | | +| $x_1$ | $f(x_1)$ | $f[x_1, x_2]$ | $f[x_1, x_2, x_3]$ | | | +| ... | | | | | | +Czyli dla przykładu: +x: $-3\ 1\ 3$ +f(x) $8\ 16\ 24$ + +| -3 | 8 | $\frac{8}{4} = 2$ | $\frac{1}{3}$ | +| ------------- | --- | ----------------- | ------------- | +| 1 | 16 | $\frac{8}{2} = 4$ | | +| 3 | 24 | | | +| | | | | +| Współczynniki | 8 | 2 | $\frac{1}{3}$ | +$W(x) = 8 + 2(x + 3) + \frac{1}{3}(x + 3)(x - 1)$ + +# [Hesjan](https://pl.wikipedia.org/wiki/Macierz_Hessego) +Elitarna macierz drugich pochodnych cząstkowych fla funkcji wielu zmiennych o wartościach rzeczywistych dwukrotnie różniczkowalnych w punkcie, w którym liczone są te pochodne. +Cała ta macierz opisuje własności krzywizny wykresu funkcji w otoczeniu punktu. Jest wyznaczana np. w punktach krytycznych przy wyszukiwaniu ekstremów i punktów przegięcia (przegioles). +Macierz ta jest kwadratowa. + +Buduje się ją w sposób następujący: +Musimy policzyć wszytkie możliwe pochodne drugiego stopnia. Po x, po xy, po xz, po z, po zx, itp, itd. Potem układamy to tak: +- po przekątnej pochodne po jednej zmiennej +- w kolumnach ustawiamy pochodne po kolejnych zmiennych +- w rzędach zmieniamy naszą pojedyńczą zmienną pochodną + +# [Wersor](https://pl.wikipedia.org/wiki/Wektor_jednostkowy) +Wektor jednostkowy o długości jeden. Wskazuje kierunek i zwrot pewnego wektora początkowego, któremu ten wersor się przypisuje. Każdy wektor można zamienić na wersor (znormalizować) dzieląc go przez jego długość. +![](../../../media/Pasted%20image%2020260617124924.png) + + +# [Aproksymacja sześcienna] +POlega na przykładaniu funkcje szcześciennej do wykresu. Jako punktów uzywa wartości funkcji i jej pochodnych. żeby to było możliwe to oczywiście musi być na krańcach przedzialu pochodna wymnożona na minusie tego typu +# Wyznaczniki macierzy trójkątnych i diagonalnych +Oblicza się je mnożąć przez siebie wszytkie liczby leżące na głównej przekątnej. +Macierz diagonalna – ma liczby tylko na przekątnej. +Macierz trójkątna górna –» ma zera pod główną przekątną. +# [Metoda falsi](https://pl.wikipedia.org/wiki/Regula_falsi) +Falsi z nazwy metody oznacza oczywiście, że jest do dupy. + +Bierzemy jeden punkt powyżej zera, drugi jest poniżej. Rysujemy prostą (falsi[znaną również pod nazwą prosta sraka]). Tam gdzie sraka przechodzi przez oś OX liczymy wartość funkcji. Zamieniamy to na ten poprzedni punkt. Powtarzamy, aż osiągniemy dokładność. +# [Metoda Powella](http://www.optymalizacja.w8.pl/Powella.html) +Metoda bezgradientowa, która wykorzystuje wersor. + +Mamy dwa wektory (kierunki) $d_i$ i $d_j$. Nazywamy je sprzężonymi względem kwadratowej dodatnio określonej macierzy H jeśli: +$$ +d_i^T \cdot H d_i = 0 \ \text {dla } i \neq j +$$ +# [Metoda najszybszego spadku](https://pl.wikipedia.org/wiki/Metoda_najszybszego_spadku) +Metoda gradientowa. +$$ +\begin{align} +P_{k + 1} &= P_K - \alpha_k \nabla_k \\ +\text{gdzie:} +\alpha_k &= \frac{\nabla_k^T \nabla_k}{\nabla_K^T H_k \nabla_K} +\end{align} +$$ +Co to dokładnie oznacza? Każdy kolejny krok to aktualny krok minus gradient razy krok. A krok $\alpha$ to transponowany gradient razy gradient podzielone przez transponowany gradient razy gradient razy hesjan. + +Problemem tej metody jest to, że potrafi zwolni i zyzgakować, jeśli funkcja ma kształt wąskiej, ukośnej doliny. Żeby temu zapobiec: +- wprowadza się nowe kierunki za pomocą $\Delta = P_{k + 1} - P_K$ +- obraca układ współrzędnych, by jedna z osi pokrywała się z dnem. + +# Warunek przy szukaniu pierwiastków równań nieliniowych +$f(a) \cdot f(b) < 0$ +W komputerach może występować zjawisko niedomiaru, które wynika z zaokrągleń do zera przy bardzo niskich wartościach. Aby temu zapobiec robi się taki myk: +```matlab +sign(f(a)) != sign(f(b)) +``` + + +# [Metoda połowienia](https://pl.wikipedia.org/wiki/Metoda_r%C3%B3wnego_podzia%C5%82u) +Algorytm poszukiwania pierwiastka funkcji. + +Metoda polega na dzieleniu przedziału na pół, a następnie odrzucaniu podprzedziału, w którym dwa jego krańce nie zmieniają znaku. +# [Metoda złotego podziału](https://pl.wikipedia.org/wiki/Metoda_z%C5%82otego_podzia%C5%82u) +Metoda optymalizacji jednowymiarowej funkcji celu. + +Polega na wybraniu dwóch punktów, które wyznaczane tak, by dzieliły aktualny przedział w proporcji złotego podziału ( $61.8\%$) Wybieramy punkt, który jest niżej, aż osiągniemy zadaną dokładność $\epsilon$. Wtedy bierzemy środek tego przedziału i *beng* mamy minimum! + +# [Simplex](https://pl.wikipedia.org/wiki/Algorytm_sympleksowy) +## Jak to liczyć? +1. Mamy układ równań, przekształcamy go w taki sposób, by po prawej stronie wszystkie wartości były dodatnie +$$ +\begin{align*} +&f(x,y) = 7x + 2y - 3z \to \max \\ +&\text{ograniczenia:} +\end{align*} +$$ +$$ +\begin{align*} +\begin{cases} +2x + y - 3z &= 12 \\ +5x - 3y + z &\geqslant 15 \\ +-x + 2y &\geqslant -1 \implies x - 2y \leqslant 1 +\end{cases} +\end{align*} +$$ +2. Zamieniamy układ w taki sposób, by wyeliminować wszystkie symbole mniejsze/większe niż. Żeby to zrobić musimy albo odjąć albo dodać jakąś liczbę. Jeżeli mam więcej niż to musimy odjąć i analogicznie. +$$ +\begin{align} +\begin{cases} +2x + y - 3z &= 12 \\ +5x - 3y + z - S_1 &= 15 \\ +-x + 2y + S_2 &= 1 +\end{cases} +\end{align} +$$ +3. Zapisujemy funkcję celu +$$ +f(x, y, z, S_1, S_2) = 7x + 2y -3z + 0S_1 + 0S_2 +$$ +4. Sprawdzamy czy mamy macierz jednostkową, czyli czy mamy kwadrat 3x3 w którym w każdej kolumnie jest jedna jedynka i jest ona w różnych wierszach. W tym przypadku pasuje nam tylko $S_2$, więc dodajemy zmienne sztuczne. +$$ +\begin{align} +\begin{cases} +2x + y -3z + a_1 &= 12 \\ +5x - 3y + z - S_1 + a_2 &= 15 \\ +x - 2y + S_2 &= 1 +\end{cases} +\end{align} +$$ + +Macierz jednostkowa +![](../../../media/Pasted%20image%2020260616200137.png) + +5. Dodajemy albo bardzo dużo zmiennych sztucznych albo bardzo mało. Jeśli szukamy max to muszą być ujemne i analogicznie. +$$ +f(x, y, z, S_1, S_2) = 7x + 2y -3z + 0S_1 + 0S_2 - 1000a_1 - 1000a_2 +$$ +6. Przechodzimy do uzupełniania tabelki +- 1. W $Cj\rightarrow$ wpisujemy ilość następnych zmiennych. +- 2.Wiersz pod tym wpisujemy co to są za zmienne. +- 3. Pod $B\downarrow$ (bazowa) wpisujemy wyrazy wolne po przekształceniach. +- 4. Pod $WB\downarrow$(wektory bazowe) wpisujemy nasze nowe sztuczne zmienne, które tworzą naszą macierz jednostkową. +- 5. Pod $CB\downarrow$ wpisujemy ile mamy tych zmiennych. +- 6. Uzupełniamy środek tabelki wartościami z układu równań dla poszczególnych zmiennych +- 7. $WW\rightarrow$(wiersz wskaźnikowy) uzupełniamy poprzez wymnożenie wartości z kolumny CB razy poszczególne kolumny zmiennych. Np. $-1000 \cdot 12 + (-1000) \cdot 15 + 0 \cdot 1 = -27000$ . Potem od wyniku odejmujemy wartość funkcji celu. Czyli np dla kolumny $x$ wykonujemy działanie $-7000 - 7 = -7007$ + +| | | Cj$\rightarrow$ | 7 | 2 | -3 | 0 | 0 | -1000 | -1000 | | +| -------------- | --------------- | --------------- | ----- | ---- | ---- | ----- | ----- | ----- | ----- | ------------------ | +| CB$\downarrow$ | WB $\downarrow$ | B $\downarrow$ | x | y | z | $S_1$ | $S_2$ | $a_1$ | $a_2$ | $\frac{B}{x}$ | +| -1000 | $a_1$ | 12 | 2 | 1 | -3 | 0 | 0 | 1 | 0 | $\frac{12}{2} = 6$ | +| -1000 | $a_2$ | 15 | 5 | -3 | 1 | -1 | 0 | 0 | 1 | $\frac{15}{5} = 3$ | +| 0 | $S_2$ | 1 | 1 | -2 | 0 | 0 | 1 | 0 | 0 | 1 | +| | $WW\rightarrow$ | -27000 | -7007 | 1998 | 2003 | 1000 | 0 | 0 | 0 | | + +- 8. Gdy już mamy taką tabelkę szukamy wartość zmiennej wiersza wskaźnikowego, która jest najmniejsza. W tym przypadku będzie to zmienna $x$. Ta kolumna będzie kluczowa. +- 9. Szukamy najmniejszej wartości w $\frac{B}{x}$. To będzie wiersz kluczowy. +- 10. Kolumna kluczowa wskazuje na wartość, która wchodzi do bazy, a wiersz kluczowy na wartość, która wylatuje. +- 11. Wartość na przecięciu się kolumny i wiersza kluczowego to pivot. Pivot musi być 1, a jeśli nie jest to mnożymy wiersz razy tyle, żeby był. +- 12. Używamy pivota, żeby zlikwidować liczbę nad nim. Ma być zerem! (Podobnie jak w metodzie eliminacji gaussa) + +- 13.(7) Ponownie liczymy wiersz wskaźnikowy. + +| | | Cj$\rightarrow$ | 7 | -2 | -3 | 0 | 0 | -1000 | -1000 | | +| -------------- | --------------- | --------------- | --- | ------ | ---- | ----- | ----- | ----- | ----- | ------------------ | +| CB$\downarrow$ | WB $\downarrow$ | B $\downarrow$ | x | y | z | $S_1$ | $S_2$ | $a_1$ | $a_2$ | $\frac{B}{y}$ | +| -1000 | $a_1$ | 10 | 0 | 5 | -3 | 0 | -2 | 1 | 0 | $\frac{10}{5} = 2$ | +| -1000 | $a_2$ | 10 | 0 | 7 | 1 | -1 | -5 | 0 | 1 | $\frac{10}{7}$ | +| 7 | $x$ | 1 | 1 | 2 | 0 | 0 | 1 | 0 | 0 | $\frac{1}{-2}$ | +| | $WW\rightarrow$ | -19993 | 0 | -12016 | 2003 | 1000 | 7007 | 0 | 0 | | +- 14.(8) szukamy najmniejszej wartości w WW. To jest kolumna kluczowa. Tutaj to będzie y +- 15. (9) wyliczamy $\frac{B}{x}$ i szukamy najmniejszej nieujemnej wartości. Wiersz kluczowy! +- 16. Powtarzamy te kroki, aż w WW nie będzie ujemnych wartości. Wtedy w kolumnie WB w WW będzie wartość optymalna funkcji celu. W kolumnie B znajdziemy wtedy wartości zmiennych funkcji, które wpływają na ten wynik. Gdybyśmy liczyli minimum to nie mogłoby być żadnych liczb dodatnich! + +**Na egzaminie nie trzeba rysować drugiej tabelki!** \ No newline at end of file