opgaver:Uge1mandag
Pia Jensen (Talk | contribs) |
Pia Jensen (Talk | contribs) |
||
| Line 126: | Line 126: | ||
# <code>0, 1/2, 2/3, 3/4, 4/5, ... </code> | # <code>0, 1/2, 2/3, 3/4, 4/5, ... </code> | ||
{{hidden begin|toggle=right|title=Løsning|titlestyle=background:#ccccff|bg2=#eeeeee}} | {{hidden begin|toggle=right|title=Løsning|titlestyle=background:#ccccff|bg2=#eeeeee}} | ||
| + | Igen skal man gøre brug af kolon operatoren, mens man i de sidste to løser opgaven nemmest ved at lave flere lister og kombinere dem: | ||
| + | # <code>2:2:10</code> | ||
| + | # <code>10:-2:-4</code> | ||
| + | # <code>x = 1:10; 1./x</code> | ||
| + | # <code>x = 0:10; y = 1:11; x./y</code> | ||
| + | |||
| + | I de uendelige lister vælger man selvfølgelig bare selv sit maksimum. Det er dog vigtigt i den sidste af opgaverne at man laver de to arrays lige lange - ellers er det ikke muligt at dividere dem med hinanden element-vist. | ||
{{hidden end}} | {{hidden end}} | ||
| Line 135: | Line 142: | ||
Summér elementerne i versionen af denne vektor med 100 elementer. | Summér elementerne i versionen af denne vektor med 100 elementer. | ||
| − | {{hidden begin|toggle=right|title= | + | {{hidden begin|toggle=right|title=Hint|titlestyle=background:#ccccff|bg2=#eeeeee}} |
| + | Start med at lave en liste <code>n</code>, som du så kan udregne <code>x</code> ved hjælp af. Hvis du vil have 100 elementer laver du f.eks. | ||
| + | |||
| + | :<code> n = 1:100</code> | ||
| + | Og husk alle de nødvendige punktummer i udtrykket når du udregner <code>x</code>. | ||
| + | {{hidden end}} | ||
| + | {{hidden begin|toggle=right|title=Løsning|titlestyle=background:#ccccff|bg2=#eeeeee}} | ||
| + | Opgaven kan løses ved de følgende tre linier: | ||
| + | :<code>n = 1:100</code> | ||
| + | :<code>x = ( (-1).^(n+1) ) ./ (2*n - 1)</code> | ||
| + | :<code>sum(x)</code> | ||
{{hidden end}} | {{hidden end}} | ||
| Line 148: | Line 165: | ||
# cot(t) | # cot(t) | ||
# sec<sup>2</sup>(t) + cot(t) - 1 | # sec<sup>2</sup>(t) + cot(t) - 1 | ||
| + | {{hidden begin|toggle=right|title=Hint|titlestyle=background:#ccccff|bg2=#eeeeee}} | ||
| + | Alle de typiske matematiske funktioner som cos, sin og exp er allerede standardfunktioner i MATLAB, og de opererer element-vist. Husk dine punktummer når de er nødvendige. | ||
| + | {{hidden end}} | ||
{{hidden begin|toggle=right|title=Løsning|titlestyle=background:#ccccff|bg2=#eeeeee}} | {{hidden begin|toggle=right|title=Løsning|titlestyle=background:#ccccff|bg2=#eeeeee}} | ||
| + | Udtrykkene er stort set som man forventer de vil være, når bare man husker sine punktummer hvor de er nødvendige. Bemærk at inverse trigonometriske funktioner starter med <code>a</code>, som f.eks. <code>acos</code> eller <code>atan</code>. | ||
| + | #<code>log(2+t+t.^2)</code> | ||
| + | #<code>exp(t).*(1 + cos(3*t))</code> | ||
| + | #<code>cos(t).^2 + sin(t).^2</code> | ||
| + | #<code>atan(t)</code> | ||
| + | #<code>cot(t)</code> | ||
| + | #<code>sec(t).^2 + cot(t) - 1</code> | ||
{{hidden end}} | {{hidden end}} | ||
Revision as of 12:23, 5 March 2012
Disse simple øvelser skal laves som en warm-up til MATLAB. Dem der allerede er trygge ved at bruge MATLAB kan nøjes med at skimme dem igennem - men det anbefales at I går igennem dem alle for at være sikre på at I ved hvad I laver!
Opgaverne er taget fra [1] (Basic syntax and command-line exercises 1-6 og 8, samt Basic array exercises 1-3).
Contents |
Simpel syntaks og kommando-linie
Opgave 1
Lav en vektor bestående af de lige heltal mellem 31 og 75.
Gå ind i MATLAB hjælp og søg på colon, og læs så
Her vil du se at man kan lave en liste ved hjælp af koden j:i:k. Denne liste vil starte i j og lægge i til for hvert element, indtil den kommer til k. Hvis den ikke rammer k præcist, så vil den stoppe før den kommer over.
Dette gøres med kolon-operatoren, der kan lave lister:
x = 32:2:75
Med denne kode starter listen ved 32 (det første af de lige tal vi gerne vil have), og stiger 2 hver gang, indtil den når til 75. Den når faktisk til 74, og ser så at 76 er over 75, og stopper derfor ved 74.
Opgave 2
Lad x = [2 5 1 6].
- Addér 16 til hvert element.
- Addér 3 til alle ulige-indeks elementerne.
- Udregn kvadratroden af hvert element.
- Udregn kvadratet af hvert element.
For at løse disse opgaver kan det være en god idé at se på Getting Started with MATLAB bogen, der ligger på Absalon. Læs om arrays fra side 3-21. Vær sikker på at du forstår hvornår man skal bruge punktum foran matematiske operatorer.
1. For at lægge et tal til alle elementerne i et array skal man simpelthen skrive
x + 16
2. For at lægge et tal til kun nogle bestemte elementer i et array skal man lægge det til array'et med indekserne man gerne vil have - skrevet her som en liste:
x(1:2:end) + 3
3. Kvadratroden er en funktion i MATLAB, og funktioner der fungerer på skalarer fungerer automatisk element-vist. Derfor fås kvadratroden af hvert element i x ved
sqrt(x)
4. Kvadratet af en vektor er ikke det samme som kvadratet af alle elementerne hver for sig, og derfor skal man huske et punktum her:
x.^2
Man kunne også vælde at bruge
x.*x
Opgave 3
Lad x = [3 2 6 8]' og y = [4 1 3 5]'. (NB. x og y bør være søjlevektorer).
- Addér summen af elementerne i
xtily. - Hæv hvert element af
xi potensen specificeret af det tilsvarende element iy. - Dividér hvert element i
ymed det tilsvarende element ix. - Gang hvert element i
xmed det tilsvarende element iy, og gem resultatet i den ny variabelz. - Addér elementerne i
zog gem resultatet i den ny variabelw. - Udregn
x'*y - wog fortolk resultatet.
Det kan være en god idé at se på MATLAB funktionen sum. Find selv information om den ved at åbne MATLAB hjælp og søge på den. Her skal du også være sikker på at du forstår hvornår man skal bruge punktum foran matematiske operatorer.
1. Her skal man bruge funktionen sum, der summer over alle elementerne i den vektor den bliver fodret med,
y + sum(x)
2. Da vektorerne x og y har samme længde, er det bare at skrive (husk at bruge et punktum, så MATLAB forstår at du vil gøre det element-vist)
x.^y
3. Tilsvarende divideres der simpelthen ved
y./x
4. Den nye variabel defineres ved
z = x.*y
5. Igen bruges sum funktionen, og
w = sum(z)
6. Hvis alt er blevet gjort korrekt, skulle man gerne se at x'*y - w er lig med nul. Udtrykket x'*y er et normalt prikprodukt mellem en rækkevektor (x der er blevet transponeret med '-operatoren) og en søjlevektor. Derimod er w det manuelt udregnede prikprodukt.
Opgave 4
Udregn de følgende MATLAB udtryk i hånden, og brug så MATLAB til at tjekke resultaterne:
-
2 / 2 * 3 -
6 - 2 / 5 + 7 ^ 2 - 1 -
10 / 2 \ 5 - 3 + 2 * 4 -
3 ^ 2 / 4 -
3 ^ 2 ^ 2 -
2 + round(6 / 9 + 3 * 2) / 2 - 3 -
2 + floor(6 / 9 + 3 * 2) / 2 - 3 -
2 + ceil(6 / 9 + 3 * 2) / 2 - 3
Det kan hjælpe at skrive udtrykkene ind i MATLAB og tilføje parenteser så de bliver nemmere at læse. Hvis resultatet stadig er det samme efter parenteserne er blevet sat, så har du højest sandsynligt ikke misforstået noget.
For at vise hvad udtrykkene betyder, er de her omskrevet med ekstra parenteser:
2 / 2 * 3 = (2/2) * 36 - 2 / 5 + 7 ^ 2 - 1 = 6 - (2/5) + (7^2) - 110 / 2 \ 5 - 3 + 2 * 4 = (10/2)/5 - 3 + (2*4)3 ^ 2 / 4 = (3^2) / 43 ^ 2 ^ 2 = (3^2) ^ 22 + round(6 / 9 + 3 * 2) / 2 - 3 = 2 + (round((6/9) + (3*2))/2) - 32 + floor(6 / 9 + 3 * 2) / 2 - 3 = 2 + (floor((6/9) + (3*2))/2) - 32 + ceil(6 / 9 + 3 * 2) / 2 - 3 = 2 + (ceil((6/9) + (3*2))/2) - 3
Bemærk at round() funktionen runder op/ned automatisk, mens floor() altid runder ned og ceil() altid runder op.
Opgave 5
Lav en vektor x med elementerne...
-
2, 4, 6, 8, ... -
10, 8, 6, 4, 2, 0, -2, -4 -
1, 1/2, 1/3, 1/4, 1/5, ... -
0, 1/2, 2/3, 3/4, 4/5, ...
Igen skal man gøre brug af kolon operatoren, mens man i de sidste to løser opgaven nemmest ved at lave flere lister og kombinere dem:
-
2:2:10 -
10:-2:-4 -
x = 1:10; 1./x -
x = 0:10; y = 1:11; x./y
I de uendelige lister vælger man selvfølgelig bare selv sit maksimum. Det er dog vigtigt i den sidste af opgaverne at man laver de to arrays lige lange - ellers er det ikke muligt at dividere dem med hinanden element-vist.
Opgave 6
Lav en vektor x med elementerne
- xn = (-1)n+1/(2n-1)
Summér elementerne i versionen af denne vektor med 100 elementer.
Start med at lave en liste n, som du så kan udregne x ved hjælp af. Hvis du vil have 100 elementer laver du f.eks.
n = 1:100
Og husk alle de nødvendige punktummer i udtrykket når du udregner x.
Opgaven kan løses ved de følgende tre linier:
n = 1:100x = ( (-1).^(n+1) ) ./ (2*n - 1)sum(x)
Opgave 8
Givet en vektor, t, med længde n, lav MATLAB udtryk der vil udregne de følgende udtryk. Test at de virker for t = 1:0.2:2.
- ln(2 + t + t2)
- et(1 + cos(3t))
- cos2(t) + sin2(t)
- tan-1(1) (dette er den inverse tangentfunktion)
- cot(t)
- sec2(t) + cot(t) - 1
Alle de typiske matematiske funktioner som cos, sin og exp er allerede standardfunktioner i MATLAB, og de opererer element-vist. Husk dine punktummer når de er nødvendige.
Udtrykkene er stort set som man forventer de vil være, når bare man husker sine punktummer hvor de er nødvendige. Bemærk at inverse trigonometriske funktioner starter med a, som f.eks. acos eller atan.
log(2+t+t.^2)exp(t).*(1 + cos(3*t))cos(t).^2 + sin(t).^2atan(t)cot(t)sec(t).^2 + cot(t) - 1
Simple arrays
Opgave 1
Givet x = [3 1 5 7 9 2 6], forklar hvad de følgende kommandoer "betyder":
-
x(3) -
x(1:7) -
x(1:end) -
x(1:end-1) -
x(6:-2:1) -
x([1 6 2 1 1]) -
sum(x)
Opgave 2
Givet array'et A = [ 2 4 1 ; 6 7 2 ; 3 5 9], find kommandoerne der behøves for at...
- definere en vektor
x1der indholder den første række afA. - definere et array
yder indeholder de sidste 2 rækker afA. - udregne summen over kolonnerne i
A. - udregne summen over rækkerne i
A. - udregne standardafvigelsen af gennemsnittet for hver kolonne i
A(NB. standardafvigelsen af gennemsnittet er defineret som standardafvigelsen divideret med kvadratroden af antallet af elementer der bruges til at udregne gennemsnittet).
Opgave 3
Givet arrays'ene x = [1 4 8], y = [2 1 5] og A = [3 1 6 ; 5 2 7], find ud af hvilke udtryk der kan udføres korrekt, og giv resultaterne. Hvis udtrykket ikke kan udføres, så find ud af hvorfor. Kommandoen whos kan være smart at bruge her.
-
x + y -
x + A -
x' + y -
A - [x' y'] -
[x ; y'] -
[x ; y] -
A - 3
