MATLAB
文章平均质量分 64
Hironono
on the way
展开
-
GNSS速度场简易MATLAB克里金插值
GNSS速度场简易MATLAB克里金插值原创 2022-01-14 21:19:46 · 2569 阅读 · 1 评论 -
MATLAB 利用SVD分解拟合三维平面
方法简述:利用SVD分解求拟合平面 拟合方程为ax+by+cz=d,约束条件为a2+b2+c2=1,目标是使得尽可能多的点在平面上。构建矩阵为AX=0,将A矩阵进行奇异值分解之后,最小奇异值对应的特征值向量即为拟合平面的系数向量。详细的原理可以查看:SVD分解平面拟合原理在这里我举的例子是利用余震的三维空间坐标拟合余震的可能发震平面。as=load('data.txt');%根据需要选择合适经纬度范围的数据点as1 = [];for ir = 1:m if (as(ir,1)>=原创 2020-09-16 17:59:10 · 3046 阅读 · 2 评论 -
MATLAB矩阵计算陌生命令总结(一)
在进行MATLAB矩阵科学计算时,会用到以下命令,由于不是很熟悉,我将其总结记录如下:Meshgrid生成二维网格,用法为:[x y]=meshgrid(a b);% a 和b是一维数组,如a=[1 2 3]; b= [2 3 4]; 则生成的 X 和 Y 都是为 3X3 维的矩阵,[x y]=meshgrid(a,b)x =1 2 31 2 31...原创 2018-11-17 09:55:47 · 292 阅读 · 0 评论 -
如何正确优雅地使用武大超算进行MATLAB计算
最近总会跑一些计算量大的matlab程序,然后PC就是灾难现场了,不仅风扇转个不停,声音巨大,而且基本上是不敢再开其他的应用了,cpu占用率惊人,加速球变红…(此时应该有捂脸的表情)这样折腾了几天之后,发现不行,对电脑的伤害太大了,虽然说,电脑就是用来用的,但是也还是得稍微珍惜一下的。之后我就想到了我们学校的超算中心,这么好的资源,不用白不用。于是我就开始了自己的探索过程。1,进入武汉大学超算...原创 2018-12-04 16:51:47 · 8099 阅读 · 5 评论 -
MATLAB中如何保存嵌套循环的中间变量值
%思想:先预置一个空矩阵,然后把每次循环的结果分别逐个加入这个空矩阵。% Purpose:Save intermediate variable values for nested loops% Programmer:YQH% Date:20181206clear all;tic;lamda1=[0:0.1:1];lamda2=1-lamda1;dip1=[70,80,90]; %...原创 2018-12-06 14:38:53 · 4798 阅读 · 0 评论 -
MATLAB中如何实现mesh三维图
当有两个自变量和一个因变量需要进行可视化时,可以使用mesh命令绘制双变量的三维图。绘制三维图以下是MATLAB官网中对于mesh命令的解读:mesh(X,Y,Z)mesh(X,Y,Z) 使用 Z 确定的颜色绘制线框网格,因此其颜色与曲面高度成比例,其中Z一定是矩阵,X,Y可以是矢量。如果 X 和 Y 为向量,length(X) = n 且 length(Y) = m,其中 [m,n]...原创 2018-12-10 15:03:27 · 50758 阅读 · 1 评论 -
如何正确优雅地使用武大超算进行MATLAB计算(续)
武大超算集群在2018年12月26日进行为期三天的停机升级,同时升级了存储系统。升级之后的超算换汤不换药,原来的配方基础上增加了性能,当然,在超算上运行MATLAB程序也要有所改变了。升级之后官方的swarm集群文件系统说明/home:用户主目录,主要存放用户的环境变量等信息,限额1G,长期保存;/project:用户项目文件、运行程序等,限额1T,长期保存;/workfs:用户数据文件,...原创 2019-01-03 16:44:30 · 2311 阅读 · 0 评论 -
模拟退火法
模拟退火法概述模拟退火法来源于固体退火原理,将固体加热至充分高,再让其慢慢冷却,加热时,固体内部粒子随温度升高变为无序状态,内能增大,而慢慢冷却时粒子逐渐有序,在每个温度都达到平衡态,最后再常温时达到基态,内能减为最小。根据Metropolis准则,粒子在温度T时趋于平衡的概率为 e-ΔE/(KT),其中E为温度T时的内能,ΔE为其改变量,K为Boltzmann常数。用固体退火模拟组合优化问题...原创 2019-01-21 20:53:02 · 530 阅读 · 0 评论