随机微分方程的MATLAB数值求解

在这里插入图片描述
在这里插入图片描述

dt=0.01;
tout=200; %总时间为2
x=zeros(1,tout);
x(1)=0.5; %初始位置
mu=0.2; sigma=1;
Wt=sqrt(dt)*randn(1,tout); %产生随机序列Wt
for t=1:tout-1x(t+1)=x(t)+mu*x(t)*dt+sigma*x(t)*Wt(t);
end
t1=1:10:tout; %对原时间序列进行抽样
xt=zeros(1,length(t1));
i=1;
for t=t1xt(i)=0.5*exp((mu-0.5*sigma^2).*dt*t+ ...sigma*sum(Wt(1:t))); %使用理论解对原序列进行抽样计算i=i+1;
end
plot(dt*(1:tout),x,LineWidth=1.5)
hold on
plot(dt*t1,xt,LineWidth=1.5)
hold off
xlabel("$t$",FontSize=20,Interpreter="latex")
dt=0.01;
tout=200;
n=10000;
tends=zeros(1,n);
x=zeros(1,tout);
x(1)=0.1;
mu=0.2; sigma=1;
for k=1:nWt=sqrt(dt)*randn(1,tout); %产生随机序列Wtfor t=1:tout-1x(t+1)=x(t)+mu*x(t)*dt+sigma*x(t)*Wt(t);endtends(k)=x(end);
end
histogram(tends,200)
dt=0.01;
tout=200;
x=zeros(1,tout);
x(1)=0; %初始位置
mu=5; sigma=2;
for i=1:5Wt=sqrt(dt)*randn(1,tout); %产生随机序列Wtfor t=1:tout-1x(t+1)=x(t)-mu*x(t)*dt+sigma*Wt(t);endplot(dt*(1:tout),x,LineWidth=1.5)hold on
end
xlabel("$t$",FontSize=20,Interpreter="latex")
hold off
dt=0.01;
tout=200;
x=zeros(1,tout);
x(1)=0; %初始位置
mu=5; sigma=2;
tends=zeros(1,100000);
for i=1:100000Wt=sqrt(dt)*randn(1,tout); %产生随机序列Wtfor t=1:tout-1x(t+1)=x(t)-mu*x(t)*dt+sigma*Wt(t); %Langevin方程endtends(i)=x(end);
end
s1=sqrt((1-exp(-2*mu*dt*tout))/(2*mu))*sigma;
histogram(tends,Normalization="pdf")
hold on 
plot(-3:0.1:3,1/sqrt(2*pi*s1^2)* ...exp(-(-3:0.1:3).^2/(2*s1^2)),LineWidth=1.5) %理论分布
hold off
dt=0.01;
tout=100;
x=zeros(1,tout); 
x(1)=0; %初始位置
y=x;for i=1:20Wtx=sqrt(dt)*randn(1,tout); %产生随机序列WtxWty=sqrt(dt)*randn(1,tout); %产生随机序列Wtyfor t=1:tout-1x(t+1)=x(t)+Wtx(t); %布朗运动y(t+1)=y(t)+Wty(t);endpatch([x,NaN],[y,NaN],[1:tout,NaN],LineWidth=1.5,EdgeColor='interp')colormap('jet')
end
xlabel("$x$",FontSize=20,Interpreter="latex")
ylabel("$y$",FontSize=20,Interpreter="latex")

^Sauer T. Computational solution of stochastic differential equations[J]. Wiley Interdisciplinary Reviews: Computational Statistics, 2013, 5(5): 362-371. https://wires.onlinelibrary.wiley.com/doi/full/10.1002/wics.1272

详见知乎:https://zhuanlan.zhihu.com/p/666480835?utm_campaign=&utm_medium=social&utm_oi=1315073218793488384&utm_psn=1710572578152996864&utm_source=qq

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/208574.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Missing file libarclite_iphoneos.a 问题解决方案

问题 在Xcode 运行项目会报以下错误 File not found: /Applications/Xcode-beta.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/usr/lib/arc/libarclite_iphoneos.a解决方案 打开URL https://github.com/kamyarelyasi/Libarclite-Files ,下载liba…

校园圈子论坛,交友,帖子内短视频,二手市场,APP小程序H5三端交付,源码交付,支持二开

校园圈子论坛,交友频道,商城,二手市场,活动专区,短视频,从校园生活的方方面面展现出了充满活力和创造力的镜头。这个频道是一个让学生们相互交流、结识新朋友的平台,不仅有交友功能,…

系列八、key是弱引用,gc垃圾回收时会影响ThreadLocal正常工作吗

一、key是弱引用,gc垃圾回收时会影响ThreadLocal正常工作吗 到这里,有些小伙伴可能有疑问,ThreadLocalMap的key既然是 弱引用,那么GC时会不会贸然地把key回收掉,进而影响ThreadLocal的正常使用呢?答案是不会…

office word 使用笔记

office word 使用笔记 1. 功能1.1 格式快捷键1.2 复选框 2 遇到过的问题2.1 表格标题和表格距离过大 1. 功能 1.1 格式快捷键 复制格式:ctrl shift c 粘贴格式:ctrl shift v 1.2 复选框 方框位置和类型:“插入——高级符号——字体”选…

阿里云发送短信

官方代码如下: // This file is auto-generated, dont edit it. Thanks. package com.aliyun.sample;import com.aliyun.tea.*;public class Sample {/*** 使用AK&SK初始化账号Client* param accessKeyId* param accessKeySecret* return Client* throws Excep…

基于C#实现树状数组

有一种数据结构是神奇的,神秘的,它展现了位运算与数组结合的神奇魅力,太牛逼的,它就是树状数组,这种数据结构不是神人是发现不了的。 一、概序 假如我现在有个需求,就是要频繁的求数组的前 n 项和&#x…

黑马React18: ReactRouter

黑马React: ReactRouter Date: November 21, 2023 Sum: React路由基础、路由导航、导航传参、嵌套路由配置 路由快速上手 1. 什么是前端路由 一个路径 path 对应一个组件 component 当我们在浏览器中访问一个 path 的时候,path 对应的组件会在页面中进行渲染 2. …

如何在Ubuntu的Linux系统中安装MySQL5.7数据库

前往MySQL数据库官网链接地址下载5.7数据库。 MySQL :: Download MySQL Community Server (Archived Versions)使用ssh的可视化工具将下载的mysql-5.7.40-linux-glibc2.12-x86_64.tar.gz文件上传到Linux服务器,并解压文件 tar -zxvf mysql-5.7.40-linux-glibc2.12-x…

lv11 嵌入式开发 FS4412开发环境搭建 8

目录 1 开发板硬件资源介绍 2 初始原理图 3 交叉开发环境搭建 3.1 在ubuntu下安装交叉编译工具链 3.2 在windows下安装SecureCRT 3.3 在windows下安装USB转串口驱动 3.4 下载测试程序 4 硬件控制原理 4.1 地址映射表 4.2 查看数据手册 1 开发板硬件资源介绍 …

外网IP和内网ip的区别

首先得先知道什么是ip地址,它就是唯一标识连接网络的设备的,即IP地址充当了设备在网络中的“住址”,使得设备能够相互通信和交换数据。 我们常听开发人员说外网内网,那么它们有什么区别呢? 外网可以理解为互联网&…

开发上门送桶装水小程序要考虑哪些业务场景

上门送水业务已经有很长一段时间了,但是最开始都是给用户发名片、贴小广告,然后客户电话订水,水站工作人员再上门去送,这种人工记单和派单效率并不高,并且电话沟通中也比较容易出现偏差,那么根据这个情况就…

二分查找——经典题目合集

文章目录 🦜69. x 的平方根🌼题目🌻算法原理🌷代码实现 🐳35. 搜索插入位置🌼题目🌻算法原理🌷代码实现 🦭852. 山脉数组的峰顶索引🌼题目🌻算法原…