matlab串联数组不一致,要串联的数组的维度不一致。

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

读入文件