Hjälp med att skriva funktioner i Matlab

Hur skriver jag de här funktionerna i Matlab?


x^(1/5)cos(x)e^(x^(1/7))



e^(sin(x))

Ska du bara skriva funktionen rakt av och inte räkna med någon vektor skriver du bara:


Vill du definiera det som en funktion f(x) kan du skriva:


Där du sen skriver y(x) för att räkna ut värdet för ett visst x.

Jag ska integrera de

spännande

Fan, i mitt nästa liv ska min Flashbackanvändare heta pMorberg!

Då får du kolla på Euler framåt eller Runge-Kutta eller något liknande.

Okej, jag skriver hela uppgiften...


Använd Matlab och Stora talens lag för att (med "Monte-Carlo metoden" ovan) med god noggrannhet bestäma

a) integraltecken från 0 till 1 e^(x^(2))dx

b) integraltecken från 0 till 1 x^(1/5)cos(x)e^(x^(1/7))dx

c) integraltecken från 0 till 2,38 e^(sin(x))dx

Som redovisning lämnas svar och Matlabkod.


Här är en bild på uppgiften ifall ni inte förstår: http://imageshack.us/a/img145/8292/sannstat2.jpg


a) Är det inte bara göra så här?

fun = @ (x) exp(x.^2) ;
q = integral (fun,0,1)
q = 1,4627


Så det jag undrar över är hur jag skriver b) och c) uppgiften i Matlab?

(1) Vet du vad Monte-Carlo är?
(2) Använder du någonstans Monte-Carlo?

Om svaret på (1) är ja då borde du inse att svaret på (2) är nej och därmed är lösningen fel.
Om svaret på (1) är nej då borde du läsa på vad Monte-Carlo är och inse att svaret på (2) är nej.

Således: Du använder inte Monte-Carlo och därmed är uppgiften inte korrekt löst.

Monte-Carlo går väl ut på att man ska skriva om integralen som ett väntevärde typ:

I = integraltecken från 0 till 1 g(x)dx = integraltecken från 0 till 1 g(x)*1 dx = E(g(U)) dvs I ≈ 1/n ∑^n g(u_k)

Men jag förstår inte hur man ska kombinera det med Matlab i sådanna fall? Hur blir funktionerna? Ngn som kan hjälpa mig?

Ja. Har du gjort så? Nej, alltså är det fel. Säg att du vill integrera funktionen f(x) = x mellan 0 till 1 med 1000 "samplingar".

N = 10^3 ;
f = @(x) x ;
1/N * sum(f(rand(1,N)))

rand genererar ett slumptal mellan 0 och 1 och rand(1,N) genererar N stycken sådana. Nu måste du göra likadant för (a),(b) fast andra funktioner och för (c) måste du skifta så den generar mellan 0 och 2.38. Hur gör du det?

Det borde väl bli

N = 10^3 ;
f = @(x) x ;
1/N * sum(f(rand(2.38,N)))

Haha tack Men var får alla ett R ifrån