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