EX:Fluent导出一个记录小球重心坐标随时间变化的文件,在后处理中需要用MATLAB导入并且绘制图形。
文件如图所示:
由于带有5个标题行,所以textscan(fid,'%f %f %f %f','headerlines',5)需要填写额外的控制参数以便跳过标题行。
代码如下:
%% 该程序用于读取重心坐标的信息,并且生成图形
fid=fopen('fluent_ball-latter','r');
%% 调用fopen函数打开小球重心坐标文件,返回文件句柄给fid
motion_xy=textscan(fid,'%f %f %f %f','headerlines',5);%其中fid只是一个变量,可以随意命名
%% textscan读取的文件储存在元胞数组motion_xy中,分别定义时间t,x坐标,y坐标向量
t=motion_xy(1,1);
x=motion_xy(1,2);
y=motion_xy(1,3);
%% 元胞数组不能直接用于plot函数,这里先用cell2mat转化成矩阵(2意味着“to”)
xx=cell2mat(x);
yy=cell2mat(y);
tt=cell2mat(t);
%% 利用plot函数绘图
plot(tt,xx