absT123321-CSDN博客

转载 MATLAB 制作gif动画

1、已有一系列时序图片的基础上，可借助matlab制作gif图片，虽然比不上专业绘图制作的精良，好在比较方便。2、假定时序图片的命名为：dens000000.tiff,dens000500,dens0010000,……3、效果展示4、代码呈现% make gifclc;clear;close all;% fileneme for gifpref...

2019-03-26 10:36:00 609

转载 MATLAB求解浸润角

1、概述　　通过后处理软件获得流场数据的密度云图，做出附着于壁面上液滴或气泡在接触点上的切线，测量可以获得浸润角度，手工做切线存在较大的自由度，所得的角度往往并不精确，可以考虑用圆拟合液滴或者气泡，通过拟合出的圆半径及圆心位置，则可以精确的计算出相应的角度，相比做切向的方法，准确度有较大的提升，但是如果还是通过手工来找准这个拟合圆，则数据准确性也受影响，可以考虑采用程序，通过图像处理...

2019-01-06 13:43:00 1015

转载 MATLAB 解代数方程组

今天用matlab给了我三重惊喜，简直打开了新世界的大门：1、虽然知道matlab有内置的符号工具箱，但以往用的很少，直到今天，需要求解一个方程组，方程本身到不是多么复杂，只不过变量众多，非常的恶心，手工求解出错可能性非常高，尝试了下用matlab来求解，求解起来是如此的顺利，这超出我的预料；2、在live Script 模式下，符号变量的显示非常自然，类手写的显示方式，help...

2018-11-05 19:22:00 3056

转载空间直线同任意形状椭球交点

空间直线同空间中三维椭球相交，其交点即为空间直线方程同椭球方程的解，对于空间直线方程，只要知道两点空间坐标即可，而欧拉角不为零的三维椭球方程则较难描述，但可以考虑对椭球进行变换，使其欧拉角为零，进而转化为标准椭球，标准椭球方程则容易描述，相应在对椭球进行变换的同时也许对空间直线进行相应的变换。由此，通过标准椭球方程同变换后的直线方程进行联立求解，即可获取交点坐标，下面简述求解过程及列出主...

2018-05-26 20:09:00 801

转载空间旋转——四元数表示

1、旋转的四元数表示空间中的旋转可用一个四元数来表述，如点 $P(\,x,y,z\,)$ 绕向量 $V(\,u,v,w\,)$旋转，此旋转过程可表示为四元数Q $[\,\,cos( \frac{\theta}{2} ),\,sin( \frac{\theta}{2} )\cdot (u,v,w)\,]$。2、基本旋转的四元数表示空间中的三维旋转可视为绕三个基本轴的旋...

2018-05-26 15:20:00 959

转载坐标系旋转 & 坐标点旋转

1、坐标系旋转坐标系逆时针旋转角度 θ后，原坐标系下坐标变换为新坐标系下坐标，示意图见图1，变换公式如下：2、坐标点旋转坐标系下的点 P 绕圆心逆时针旋转 θ后，变换后坐标点 P’，示意图见图2，坐标点变换公式如下：3、总结坐标系旋转贺坐标点旋转两者含义不同，但变换公式非常接近，仅仅相差一个负号，如不注意，极容易搞混。转载于:https://ww...

2018-05-24 19:11:00 3077

转载展开循环来加速计算

对于循环次数较为有限的for循环，可通过手工展开该循环达到消除此循环的目的，达到加速运行的目的，下面给出一个代码段实例。1、for循环体1 // Q为常数2 for ( k = 0; k < Q; k++ ) {3 for ( int j = 0; j < Q; j++ ) {4 m[k] += fp[z][y][x][j] *...

2018-05-12 10:54:00 491

转载空间直线同球体交点求解

一、问题求解空间直线同球体交点问题求解，为求解二元二次方程问题，下面是求解过程二、代码实现1、代码实现 1 void cal_q ( double par_x, double par_y, double par_z, double par_r, double bpx, double bpy, double bpz, double spx, doubl...

2018-05-08 10:06:00 3665

转载变量初始化问题

代码段：1 double delta_fx;2 if ( distance <= ksi ) {3 delta_fx = eta_p / C_eta * ( x1 - x2 ) * ( ksi - distance ) * ( ksi - distance );4 }5 tor += - delta_fx * ( y1 - par_y[num] );...

2018-04-09 22:38:00 155

转载自转、公转模型

原以为硕士期间做的叶片造型和现在研究的颗粒什么的完全不相关，没想到还是有些共同的的地方，就比如颗粒的运动过程中的坐标跟踪，而且比较起来，叶片造型反而更复杂些。1、椭圆粒子自转2、椭圆粒子公转3、椭圆粒子自转 + 公转转载于:https://www.cnblogs.com/kljfdsa/p/8662405.html...

2018-03-28 10:40:00 585

转载 Centos 下配置 SSH

一、设置步骤1、安装ssh，不过系统默认自带ssh服务：yum install-y openssl openssh-server2、修改ssh 配置文件，具体配置见（二、配置“/etc/ssh/sshd_config”文件）vim /etc/ssh/sshd_config3、启动ssh的服务：systemctl start sshd.service4、设置ssh开机自启sys...

2018-03-20 21:58:00 148

转载不可压格子波尔兹曼模型-D2G9

D2Q9模型是一种弱可压模型，郭老师在此基础上提出了D2G9模型，是一种不可压模型，两个模型主要有以下不同点：控制方程：1、平衡态方程定义不同对于D2Q9：，对于D2G9：，，，其中Si(u)定义为：2、演化方式一个基于密度演化（D2Q9），一个基于压力演化（D2G9），体现在宏观量求解量不同，且定义不同对于D2Q9：，对于D2G9...

2017-12-04 20:21:00 869

转载 MATLAB读取TECPLOT笛卡尔网格三维流场数据

基于笛卡尔网格的三维流动数值模拟，其流场信息可以通过tecplot格式进行输出，方便tecplot对流场进行可视化处理，但对数据进行在加工时，还是导入matlab中比较方便，那么对于一个tecplot数据文件，matlab是不能直接读取的，必须有一个函数将tecplot数据文件中数据转换成matlab便于操作的数据格式。tecplot数据文件前三行是文件头，其中第一行是数据文件说明，...

2017-11-16 18:47:00 3989

转载 c++链表简单功能实现

这周主要学习了链表这种数据结构，基本实现了链表的各项功能，也算是对类的创建以及指针的操作有了进一步的了解，下面给出实现的链表的代码。1、链表节点#ifndef NODE_H_INCLUDED#define NODE_H_INCLUDED#include <iostream>#include <iomanip>using name...

2017-11-03 22:38:00 184

转载 vector 学习

#include <iostream>#include <iomanip>#include <vector>#include <algorithm>using namespace std;void printarray ( const vector <int> & );int ma...

2017-10-28 21:22:00 157

转载格子方法曲边处理

　　应用格子玻尔兹曼方法做气固两相流时，对于颗粒相，借助相关作用力模型，可视为点源处理，也可以直接解析颗粒，精确刻画颗粒气体间相互作用力。全解析颗粒时，一般可将颗粒视为球体（三维）或者圆（二维），气体颗粒间相互作用体现在颗粒边界的的处理上，对于这样的曲边界，其处理方法的好坏直接影响计算的精度和稳定性，对于无滑移的曲边界，有将曲边近似为阶梯状的直接反弹法、在其基础上改进的插值反弹法，及基于...

2017-08-27 21:47:00 360

absT123321的博客