- 博客(16)
- 收藏
- 关注
RSS订阅转载 Matlab去除图像右边和上边的刻度线
一般情况下，使用matlab画图出来的效果是这样的：x = -pi:pi/100:pi;y = tan(sin(x)) - sin(tan(x));plot(x,y,'rs-','LineWidth',1,'MarkerEdgeColor','k',... 'MarkerFaceColor','g','MarkerSize',10);图像的右边和上边的刻度线很扰人，想要去掉刻度线，但是保留坐标线，可以采用新建坐标轴的方法：box offax2 = axes('Position',ge
2022-02-24 09:09:18
14483
原创 线性波传播至垂直壁面反射形成驻波的动画MATLAB程序
最近在写博士论文非常繁忙，很久没有过啦更新，这次带来个之前写的好玩的程序，可以帮助大家理解线性波传递以及驻波的形成。%% linear waves add%% code by JN_CUI%% Variables% eta - water surface profile% u - horizontal velocity of water particle% w - vortical velocity of water particle% ze
2021-11-22 21:00:01
1024
3
原创 Linux 操作常用命令
快捷打开终端ctrl + alt + t打开文件sudo /path/filename修改文件夹及其子文件夹权限为777chmod -R 777 /path/foldername强制删除文件夹及子文件夹rm -rf /path/foldername查看cpu使用量top持续更新中
2021-10-15 15:30:42
117
原创 MATLAB 快速傅里叶变换分析
快速傅里叶变换快速傅里叶变换原理有知乎大神的详细讲解，详细讲解传送门这里我主要介绍快速傅里叶变换的用法以及输出结果的含义，以随机海浪谱为例：随机波浪序列的生成：随机波浪模拟传送门当我们拥有一个频率为100Hz的随机海浪序列wsw_sws之后，想要通过快速傅里叶变换分析其频率组成：%% FFT analysis% code by JN_Cuifunction [S,F,Mag] = fft_analysis(x,frequency)Y = fft(x);N = length(Y)
2021-06-19 09:55:41
2858
11
原创 Matlab 四阶龙格库塔法求解二元常微分方程组
龙格库塔法是一种求解高阶常微分方程的常用方法，在工程当中应用广泛，例如求解物体的运动方程等。这里我们通过matlab程序编写龙格库塔算法求解二元常微分方程组，假设有常微分方程组：{x¨−x˙+2y¨+y˙=−2sint−3costx¨+y¨=−sint−costx(0)=0y(0)=1x˙(0)=1y˙(0)=0\left\{\begin{array}{lr}\ddot{x}-\dot{x}+2\ddot{y}+\dot{y}=-2\rm sin \it t - \rm 3\rm cos \it t
2021-03-04 15:40:27
13627
15
原创 MATLAB将图片输出视频及GIF文件
将图片输出为视频及GIFclear;clcPath = 'I:Path\';Folder = 'Foldername\';Data = imageDatastore([Path,Folder]);V = VideoWriter([Path,Folder,'Movie.avi']); % 输出AVI文件名字V.FrameRate = 50; % 帧率，40fpsV.Quality = 90; % 质量，百分比，90%GIF = [Path,Folder,'GIF.gif'];open(V
2020-12-30 19:42:05
1464
原创 Matlab 读取欧空局*.N1卫星遥感图像
卫星遥感图像往往具有复杂的数据格式，其文件当中既包含了字符类型，也包含了二进制，并且编码类型也各不相同，如果采用直接打开文件的方式，很可能读取到的都是乱码，因此要分步读取文件内容，并对数据进行加工处理。作者目前了解到遥感数据根据数据的处理程度分为了L0，L1，L2，L3等级别，等级越高的数据处理程度越强，越容易直接读取。（L0级数据就是卫星得到的原始遥感数据）。这里以欧空局L1级IMP.N1数据为例来展示读取遥感数据的步骤，这个步骤对于读取其他复杂文件同样适用。clear;clc;fid = fope
2020-10-19 21:06:46
584
3
原创 Matlab 打开文件夹中的多个文件
处理数据时，往往会遇到大量的同类型数据需要读取，在不确定数据文件个数和名称的情况下如何利用循环逐个打开并分析处理数据呢？这里就需要用到dir函数，以下是示例：clear;clc;Path = 'c:\filepath\';Dir = dir(strcat([Path,'*.dat']));for i = 1:length(Dir) fid = fopen([Path,Dir(i).name],'rb'); Data = fscanf(fid,'%f');end...
2020-09-18 14:13:20
2386
原创 随机波浪谱分析——自相关函数法 （MATLAB）
基本原理自相关函数假设采样时距为Δt\Delta tΔt，则海浪随机序列自相关函数可以写成R(υΔt)=1N−υ∑n=1N−υx(tn+υΔt)x(tn)R(\upsilon\Delta t)=\frac{1}{N-\upsilon}\sum\limits_{n=1}^{N-\upsilon}x(t_n+\upsilon\Delta t)x(t_n)R(υΔt)=N−υ1n=1∑N−υx(tn+υΔt)x(tn)τ=υΔt,υ=0,1,2,...,m\tau=\upsilon\Delta t
2020-08-12 09:19:51
5044
15
原创 散点数据的包络线获取（MATLAB）
在科研数据处理当中，往往会遇到许多不规则的散点数据，通过程序将这些数据的包络线（Envelope）进行描述可以更好地展现数据的取值范围。在进行数据统计当中，对于样本数据包络线的描述一般不可或缺。这里采用MATLAB软件，自己编写了一段描述散点数据包络线的代码，也是我平时使用的代码，与大家交流协同改进。程序当中变量定义如下：X：散点数据的横坐标；Y：散点数据的纵坐标；n：需要将包络线进行多项式拟合的阶数；x_min, x_max, x_int：需要进行多项式拟合的坐标下限、上限以及坐标步长；wi
2020-08-10 17:36:41
11662
9
原创 等分频率法模拟随机波列（线性波叠加原理）
线性叠加法海浪可看做一系列不同周期不同初相位的线性波叠加而成的：η(t)=∑i=1Maicos(kix−ωit+ϵi)\eta(t)=\sum\limits_{i=1}^{M}a_i\cos(k_ix-\omega_it+\epsilon_i)η(t)=i=1∑Maicos(kix−ωit+ϵi),aia_iai为第iii个组成波的振幅，ki和ωik_i和\omega_iki...
2019-12-09 20:53:57
2597
30
原创 随机波浪P-M谱
随机海浪往往具有统计特征，组成频率会呈现出某一频率集中的特征。由此而衍生出的海浪谱多种多样。P-M谱由于其结构简单，便于计算，运用于许多领域当中。用圆频率表示的P-M谱有如下形式：S(ω)=0.78ω5exp[−3.11ω4Hs2]S(\omega)=\frac{0.78}{\omega^5}\exp{[-\frac{3.11}{\omega^4H_s^2}]}S(ω)=ω50.78exp...
2019-12-09 19:57:30
8267
3
原创 随机波浪Jonswap谱
随机海浪往往具有统计特征，组成频率会呈现出某一频率集中的特征。由此而衍生出的海浪谱多种多样。其中较为著名的一种海浪谱Jonswap被广泛应用在海洋科学、海洋工程领域。以合田改进的Jonswap谱（1999）为例：S(f)=βjH1/32TP−4f−5exp[−54(TPf)−4]γexp[−(ffP−1)2/2σ2]S(f)=\beta_jH_{1/3}^2T_P^{-4}f^{-5}\e...
2019-12-09 19:38:54
16421
14
原创 PIV实验流场流速云图（MALAB contour函数用法）
clear;clc;for k=1:5 str1=['filename1']; str2=['filename2']; str3=['filename3']; str=eval(['str',num2str(k)]); load(str); [m n]=size(y); y=flipud(y); v=-v; x=x-x(1,...
2019-07-20 10:51:24
2550
3
原创 PIV流场流速矢量图、流线图（MATLAB quiver函数的用法）
PIV实验结果流速矢量图、流线图clear;clc;for k=1:n str1=['file position']; str2=['file2 position']; str3=['file3 position']; ... str=eval(['str',num2str(k)]); load([str,'filename']); [m...
2019-07-19 16:53:31
9876
5
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人