Matlab批量计算CP2K的差分电荷的代码

```
clear;
clc;
close all;

% 初始位置
st = 1;
% 步长
h = 60;
% 终位
en = 1000;
%% 文件读取节点
for count = st:h:en
    disp('-------------------'+string(count)+'---------------------------');
    str_file1 = string(count)+".cube";  % 前文件
    % str_file1 = "2100.cube";
    str_file2 = "C:\Users\82399\Desktop\BaZrS_aimd-cube-ELECTRON_DENSITY-1_"+string(count)+".cube";  % 后文件
    str_out = "C:\Users\82399\Desktop\out-"+string(count)+".cube"; % 存储文件
file_1 = readlines(str_file1);
file_2 = readlines(str_file2);
 
%% 得到分界点及其索引cont
TIC = tic();
for i = 1:length(file_1)
    if length(str2num(file_1(i)))>5
        cont = i;
        break;
    end
    a(i) = file_1(i);
end
% 把开头文件写入结果文件
file_out = fopen(str_out,'w');
[R,col] = size(a);
for i = 1:col
    fprintf(file_out,a(i)+'\n');
end
data1 = str2num(char(file_1(cont:end)));
data2 = str2num(char(file_2(cont:end)));
data_out = data1-data2;
[row,C] = size(data_out);
for j = 1:row
    da = '  ' + convertCharsToStrings(num2str(data_out(j,:))) + '\n'; 
    fprintf(file_out,da);
end

fclose(file_out);
toc(TIC);
end
```

Matlab批量计算CP2K的差分电荷的代码

【温馨提示：本站文档可配置可见范围，如登录后可见、对特定群组可见等，看不到就是没权限】

注册码获取邮箱
work@liuyaoze.com

Matlab批量计算CP2K的差分电荷的代码

【温馨提示：本站文档可配置可见范围，如登录后可见、对特定群组可见等，看不到就是没权限】

注册码获取邮箱 work@liuyaoze.com

注册码获取邮箱
work@liuyaoze.com