1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
|
!-------------------------------------------------------------------------------
!
! This file is part of the WSPR application, Weak Signal Propagation Reporter
!
! File Name: wsprcode.f90
! Description: This program provides examples of the source encoding,
! convulsional error-control coding, bit and symbol ordering,
! and synchronizing information contained in WSPR messages.
!
! Copyright (C) 2001-2014 Joseph Taylor, K1JT
! License: GPL-3
!
! This program is free software; you can redistribute it and/or modify it under
! the terms of the GNU General Public License as published by the Free Software
! Foundation; either version 3 of the License, or (at your option) any later
! version.
!
! This program is distributed in the hope that it will be useful, but WITHOUT
! ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
! FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
! details.
!
! You should have received a copy of the GNU General Public License along with
! this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
! Street, Fifth Floor, Boston, MA 02110-1301, USA.
!
!-------------------------------------------------------------------------------
program wsprcode
parameter (NSYM=162)
parameter (MAXSYM=176)
character*22 msg,msg2
integer*1 data0(13)
integer*1 data1(13)
integer*1 dat(206)
integer*1 softsym(206)
! Define the sync vector:
integer*1 sync(NSYM)
data sync/ &
1,1,0,0,0,0,0,0,1,0,0,0,1,1,1,0,0,0,1,0, &
0,1,0,1,1,1,1,0,0,0,0,0,0,0,1,0,0,1,0,1, &
0,0,0,0,0,0,1,0,1,1,0,0,1,1,0,1,0,0,0,1, &
1,0,1,0,0,0,0,1,1,0,1,0,1,0,1,0,1,0,0,1, &
0,0,1,0,1,1,0,0,0,1,1,0,1,0,1,0,0,0,1,0, &
0,0,0,0,1,0,0,1,0,0,1,1,1,0,1,1,0,0,1,1, &
0,1,0,0,0,1,1,1,0,0,0,0,0,1,0,1,0,0,1,1, &
0,0,0,0,0,0,0,1,1,0,1,0,1,1,0,0,0,1,1,0, &
0,0/
! Metric table for decoding from soft symbols
integer mettab(0:255,0:1)
data mettab/ &
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, &
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, &
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, &
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, &
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, &
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, &
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, &
5, 5, 5, 5, 5, 5, 5, 5, 5, 4, &
4, 4, 4, 4, 4, 4, 4, 4, 4, 4, &
4, 4, 4, 4, 4, 4, 4, 4, 4, 4, &
3, 3, 3, 3, 3, 3, 3, 3, 3, 2, &
2, 2, 2, 2, 1, 1, 1, 1, 0, 0, &
-1, -1, -1, -2, -2, -3, -4, -4, -5, -6, &
-7, -7, -8, -9, -10, -11, -12, -12, -13, -14, &
-15, -16, -17, -17, -18, -19, -20, -21, -22, -22, &
-23, -24, -25, -26, -26, -27, -28, -29, -30, -30, &
-31, -32, -33, -33, -34, -35, -36, -36, -37, -38, &
-38, -39, -40, -41, -41, -42, -43, -43, -44, -45, &
-45, -46, -47, -47, -48, -49, -49, -50, -51, -51, &
-52, -53, -53, -54, -54, -55, -56, -56, -57, -57, &
-58, -59, -59, -60, -60, -61, -62, -62, -62, -63, &
-64, -64, -65, -65, -66, -67, -67, -67, -68, -69, &
-69, -70, -70, -71, -72, -72, -72, -72, -73, -74, &
-75, -75, -75, -77, -76, -76, -78, -78, -80, -81, &
-80, -79, -83, -82, -81, -82, -82, -83, -84, -84, &
-84, -87, -86, -87, -88, -89, -89, -89, -88, -87, &
-86, -87, -84, -84, -84, -83, -82, -82, -81, -82, &
-83, -79, -80, -81, -80, -78, -78, -76, -76, -77, &
-75, -75, -75, -74, -73, -72, -72, -72, -72, -71, &
-70, -70, -69, -69, -68, -67, -67, -67, -66, -65, &
-65, -64, -64, -63, -62, -62, -62, -61, -60, -60, &
-59, -59, -58, -57, -57, -56, -56, -55, -54, -54, &
-53, -53, -52, -51, -51, -50, -49, -49, -48, -47, &
-47, -46, -45, -45, -44, -43, -43, -42, -41, -41, &
-40, -39, -38, -38, -37, -36, -36, -35, -34, -33, &
-33, -32, -31, -30, -30, -29, -28, -27, -26, -26, &
-25, -24, -23, -22, -22, -21, -20, -19, -18, -17, &
-17, -16, -15, -14, -13, -12, -12, -11, -10, -9, &
-8, -7, -7, -6, -5, -4, -4, -3, -2, -2, &
-1, -1, -1, 0, 0, 1, 1, 1, 1, 2, &
2, 2, 2, 2, 3, 3, 3, 3, 3, 3, &
3, 3, 3, 4, 4, 4, 4, 4, 4, 4, &
4, 4, 4, 4, 4, 4, 4, 4, 4, 4, &
4, 4, 4, 4, 5, 5, 5, 5, 5, 5, &
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, &
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, &
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, &
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, &
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, &
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, &
5, 5, 5, 5, 5, 5, 5, 5, 5, 5, &
5, 5/
! Get command-line argument(s)
nargs=iargc()
if(nargs.ne.1) then
print*,'Usage: WSPRcode "message"'
go to 999
endif
call getarg(1,msg) !Get message from command line
write(*,1000) msg
1000 format('Message: ',a22)
nbits=50+31 !User bits=50, constraint length=32
nbytes=(nbits+7)/8
ndelta=50
limit=20000
data0=0
call wqencode(msg,ntype0,data0) !Source encoding
write(*,1002) data0(1:7),data0(1:6),data0(7)/64
1002 format(/'Source-encoded message, 50 bits:'/'Hex: ',7z3.2/ &
'Binary: ',6b9.8,b3.2)
call encode232(data0,nbytes,dat,MAXSYM) !Convolutional encoding
call inter_mept(dat,1) !Interleaving
write(*,1004)
1004 format(/'Data symbols:')
write(*,1006) (dat(i),i=1,NSYM)
1006 format(5x,30i2)
write(*,1008)
1008 format(/'Sync symbols:')
write(*,1006) (sync(i),i=1,NSYM)
write(*,1010)
1010 format(/'Channel symbols:')
write(*,1006) (2*dat(i)+sync(i),i=1,NSYM)
call inter_mept(dat,-1) !Remove interleaving
softsym=-dat !Simulate soft symbols
! Call the sequential (Fano algorithm) decoder
call fano232(softsym,nbits,mettab,ndelta,limit,data1,ncycles,metric,nerr)
call wqdecode(data1,msg2,ntype1)
write(*,1020) msg2,ntype1
1020 format(/'Decoded message: ',a22,' ntype:',i3)
999 end program wsprcode
include 'wspr_old_subs.f90'
|
