clc
clear
for i=1:5
fileName = [num2str(i,'%04d') '.txt']; %读取瞬时流速
a=load(fileName);
a(1:80,3)=0;
a(:,6)=-a(:,6);
E=[];
for j=1:238
for k=((j-1)*79+j):1:(j*79+(j-1))
B=(a(k+1,5)-a(k,5))/0.07678;
A=(a(k+80,5)-a(k,5))/0.07678;
D=(a(k+1,6)-a(k,6))/0.07678;
C=(a(k+80,6)-a(k,6))/0.07678;
M=([A,B;C,D]); %求速度梯度矩阵M
if ((-A-D)^2-4*(A*D-B*C)<0)
TZZ=eig(M);
TZZ1=[TZZ(1,1) TZZ(2,1)]; %将两个特征值放在一行
lamedaci=abs(imag(TZZ1)); %两个特征值取绝对值
else
lamedaci=[0 0];
end
E=[E;lamedaci];
% E1=[0.09597+(k-1)*0.7678 (k-1)*0.7678 E];
end
E1=[a(((j-1)*79+j):(j*79+(j-1)),1:4) E];
end
str=['gk0201_xzqdf' num2str(i),'.dat'];
fid=fopen(str,'w');
fprintf(fid,'%d\t %d\t %f\t %f\t %f\t %f\r\n',E1(:,1:6));
%fprintf(fid,'%f\t %f\t %f\t %f\r\n',E1(:,1:4));
fprintf(fid,'\n');
fclose(fid);
end
2020-4-4 21:58 上传
点击文件名下载附件
1.19 MB, 下载次数: 8
读入文件