#!/bin/csh -f

grep PA1EJO wsprspots*.csv |awk -F'[,:]’ ‘{ print $3″ “$6” “$10” “$12” “$13” % T=”$8 }’ |grep ‘T=PA1EJO’ > filter.d

and next you put the filter.d file in matlab:

clear;

data = load(‘filter.d’);

hour = data(:,1)/3600;

dis = data(:,4);

azi = data(:,5);

snr = data(:,2);

pwr = data(:,3);

y = dis.*cos(azi/180*pi);

x = dis.*sin(azi/180*pi);

t0 = floor(min(hour));

t1 = ceil(max(hour));

t1-t0;

step = 24;

k=0;

for tt=t0:step:t1,

idx = find( (hour >= (tt-0.5*step)) & (hour 0),

scatter(x(idx),y(idx),80,snr(idx),’filled’); title(sprintf(‘t=%10.1f h=%.0f’,floor(tt/24),mod(tt,24)));

axis([-9000 9000 -9000 9000]); pause(1);

k = k + 1;

xk(k) = tt;

yk(k,1) = mean(dis(idx));

yk(k,2) = median(dis(idx));

yk(k,3) = max(dis(idx));

end

end

figure(2);

plot((xk-t0)/24,yk(:,1),’o-‘,(xk-t0)/24,yk(:,2),’.-‘,(xk-t0)/24,yk(:,3),’.-‘);

I hope this sort of helps you to get started.

]]>I wonder how you pick up the data you need from the files ? ]]>

If all say the Earth is flat, that doesn’t make it flat. If all say a repeater(manager) has (to uphold the) houserules, that doesn’t make it so. If you follow the herd, it is unlikely you get any further than the herd.

]]>