昨晚在看Lammps前辈 Michael的博客，看到他《MATLAB计算平衡晶格常数》的博文，今早拿着Lammps的脚本算出来，后用他给的脚本在Matlab中运行，结果提示：

Warning: Polynomial is badly conditioned. Add points with distinct X values, reduce the degree of the polynomial, or try centering and scaling as described in HELP POLYFIT.

他博文中采用了8阶多项式插值，我们可以看到Lammps循环输出的数据点>9，因此多项式插值的方程数条件是满足的，那么问题出在哪呢？随后我将数据导入到cftool中，在没用smooth的前提下，发现高阶多项式插值效果都还不错。也就是说cftool的polyN和polyfit的算法还是不同的，同时高阶多项式插值是可以得到插值系数的。

多番尝试后，我将x轴的数据做了线性处理：x= (x-mean)./std(x)后，发现polyfit居然可以工作了。。。究其原因应该是在这套数据体系下，x的差别特别小，使得求解插值系数的矩阵病态。

这是修改后的代码：

今天用G95编译器编译一些老的Fortran代码，额，之所以说老就是因为gfortran编译不了啊。。。。

结果报错：

​​ld: cannot find crt1.o: No such file or directory

ld: cannot find crti.o: No such file or directory

检查了 glibc 什么的都是正确的，只能说。。。 我勒个去了。。。其实这就是lib路径的引用问题，所以临时解决方法：

​​​ export LIBRARY_PATH=/usr/lib/x86_64-linux-gnu

反正不常用，编译完了就好~~

mxio

2014.10.5​​

Origin自带的差分工具只能进行两列差分运算，当有多列Y值时，手动处理起来只能呵呵了~~

这个时候采用Origin的command窗口，进行循环处理是不二的选择。

输入：

loop(idx,2,wks.ncols)
{differentiate (WCol(1), WCol(idx));}

敲回车，效果如图：

adaptmesh函数

ATh = adaptmesh(Th, func1,func2,...)

ATh: 生成的自适应网格； Th: 原始网格； funcN: 函数。

功能：当func的值梯度很大使，可使用adaptmesh进行网格Hessian自适应划分。

其他参数：

hmin,hmax: 最小,大网格尺寸;     err: P1插值误差默认值0.01; 

errg: 相对几何误差，默认值0.01，在默些情况下必须小于1/sqrt(2)，由于几何的约束这一项可能使得某些网格尺寸小于-hmin。

nbvx: 网格划分生成的最大节点数，默认值9000;

nbsmooth: 平滑处理的次数，默认值5;

nbjacoby: 在一次平滑步中平滑计算的迭代步数，默认值6，设置为0为不平滑。

ratio: 平滑比例，当比例 0<=ratio<=1.1时不平滑，默认值1.8，当ratio>1.1时，网格变化速度受制于log(ratio)。这样有助于控制加密区域如激波区域或边界层的厚度。

omega: 平滑过程的松弛因子，默认值1;

iso: 加载到metric上的力为各项同性的，默认值：false;

abserror: 绝对误差，默认值false；

cutoff: 截断误差，默认1e-6;

verbosity: 信息标记；

inquire: 默认值f

本人最近身体状态较差，受邀做一期Comsol with Matlab的专题视频，PPT早已做好，录视频估计还要给自己身体一段时间。见期望较高，就提前流出PPT吧。届时视频也在本页更新。

转载或套用请标明出处，谢谢！本页不做答疑。