16 svar
253 visningar
naytte Online 5284 – Moderator
Postad: 9 jan 16:53 Redigerad: 9 jan 17:11

Har jag beräknat kroppens area rätt?

Halloj!

Jag sitter med följande uppgift:

Jag har ett försök men jag är osäker på om jag har gjort rätt och vill helst undvika att tjuvkika i facit. Så här har jag gjort:

Jag har ritat upp kurvan enligt följande skiss:

Jag har även skapat ett infinitesimalt tunnt rektangulärt område med bredd dy\mathrm{d}y. Min tanke är att om man roterar detta område runt yy-axeln så borde man få en typ av skiva. Denna skivan kommer ju ha höjden dy\mathrm{d}y och omkretsen 2πx2\pi x, så arean av skivans sida borde ju bli dA=2πxdy \mathrm{d} A = 2\pi x \mathrm{d}y.

Sedan tänker jag att man helt enkelt gör en oändlig summering av alla dessa infinitesimalt små ytor längs yy-axeln enligt:

A=2πy=3y=5xdy=[dy=2xdx]=4πx=0x=2x2dx=8π3·2  a.e.\displaystyle A= 2\pi\int_{y=3}^{y=5}x\mathrm{d}y=[\mathrm{d}y=2x\mathrm{d}x]=4\pi\int_{x=0}^{x=\sqrt2}x^2\mathrm{d}x=\frac{8\pi}{3}\cdot\sqrt2 \;\;\text{a.e.}

Är detta ett OK resonemang eller är det helt fel?


EDIT: vänta nu, är det kanske ds\mathrm{d}s man ska använda, alltså förändringen i kurvans längd? Det kanske spelar roll här eftersom vi tittar på area och inte volym.

Laguna Online 30830
Postad: 9 jan 19:09

Ja, ditt dA är i själva verket en volym. Arean beror på hur mycket kurvan lutar.

naytte Online 5284 – Moderator
Postad: 10 jan 08:35 Redigerad: 10 jan 08:35

Varför beror endast arean på lutningen men inte volymen? Råkar det bara vara så att approximeringen dy är tillräckligt bra i volymberäkningen att integralen konvergerar mot volymen, medan samma approximering inte fungerar här eftersom felet blir för stort?

Trinity2 Online 2091
Postad: 10 jan 09:04 Redigerad: 10 jan 09:15
naytte skrev:

Varför beror endast arean på lutningen men inte volymen? Råkar det bara vara så att approximeringen dy är tillräckligt bra i volymberäkningen att integralen konvergerar mot volymen, medan samma approximering inte fungerar här eftersom felet blir för stort?

Så här?

Jag kanske var otydlig i min fråga.

Det jag undrar är varför man måste räkna med längdelementet:

ds=dx1+dydx2\displaystyle \mathrm{d}s=\mathrm{d}x\sqrt{1+\left(\frac{\mathrm{d}y}{\mathrm{d}x}\right)^2}

då man beräknar areor, men inte behöver göra detta då man beräknar volymer. När man vill beräkna rotationsvolymer och delar in det aktuella området i rektanglar med bredd dx\mathrm{d}x eller dy\mathrm{d}y så är ju detta också egentligen en approximering som inte tar hänsyn till att kurvan böjs, men denna approximering ger ändå rätt svar i slutändan. Men en sådan approximering verkar ju inte fungera i fallet då man beräknar en area, utan där måste man ta hänsyn till att kurvan lutar. Varför är det så?

JohanF Online 5784 – Moderator
Postad: 10 jan 22:31 Redigerad: 10 jan 22:35

Jag har inget svar, men eftersom det känns som en väldigt fundamental grej att förstå så bumpar jag till den lite genom att ge ett annat exempel som jag tror kanske är relaterat.

Om vi utgår från din funktion f(x)=3+x2, varför kan man skriva arean mellan kurvan och x-axeln som

x=0x=2f(x)dx (summera rektanglar med bredden dx) 

men man kan ju uppenbarligen inte skriva längden på kurvan som

x=0x=2dx

?

Anto 329
Postad: 10 jan 23:00

Jag tror att felet inte blir försumbart. Förhållandet mellan dx och ds på ett kurvstycke är oförändrat då de blir mindre och mindre. Rita en kurva och rita ut ds och dx som en triangel ser du nog det. Det gäller dock inte felet som skapas av volymen. Förhållandet mellan felet och det verkliga värdet förändras och felet blir (relativt) mindre. Felet i detta fall är triangeln som du ritade med ds och dx!

JohanF Online 5784 – Moderator
Postad: 10 jan 23:39 Redigerad: 10 jan 23:40

Nä just det, felet blir inte försumbart. Men hur vet man när felet konvergerar mot noll och när det inte gör det?

Anto 329
Postad: 11 jan 00:00

Jag är inte säker men lite slarvigt kan man kanske säga att om felet går mot 0 snabbare än det man försöker approximera blir det försumbart? 

Jan Ragnar 1960
Postad: 11 jan 11:43

D4NIEL 3009
Postad: 11 jan 15:19 Redigerad: 11 jan 15:21
JohanF skrev:

Nä just det, felet blir inte försumbart. Men hur vet man när felet konvergerar mot noll och när det inte gör det?

Ett sätt är att uppskatta felet med en Taylorutveckling.

För en strängt växande funktion är arean av en infinitesimal stapel större än a1=f(x)Δxa_1=f(x)\Delta x och mindre än a2=f(x+Δx)Δxa_2=f(x+\Delta x)\Delta x. Uttrycket för a2a_2 kan vi taylorutveckla:

a2=(f(x)+f'(x)Δx+)Δxa_2=(f(x)+f^\prime(x)\Delta x+\dots)\Delta x

Nu ser vi att det största möjliga felet a2-a1f'(Δx)2+O(Δx)3a_2-a_1\leq f^\prime(\Delta x)^2+O(\Delta x)^3

Låter vi Δx\Delta x bli väldigt liten och använder snälla funktioner kommer bidraget från stapeln själv, dvs f(x)Δxf(x)\Delta x alltså vara mycket större än det maximala felet som är av storleken f'(x)(Δx)2f^\prime(x)(\Delta x)^2.

På samma sätt skulle varje skivas volym vara mycket större än felet vid en volymberäkning enligt skivmetoden. Men när vi räknar areaelement på det sätt naytte ville göra i den ursprungliga uppgiften visar det sig att felet och bidraget är av samma storleksordning. Vi behöver hitta en mer precis metod för att bestämma ett lämpligt areaelement. 

D4NIEL skrev:
JohanF skrev:

Nä just det, felet blir inte försumbart. Men hur vet man när felet konvergerar mot noll och när det inte gör det?

Ett sätt är att uppskatta felet med en Taylorutveckling.

För en strängt växande funktion är arean av en infinitesimal stapel större än a1=f(x)Δxa_1=f(x)\Delta x och mindre än a2=f(x+Δx)Δxa_2=f(x+\Delta x)\Delta x. Uttrycket för a2a_2 kan vi taylorutveckla:

a2=(f(x)+f'(x)Δx+)Δxa_2=(f(x)+f^\prime(x)\Delta x+\dots)\Delta x

Nu ser vi att det största möjliga felet a2-a1f'(Δx)2+O(Δx)3a_2-a_1\leq f^\prime(\Delta x)^2+O(\Delta x)^3

Låter vi Δx\Delta x bli väldigt liten och använder snälla funktioner kommer bidraget från stapeln själv, dvs f(x)Δxf(x)\Delta x alltså vara mycket större än det maximala felet som är av storleken f'(x)(Δx)2f^\prime(x)(\Delta x)^2.

På samma sätt skulle varje skivas volym vara mycket större än felet vid en volymberäkning enligt skivmetoden. Men när vi räknar areaelement på det sätt naytte ville göra i den ursprungliga uppgiften visar det sig att felet och bidraget är av samma storleksordning. Vi behöver hitta en mer precis metod för att bestämma ett lämpligt areaelement. 

Tack D4NiEL! (det smög sig in ett litet typo där eller hur?  Största möjliga fel a2-a1f'(x)(x)2+O(x)3)

Det som jag upplever lite frustrerande med det som naytte's exempel visar, är att den analys som du just gjorde verkar ska tas som som en självklarhet och trivialitet. I matematik på universitetsnivå brukar inte litteraturen nöja sig med att presentera lösningsmetoder likt "recept", utan att samtidigt argumentera för varför just den och ingen annan metod är riktig.

naytte Online 5284 – Moderator
Postad: 11 jan 23:39 Redigerad: 11 jan 23:43

Jag gillar resonemanget med Taylorutvecklingar mycket! Era kommentarer fick mig att fundera lite (särskilt din, JohanF) och just angående varför integralen kan uttryckas som den gör kan man nog resonera så här också (lite heuristiskt):

Låt ff vara en icke-avtagande funktion. I så fall har vi, som D4NIEL påpekade, att:

fx-dxdxfxdxfx+dxdx\displaystyle f\left(x-\mathrm{d}x\right)\mathrm{d}x\le f\left(x\right)\mathrm{d}x\le f\left(x+\mathrm{d}x\right)\mathrm{d}x, för någon infinitesimal dx\mathrm{d}x

Om vi använder Taylorutvecklingar här igen har vi alltså:

(fx-f'xdx+12f''xdx2-...)dxfxdx(fx+f'xdx+12f''xdx2+...)dx\displaystyle (f\left(x\right)-f'\left(x\right)\mathrm{d}x+\frac{1}{2}f''\left(x\right)\mathrm{d}x^2-...)\mathrm{d}x \le f\left(x\right)\mathrm{d}x\le(f\left(x\right)+f'\left(x\right)\mathrm{d}x+\frac{1}{2}f''\left(x\right)\mathrm{d}x^2+...)\mathrm{d}x

Här kan vi alltså stänga in den sökta integralen av fxdxf\left(x\right)\mathrm{d}x mellan integralerna av båda våra taylorserier. Men eftersom alla derivator av grad nn multipliceras med en infinitesimal av grad n+1n+1 (förutom termen fxf\left(x\right)) måste det betyda att både VL och HL blir den sökta integralen vid integrering.

JohanF Online 5784 – Moderator
Postad: 12 jan 12:58 Redigerad: 12 jan 12:59

Jag tänker ungefär lika som du.

(Puh! Så skönt att slippa behöva börja tvivla på att integralkalkylens fundamentalsats verkligen stämmer 😃)

naytte Online 5284 – Moderator
Postad: 12 jan 15:03 Redigerad: 12 jan 15:03

Men jag undrar hur Leibniz kom fram till det här. Han kan ju knappast ha använt taylorpolynom, för dessa fanns ju inte när han uppfann integralen.

D4NIEL 3009
Postad: 12 jan 15:11 Redigerad: 12 jan 15:19

Jag kan hålla med om att grundkurserna i många fall hoppar över integrationsmetoder i en dimension, ofta får det räcka med det man lärde sig i gymnasiet. Detta kan delvis förklaras av att kursplanen måste prioritera ett brett spektrum av viktiga moment på väldigt kort tid, särskilt inom civilingenjörsprogrammen. Samtidigt verkar tanken vara att dessa frågor främst ses som en del av flerdimensionella problem, där lösningsmetoderna behandlas mer ingående i senare kurser.

Apropå flervarre tänkte jag visa hur man löser ytintegralen med just sådana metoder. Det första vi behöver göra är att parametrisera ytan. Eftersom vi snurrar runt y-axeln är problemet symmetriskt i radiell led (i xzxz-planet). Vidare är höjden över xzxz-planet och radien sammankopplade, för en given radie är höjden över xzxz-planet h=3+r2h=3+r^2. Om vi lånar standarduttrycket för cylinderkoordinater och skyfflar om så att yy är "höjden" får vi

S:r(r,θ)=rcos(θ),rsin(θ),3+r2,  r[0,2],θ[0,2π]S:\,\, \mathbf{r}(r,\theta)=\left(r\cos(\theta), r\sin(\theta), 3+r^2\right),\quad r\in[0,\sqrt2],\, \theta \in [0,2\pi]

Kryssprodukten mellan två vektorer ger oss arean av en parallellogram, och som man visar i flervarren kan vi nu bilda ett infinitesimalt (vektoriellt) ytelement enligt

dS=rr×rθdrdθ=(ry^-2r2r^)drdθd\mathbf{S}=\left(\frac{\partial\mathbf{r}}{\partial r}\times \frac{\partial\mathbf{r}}{\partial \theta}\right)\,\mathrm{d}r \mathrm{d}\theta=(r\hat{y}-2r^2\hat{r})\,\mathrm{d}r\mathrm{d}\theta

Arean bestäms genom att integrera absolutbeloppet av dSd\mathbf{S}:

sdS=θ=02πr=02r2+4r4drdθ=13π3\displaystyle \int_s \left|\mathrm{d}\mathbf{S}\right|=\int_{\theta=0}^{2\pi} \int_{r=0}^{\sqrt{2}} \sqrt{r^2+4r^4}\,\mathrm{d}r\mathrm{d}\theta=\frac{13\pi}{3}

JohanF Online 5784 – Moderator
Postad: 12 jan 20:11 Redigerad: 12 jan 20:30

Jag blev lite nyfiken på detta.

Några övningar kring volym och mantelarea på en sfär (V=4πr33, A=4πr2)

Volymen av rotationskroppen kring x-axeln av funktionen f(x)=r2-x2 i intervallet -r<x<r,

med volymdifferentialen dV=πf(x)2dx=πr2-x22dx=πr2-x2dx (dvs med tjockleken dx på "skivan")

V=dV=-rrπr2-x2dx=πr2x-x33r-r=4πr33

Mantelarean av rotationskroppen kring x-axeln, med ytdifferentialen 

dA=2πf(x)1+f'(x)2dx=2πr2-x21+-xr2-x22dx=2πr2-x2rr2-x2dx=2πrdx (dvs med bredden dl på "skivan")

A=dA=-rr2πrdx=2πrxr-r=4πr2

"Mantelarean" av rotationskroppen kring x-axeln, med ytdifferentialen 

dAfel=2πf(x)dx=2πr2-x2dx (dvs med bredden dx på "skivan")

Afel=dAfel=-rr2πr2-x2dx=2π-rrr2-x2dx=2π·(arean av en halcirkel med radien r)=2π·πr22=π2r2


Tillägg: 12 jan 2025 20:27

(Man skulle också kunna lösa sista integralen med en tabellerad primitiv funktion, tex nr30 https://www.integral-table.com/)

Svara
Close