已结贴√
问题点数:20 回复次数:2
这个错误可能的原因是你忘了初始化变量
function EMDGETK
%计算Argon热传导系数
%采用快速Fouie变换求解EMD算法得到的热传导系数
%ref: Molecular simulation of liquids.p190
%
%为了节约计算时间和所需内存量,这里将热流文件分开成几个文件,如每1000000个数据放在一个文件中.
%
%计算用到的文件:initpara.dat, datafile.in, flux*.dat,avg.dat;
%flux1.dat, flux2.dat,类似文件的个数
%*********************************************
%依据所拥有的flux*.dat的个数,修改filenum的值*
%*********************************************
filenum=1;
%最大相关长度
maxcor=100000;
%通过一下相关长度计算热传导系数
cor=10000:10000:maxcor;
Kb=1.38E-23;
[filename, filepath] = uigetfile('initpara.dat','载入数据');
if filepath==0
return;
end