資源簡(jiǎn)介
matlab開發(fā)-隨機(jī)微分方程解算。用于求LSDE前兩個(gè)矩的函數(shù)

代碼片段和文件信息
%?This?function?finds?out?the?first?two?moments?
%?of?a?linear?Stochastic?Differential?Equation?(LSDE)
%?LSDE?:?dy(t)?=?ay(t)+?b*sigma*dW(t)?where?W(t)?is?Wiener?Process
%?Also?dW(t)=n(t)dt?where?n(t)?is?Gaussian?White?noise
%?Inputs?:?a?b?sigma(?magnitude?of?sqrt?of?PSD)?and?
%??????????c?(?=y(0)?initial?conditions)
%??????????t?(Time?Vector)?
%?Output?:m(t)(?First?Moment?-?Mean)
%?????????M(t)?(Second?Moment)
%?Note?the?solution?statistics?shown?below?have?been?derived?using?Ito‘s
%?Calculus
%?Author:?abhiruplahiri@yahoo.com
function?[mM]?=?LSDE(absigmact)
m=c*exp(a*t);
q=b^2*sigma^2/(2*a)
p=c^2+q;;
M=-q?+?p*exp(2*a*t);??????%?derived?using?Ito‘s?Formula
%plot(tmtM)
%grid
???????
?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????文件?????????767??2018-07-11?20:20??LSDE.m
評(píng)論
共有 條評(píng)論