Question
Say I have a very complicated probability distribution function: x^(-2)*y^(-3)*exp(-z/t), where x y z t are all continuous positive numbers. The integration in x,y,z,t space is not 1, so there has to be some scaling factor. How can I devise a function, or is there a convenient function with some symbolic that allows me to sample a random, continuous set of numbers from this distribution? What if I want to sample a full conditional distribution, say x, with y,z,t as known?
Expert Answer
John Michell
PhD Expert
Answered Aug 31, 2026
There are two key elements to your problem:
- Sampling from an arbitrary distribution
- Sampling from a multivariate distribution (and are those variables correlated?)
To do the first part -- for example just getting samples from x^(-2) -- there is a method using the inverse of the cdf. I suggest you take a look at this thread in the old MATLAB forum to get started.
If you have independent variables, then you can just multiply the distributions of each variable. If you have dependent variables, then your best bet is to use the copula approach.
There are examples in that documentation that will help you get started. I suggest that you try some baby steps for each of these two parts, so that you build your understanding, rather than trying to solve your whole problem at once.
100% Run Guarantee
3-Hour Fast-Track Delivery
Need a Custom Version or Complete Simulation for This Problem?
Our 500+ PhD engineers build, debug, and optimize working MATLAB scripts and Simulink (.slx) models tailored to your exact assignment rubrics with zero plagiarism.
Tested on MATLAB R2024b / R2026a
Turnitin 0% Plagiarism Report
Free 7-Day Revisions Guarantee
Have a different question? Ask here
Related probability distri... Questions & Solutions
Browse All →
Explore similar technical troubleshooting questions and verified MATLAB solutions:
Ready-to-Run MATLAB & Simulink Projects
Browse All Projects →