Matematika/Gauso formulė

Iš testwiki.
18:36, 20 gruodžio 2021 versija, sukurta imported>Paraboloid (Pavyzdžiai)
(skirt) ← Ankstesnė versija | Dabartinė versija (skirt) | Vėlesnė versija → (skirt)
Pereiti į navigaciją Jump to search

Formulė Gauso - Ostrogradskio yra analogas Gryno formulės. Tada kai formulė Gryno - Ostrogradskio suriša kreivinį integralą antrojo tipo uždara kreive su dvilypiu integralu plokščia sritimi, apribota šia kreive, tai formulė Gauso - Ostragradskio nustato ryšį tarp paviršinio integralo (antrojjo tipo) uždaru paviršiumi ir trilypiu integralu palei erdvinę sritį, apribotą šiuo paviršiumi.


Formulė Gauso - Ostrogradskio:
V(Px+Qy+Rz)dxdydz=SPdydz+Qdzdx+Rdxdy.
Ji išreiškia paviršinį integralą bendro pavidalo palei išorinę pusę uždaro paviršiaus S per trilypį integralą palei trimatę sritį V, apribotą šiuo paviršiumi.
Formulę Gauso - Ostrogradksio galimą naudoti apskaičiavimui paviršinių integralų uždaru paviršiumi.


Vienam iš pritaikymų formulės Gauso - Ostrogradskio, paimkime uždavinį apie apskaičiavimą kūno tūrio su paviršiniu integralu, palei išorinę pusę paviršiaus, apribojantį šitą kūną.
Tikrai, jei sritis V turi nurodyta aukščiau formą, tai
P(x;y;z)=x,Q(x;y;z)=0,R(x;y;z)=0,
pagal formulę Gauso - Ostrogradksio randame:
Sxdydz=Vdxdydz=v.
Keisdami rolėmis x, y, z, gausime taip pat, kad
Sydzdx=v,Szdxdy=v.
tokiu budu yra formulės:
Sxdydz=v,Sydzdx=v,Szdxdy=v,
išreiškiančios kūno tūrį v per integralą palei išorinę pusę jo paviršiaus.
Paėmę funkciją P(x;y;z)=13x,Q(x;y;z)=13y,R(x;y;z)=13z,
gausime formulę, išreiškiančią kūno tūrį per paviršinį integralą bendro pavidalo:
v=13Sxdydz+ydzdx+zdxdy.
Nes tada Px+Qy+Rz=13+13+13=1, kai P=x3,Q=y3,R=z3 ir tada
V(Px+Qy+Rz)dxdydz=Vdxdydz=v.
Todėl kūno V, apriboto paviršiumi S tūris v lygūs:
v=Vdxdydz=13Sxdydz+ydzdx+zdxdy.


Pavyzdžiai

  • Pavyzdis. Apskaičiuoti integralą Sx3dydz+y3dzdx+z3dxdy pagal išorinę pusę sferos x2+y2+z2=R2.
Taikydami formulę Gauso - Ostragradksio, gauname:
Sx3dydz+y3dzdx+z3dxdy=V(3x2+3y2+3z2)dxdydz=302πdϕπ2π2cosθdθ0Rρ4dρ=
=302πdϕπ2π2cosθdθρ55|0R=3R5502πdϕπ2π2cosθdθ=3R5502πdϕsinθ|π2π2=3R5502π(sinπ2sinπ2)dϕ=
=3R5502π(1(1))dϕ=6R55ϕ|02π=6R552π=12πR55.
Teisingiau skaičiuoti taip (x2+y2+z2=ρ2):
Sx3dydz+y3dzdx+z3dxdy=V(3x2+3y2+3z2)dxdydz=302πdϕ0πsinθdθ0Rρ2ρ2dρ=
=302πdϕ0πsinθdθρ55|0R=3R5502πdϕ0πsinθdθ=3R5502πdϕ(cosθ)|0π=
=3R5502πdϕ([11])=6R5502πdϕ=6R55(2π0)=12πR55.
Patikrinsime apskaičiuodami Sx3𝐝y𝐝z,Sy3𝐝z𝐝x ir Sz3𝐝x𝐝y sumą.
x2=R2y2z2,
x=R2y2z2;
xy=(R2y2z2)y=2y2R2y2z2=yR2y2z2;
xz=(R2y2z2)z=2z2R2y2z2=zR2y2z2.
Vx=Sx31+(xy)2+(xz)2𝐝y𝐝z=
=Sx31+(yR2y2z2)2+(zR2y2z2)2𝐝y𝐝z=
=Sx31+y2R2y2z2+z2R2y2z2𝐝y𝐝z=
=Sx3R2y2z2+y2+z2R2y2z2𝐝y𝐝z=
=Sx3R2R2y2z2𝐝y𝐝z=S(R2y2z2)3R2R2y2z2𝐝y𝐝z=S(R2y2z2)R𝐝y𝐝z=
=S(R2ρ2)Rρ𝐝ρ𝐝ϕ=02π(0R(R3ρRρ3)𝐝ρ)𝐝ϕ=
=02π(R3ρ22Rρ44)|0R𝐝ϕ=02π(R3R22RR44)𝐝ϕ=
=02π(R52R54)𝐝ϕ=02πR54𝐝ϕ=R542π=πR52.
Kadangi ir teigiama ir neigiama kryptimi reikia apskaičiuoti, tai VX=2Vx=2πR52=πR5. Ir kadangi funkcijų x3, y3 ir z3 laipsniai vienodi, tai
M=3VX=3πR5.
Gauta masė (skardinės sferos tankis priklauso nuo x reikšmės trečiame laipsnyje) nesutampa su Gauso formulės logika.
Update. Pagal oficialią teoriją taip ir neturi sutapti atsakymai ir kas ką tik buvo apskaičiuota (M=3VX=3πR5) neturi nieko bendro su Gauso formule.
Pastaba. Rutulio paviršiaus plotą įmanoma apskaičiuoti cilindinėse koordinatėse. Bandydami, gauname:
R2=x2+y2+z2,
z=R2x2y2;
zx=2x2R2x2y2=xR2x2y2;
zy=2y2R2x2y2=yR2x2y2;
1+(zx)2+(zy)2=1+(xR2x2y2)2+(yR2x2y2)2=
=1+x2R2x2y2+y2R2x2y2=R2x2y2+x2+y2R2x2y2=R2R2x2y2;
x2+y2=ρ2 cilindinėse ir polinėse koordinatėse;
internetinio integratoriaus 0RρR2R2ρ2dρ=(ρ2R2)R2R2ρ2|0R=RR2ρ2|0R;
Spav.=s1+(zx)2+(zy)2dxdy=sR2R2x2y2dxdy=02πdϕ0RρR2R2ρ2dρ=
=02πRR2ρ2|0Rdϕ=R02π(R2R2R202)dϕ=R02π(0R2)dϕ=R202πdϕ=2πR2.
Taigi, visas rutulio paviršiaus plotas susideda iš dviejų pusrutulių, todėl
Svisas=2Spav.=22πR2=4πR2.
Kitaip patikrinsime apskaičiuodami Sx3𝐝y𝐝z,Sy3𝐝z𝐝x ir Sz3𝐝x𝐝y sumą.
x2=R2y2z2,
x=R2y2z2.
Vx=Sx3𝐝y𝐝z=S(R2y2z2)3𝐝y𝐝z=S(R2ρ2)3ρ𝐝ρ𝐝ϕ=
=02π(0Rρ(R2ρ2)3𝐝ρ)𝐝ϕ.
Pasinaudodami internetiniu integratoriumi, gauname, kad
0Rρ(R2ρ2)3𝐝ρ=15((Rρ)(R+ρ))5/2|0R=
=15((RR)(R+R))5/2(15((R0)(R+0))5/2)=
=0+15((R0)(R+0))5/2=15(R2)5/2=R55;
Vx=02π(0Rρ(R2ρ2)3𝐝ρ)𝐝ϕ=02πR55𝐝ϕ=2πR55.
Kadangi reikia dviejų rutulio pusrutulių (teigiama ir neigiama Ox kryptimi), tai
VX=2Vx=22πR55=4πR55.
Kadangi, pagal sąlyga bus VX=VY=VZ, tai
V=3VX=34πR55=12πR55.
Ką norėta rasti ir kas rasta. Norėta apskaičiuoti (kaip supranta redaguotojas) sferos iš skardos masę. Skardos tankis vienu skaičiavimu kinta tik priklausomai nuo Ox koordinatės pagal funkcija γ(x)=x3. Antru atveju skardos tankis kinta tik priklausomai nuo y koordinatės pagal funkciją γ(y)=y3. Trečiu atveju skardos tankis kinta tik priklausomai nuo z koordinatės pagal funkciją γ(z)=z3. Kadangi rutulys simetriškas, tai užtenka apskaičiuoti, tarkime, Sz3𝐝x𝐝y ir padauginti iš 3, o paskui dar padauginti iš 2, nes ir teigiama ir neigiama kryptimi tankis didėja vienodai. Skaičiuojant analogiškai kreiviniam integralui (pirmojo tipo) gauname atsakyma M=6πR3. Atsakymas M=6πR3 ir yra skardinės sferos masė išintegruota trimis ašimis M=Sx3dydz+Sy3dzdx+Sz3dxdy. Pagal Gauso samprotavimo formulę skardinės sferos masė yra M=12πR55. Kiek suprantu, Gauso formulė iškraipo prasmę integravimas paviršiumi, nes tik su iškraipyta prasme integravimas paviršiumi Gauso formulė gali egzistuoti kaip teisinga formulė.
Gauto rezultato Sz3𝐝x𝐝y=1312πR55=4πR55 prasmę galima išaiškinti taip: sfera spindulio R=10 projektuojasi į plokštumą xOy; sferos centras yra taškas O, sferos projekcija į xOy plokštumą yra skritulys, kurio centras koordinačių pradžios taškas O; skritulio formulė yra x2+y2R2; skritulio plotas telpa į apskritimą kurio formulė x2+y2=R2; kadangi apskritimo plokštumoje xOy spindulys r=10 kaip ir sferos spindulys, tai į plotą πr2=π102=100π=314.1592653589793 telpa 314 strypų lygiagrečių Oz ašiai ir atstumas tarp strypų ant xOy plokštumos yra vienodas; kiekvienas strypas susikerta su sferos paviršiumi ir kiekvieno strypo ilgis yra nuo plokštumos xOy iki susikirtimo su sferos paviršiumi (mes skaičiuojame tik vienam pusrutuliui, tik teigiama Oz ašimi); trumpiausias strypo ilgis yra 0, o ilgiausias strypo ilgis yra iš centro O ir lygus R=10; dabar kiekvieną strypo ilgį reikia pakelti kubu, nes z3; Todėl strypo iš centro O (sutampančiam su ašimi Oz) ilgis yra z3=R3=103=1000; tolstančių nuo centro strypų ilgis trumpėja, o ant apskritimo kraštų strypų ilgis artimas arba lygus nuliui; rezultatas 4πR55=4π1055=80000π=251327.412287 yra visų strypų ant plokštumos xOy ir lygiagrečių ašiai Oz ilgių suma (dviejų pusrutulių).


  • Pavyzdis. Apskaičiuoti integralą Sx2dydz+0dzdx+0dxdy pagal išorinę pusę sferos x2+y2+z2=R2.
Taikydami formulę Gauso, gauname:
Sx2dydz=V(2x+0+0)dxdydz=202πdϕπ2π2cosθdθ0Rρ3dρ=
=202πdϕπ2π2cosθdθρ44|0R=2R4402πdϕπ2π2cosθdθ=R4202πdϕsinθ|π2π2=R4202π(sinπ2sinπ2)dϕ=
=R4202π(1(1))dϕ=R4ϕ|02π=R42π=2πR4.
Neteisingai skaičiuota, nes x polinėse ir sferinėse koordinatėse nėra ρ, bet yra x=R2y2z2 ir per ρ sferinėse koordinatėse net išreikšti negalima (nes sferinėje koordinačių sistemoje ρ2=R2=x2+y2+z2, o x išreikšti galima taip: x=Rsinθcosϕ (o gal taip: x=ρsinθcosϕ)). Žemiau turėtų būti teisingai paskaičiuota. Bet šansai, kad jei x pakeisti R2y2z2=R2ρ2 ir bus gautas teisingas atsakymas (skaičiuojant kaip aukščiau sferinėse koordinatėse) yra labai maži (bandžiau integruot su internetiniu integratoriumi ir gaunasi dalyba iš nulio, šaknyje minusas ir/ar panašiai).
Kitaip patikrinsime apskaičiuodami Sx2𝐝y𝐝z.
x2=R2y2z2,
x=R2y2z2.
Vx=Sx2𝐝y𝐝z=S(R2y2z2)2𝐝y𝐝z=S(R2ρ2)2ρ𝐝ρ𝐝ϕ=
=02π(0Rρ(R2ρ2)𝐝ρ)𝐝ϕ.
Pasinaudodami internetiniu integratoriumi, gauname, kad
0Rρ(R2ρ2)𝐝ρ=0R(ρR2ρ3)𝐝ρ=(ρ2R22ρ44)|0R=
=R2R22R44=R44.
Vx=02π(0Rρ(R2ρ2)2𝐝ρ)𝐝ϕ=02πR44𝐝ϕ=2πR44=πR42.
Kadangi reikia dviejų rutulio pusrutulių (teigiama ir neigiama Ox kryptimi), tai
VX=2Vx=2πR42=πR4.
Atsakymai 2πR4 ir πR4 nesutampa (gal atsakymas πR42 yra tik ketvirtadalis sferos, bet tada kyla klausimas: kodėl ne aštuntadalis sferos?).
Šiaip, atsakymas πR4 turėtų būti teisingas.
Teisingas skaičiavimas pagal tūrio formulę (skaičiavimas sferinėje koordinačių sistemoje). Mums prireiks šitos trigonometrinės formulės cos(2θ)=12sin2θ, sin2θ=1cos(2θ)2. Pagal Gauso formulę (sferinėse koordinatėse turime, kad x=ρsinθcosϕ):
Sx2dydz=V(2x+0+0)dxdydz=202πcosϕdϕ0πsin2θdθ0Rρρ2dρ=
=202πcosϕdϕ0πsin2θdθρ44|0R=2R4402πcosϕdϕ0πsin2θdθ=
=R421202πcosϕdϕ0π(1cos(2θ))dθ=R421202πcosϕdϕ(θsin(2θ)2)|0π=
=R4402πcosϕdϕ([π0][00])=πR4402πcosϕdϕ=πR44sinϕ|02π=πR44(sin(2π)sin0)=0.
Kažkas nesiintegruoja. Galima pabandyti integruoti per ϕ nuo 0 iki π/2 ir paskui viską padauginti iš 4:
4πR44sinϕ|0π/2=4πR44(sin(π/2)sin0)=πR4(10)=πR4.

Taip pat skaitykite