%不考虑源项,材料组合与厚度值
clc;clear;
vardata=load('20210506_174438.log');
var1=vardata(2:10801,2:6);
s1=zeros(2700,5);
for i1=1:27:2700
for j1=1:27
s1(fix(i1/27)*27+j1,:)=var1(fix(i1/27)*108+j1,:);
end
end
s1(:,2)=s1(:,2)-s1(:,1);
fid=fopen('inpforsec.txt','wt');
[m,n]=size(s1);
for ii=1:m
for jj=1:n
if jj==n
fprintf(fid,'%e\n',s1(ii,jj));
else
fprintf(fid,'%e\t',s1(ii,jj));
end
end
end
fclose(fid);
%%
%将4种源项的计算结果分开
clc;clear;
outt=importdata('outall.txt');
out=outt.data;
s1=zeros(2700,2);s2=zeros(2700,2);s3=zeros(2700,2);s4=zeros(2700,2);
for ii=1:27:2700
for jj=1:27
s1(fix(ii/27)*27+jj,:)=out(fix(ii/27)*108+jj,:);
s2(fix(ii/27)*27+jj,:)=out(fix(ii/27)*108+jj+27,:);
s3(fix(ii/27)*27+jj,:)=out(fix(ii/27)*108+jj+54,:);
处理数据样本及优化后数据
最新推荐文章于 2023-11-08 16:57:18 发布