impdata=importdata('C:\Users\12345\Desktop\User\GenLinWin\Cavity\Ez_field.txt');
[m,n]=size(impdata.data); %求数组的大小
%disp([m,n]);
figure;
plot(impdata.data(:,1)/1000,impdata.data(:,2)/1e6,'-r','LineWidth',2);
%axis tight; %设置坐标范围等于数值范围
xlim([0,220/1000]);
grid on; %显示网格线
xlabel('z (m)','fontsize',16)
ylabel('Ez (MV/m)','fontsize',16)
% fid=fopen('C:\Users\12345\Desktop\User\GenLinWin\Cavity\Ez_matlab.txt','wt');
% data=[impdata.data(:,1)/1000,impdata.data(:,2)/1e6];
% for i=1:1:m
% for j=1:1:n
% if j==n
% fprintf(fid,'%g\n',data(i,j));
% else
% fprintf(fid,'%g\t\t',data(i,j));
% end
% end
% end
% fclose(fid);
norm=1;
fid=fopen('C:\Users\12345\Desktop\User\GenLinWin\Cavity\Ez_1.edz','wt');
data=[impdata.data(:,1)/1000,1.82*impdata.data(:,2)/1e6];
fprintf(fid,'%d %.6f\t\n',m-1,data(m,1));
fprintf(fid,'%.6f\n',norm);
for i=1:1:m
fprintf(fid,'%.6f\n',data(i,2));
end
fclose(fid);
fid=fopen('C:\Users\12345\Desktop\User\GenLinWin\Cavity\Ez_2.edz','wt');
data=[impdata.data(:,1)/1000,2.95*impdata.data(:,2)/1e6];
fprintf(fid,'%d %.6f\t\n',m-1,data(m,1));
fprintf(fid,'%.6f\n',norm);
for i=1:1:m
fprintf(fid,'%.6f\n',data(i,2));
end
fclose(fid);