Inlägg

Inlägg som raol har skrivit i forumet
Av raol
Citat:

Ursprungligen inskrivet av Satsujin
raol:

Du råkar inte kunna Mathematica också? Vi har inte matlab.. =(

Vad är ode-rutinen för något?

Tack som fan för hjälpen iaf!

aldrig använt det programmet, tyvärr..
ODE23, ODE45 osv är rutiner i matlab för att numeriskt lösa ordinära differentialekvationer

Av raol
Citat:

Ursprungligen inskrivet av Gibson
Jag har slagit vad om 50 spänn...

kompis:

-5--5 blir 10 alltså minus 5, minus minus 5 blir tio

jag: nej, det blir 0

Jag slår in på miniräknaren och visar att det blir 0. Han ger mig 50 spänn men vill ha lite mer bevis, så jaaa, vad väntar ni på

-5-(-5)=0
-(5-(-5))=-10

a-b är definierat som det tal x som satisfierar b+x=a
Vill vi visa att -5-(-5)=0 måste vi därför visa att
-5+0=-5
Men a+0=a för varje heltal a, detta är själva definitionen av nollelementet.
Det torde väl övertyga din kompis...

Av raol
Citat:

Ursprungligen inskrivet av Sweclocker
1+1=?

1+1=2

Utifrån ett axiomsystem för de naturliga talen, t ex Peanos, kan bevisa påståendet. Där byggs de naturliga talen upp utifrån ett tal 0 som inte är efterföljare till något tal. Om n är ett tal är också efterföljaren S(n) ett tal. Man kan då definiera 1 som S(0) och 2 som S(S(0)) = S(1). Addition definieras rekursivt genom att 0 + n = n och S(m) + n = S(m + n). Enligt definitionen är

1 + 1 = S(0) + 1 = S(0 + 1) = S(1) = 2

Av raol

Re: Ny uppgift..

Citat:

Ursprungligen inskrivet av Satsujin
Ok, här kommer en till uppgift för er som inte har något bättre för er..

Då ett föremål faller mot marken utsätts det dels för tyngdkraftsaccelerationen g=9.82m/s^2, dels för en bromsande kraft, luftmotståndet, som ökar med fallhastigheten. En vanlig modell är att luftmotståndet är proportionellt mot m*(v^a), där m är massan och v fallhastigheten vid en viss tidpunkt. Eftersom accelerationen är derivatan av hastigheten kan vi formulera detta med en differentialekvation:

v´=9.82 - k*(v^a)

Vi antar i det följande att utgångshastigheten v[0] = 0.

a) Lös denna differentialekvation numeriskt med k = 0.19 och a = 1.4. Vilken blir gränshastigheten?

b) Antag att a = 1.4. Försök att bestämma k så att gränshastigheten blir 5 m/s.

Finns det nån som kan hjälpa till med den här?? Och vad menas med gränshastighet??

löste den i matlab...
hastighet.m:

function vdot=hastighet(t,v,flag,k,a) vdot=9.82 - k*(v^a);

exekvera:

options = odeset('RelTol',1e-8,'AbsTol',1e-8); [t,y]=ode45('hastighet',[0 20],0,options,0.19,1.4); plot(t,y)

ungefär 16,743 blir gränshastigheten, dvs hastigheten när t går mot oändligheten
Dels ser man det i resultatet, men man kunde även ha räknat fram det genom att lösa ekvationen 9.82-0.19*v^1.4=0

för att lösa b) vill man hitta en lösning till ekvationen
9.82-k*5^1.4=0
Det ger k=1.0317
Verifieras genom att köra

[t,y]=ode45('hastighet',[0 20],0,options,1.0317,1.4); plot(t,y)

[0 20] avser tidsintervallet, det går givetvis bra att ändra
options anger toleranser för ode-rutinen i matlab, de kan man ju också ända om man vill ha hög noggrannhet

Av raol
Citat:

Ursprungligen inskrivet av Ulvenstein

Vad pluggar du? Du verkar väl insatt i detta.

Jag läser min fjärde termin till civ.ing. i teknisk fysik och elektroteknik (y-programmet kallat).

Av raol
Citat:

Ursprungligen inskrivet av Ondjultomte

skulle du kunna ge ett exempel?

http://mathworld.wolfram.com/LebesgueIntegral.html
här står det om Lebesgueintegralen

Ett exempel är Dirichlets funktion:
f(x)=1, x irrationellt
f(x)=0, x rationellt
Den är lebesgueintegrerbar men inte riemannintegrerbar.

Mängden rationella tal är uppräkningsbart och har därför lebesguemåttet noll.
Lebesgueintegralen från noll till ett av Dirichlets funktion är ett, då lebesguemåttet av alla irrationella tal på intervallet [0,1] är 1.

Av raol
Citat:

Ursprungligen inskrivet av Ulvenstein

När man beräknar en integral använder man bara just en metod för att beräkna gränsvärdet där översumman går mot undersumman. Metoden att använda primitiv funktion kommer från Leibniz, och Riemann visade att denna metod genererar en integral som uppfyller de aktuella villkoren. Det finns en grupp av bevis för detta, som jag visserligen tentat på men inte kommer att gå in på här.

"men inte kommer att gå in på här".. hehe, låter lite som på en föreläsning.
Men det är bevis för analysens huvudsats du avser..?

Vore kul att höra vad du har för examen Ulvenstein och hur många poäng matematik har du läst? (själv har jag nått upp i ca 35..)

Av raol
Citat:

Ursprungligen inskrivet av magnifique

Bevis: Givet epsilon > 0 så gäller det att hitta F,Y så att
I(Y) - I(F) = Int(Y(x)dx) - Int(F(x)dx) (med gränserna från a till b) < epsilon.

Antag att f är växande och dela in intervallet
i n lika delar a = x0 < x1 < x2 <K< xn-1 < xn = b
där xk - xk-1 = (b - a)/n. Vi definierar nu F,Y genom
att på [xk-1, xk] sätta F(x) = f (xk-1) och Y(x) = f (xk).

(med xk menar jag x index k)

I(Y) - I(F) = Summa (k=1 till n)(f(xk)(xk - xk-1 )) - Summa (k=1 till n)(f(xk-1)(xk - xk-1)) = (b - a)/n * Summa (k=1 till n)(f(xk) - f(xk-1)) = (b - a)/n*(f(b) - f(a)) < epsilon
för n tillräckligt stort. VSB.

Jovisst, gjorde ju också precis så för att visa det.
Om man ska göra nåt påpekande så har du ju inte motiverat att Y(x)>=F(x), a<=x<=b, men det följer ju av att f är växande.
Dessutom så måste man ju sätta F(x) = f (xk-1) och Y(x) = f (xk) på ]xk-1, xk[ och inte på ett slutet intervall, för annars blir ju inte F och Y så väldefinierade i xk.

Av raol
Citat:

Ursprungligen inskrivet av Ondjultomte
är det nån som kan integrera på andra sätt utan riemann summor ?

Hur menar du...? (Riemann-)Integralbegreppet är ju definierat utifrån över- och undersummor, men när man beräknar en integral så använder man ju sällan sådana summor utan tar fram en primitiv funktion och utnyttjar analysens huvudsats.

Av raol
Citat:

Ursprungligen inskrivet av Rambo

LOL!

skratta på du... nu har jag lyckats visa det.

Av raol

Nån som knäcker den här?
Visa att en monoton funktion definierad på ett slutet intervall är (Riemann)integrerbar.

Av raol

Re: Problemstimulering

Citat:

Ursprungligen inskrivet av magnifique
Det är inget riktigt drag här nu. Känner att många väntar på ett litet problem att knäcka.

En gumma har n st. ägg i en korg. Om hon tar upp 3 ägg åt gången får hon till slut 2 ägg kvar. Om hon tar 5 ägg åt gången får hon till slut 3 ägg kvar och tar hon 7 ägg åt gången får hon till slut 5 ägg kvar. Vad är det minsta antalet ägg hon kan ha i sin korg?

n=3x+2
n=5y+3
n=7z+5

Vi har att:
3x+2=5y+3 <=> 3x-5y=1
(2,1) är en lösning
Skriver om ekvationen:
3(x-2)=5(y-1)
vi måste ha heltalslösningar så detta ger att y-1=3*k, x-2=5*k
Alla lösningar:
x=2+5*k
y=1+3*k

Vi har också:
5y+3=7z+5 <=> 5y-7z=2
(6,4) är en lösning.
5(y-6)=7(z-4)
således y-6=7*m, z-4=5*m
y=6+7*m
z=4+5*m

Har nu att 1+3k=6+7m <=> 3k-7m=5
(4,1) är en lösning
3(k-4)=7(m-1)
k=4+7*l
m=1+3*l

Således:
x=2+5*k=2+20+35*l=22+35*l
y=1+3*k=1+12+21*l=13+21*l
z=4+5*m=4+5+15*l=9+15*l

n=3x+2=3*(22+35*l)+2=68+105*l

Således har hon minst 68 ägg!

(på köpet fick du också alla lösningar )

Av raol
Citat:

Ursprungligen inskrivet av Kevlar
"Ett tåg lämnar stationen och färdas norrut med hastigheten 75 km/h. Två timmar senare lämnar ett annant tåg med hastigheten 125 km/h som kör i samma riktning. När hinner det senare upp det första? Hur längt från stationen befinner de sig då?"
Vet inte riktigt hur jag ska ställa upp ekvationssystemet, kan nån hjälpa mig?

s=75*t
s=125*(t-2)

Av raol
Citat:

Ursprungligen inskrivet av Kevlar
Men... (4x/3) + 5x = 19x/3 ?

ja just det..

Av raol
Citat:

Ursprungligen inskrivet av Kevlar

Edit2: -... som då blir 5x + (2x/6) = -19 ...

Men hur ska jag göra nu för att lösa ut de 2x, bara multiplicera allt med 6 ?

nä.. det ska stå 4x/3 där det står 2x/6, och nästa steg jag gör är att skriva det på gemensamtbråkstreck (utnyttjar 5x=15x/3)

Av raol
Citat:

Ursprungligen inskrivet av Kevlar
Hmmm.... jag är lite trött nu , ska det vara:

5x-2(-2x/3)=-19 => 5x - (2x/6) =-19

Om det är rätt, hur fortsätter jag då? hehe

Resten är ju trivialiteter...
Nej 5x-2(-2x/3)=-19 <=> 5x + (4x/3) =-19 <=> 19x/3=-19 <=> x=-3

Av raol
Citat:

Ursprungligen inskrivet av Kevlar

Hmm... det ska väl vara 5x-2(-2x/3)=-19

stämmer, slarvigt av mig

att kontrollera sitt svar i de ursprungliga ekvationerna när man är klar är alltid bra ...

Av raol
Citat:

Ursprungligen inskrivet av Kevlar
Problem med ekvationssystem..
__________________________________

{ 5x - 2y = -19
{ 2x + 3y = 0

Fast klammerna ska vara en klammer

Jag kommer en bit men sen kraschar allt o jag blir sur

Kan nån förklara hur man löser denna med substitutionsmetoden

5x - 2y = -19 (1)
2x + 3y = 0 (2)

(2) <=> y=-2x/3, sätt in i (1):
5x-2(2x/3)=-19
lös detta och sätt in värdet på x i y=-2x/3

Av raol
Citat:

Ursprungligen inskrivet av Ondjultomte
nån som vet hur man använder commandot roots i matlab?

roots hittar nollställena till ett valfritt polynom... du skriver in polynomets koeffecienter i vektorform
t.ex. roots([1,2,3])
finner nollställena till x^2+2*x+3

Av raol
Citat:

Ursprungligen inskrivet av pavs
Nån som vet om man kan integrera e^x^2 analytiskt?

edit: jag är dum

om du menar (e^x)^2 så är det inga problem
En primitiv funktion är 1/2*e^(2x)

Däremot kan inte en primitiv funktion till e^(x^2) utryckas i elementära funktioner...