matlab中提取微分方程某点的值,从随机微分方程解中提取某些参数的值

我在matlab中求解随机微分方程。

例如：

考虑随机微分方程

dx=k A(x,t)dt+ B(x,t)dW(t)

其中k是常数，A和B是函数，dW(t)是Wiener过程。

我在[0,20]中为所有t绘制解决方案。我们知道dW(t)是随机生成的。我的问题是：我想知道特定的t值和特定子区间的A(x，t)，B(x，t)，dW(t)的值，比如说[3,6]。我可以使用Matlab中的哪个命令？

以下是我根据D.Higham的论文使用的代码：

clear all

close all

t0 = 0; % start time of simulation

tend = 20; % end time

m=2^9; %number of steps in each Brownian path

deltat= tend/m; % time increment for each Brownian path

D=0.1; %diffsuion

R=4;

dt = R*deltat;

dW=sqrt( deltat)*randn(2,m);

theta0=pi*rand(1);

phi0=2*pi*rand(1);

P_initial=[ theta0; phi0];

L = m/ R;

pem=zeros(2,L);

EM_rescale=zeros(2,L);

ptemp=P_initial;

for j=1:L

Winc = sum(dW(:,[ R*(j-1)+1: R*j]),2);

theta=ptemp(1);% updating theta

phi=ptemp(2); % updating phi

%psi=ptemp(3); % updating psi

A=[ D.*cot(theta);...

0];% updating the drift

B=[sqrt(D) 0 ;...

0 sqrt(D)./sin(theta) ]; %% updating the diffusion function

ptemp=ptemp+ dt*A+B*Winc;

pem(1,j)=ptemp(1);%store theta

pem(2,j)=ptemp(2);%store phi

EM_rescale(1,j)=mod(pem(1,j),pi); % re-scale theta

EM_rescale(2,j)=mod(pem(2,j),2*pi); % re-scale phi

end

plot([0:dt:tend],[P_initial,EM_rescale],'--*')

假设我想知道每个特定时间点或任何时间间隔的所有参数(包括随机：布朗)。怎么做？

matlab中提取微分方程某点的值,从随机微分方程解中提取某些参数的值相关推荐

python是值传递还是引用传递_Python里参数是值传递还是引用传递？
Python就是把一些参数从一个函数传递到另一个函数,从而使其执行相应的任务.但是你有没有想过,参数传递的底层是如何工作的,原理又是怎样的呢? 实际工作中,很多人会遇到这样的场景:写完了代码,一测试, ...
python中序列和列表区别细菌真菌病毒_生物信息中的Python 02 | 用biopython解析序列...
上一篇文章生物信息中的Python 01 | 从零开始处理基因序列自己造轮子实现了序列的基础操作,但是在Python的世界里,一项工作只要重复的次数多了,那么一定就会有大神来开发相应的包来解决,这个包 ...
ES6：Rest 参数和参数默认值【转】
Rest 参数通常,我们需要创建一个可变参数的函数,可变参数是指函数可以接受任意数量的参数.例如,String.prototype.concat 可以接受任何数量的字符串作为参数.使用 Rest 参 ...
matlab中函数或变量无法识别怎么办_用MATLAB巧解微分方程实例分析
点"考研竞赛数学"↑可每天"涨姿势"哦! MATLAB巧解微分方程实例分析王少华西安电子科技大学微分方程求解难, 字母一堆看着烦. 写错数字一时爽, 一直 ...
组数(【问题描述】输入一行字符串（设字符数不大于80），提取该字符串中的数字字符并组成一个整数，输出该整数及其两倍的值。)
组数 [问题描述]输入一行字符串(设字符数不大于80),提取该字符串中的数字字符并组成一个整数,输出该整数及其两倍的值.要求在主函数中输入字符串,并输出结果.在子函数中提取该字符串中的数字字符并组成一 ...
二自由度振动仿真：matlab直接解微分方程virtual.lab motion仿真
一.二自由度振动解析法与多种数值算法见:https://blog.csdn.net/lijil168/article/details/67647924 二.用matlab符号建立微分方程,并用龙格库塔 ...
用MATLAB求解非线性微分方程
总结一下MATLAB中求解微分方程的思路和步骤.固然,网上很多关于此类的技术型文章,但往往一看下来发现,文章中的友情链接比文章字数还多,要了解这一篇文章,你要先了解那个:要了解那个,你又要了解那个那个 ...
Matlab进行gnss用户坐标计算,论文推荐 | 嵇昆浦，沈云中：含缺值GNSS基准站坐标序列的非插值小波分析与信号提取...
关键词:GNSS坐标时间序列缺失数据小波变换信号提取 Dyadic wavelet transform and signal extraction of GNSS coordinate time ...
【Matlab 控制】微分方程 ode45() 求解并绘制曲线
Matlab 微分方程 ode45 求解并绘制曲线 2. 用 ode45() 求解 2.1 ode45() 函数用法 2.2 示例:求解一阶微分方程 2.2.1 Matlab 代码如下 2.2.2 代 ...
matlab 三维全息谱,基于MATLAB软件平台的全息谱图及在风机振动故障诊断中应用的研究...
一.引言风机是化工机械中的关键设备.由于机器在运转过程中,受到外界的作用和机器本身零件之间的相互作用,不可避免地要产生振动.在产生振动的诸多因素中,机械部件本身的缺陷占有相当大的比重.由于分析方法的 ...

matlab中提取微分方程某点的值,从随机微分方程解中提取某些参数的值

matlab中提取微分方程某点的值,从随机微分方程解中提取某些参数的值相关推荐

最新文章

热门文章