这是生物学家matlab入门课程的一部分.我有一个矩阵的数据点(对于单个粒子!)有4列(时间,x,y,z)和几千行.我想要做的是使用所有时间步长的xyz坐标计算粒子的均方位移. MSD定义为MSD = average(r(t)-r(0))^ 2其中r(t)是粒子在时间t的位置,r(0)是初始位置,所以在某种意义上它是粒子在时间间隔t内行进的距离.这就是我到目前为止所拥有的.我显然是一个MATLAB新手,我非常感谢som帮助/输入.

dx=zeros(length(data),1); %Create space for displacements.

dy=zeros(length(data),1);

dz=zeros(length(data),1);

for i=1:length(pos)

dx(i)=data(i+1,2)-data(1,2); %Calculate the distance at each time step

dy(i)=data(i+1,3)-data(1,3); %back to the origin.

dz(i)=data(i+1,4)-data(1,4);

end

我接下来需要做的是计算这些值的均方.但令我担心的一件事是,如果粒子的初始位置保持不变,那么数据矩阵中的大部分信息会被忽略.我认为沿着这些方向的算法可以更好地使用(更准确,更好的统计数据)数据,但我不是一个经验丰富的程序员来实现它…

注意:dt = delta t =时间步长,Pos =位置

>计算Pos 2和Pos 1之间的位移(dt = 1)

>计算Pos 3和Pos 1之间的位移(dt = 2)

>计算Pos 4和Pos 1之间的位移(dt = 3),依此类推……

现在将每个位置与Pos 2进行比较.

>计算Pos 3和Pos 2之间的位移(dt = 1)

>计算Pos 4和Pos 2之间的位移(dt = 2)

>计算Pos 5和Pos 2之间的位移(dt = 3)等等……

现在将每个位置与Pos 3进行比较

>计算Pos 4和Pos 3之间的位移(dt = 1)

>计算Pos 5和Pos 3之间的位移(dt = 2)

>计算Pos 6和Pos 4之间的位移(dt = 3),依此类推……

位置4 ……

>计算Pos 5和Pos 4之间的位移(dt = 1)

>计算Pos 6和Pos 4(dt = 2)…之间的位移数千行.

这样,每个dt会有更多的数据点.对其他人这样做是否有意义?这会使计算减少“嘈杂”(即使msd与时间的关系曲线更平滑吗?)也许有一些方便的matlab函数可以使这个问题更容易?我真的很感激一些意见.

matlab怎么算均方位移,用MATLAB计算均方位移(msd)相关推荐

  1. matlab状态空间法算反馈阵,matlab中已知系统的状态方程怎样绘制系统阶跃响应曲线...

    matlab中,怎样设置坐标系统的原点在左上角? axisijx=-pi:0.1:pi;y=sin(x);plot(x,y,'g:'); 谁能帮我做Matlab的作业:分析系统能否通过状态反馈实现镇定 ...

  2. matlab能否算高中数学题,Matlab软件在高中数学学习中的应用

    强雨筱 [摘要]运用Matlab软件强大的计算及图像处理功能求解高中数学教材中2种典型例题.首先基于二分法,求解了一般非线性方程;其次基于随机模拟方法,估计了圆周率π的值.通过例题强化学生对Matla ...

  3. matlab解算平差实例,MATLAB软件在测量平差解算中的应用

    MAT LAB 软件在测量平差解算中的应用 胡远新,赵奋军 (浙江省第七地质大队, 浙江丽水市 323000) 摘 要:阐述了4种经典的测量平差的解算过程,并结合平差实例,说明了MAT LAB 在测 ...

  4. 轨迹分析—Matlab计算均方位移

    关注 Mr.material,\color{Violet} \rm Mr.material\ ,Mr.material , 更\color{red}{更}更多\color{blue}{多}多精\col ...

  5. matlab 简单算例,(简单算例)基于Matlab的电力系统潮流编程计算.pdf

    (简单算例)基于Matlab的电力系统潮流编程计算 基于Matlab的电力系统潮流编程计算 口黄扬威吴喜春郭志峰张斯翔 (三峡大学电气与新能源学院湖北·宜昌443002) 摘要:通过介绍电力系统的实际 ...

  6. matlab算线性方程解,MATLAB计算方法3解线性方程组计算解法.pptx

    第三章线性方程组数值解法解线性方程组 §3.1 直接法一. Gauss 消去法设 有消 元: 用Matlab实现顺序Gauss消去法在Matlab程序编辑器中输入:function x=nagauss ...

  7. 5种方法计算均方位移MSD

    5种方法计算均方位移MSD (1)lammps计算 (2)VMD计算 (3)OVITO计算 (4)ISAACS软件计算 (5)自编程计算-附MATLAB代码 总结 在统计力学中,均方位移(MSD,均方 ...

  8. Matlab解算空间后方交会外方位元素

    目录 1 问题描述 2 实现部分 参考文献   本问题为武汉大学摄影测量学教材课后习题,现在用MATLAB实现,供大家学习参考. 1 问题描述   已知摄像机主距f=153.24mm,四对点的像点坐标 ...

  9. 【椭球大地测量学】Python及MATLAB实现贝塞尔大地问题正反解计算编程(含流程图)

    目录 一.功能 二.使用方法 三.算法流程图 (1)贝塞尔大地问题正算 (2)贝塞尔大地问题反算 四.Python程序 (1)计算程序 (2)DOS接口程序 五.MATLAB程序 (1)贝塞尔大地问题 ...

  10. 有限元悬臂梁matlab,悬臂梁MATLAB有限元算例注释.doc

    悬臂梁MATLAB有限元算例注释 用有限元法对悬臂梁分析的算例 算例: 如下图所示的悬臂梁,受均布载荷q=1N/mm2作用.E=2.1×105N/mm2, μ=0.3厚度h=10mm.现用有限元法分析 ...

最新文章

  1. C#创建Windows服务
  2. 【机器学习基础】深入讨论机器学习 8 大回归模型的基本原理以及差异!
  3. C# GDI+ 画坐标(x,y)
  4. 解决iPhone网络软件在睡眠情况断线问题
  5. HEC-RAS二维溃坝洪水模拟(尾矿库)
  6. 十五天数据分析-笔记
  7. C++结构体,联合体与枚举类
  8. 苹果电脑如何查看本机IP地址和DNS?
  9. Hadoop环境搭建 MYSQL环境配置
  10. 离散数学编程实现真值表输出
  11. PHP编程学习之路 2
  12. Java学习决心计划书
  13. Voldemort启动失败案例一则
  14. cmd介绍及常用命令
  15. 支持在线大数据SQL查询平台开源项目
  16. python实现SVG文件解析
  17. 打破FTTR布线标准瓶颈,让千兆光纤网络走进千家万户
  18. qt 添加Q_OBJECT报错
  19. 又一波网红餐饮店惨死,餐饮+新零售或将产生新格局
  20. werkzeug实现WSGI Application

热门文章

  1. HUAWEI MateBook X Pro 2022款 评测
  2. 计算机重启动桌面干净,电脑不断刷新重启桌面
  3. 小狼毫输入法标点符号直接上屏,反斜杠(\)输入顿号(、) BY 冷家锋、刘亭亭
  4. 如何制作抖音照片图集 照片视频怎么制作
  5. python3.7安装pycrypto
  6. Python-Pandas-Excel/CSV 数据处理大全整理 (一)
  7. 微信小程序的云开发以及与传统开发的比较
  8. python花猫_Python猫荐书系列之五:Python高性能编程
  9. scrcpy手机投屏电脑教程,可操控手机,低延迟
  10. 树莓派文件系统查看与清理垃圾箱