Afsnit ML.5: Gentagne beregninger, betingede beregninger

Indholdet i dette afsnit skal ikke bruges i kurset Introducerende Statistik og Dataanalyse med MATLAB, men er medtaget for fuldstændighedens skyld.
Ofte kommer man ud for, at man ønsker at gentage et stykke kode mange gange. Dette kan gøres ved en såkaldt for-løkke. Denne starter med for og angivelse af, at en variabel skal løbe i en bestemt mængde, og slutter med end.
Følgende program finder for et givet heltal hvilke tal mellem 1 og , der går op i tallet. Dette gøres ved en for-løkke med et index , der løber fra 1 til , og i hvert trin udregnes resten ved division af med . Efter løkken angives de , for hvilke resten er nul. Undervejs bruges funktionen floor, der giver det største hele tal mindre end eller lig med input til funktionen.

>> n=12;

>> divisor=[1:n];

>> rest=zeros(1,n);

>> for j=1:n
 rest(j)=n-floor(n/j)*j;
 end

>> divisor(rest==0)
 ans =     1     2     3     4     6    12

En for-løkke er velegnet, når man på forhånd ved, hvor mange gange man skal løbe gennem løkken. Hvis man ikke ved dette, og det afgøres af de beregninger, der laves undervejs, bruges en while-løkke. Denne starter med while og et sandshedsudsagn. Så længde sandshedsudsagnet er sandt, løber man gennem løkken (her skal man altså passe på ikke at lave en løkke, der aldrig stopper).
Det følgende program beregner alle fibonacci tal indtil første gang, man får et tal, der er større end eller lig med 1000.

>> fib=[1,1];

>> while fib(end)<1000
  nye=fib(end)+fib(end-1);
  fib=[fib,nye]; % vektoren fib forlænges med et tal
  end

>> fib
 fib =
   Columns 1 through 9
     1   1   2   3   5   8   13   21   34
   Columns 10 through 17
     55  89  144  233  377  610  987  1597

En while-løkke bruger et sandhedsudsagn til at bestemme, om løkken skal fortsætte. Mere generelt kan man kontrollere, at et stykke kode udføres, alt efter om en betingelse er opfyldt ved hjælp af en if-konstruktion. En if-konstruktion er enten af typen if ... end eller af typen if ... else ... end. I den første konstruktion udføres en beregning, alt efter om en betingelse er opfyldt. I den anden konstruktion udføres en beregning, hvis en betingelse er opfyldt, og en anden beregning udføres, hvis betingelsen ikke er opfyldt.
I den første del af koden nedenfor findes alle lige tal mellem 1 og 10, og i den anden del findes separat de lige og de ulige tal.

% første del

>> lige=[];

>> for j=1:10
  rest=j-floor(j/2)*2;
  if rest==0
   lige=[lige,j];
  end
 end

>> lige
 lige =     2     4     6     8    10

% anden del

>> lige=[];

>> ulige=[];

>> for j=1:10
  rest=j-floor(j/2)*2;
  if rest==0
   lige=[lige,j];
  else
   ulige=[ulige,j];
  end
 end

>> lige
 lige =     2     4     6     8    10

>> ulige
 ulige =     1     3     5     7     9

ForegåendeNæste