cemrd4.m

来自「这是一个 York大学开发的FDTD代码。很好用的工具」· M 代码 · 共 29 行

function [f,db,a,p]=cemrd4(fn)
% [f,db,a,p]=cemrdfd('filename')
% Read CEM frequency domain file with phase
% You need to run:
% xtransall phase
% AND
% xfreq phase
% to get phase data in fd??.dat
% f=frequency MHz
% db= dB (V/m)
% a = amplitude (V/m)
% p = Phase (degrees)
tdf=fopen(fn);
if tdf==-1 
disp('Can not open file!');

else
tdid=fscanf(tdf,'%*s %d %d %d',3);
td=fscanf(tdf,'%g %g %g %g',inf);
f=[td(1:4:size(td,1))];
db=[td(2:4:size(td,1))];
a=[td(3:4:size(td,1))];
p=[td(4:4:size(td,1))];
fclose(tdf);
%size(td,1)
%size(f,1)
end

end