See my post on Discourse:

https://scilab.discourse.group/t/exponentially-modified-gaussian-distribution

S.

On 4/30/24 08:46, Stéphane Mottelet wrote:
Hello,

As I already said, there is no need for using this technique, as an EMG
random variable is the sum of a Gaussian and an exponential. If your CDF
reads as

f = (a/2)*exp((a/2)*((a*σ^2)-2*h)) .* erfc((a*σ^2-h)/σ/sqrt(2))

then the underlying variable is the sum of Gaussian of mean 0, variance
σ^2, and an exponential of rate a. Then, you can simulate n draws of
this variable in Scilab with

grand(n,1,"nor",0,σ) + grand(n,1,"exp",a)

S.

On 4/30/24 01:44, Heinz Nabielek wrote:
Very useful. I had done a similar thing with the EMG, but much more
clumsy....
Heinz

On 29.04.2024, at 16:41, Federico Miyara <fmiy...@fceia.unr.edu.ar>
wrote:


Heinz,

I don't know if this might be useful. The function I'm attaching
allows to generate random numbers according to any distribution,
either empirical or analytical. The results may be scaled if necessary.

Regards,

Federico Miyara


On 29/4/2024 00:11, Heinz Nabielek wrote:
Colleagues:

bird flight altitude probabilities are given by the exponentially
modified Gaussian distribution EMG f=f(h), in my case

f = (a/2)*exp((a/2)*((a*σ^2)-2*h)) .* erfc((a*σ^2-h)/σ/sqrt(2))

with h=hmeasured-85 in meters, a=1/60m, σ=30m. See:
<https://antiphishing.vadesecure.com/v4?f=SDlhUlZiaXFEVHROSUhTb0qrO2AR7jfbvbNtVdf1_aC917oIdR0Gc2d8Rd6FD2Bc&i=TFpzSWRVdmRmMGQwa1BoSiimVxxfUPBdU67jVuwVi4s&k=4IV5&r=bzNQRmtWdWVld3VrdkQzNUBU-sIUgBXJ8jKxvWcVsg0fKphq97waH97_kCu1wpxV&s=e2aae9ad20320b735f15a05e15e449b31fb9fc2727a89289de82079d7c61fc15&u=https%3A%2F%2Fen.wikipedia.org%2Fwiki%2FExponentially_modified_Gaussian_distribution>.

For bird flight Monte-Carlo simulations, I need random deviates
representing this distribution.
Normally, I generate the cumulative function F of f. Then, I
project uniform random numbers U(0,1) on F^-1, the inverse of the
cumulative.

With the EMG, closed solutions are not available and I need a
numerial procedure to generate EMG deviates.
What is the best way to do this?
Best greetings
Heinz
This email and any attachments are intended solely for the use of
the individual or entity to whom it is addressed and may be
confidential and/or privileged.

If you are not one of the named recipients or have received this
email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete
this email and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according
to our data privacy policy as described on our website. Should you
have any questions related to personal data protection, please
contact 3DS Data Protection Officer
https://antiphishing.vadesecure.com/v4?f=SDlhUlZiaXFEVHROSUhTb0qrO2AR7jfbvbNtVdf1_aC917oIdR0Gc2d8Rd6FD2Bc&i=TFpzSWRVdmRmMGQwa1BoSiimVxxfUPBdU67jVuwVi4s&k=4IV5&r=bzNQRmtWdWVld3VrdkQzNUBU-sIUgBXJ8jKxvWcVsg0fKphq97waH97_kCu1wpxV&s=b5d6e33a0a14988d2e95f8136aae533cede18997ddc7ac3015294495eeebdb35&u=https%3A%2F%2Fwww.3ds.com%2Fprivacy-policy%2Fcontact%2F






_______________________________________________
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://antiphishing.vadesecure.com/v4?f=SDlhUlZiaXFEVHROSUhTb0qrO2AR7jfbvbNtVdf1_aC917oIdR0Gc2d8Rd6FD2Bc&i=TFpzSWRVdmRmMGQwa1BoSiimVxxfUPBdU67jVuwVi4s&k=4IV5&r=bzNQRmtWdWVld3VrdkQzNUBU-sIUgBXJ8jKxvWcVsg0fKphq97waH97_kCu1wpxV&s=3c9b8e836fc56309ec6648e54301842439f5f2e8a78a44013cb50d3010a031b6&u=https%3A%2F%2Flists.scilab.org%2Fmailman%2Flistinfo%2Fusers


--
Este correo electrónico ha sido analizado en busca de virus por el
software antivirus de Avast.
www.avast.com
This email and any attachments are intended solely for the use of the
individual or entity to whom it is addressed and may be confidential
and/or privileged.

If you are not one of the named recipients or have received this
email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete
this email and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to
our data privacy policy as described on our website. Should you have
any questions related to personal data protection, please contact 3DS
Data Protection Officer
https://antiphishing.vadesecure.com/v4?f=SDlhUlZiaXFEVHROSUhTb0qrO2AR7jfbvbNtVdf1_aC917oIdR0Gc2d8Rd6FD2Bc&i=TFpzSWRVdmRmMGQwa1BoSiimVxxfUPBdU67jVuwVi4s&k=4IV5&r=bzNQRmtWdWVld3VrdkQzNUBU-sIUgBXJ8jKxvWcVsg0fKphq97waH97_kCu1wpxV&s=b5d6e33a0a14988d2e95f8136aae533cede18997ddc7ac3015294495eeebdb35&u=https%3A%2F%2Fwww.3ds.com%2Fprivacy-policy%2Fcontact%2F




_______________________________________________
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://antiphishing.vadesecure.com/v4?f=SDlhUlZiaXFEVHROSUhTb0qrO2AR7jfbvbNtVdf1_aC917oIdR0Gc2d8Rd6FD2Bc&i=TFpzSWRVdmRmMGQwa1BoSiimVxxfUPBdU67jVuwVi4s&k=4IV5&r=bzNQRmtWdWVld3VrdkQzNUBU-sIUgBXJ8jKxvWcVsg0fKphq97waH97_kCu1wpxV&s=3c9b8e836fc56309ec6648e54301842439f5f2e8a78a44013cb50d3010a031b6&u=https%3A%2F%2Flists.scilab.org%2Fmailman%2Flistinfo%2Fusers


--
Stéphane Mottelet
Ingénieur de recherche HDR
EA 4297 Transformations Intégrées de la Matière Renouvelable
Département Génie des Procédés Industriels
Sorbonne Universités - Université de Technologie de Compiègne
CS 60319, 60203 Compiègne cedex
Tel : +33(0)344234688
https://antiphishing.vadesecure.com/v4?f=UWNQbzIzM2JENXZKZ05qY4zc4ZhHNJwVkRuQ3XxY754AgjIp-lIF9dDd_yicYx2w&i=cnlNUk9nak5NbVU2b0IwbrA-UKliAAYtcNAOAp_it7o&k=I0Xg&r=enZpZVAxT0hiZkROWVVaWLFrrKZpuQOhzxBoJqwbZQek__lI9PCLLRgnErXqdjCr&s=cb39676f427009e8577a40d41978922b512777b8ba6d8a93e43acc7c96cd3672&u=http%3A%2F%2Fwww.utc.fr%2F~mottelet


_______________________________________________
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://antiphishing.vadesecure.com/v4?f=UWNQbzIzM2JENXZKZ05qY4zc4ZhHNJwVkRuQ3XxY754AgjIp-lIF9dDd_yicYx2w&i=cnlNUk9nak5NbVU2b0IwbrA-UKliAAYtcNAOAp_it7o&k=I0Xg&r=enZpZVAxT0hiZkROWVVaWLFrrKZpuQOhzxBoJqwbZQek__lI9PCLLRgnErXqdjCr&s=2170f077e28101444bcdebb7db61222e7cb9cfcfac669aec1d62e1d591f8c735&u=https%3A%2F%2Flists.scilab.org%2Fmailman%2Flistinfo%2Fusers

This email and any attachments are intended solely for the use of the
individual or entity to whom it is addressed and may be confidential
and/or privileged.

If you are not one of the named recipients or have received this email
in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete
this email and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to
our data privacy policy as described on our website. Should you have
any questions related to personal data protection, please contact 3DS
Data Protection Officer
https://antiphishing.vadesecure.com/v4?f=UWNQbzIzM2JENXZKZ05qY4zc4ZhHNJwVkRuQ3XxY754AgjIp-lIF9dDd_yicYx2w&i=cnlNUk9nak5NbVU2b0IwbrA-UKliAAYtcNAOAp_it7o&k=I0Xg&r=enZpZVAxT0hiZkROWVVaWLFrrKZpuQOhzxBoJqwbZQek__lI9PCLLRgnErXqdjCr&s=a4a51055eeff8ec73b211a5b5f7ac82627b452f26ba3eb82bc0912ec73bdedb9&u=https%3A%2F%2Fwww.3ds.com%2Fprivacy-policy%2Fcontact%2F

--
Stéphane Mottelet
Ingénieur de recherche HDR
EA 4297 Transformations Intégrées de la Matière Renouvelable
Département Génie des Procédés Industriels
Sorbonne Universités - Université de Technologie de Compiègne
CS 60319, 60203 Compiègne cedex
Tel : +33(0)344234688
http://www.utc.fr/~mottelet

_______________________________________________
users mailing list - users@lists.scilab.org
Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org>
https://lists.scilab.org/mailman/listinfo/users
This email and any attachments are intended solely for the use of the 
individual or entity to whom it is addressed and may be confidential and/or 
privileged.

If you are not one of the named recipients or have received this email in error,

(i) you should not read, disclose, or copy it,

(ii) please notify sender of your receipt by reply email and delete this email 
and all attachments,

(iii) Dassault Systèmes does not accept or assume any liability or 
responsibility for any use of or reliance on this email.


Please be informed that your personal data are processed according to our data 
privacy policy as described on our website. Should you have any questions 
related to personal data protection, please contact 3DS Data Protection Officer 
https://www.3ds.com/privacy-policy/contact/

Reply via email to