Simo K. Kivelä / 09.05.2003

Ympyrä ja pyörähdyspinta

•Tehtävä 1:

Ympyrän keskipiste on origo ja säde on 2. Muodosta ympyrälle parametriesitys ja piirrä ympyrä tämän avulla.

•Ratkaisu 1:

In[1]:=

Remove["Global`*"]

Trigonometristen funktioiden määritelmistä seuraa, että origoon asetettua suuntakulmaa t vastaava piste (2 cos(t), 2 sin(t)) on etäisyydellä 2 origosta, ts. se
sijaitsee halutulla ympyrällä. Jos toisaalta t saa kaikki arvot väliltä [0, 2 π[, saadaan kaikki ympyrän pisteet.

Mallikuvio voidaan piirtää Mathematicalla:

In[2]:=

kuvio = Graphics[{Line[{{-2.5, 0}, {2.5, 0}}], Line[{{0, -2.5}, {0, 2.5}}], Line[{{0, 0}, {2 C ... }, {-1, 1.5}], Text["t", {Cos[1/2]/3, Sin[1/2]/3}]}] ; Show[kuvio, AspectRatio -> 1]

[Graphics:HTMLFiles/pyorpinta_7.gif]

Out[3]=

-Graphics -

Ympyrän parametriesitykseksi siis kelpaa (2 cos(t), 2 sin(t)) ja käyrä voidaan tähän perustuen piirtää:FormBox[Cell[], TraditionalForm]

In[4]:=

ParametricPlot[2 {Cos[t], Sin[t]}, {t, 0, 2 Pi}]

[Graphics:HTMLFiles/pyorpinta_12.gif]

Out[4]=

-Graphics -

Kuvio sinänsä on oikea, mutta akseleiden yksiköt ovat eri pitkät, minkä takia tulos näyttää ellipsiltä. Asia voidaan korjata:

In[5]:=

ParametricPlot[2 {Cos[t], Sin[t]}, {t, 0, 2 Pi}, AspectRatio -> Automatic]

[Graphics:HTMLFiles/pyorpinta_15.gif]

Out[5]=

-Graphics -

•Tehtävä 2:

Käyrä y = sin(x)^2, 0 <= x <= π, pyörähtää x-akselin ympäri, jolloin syntyy pyörähdyspinta. Laske pinnan ala, sen rajaaman alueen tilavuus ja piirrä pinta.

•Ratkaisu 2:

In[6]:=

Remove["Global`*"]

Olkoon

In[7]:=

y = Sin[x]^2

Out[7]=

Sin[x]^2

Tilavuus saadaan integroimalla:

In[8]:=

v = Integrate[Pi y^2, {x, 0, Pi}]

Out[8]=

(3 π^2)/8

Huomaa: Sen enempää tilavuuden laskeminen kuin seuraava pinnan alan laskeminenkaan ei edellytä kaavojen ulkoa osaamista, vaan niiden ymmärtämistä. Tämän jälkeen kaavan mieleen palauttaminen ei ole minkäänlainen ongelma!

Pinta-alan integrointi:

In[9]:=

integrandi = 2 Pi y Sqrt[1 + D[y, x]^2]

Out[9]=

2 π Sin[x]^2 (1 + 4 Cos[x]^2 Sin[x]^2)^(1/2)

In[10]:=

a = Integrate[integrandi, {x, 0, Pi}]

Out[10]=

2 π EllipticE[-1]

Tulosta ei voida esittää alkeisfunktioiden avulla, mutta kylläkin Mathematican tunteman erikoisfunktion avulla. Tälle voidaan laskea numeerinen arvo tavalliseen tapaan:

In[11]:=

a // N

Out[11]=

12.001505309269428`

Vaihtoehtona on integroida suoraan numeerisesti:

In[12]:=

NIntegrate[integrandi, {x, 0, Pi}]

Out[12]=

12.001505309255592`

Pinta voidaan parametrisoida sen perusteella, että sen leikkaus x-akselia vastaan kohtisuoran tason kanssa on ympyrä. Tälle puolestaan voidaan käyttää edellä muodostettua parametriesitystä:

In[13]:=

pinta = {x, y Cos[t], y Sin[t]}

Out[13]=

{x, Cos[t] Sin[x]^2, Sin[t] Sin[x]^2}

In[14]:=

ParametricPlot3D[pinta, {x, 0, Pi}, {t, 0, 2 Pi}]

ParametricPlot3D :: ppcom :  Function  pinta  cannot be compiled; plotting will proceed with the uncompiled function.

[Graphics:HTMLFiles/pyorpinta_36.gif]

Out[14]=

-Graphics3D -


Converted by Mathematica  (May 9, 2003)