Difficult Probability Density Function Question

Click For Summary
SUMMARY

The discussion focuses on solving a complex probability density function (PDF) problem involving two functions, ##f(x)## and ##g(x)##, defined as ##f(x)=|a^{\frac{\sin(x)}{\ln(ax)}}-\frac{x}{a}|## and ##g(x)=|\frac{\sin(\ln(\sqrt{x}-\sqrt{a}))}{x^2-a^2}|##. Participants aim to determine the constant ##a##, the coordinates of intersection points ##Q(b,f(b))## and ##R(c,f(c))##, and to calculate the mean, variance, and standard deviation of the continuous random variable ##X##. The discussion emphasizes the need for numerical methods to find suitable values for ##a## within the constraints of the problem.

PREREQUISITES
  • Understanding of probability density functions (PDFs)
  • Familiarity with calculus, specifically integration and differentiation
  • Knowledge of numerical methods for solving equations
  • Proficiency in a programming language for numerical analysis (e.g., Python, MATLAB)
NEXT STEPS
  • Learn numerical methods for root-finding, such as the Newton-Raphson method
  • Study the properties of probability density functions and their applications
  • Explore integration techniques for calculating expected values and variances
  • Implement a numerical search algorithm in Python to find suitable values for ##a##
USEFUL FOR

Students studying probability and statistics, mathematicians working on continuous random variables, and data analysts involved in statistical modeling.

Saracen Rue
Messages
150
Reaction score
10

Homework Statement


A function, ##f\left(x\right)=\left|a^{\frac{\sin \left(x\right)}{\ln \left(ax\right)}}-\frac{x}{a}\right|##, intersects with another function, ##g\left(x\right)=\left|\frac{sin(ln(\sqrt{x}-\sqrt{a}))}{x^2-a^2}\right|##, at point ##Q(b,f(b))## and point ##R(c,f(c))##. A probability density function, ##p(x)=g(x)-f(x)##, can be formed over the domain ##[b,c]##, where ##0<b<c<2##.

Q1. Determine, correct to 16 decimal places;
(a) The value of the real-valued constant, ##a##
(b) The coordinates of the points ##Q## and ##R##
Q2. Calculate the mean, variance and standard deviation correct to 4 decimal places

Q3. Find correct to 1 decimal place; the percentage probability of the continuous random variable, ##X##, being within ##a^2## variances either side of the mean (i.e. ##Pr(μ-a^2⋅Var(X)≤X≤μ+a^2⋅Var(X))##)

Homework Equations


For PDfs: ##∫_{b}^{c}f(x)dx=1##

The Attempt at a Solution


I know that to find points ##b## and ##c## I need to solve ##f(x)=g(x)## for ##x##, however when I try to do this on my calculator it gives me an error message. I know that I should be getting ##b## and ##c## in terms of ##a##, thus ##∫_{b}^{c}p(x)dx=1## should have ##a## as the only unknown and it should be able to be solved for. After getting ##a##, I'd substitute the value into ##b## and ##c## to get those values as well, and then I'd be able to find ##f(b)## and ##f(c)## as well (giving me points ##Q## and ##R##).

To find the mean, variance and standard deviation I'd just solve the following:
##E(X)=∫_{b}^{c}x⋅p(x)dx##, ##Var(X)=∫_{b}^{c}x^2⋅p(x)dx## and ##σ=\sqrt{Var(X)}##

And for ##Q3## I'd solve this: ##Pr(μ-a^2⋅Var(X)≤X≤μ+a^2⋅Var(X))=∫_{μ-a^2⋅Var(X)}^{μ+a^2⋅Var(X)}p(x)dx##

I am fairly confident with my method for the most part (although I'd still appreciate somebody telling me if there's any errors in it), the main issue I'm having is simply not being able to solve the equations on my calculator. I'd be very thankful for anybody who has an alternative way for me to solve this problem :)
 
Physics news on Phys.org
If we can identify a nice search space then it shouldn't be too hard to numerically find a suitable value of ##a##.

To start that process, we can observe that, if we want ##f(x)-g(x)## to be defined and continuous on ##b,c##, we must have ##1\leq a<b##. Can you see why that must be the case?

Given that, we can differentiate ##h\triangleq f-g## and set it to zero to try to find extrema, of which there will have to be at least one between ##b## and ##c##.

Or - quicker and easier if it's allowed, just write a short bit of code in your favourite number-crunching language to search for a value of ##a\in[1,2]## for which ##h(x)## attains a maximum or minimum in the interval ##[a,2]## and then for such an ##a## see if the curve of ##h## crosses the ##x## axis on either side of that extremum.

By the way, the question is worded in a way that seems to imply that there is only one possible value of ##a##. In fact there is a range of possible values. So you only need to find a value of ##a## that will work.
 
andrewkirk said:
If we can identify a nice search space then it shouldn't be too hard to numerically find a suitable value of ##a##.

To start that process, we can observe that, if we want ##f(x)-g(x)## to be defined and continuous on ##b,c##, we must have ##1\leq a<b##. Can you see why that must be the case?

Given that, we can differentiate ##h\triangleq f-g## and set it to zero to try to find extrema, of which there will have to be at least one between ##b## and ##c##.

Or - quicker and easier if it's allowed, just write a short bit of code in your favourite number-crunching language to search for a value of ##a\in[1,2]## for which ##h(x)## attains a maximum or minimum in the interval ##[a,2]## and then for such an ##a## see if the curve of ##h## crosses the ##x## axis on either side of that extremum.

By the way, the question is worded in a way that seems to imply that there is only one possible value of ##a##. In fact there is a range of possible values. So you only need to find a value of ##a## that will work.

I think the question is implying that there will be only one value of ##a## for which the equation ##1 = \int_{b(a)}^{c(a)} h(x,a) \, dx## will work.
 

Similar threads

  • · Replies 19 ·
Replies
19
Views
4K
Replies
6
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
Replies
8
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
Replies
2
Views
2K
  • · Replies 9 ·
Replies
9
Views
2K
Replies
17
Views
2K