10 svar
164 visningar
AstridG behöver inte mer hjälp
AstridG 66
Postad: 30 okt 2020 17:13 Redigerad: 30 okt 2020 22:16

Amplitude probability density function i Matlab

Hej 

Jag behöver hjälp med hur man ska beräkna sannolikheten från en bar graf jag har fått fram grafen i 20 st fördelad. Men vet inte riktigt hur jag bör kunna få i sannolikhet.

A = importdata('Bearing.mat')
t = [-Inf:0.1:Inf];
[f, x] = hist(A.Bearing1, 20)
sigma = std(f)
p = (1/(sigma*sqrt(2*pi)))*exp(-((x*2)/(2*sigma*sigma)));
I = trapz(t, p(x))
bar(x,f/trapz(x,f))

vet inte riktigt om den ger mig sannolikheten för a till b, 

SaintVenant 3956
Postad: 30 okt 2020 19:02

Sannolikheten mellan a och b får du om du integrerar din funktion pxp\left(x\right) från a till b enligt:

Pa<Xb=abpxdx\displaystyle P\left(a<X\leq b\right)=\int_{a}^{b}p\left(x\right) dx

Du vet om att det är en normalfördelning? Du har integrerat täthetsfunktionen du konstruerade med hjälp av den inbyggda trapetsregeln trapz i Matlab. Det enda spåret av ursprungsdatan i det läget är standardavvikelsen du beräknade med std-funktionen.

Är datan i "Bearing.mat" verkligen kontinuerlig? Kan du ge mer information? Kanske lägga upp bilder på några plots...

AstridG 66
Postad: 30 okt 2020 21:07
Ebola skrev:

Sannolikheten mellan a och b får du om du integrerar din funktion pxp\left(x\right) från a till b enligt:

Pa<Xb=abpxdx\displaystyle P\left(a<X\leq b\right)=\int_{a}^{b}p\left(x\right) dx

Du vet om att det är en normalfördelning? Du har integrerat täthetsfunktionen du konstruerade med hjälp av den inbyggda trapetsregeln trapz i Matlab. Det enda spåret av ursprungsdatan i det läget är standardavvikelsen du beräknade med std-funktionen.

Är datan i "Bearing.mat" verkligen kontinuerlig? Kan du ge mer information? Kanske lägga upp bilder på några plots...

Båda första bilderna är från Bearing.mat som har två olika info. Och den är fördelad i 20 eftersom uppgifter ville ha i 20. För standardavvikelse är lite osäkert ifall jag har tagit från rätt information. Kan kopiera uppgiften så du förstår vad de handlar om.

 

1.3. Amplitude probability density functions
Load the file “Bearing.mat” that contains two signals of rotating machinery. One of the
signals contains a faulty bearing and one is from a healthy machine.


Assignments:


1. Calculate the probability density functions for the signals in MATLAB using the
histogram “hist” command. Plot the probability density functions as bar plots
with at least 20 segments. What is the value of the peak in the plot? Verify that the
total probability is one.


2. Calculate the skewness and kurtosis of the signals. From these quantities determine
which of the bearings you suspect is faulty. Describe how you reached your
conclusion.

SaintVenant 3956
Postad: 30 okt 2020 21:44

Oj, wow, vilken spännande uppgift.

Jag undrar varför de vill att man ska använda hist-funktionen egentligen, den rekommenderas inte längre i MATLAB efter rev. 2014b. Du kan hursomhelst läsa denna tråd för tips:

How can we calculate probability density of a histogram?

Jag hoppas det hjälper dig. Jag har inte tid att titta mer specifikt nu och får återkomma imorgon om du behöver mer hjälp.

AstridG 66
Postad: 30 okt 2020 22:00 Redigerad: 30 okt 2020 22:16
Ebola skrev:

Oj, wow, vilken spännande uppgift.

Jag undrar varför de vill att man ska använda hist-funktionen egentligen, den rekommenderas inte längre i MATLAB efter rev. 2014b. Du kan hursomhelst läsa denna tråd för tips:

How can we calculate probability density of a histogram?

Jag hoppas det hjälper dig. Jag har inte tid att titta mer specifikt nu och får återkomma imorgon om du behöver mer hjälp.

Jag kan vänta tills imorgon ifall jag löser inte idag. Jag tror de har bara missat uppdatera funktionen, de bruka använda samma uppgift för varje kursomgång  och förmodligen har dom missat ändra matlab koden.

 

Btw, de blev så här nu de länken du har skickat.

Smaragdalena 80504 – Avstängd
Postad: 30 okt 2020 22:23

Tog bort dina skrikiga versaler. Det står i Pluggakutens regler (och i rutan där du skriver in rubriken) att du skall undvika att skriva rubriken med endast stora bokstäver. /moderator

Dr. G 9500
Postad: 30 okt 2020 22:32
AstridG skrev:

Har det här histogrammet något att göra med de tidigare två?

AstridG 66
Postad: 30 okt 2020 22:41
Dr. G skrev:
AstridG skrev:

Har det här histogrammet något att göra med de tidigare två?

När du säger det, så glömde jag göra det separat. Eftersom det är två olika signaler som filen ger så fick jag byta variabeler så de bilda två separata resultat. 

Dr. G 9500
Postad: 30 okt 2020 22:50

Menar du att de två senaste histogrammen och de två första i tråden beskriver samma data?

Vad är det på x-axeln i de två fallen? 

De senare histogrammen ser ut att komma från hyfsat normalfördelad data, men inte de två första. 

AstridG 66
Postad: 30 okt 2020 23:17
Dr. G skrev:

Menar du att de två senaste histogrammen och de två första i tråden beskriver samma data?

Vad är det på x-axeln i de två fallen? 

De senare histogrammen ser ut att komma från hyfsat normalfördelad data, men inte de två första. 

dom två första fallet är från

histogram(A.Bearing1, 20)
histogram(A.Bearing2, 20)

och mätte upp the peak, problemet är att jag vet inte hur jag ska bevisa att båda integral av p(x) med intervall oändlighet negativ och positiv är lika med 1. (är inte alls duktigt med matlab, så du vet). x-axis borde va duration om jag har förstått rätt.  

SaintVenant 3956
Postad: 31 okt 2020 01:56

Du kan läsa här för att ta reda på hur du bestämmer arean i din plot:

Calculate area under histogram

Om man kikar på dina bilder övergripande så kan jag snabbt säga att det ser ut som om arean är 1 för båda plots.

Svara
Close