Download 31295013691000

%
%
%
%
Filename: WN3DPL0T.m
script to plot NA data in 3D form (pctloss, freq, G or K)
last updated: Nov 22, 1998
Author: Paulus S. Adisoemarta ([email protected])
sample = [ 1 2 4 5 ] ;
seteode = [ 'r i d g k' ] ;
params = [ 'a' 'r' ] ;
%
% Wellington sample codes
% real, imag, diss, G and K
% absolute and relative values
DO NOT CHANGE ANY PARAMETER AFTER THIS LINE
ial=l;ia2=35;nal=36;na2=236;
close
numsam = length(sample);
numset = length (seteode) ;
numpar = length(params);
%
%
%
%
%
index number of lA and NA files
close the current plot window
number of samples
number of parameters to plot
modes (absolute or relative)
% numsam = 1;: numset = 1;
% uncomment for debugging purposes
%Note that frequency axis, Y-axis, is now read from the data file
% not an individual file anymore.
% Loop for all samples
for ns=l:numsam
scimplecode = ['wn' int2str (sample (ns) )] ;
eval(['cd ' samplecode])
% enter the sample folder
pctloss = [samplecode 'pet' ] ;
% X-axis data filename
eval (['load ' pctloss '.dat'])
% load the X-axis data from file
xaxis = eval([pctloss]);% X-axis (the water loss in %)
% Loop to read all files of this particular sample
for nt=l:numset
setfile = [samplecode 'm' seteode(nt)];
% for absolute and relatif
for np=l:numpar
paramfile = [setfile params(np) ] ;
eval(['load ' paramfile '.dat'])
% load Z-variable
% first column is frequency, thus we have to skip it
mn = size(eval([paramfile]));
% matrix size
ztmp = eval ( [paramfile] ) ;
% save to a temp array
zaxis = ztmp(:,2:mn(2));
% after skipping first column
freq = ztmp(:,l);
* for Y-axis (Frequency)
% Surf Plot
hold off
surf(xaxis,freq(ial:ia2), zaxis(ial:ia2 ,:) )
, , ,
hold on
surf(xaxis,freq(nal:na2), zaxis(nal:na2,:) )
% clear graphic
% the lA freq
% combine with
% the NA freq
1 . .T ^^ M
% freq axis in log
set (gea, 'yseale', ' log )
,, .,
xlabeK'Moveable Water Loss, pet' ) , ylabel ('Frequency, Hz )
zlabel('unknown')
if setcode(nt) == 'g'
zlabel('Conductivity, S/m')
end
if setcode(nt) == 'k'
zlabel('Dielectric Constant')
set(gca,'zseale','log')
end
if setcode(nt) == 'd'
zlabel('Dissipation Factor')
end
% z-label for conductivity
% z-label for K
% z-axis in log
Figure C.6 Program Listing of WN3dPL0T.M
156