MatLab实现布朗运动
写MATLAB程序,画出动画来模拟粒子的布朗运动。
在二维坐标系下,假设有100个点,起始坐标都为(0,0)
,在坐标系下用.
画出这100个点。模拟这些点的布朗运动,试写出含有1000步的循环,每一步循环将计算出每个点下一步的新的x和y坐标,然后画出新位置下的这100个点,这样循环1000步,就会画出有1000个图片的动画。每个点下一步的新坐标位置,通过在当前的位置(x,y)
,对x和y分别加上服从正态分布且标准差为0.005的变量来完成。
%Brownian motion 布朗运动%清除命令,清空工作区,关闭所有窗口
clc,clear,close all; % 设置正态分布的均值
mean = 0; % 点的个数
n = 100;% 速度
s = 0.01;% x, y的初始值
x = zeros([n, 1]);
y = zeros([n, 1]);
h = plot(x, y, '.');title('模拟布朗运动');
for i = 1:1000x = x + normrnd(mean, 0.005, [n,1]) * s;y = y + normrnd(mean, 0.005, [n,1]) * s;set(h,'XData',x,'YData',y) % 重新绘点drawnow;
end
动态图的实现思路:
- 生成新的
xy
坐标值 - 暂时保存
- 重新绘点
视频:
matlab实现布朗运动
参考文章:https://blog.csdn.net/mygod2008ok/article/details/100855010
MatLab实现布朗运动相关推荐
- Matlab:利用Matlab实现布朗运动模拟
Matlab:利用Matlab实现布朗运动模拟 目录 输出结果 实现代码 输出结果 实现代码 n=40 s=0.04 x=rand(n,1)-0.5; y=rand(n,1)-0.5; h=plot( ...
- 【MATLAB】布朗运动动画仿真
目录 1.题目要求 2.代码实现 3.仿真示例 1.题目要求 先确定布朗运动的点数 n 和一个温度 s(或速度).比如 n=20, s=0.002.达到最好的动画效果的 n. s 与使用的计算机的性能 ...
- Matlab:利用Matlab编程实现模拟分子布朗运动的动画展示
Matlab:利用Matlab编程实现模拟分子布朗运动的动画展示 目录 输出结果 实现代码 输出结果 实现代码 %Brownian motion clf; n=20; s=0.02; x = rand ...
- Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动
Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 目录 输出结果 设计代码 输出结果 基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 设计代码 %基于Matlab实 ...
- Matlab动画模拟分子布朗运动的示例
Matlab动画模拟分子布朗运动的示例 %Brownian motion clf; n=20; s=0.02; x = rand(n,1)-0.5; y = rand(n,1)-0.5; h = pl ...
- matlab模拟分数布朗运动,分数布朗运动及其模拟FractionalBrownianMotion科研菜鸟http.PDF...
分数布朗运动及其模拟FractionalBrownianMotion科研菜鸟http 分数布朗运动及其模拟(Fractional Brownian Motion ) 科研菜鸟 /u/sanshiphy ...
- Matlab模拟分子的布朗运动
文章目录 1. 效果(动图) 2. 代码 1. 效果(动图) 2. 代码 %Brownian motion 布朗运动 clc,clear,close all; %清除命令,清空工作区,关闭所有窗口n= ...
- 【Matlab】构造布朗运动
1.布朗运动B(t)的定义 1.定义: {B(t)}布朗运动(brownian motion)也称为维纳过程,是一个随机过程,如果满足以下性质: B(0) = 0: 独立平稳增量:对任意t>s, ...
- 欧式期权matlab编码蒙特卡罗方法,基于MATLAB的金融工程方法与实践第七章 基于蒙特卡洛方法的期权定价.ppt...
蒙特卡洛的优缺点 MCMC方法的优点 分布假设更一般,描述市场因素可能变化的统计分布既可以是正态.对数正态的,也可以是带跳的扩散分布.t分布等. 随机生成风险因素的各种各样的未来假想情景,可在模型中融 ...
最新文章
- 最新电子皮肤的触觉有多灵?连空气流动都能感受到
- XAML Namespace http://schemas.microsoft.com/expression/blend/2008 is not resolved
- 前端iframe 能指定本地网页吗_微前端的技术拆分方式
- python分词和词频统计
- display:flex弹性布局
- 平板电脑可以插u盘吗_电视TV盒子安装app的六个方法,u盘/电脑/手机都可以安装...
- Spark 调用 hive使用动态分区插入数据
- 莫烦python博客_《莫烦Python》笔记 -- numpy部分
- Python 实现针对时间序列预测的特征选择
- python 路径拼接字符串_字符串游戏之拼接字符串
- .net Remoting的工作原理
- 二或四 通道USB数据采集卡如合通过CMI耦合,构成USB高速多通道数据采集卡
- [源码]UnicodeTOGB,能够将Unicode串转换成GB码,方便开发。
- 《内功修炼系列》之-递归从入门到入土(收藏版)
- Logback 学习笔记
- ASUS eeepc-1005HA安装archlinux后耳机无声解决方案
- 基于arduino uno 开发板制作的 智能书包项目
- QuickFIX 源码分析
- TRIZ创新方法——物质-场模型与标准解系统
- 创新领导力|学习华为蓝军思维,价值100万的“小报告”如何修正老板的决议?
热门文章
- 算法提高课-动态规划-树形DP-AcWing 1072. 树的最长路径:dfs写法
- PAT甲级1081 Rational Sum:[C++题解]分数求和、辗转相除法求最大公约数、long long有一个数据溢出
- PAT甲级1123 Is It a Complete AVL Tree (30分):[C++题解]建立平衡树、bfs,判完全二叉树
- 塔菲尔曲线斜率的大小_分段函数函数值(自变量)大小(范围)的六种题型
- 零基础学习pythonrequests_requests库初体验
- tikhonov正则化 matlab_4 L1和l2正则化详解(花书7.1 参数范数惩罚)
- kettle使用数据库来生成序列_Kettle的安装和使用
- 华为p8升级android8系统资源,华为鸿蒙系统已陆续推送! 安卓系统可无缝升级: 升级包容量高达6GB...
- 英特尔cpu发布时间表_10纳米来了:英特尔十一代酷睿上市,性能提升20%,AI算力乘5倍...
- php自动生成新闻页,自动发布新闻页面的php代码