Tråden för oss som använder Matlab

Hej, jag har under ett antal timmar försökt med denna uppgift som lyder:
Antag att kvantiteten q(t) av en vara på en marknad tid tiden t (mätt i dagar) är skild
från jämviktskvantiteten q∗. Det är rimligt att anta att den takt med vilken q(t) då
ändras i tiden är proportionell mot skillnaden mellan det pris konsumenter är villiga att
betala för q(t) och det pris producenter är villiga att leverera den för:
q'(t) = k(D(q(t) − S(q(t)))
där D är efterfrågefunktionen och S är utbudsfunktionen.
D = @(q) 300*exp(-0.02*q.^2);
S = @(q) 20*q.^(1.2) +2*q.^(0.4);
a) Dag 0 är kvantiteten på marknaden 6.53 *10^4 och dag 1 är den 6.43*10^4. Använd den informationen för att bestämma en approximation till q'(0).
b) Använd approximationen av q'(0) för att bestämma en approximation till k.
c) Lös begynnelsevärdesproblemet
q'(0) = k(D(q) − S(q))
q(0) = 6.53 *10^44
numeriskt.

Hur många dagar tar det innan kvantiteten på marknaden är mindre än 1000
enheter från jämviktskvantiteten q∗ ?
(d) Plotta D(q(t)) och S(q(t)) i samma figur. Dekorerar med lämplig text.

jag bifogar en m-fil med en funktion för riktningsfält om det är användbart (https://drive.google.com/open?id=1-c...MuNPPZAr7xnkP8)

Ingen fråga om något, eller är det meningen att vi ska lösa hela skiten?

Oj, såg nu att jag glömde ställa den frågan jag hade. Jag har problem med c uppgiften, jag använder mig av matlabs inbyggda ode45 men lyckas inte komma fram till ett rimligt svar. Jag undrar därför om det är någon som kan hjälpa mig, jag är nybörjare på matlab.

Ska du verkligen börja med 10^44 då? Det är ju ett jobbigt stort tal.

oj, det ska stå 10^4

Kollar man Demandfunktionen kan man se att efterfrågan dör av väldigt snabbt. Verkar vara något konstigt med problemformuleringen där. Åtminstone lite oväntat att det ska vara så få produkter, med tanke på hur problemet är formulerat.
Wolfram

 Ursprungligen postat av Swordplay

När man beräknar integraler med trapz, som exemplet i
https://se.mathworks.com/help/matlab/ref/trapz.htm

I = trapz(y,trapz(x,F,2))

om man vill bestämma antalet intervaller till ex 3 eller 4, hur skrivs det?


Om man vill bestämma intervallet själv?

Du ska tillhandahålla de punkter som används i beräkningen, så det är du som bestämmer.

Hej! Jag håller för tillfället på att lära mig om Fourier Series och snubblade över ett övningsproblem som lyder följande:
Create a Matlab Script, where you can plot the Fourier series of the square wave with 1 to 10 coefficients. Explain what is happening and why. Sketch the absolute value of the frequency spectrum of the signal. Add a plot of the frequency spectrum of the signal. You should use the fft function in matlab for this. Look at the help in matlab, you can find examples of how to do there.

Square Wave finns som bild på länken nedan:

https://imgur.com/a/TZFXtiA

Jag kan beräkna det hela teoretiskt med differentialekvationer, integraler, fourier formler osv. Men förstår verkligen inte hur man gör det i Matlab, vilket egentligen är vad jag behöver hjälp med här. Kollade igenom fft function dokumentationen som fanns i Matlab, men har fortfarande ingen direkt aning om hur jag ska gå till väga för att uppnå detta. Finns det någon här som vet hur något sådant script hade kunnat se ut i matlab, specifikt för en sådan square wave?

Sampla fyrkantvågen, så du har en vektor med fyrkantsvågens värden vid samplingstiderna. Tryck in vektorn i fft, ut kommer fouriertransformen.

Skall lösa Väntevärde μ med integraler

μ = E(X) = ... x f(x) dx
se form för u Kontinuerlig slumpvariabel på https://sv.wikipedia.org/wiki/Standardavvikelse

där f(x) frekvensfunktion = 0.2e^-0.2x
, där x >= 0

Funktionen som skall ge u tolkar jag det skall då se ut så här
f = x*(0.2e^-0.2x);

men hur gör man med intervallet då man bara har x >= 0?

Jag tolkar det som att man kan sätta in ett godtyckligt ändvärde utifrån grafen, tex 100.

syms x
f = x*(0.2e^-0.2x);
A = int(f, x, 0, 100)

eller?

Skall även hitta Standardavvikelse σ med integraler för stokastiska variabeln ovan,
Standardavvikelse σ = sqrt( .... (x - μ)^2 f(x) dx)
Se länken om Kontinuerlig slumpvariabel https://sv.wikipedia.org/wiki/Standardavvikelse

där f(x) frekvensfunktion = 0.2e^-0.2x
, där x >= 0

Hur skall detta se ut i Matlab?
Funderar vidare på intervallet som är samma