MATLAB 5.0 MAT-file, Platform: PCWIN, Created on: Tue Feb 17 15:07:10 2004                                                  IM0jn°ykomy8ipndp
Bsč<ysŔGŽ% ******************* MATEMATICKÝ SEMINÁŘ č. 3 ************************h8% 3.1 SUBMATICE A INDEXOVÁNÍ4h% 3.2 SUBMATICE NA LEVÉ STRANĚ PŘIŘAZOVACÍHO PŘÍKAZU0%ĐP % Copyright (C) 1996-2003 by V. Veselý (Masaryk University Brno, Czech Republic)0%.\% Dávka je kompletně upravena pro menu-systém.0@echo on@more onPformat compact8clcŔGŽ%======================================================================h8% 3.1 SUBMATICE A INDEXOVÁNÍŔGŽ%----------------------------------------------------------------------ŔH%% v(r)   ... subvektor vektoru v určený včetně pořadí vektorem indexů rh6%% A(r,s) ... submatice v A¸D%%            r=vektor řádkových indexů, s=vektor sloupcových indexůŔFŚ%%            za r nebo s lze dosadit také indikátorový vektor z 0 a 1€'N%% A(r,:) ... řádky r a všechny sloupce€'N%% A(:,s) ... všechny řádky a sloupce sH%% A(:,:)=A*T  A=(1:10)'*(1:8), %  tedy máme A(r,s)=r*s1b% vyzkoušíme vytvořit novou matici pomocí indexů:`0  S=A([1,3,5,7],[3,7,8])@pause;H% nebo takéX(  S=A(1:2:7,[3,7,8])@pause;H% nebo takéŔEŠ  S=A(logical([1 0 1 0 1 0 1 0 0 0]),logical([0 0 1 0 0 0 1 1]))     €%J% speciálně r nebo s může být skalár:¸C†  A(3,7),          %  jeden prvek matice A ve 3. řádku a 7. sloupcix"D  r4=A(4,:)        %  4. řádek   Ax"D  s6=A(:,6)        %  6. sloupec AČL  Aperm=A(10:-1:1,:) % Řádkově permutovaná matice (řádky v obráceném pořadí)H	keyboard;0ŔH%% A(:)   ... sloupce matice A naskládané pořadě pod sebe do sl. vektoru°>|%%            (matematická operace obvykle značená jako vec A)4h%%            odpovídá uložení prvku matice v paměti€&L%% A(r)   ... totéž jako: v=A(:); v(r)ČJ”%%            A řád./sl. vektor => výběr řád./sl. subvektoru na pozicích rX&%% A=reshape(v,m,n)ČK–%%   reshape(v,[m,n]) ... opak k A(:) - z řád. nebo sl. vektoru v délky m*nŔEŠ%%                        vytvoří matici A typu m x n, pro níž v=A(:)p:  u=1:30 % vytvoříme vektor uŔGŽ  v=reshape(u,length(u),1)  % transponuje řádkový vektor u na sloupcový@pause; 8p  u(1:2:length(u))          % subvektor s lichými indexy 8p  u(2:2:length(u))          % subvektor se sudými indexy@pause;€%J% vektor v přeuspořádáme do matice   `*  A1=reshape(v,2,15) @% a zpět@  A1(:)@pause;`,  A2=reshape(v,[3,10])`0  A2=reshape(v,size(A2))@  A2(4)@pause;X$  A=reshape(u,5,6)H
  A([4,8])H	keyboard;0ŔGŽ%======================================================================4h% 3.2 SUBMATICE NA LEVÉ STRANĚ PŘIŘAZOVACÍHO PŘÍKAZUŔGŽ%----------------------------------------------------------------------¸A‚%% v(r)=b   ... přepis vektoru v vektorem b na pozicích indexu r ŔEŠ%% A(r,s)=B ... Submatice A na pozicích určených vektory indexu r a sČK–%%              se přepíše prvky matice B (B je typu length(r) x length(s))ŔH%%              Jestliže některé pozice prvků jdou mimo stávající rozměr¸B„%%              matice A, pak A se automaticky zvětší na potřebnou)R%%              velikost a doplní nulami.P
%% Analogicky 8p%% A(r,:)=B ... přepis na řádcích r a ve všech sloupcích€(P%%              B=[] => odstraní řádky r 8p%% A(:,s)=B ... přepis ve všech řádcích a ve sloupcích s*T%%              B=[] => odstraní sloupce sŔFŚ%% A(r)=b   ... přepis prvku na libovolných pozicích v uspořádání A(:)ĐMš%%              A(r)=b, b neprázdná <=> v=A(:); v(r)=b; A=reshape(v,size(A));¨<x%%              A(r)=[]             <=> v=A(:); v(r)=[]; A=v0@  u=1:30X$  A=reshape(u,5,6)x!B  A(5,:)=[]    % Smažeme 5. řádekp:  A(:,6)=[]    % a 6. sloupec@pause;h2  A([1,3],[2,3])=zeros(2)°?~% nebo také jednodušeji lze uvést na pravé straně jen skalár 0:X$  A([1,3],[2,3])=0@pause;H  A(5,6)=-1@pause;°@€  A(logical([1 0 1 0 0]),logical([0 1 1 0 0 0]))=-ones(2)       €%J% speciálně r nebo s může být skalár:ČL  A(3,5)=1000,             %  přepis prvku matice A ve 3. řádku a 5. sloupci@pause;3f  A(4,:)=100*ones(1,size(A,2)) %  přepis 4. řádku A@pause;ČJ”  A(:,6)=100                   %  přepis 6. sloupce A jednodušeji skalárem@pause;0ŔGŽ%----------------------------------------------------------------------1b% Příklad na přepis prvků na libovolných pozicích0`% Chceme obrátit znaménka prvků A(2,1) a A(1,2):X(  A=reshape(1:6,2,3)p<  P(2,1)=1, P(1,2)=1, P(1,3)=0°?~% Předchozí 3 příkazy třikrát realokují matici na větší rozměr,`0% proto lepší postup je:ŔFŚ  P=zeros(size(A)); P(2,1)=1, P(1,2)=1, P=logical(P)                  ¸A‚  A(P(:))=-A(P(:))      % Obrácení znaménka prvků A(2,1) a A(1,2)H	keyboard;0ŔGŽ%================================= KONEC ==============================0%0%0%%@more off@echo off0˙˙0cpk