写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

动态图的实现思路:

  1. 生成新的xy坐标值
  2. 暂时保存
  3. 重新绘点

视频:

matlab实现布朗运动


参考文章:https://blog.csdn.net/mygod2008ok/article/details/100855010

MatLab实现布朗运动相关推荐

  1. Matlab:利用Matlab实现布朗运动模拟

    Matlab:利用Matlab实现布朗运动模拟 目录 输出结果 实现代码 输出结果 实现代码 n=40 s=0.04 x=rand(n,1)-0.5; y=rand(n,1)-0.5; h=plot( ...

  2. 【MATLAB】布朗运动动画仿真

    目录 1.题目要求 2.代码实现 3.仿真示例 1.题目要求 先确定布朗运动的点数 n 和一个温度 s(或速度).比如 n=20, s=0.002.达到最好的动画效果的 n. s 与使用的计算机的性能 ...

  3. Matlab:利用Matlab编程实现模拟分子布朗运动的动画展示

    Matlab:利用Matlab编程实现模拟分子布朗运动的动画展示 目录 输出结果 实现代码 输出结果 实现代码 %Brownian motion clf; n=20; s=0.02; x = rand ...

  4. Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动

    Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 目录 输出结果 设计代码 输出结果 基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 设计代码 %基于Matlab实 ...

  5. Matlab动画模拟分子布朗运动的示例

    Matlab动画模拟分子布朗运动的示例 %Brownian motion clf; n=20; s=0.02; x = rand(n,1)-0.5; y = rand(n,1)-0.5; h = pl ...

  6. matlab模拟分数布朗运动,分数布朗运动及其模拟FractionalBrownianMotion科研菜鸟http.PDF...

    分数布朗运动及其模拟FractionalBrownianMotion科研菜鸟http 分数布朗运动及其模拟(Fractional Brownian Motion ) 科研菜鸟 /u/sanshiphy ...

  7. Matlab模拟分子的布朗运动

    文章目录 1. 效果(动图) 2. 代码 1. 效果(动图) 2. 代码 %Brownian motion 布朗运动 clc,clear,close all; %清除命令,清空工作区,关闭所有窗口n= ...

  8. 【Matlab】构造布朗运动

    1.布朗运动B(t)的定义 1.定义: {B(t)}布朗运动(brownian motion)也称为维纳过程,是一个随机过程,如果满足以下性质: B(0) = 0: 独立平稳增量:对任意t>s, ...

  9. 欧式期权matlab编码蒙特卡罗方法,基于MATLAB的金融工程方法与实践第七章 基于蒙特卡洛方法的期权定价.ppt...

    蒙特卡洛的优缺点 MCMC方法的优点 分布假设更一般,描述市场因素可能变化的统计分布既可以是正态.对数正态的,也可以是带跳的扩散分布.t分布等. 随机生成风险因素的各种各样的未来假想情景,可在模型中融 ...

最新文章

  1. 最新电子皮肤的触觉有多灵?连空气流动都能感受到
  2. XAML Namespace http://schemas.microsoft.com/expression/blend/2008 is not resolved
  3. 前端iframe 能指定本地网页吗_微前端的技术拆分方式
  4. python分词和词频统计
  5. display:flex弹性布局
  6. 平板电脑可以插u盘吗_电视TV盒子安装app的六个方法,u盘/电脑/手机都可以安装...
  7. Spark 调用 hive使用动态分区插入数据
  8. 莫烦python博客_《莫烦Python》笔记 -- numpy部分
  9. Python 实现针对时间序列预测的特征选择
  10. python 路径拼接字符串_字符串游戏之拼接字符串
  11. .net Remoting的工作原理
  12. 二或四 通道USB数据采集卡如合通过CMI耦合,构成USB高速多通道数据采集卡
  13. [源码]UnicodeTOGB,能够将Unicode串转换成GB码,方便开发。
  14. 《内功修炼系列》之-递归从入门到入土(收藏版)
  15. Logback 学习笔记
  16. ASUS eeepc-1005HA安装archlinux后耳机无声解决方案
  17. 基于arduino uno 开发板制作的 智能书包项目
  18. QuickFIX 源码分析
  19. TRIZ创新方法——物质-场模型与标准解系统
  20. 创新领导力|学习华为蓝军思维,价值100万的“小报告”如何修正老板的决议?

热门文章

  1. 算法提高课-动态规划-树形DP-AcWing 1072. 树的最长路径:dfs写法
  2. PAT甲级1081 Rational Sum:[C++题解]分数求和、辗转相除法求最大公约数、long long有一个数据溢出
  3. PAT甲级1123 Is It a Complete AVL Tree (30分):[C++题解]建立平衡树、bfs,判完全二叉树
  4. 塔菲尔曲线斜率的大小_分段函数函数值(自变量)大小(范围)的六种题型
  5. 零基础学习pythonrequests_requests库初体验
  6. tikhonov正则化 matlab_4 L1和l2正则化详解(花书7.1 参数范数惩罚)
  7. kettle使用数据库来生成序列_Kettle的安装和使用
  8. 华为p8升级android8系统资源,华为鸿蒙系统已陆续推送! 安卓系统可无缝升级: 升级包容量高达6GB...
  9. 英特尔cpu发布时间表_10纳米来了:英特尔十一代酷睿上市,性能提升20%,AI算力乘5倍...
  10. php自动生成新闻页,自动发布新闻页面的php代码