Nick_Cui的博客

龙格库塔法是一种求解高阶常微分方程的常用方法，在工程当中应用广泛，例如求解物体的运动方程等。这里我们通过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

处理数据时，往往会遇到大量的同类型数据需要读取，在不确定数据文件个数和名称的情况下如何利用循环逐个打开并分析处理数据呢？这里就需要用到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...

在科研数据处理当中，往往会遇到许多不规则的散点数据，通过程序将这些数据的包络线（Envelope）进行描述可以更好地展现数据的取值范围。在进行数据统计当中，对于样本数据包络线的描述一般不可或缺。这里采用MATLAB软件，自己编写了一段描述散点数据包络线的代码，也是我平时使用的代码，与大家交流协同改进。程序当中变量定义如下：X：散点数据的横坐标；Y：散点数据的纵坐标；n：需要将包络线进行多项式拟合的阶数；x_min, x_max, x_int：需要进行多项式拟合的坐标下限、上限以及坐标步长；wi