CN110311759A - A kind of magnetic induction Communication System Design method and communication system based on quasi-cyclic LDPC code - Google Patents
A kind of magnetic induction Communication System Design method and communication system based on quasi-cyclic LDPC code Download PDFInfo
- Publication number
- CN110311759A CN110311759A CN201910709725.0A CN201910709725A CN110311759A CN 110311759 A CN110311759 A CN 110311759A CN 201910709725 A CN201910709725 A CN 201910709725A CN 110311759 A CN110311759 A CN 110311759A
- Authority
- CN
- China
- Prior art keywords
- matrix
- magnetic induction
- communication system
- num
- ldpc code
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004891 communication Methods 0.000 title claims abstract description 56
- 230000006698 induction Effects 0.000 title claims abstract description 44
- 238000013461 design Methods 0.000 title claims abstract description 26
- 238000000034 method Methods 0.000 title claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims abstract description 55
- 239000002245 particle Substances 0.000 claims abstract description 22
- 238000005457 optimization Methods 0.000 claims abstract description 5
- 230000009977 dual effect Effects 0.000 claims description 19
- 230000005540 biological transmission Effects 0.000 claims description 14
- 230000005284 excitation Effects 0.000 claims description 3
- 230000008054 signal transmission Effects 0.000 abstract description 3
- 230000007812 deficiency Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 3
- 238000005265 energy consumption Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 3
- 241000209094 Oryza Species 0.000 description 2
- 235000007164 Oryza sativa Nutrition 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 235000009566 rice Nutrition 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L1/00—Arrangements for detecting or preventing errors in the information received
- H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
- H04L1/0056—Systems characterized by the type of code used
- H04L1/0057—Block codes
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L1/00—Arrangements for detecting or preventing errors in the information received
- H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
- H04L1/0056—Systems characterized by the type of code used
- H04L1/0061—Error detection codes
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B5/00—Near-field transmission systems, e.g. inductive or capacitive transmission systems
- H04B5/20—Near-field transmission systems, e.g. inductive or capacitive transmission systems characterised by the transmission technique; characterised by the transmission medium
- H04B5/24—Inductive coupling
- H04B5/26—Inductive coupling using coils
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Error Detection And Correction (AREA)
Abstract
Poor anti jamming capability of the present invention for signal transmission in existing magnetic induction communication, the more serious technological deficiency of error propagation when wireless communication, propose magnetic induction Communication System Design method and communication system based on quasi-cyclic LDPC code, the present invention generates the basic matrix of quasi-cyclic LDPC code using Discrete Particle Swarm Optimization Algorithm and PEXIT algorithm and calculates the setting distance of transmitting coil and receiving coil in magnetic induction communication system, and the magnetic induction communication system based on design method design can be effectively improved the performance of magnetic induction communication system under the Complex Channel environment of underground.
Description
Technical field
The present invention relates to magnetic induction field of communication technology, especially a kind of magnetic induction communication system based on quasi-cyclic LDPC code
Design method of uniting and communication system.
Background technique
Under the complex environment of underground, the magnetic induction communication technology has more significant excellent compared with traditional electromagnetic communication
Gesture.The reason is that electromagnetic communication is in for example different types of rock of various underground mediums, the factors shadow such as soil of different moisture content
Under sound, decay larger, seriously affect the quality of signal propagation, is not suitable for the signal transmission of underground complex environment；And magnetic strength
The communication technology is answered, signal is transmitted by the medium influence very little in underground environment, there is more apparent advantage.
Magnetic induction communication, can be with according to the deployment way of coil using the wireless transmission of the coupling progress signal between coil
It is divided into direct magnetic induction communication and magnetic induction wave guide communication.Channel coding occupies an important position in wireless communications, by biography
Defeated signal carries out channel coding, can effectively improve the reliability of signal wireless transmission, reduce the bit error rate of wireless communication system
(BER:Bit error rate).
LDPC (Low Density Parity Check, low density parity check code) code is as important in channel coding
A member is increasingly taken seriously, and recent LDPC code is selected in 5G mobile communication standard, quasi- circulation as one of channel coding schemes
LDPC code (QC-LDPC code), as submatrix, can be facilitated using shift register using cyclic shift matrices and null matrix
It realizes, is a kind of mostly important LDPC code, has been used for IEEE 802.16e and IEEE 802.11n and DVB- at present
The international standards such as S2.The research of channel coding involved in existing magnetic induction Communication Studies is seldom, and not yet discovery is in magnetic induction
The application study of LDPC code in communication.And in existing magnetic induction Communication Studies, the letter communicated dedicated for magnetic induction is not found
Road code optimization design studies this is because traditional LDPC code coding needs to be converted into generator matrix G from H-matrix, then carry out
Coding, it is a bottleneck of LDPC code application that encoder complexity is relatively large.
Summary of the invention
Goal of the invention: for the blank for making up the prior art, and for the magnetic induction communication requirement in the complex environment of underground, this
Invention using quasi-cyclic LDPC code code structure it is relatively easy, matrix memory space is relatively small, is well suited for underground wireless sensing
The characteristics of net applies the low energy consumption for sensing node to require proposes a kind of based on quasi-cyclic LDPC for underground complex environment
The magnetic induction communication system and its design method of code, can be effectively improved the property of magnetic induction communication system under the Complex Channel environment of underground
Energy.
Technical solution: in order to realize the above technical effect, the present invention proposes following technical scheme:
A kind of magnetic induction Communication System Design method based on quasi-cyclic LDPC code, this method are optimized using discrete particle cluster
Algorithm and PEXIT algorithm generate the basic matrix of quasi-cyclic LDPC code and calculate transmitting coil and reception in magnetic induction communication system
The setting distance of coil, comprising steps of
(1) setting the number of iterations is r, initializes r=1, initializes local optimum pbest at this timerFor a very little
Value；Maximum number of iterations I is setmaxWith particle number Pnum；
(2) P is generated at randomnumA dimension is the dual vector of M × KEach dual vectorThe letter of corresponding basic matrix B
Cease part, i.e. B (HI), B (HI) size be M × K, B (HI) with size be K × K basic matrix check part B (Hp) group symphysis
The binary basic matrix for being M × N at size: B (H)=[B (HI)|B(HP)], wherein N=M+K；B(Hp) it is shown below:
Wherein, the position of " 1 " of first row and most next two columns is fixed, and first row column weight is 3, remaining column column weight is 2, upper right
Diagonal line is " 1 ", and in addition to first row and most next two columns, the position of another " 1 " of remaining each column is random, but is to ensure that column weight
It is 2；
PnumA dual vectorP is constructed altogethernumA basic matrix after obtaining each basic matrix, reconstructs its corresponding H square
Battle array, that is, have PnumA H-matrix；
With PnumA dual vector is PnumA particle initializes population；NoteFor the corresponding binary arrow of p-th of particle
AmountT-th of bit；
(3) it to the H-matrix of each quasi-cyclic LDPC code, is calculated in the magnetic induction communication system and is sent out by PEXIT algorithm
It penetrates, the distance between receiving coil, specific steps include:
(31) the distance between initialization transmitting, receiving coil d=d0；Initialize j-th of variable point vjIt is verified with i-th
Point ciAssociated each edge passes to vjLikelihood information and vjPriori mutual information I beforeAv(i, j) is 0；Calculate vjIt is initial
Likelihood information are as follows:
Wherein,
σ*=1.6363
aJ, 1=-0.0421061, bJ, 1=0.209252
cJ, 1=-0.00640081,
aJ, 2=0.00181491, bJ, 2=-0.142675
cJ, 2=-0.0822054, dJ, 2=0.0549608
Wherein,R indicates the code rate of quasi-cyclic LDPC code, f (d)=Pt-LMI-Pn, PtTo emit function
Rate, PnFor noise power, LMIFor the path loss of magnetic induction communication, LMIExpression formula are as follows:
Nt、NrRespectively the number of turns of transmitting coil and receiving coil, at、arRespectively the half of transmitting coil and receiving coil
Diameter；
(32) it updates:
More new variables point vjPass to checkpoint ciLikelihood information and vjBetween related mutual information IEv(i, j):σEvThe expression formula of (i, j) are as follows:
I*=0.3646
aσ, 1=1.09542, bσ, 1=0.214217, cσ, 1=2.33727
aσ, 2=0.706692, bσ, 2=0.386013, cσ, 2=-1.75017
Wherein bI, jFor in the corresponding bipartite graph of H-matrix of quasi-cyclic LDPC code, link variable node vjWith check-node ci
Between side, be 1 when position is the element of (i, j) in H-matrix, show vjAnd ciBetween have Bian Xianglian, then its corresponding bI, j=
1, otherwise bI, j=0,
Update checkpoint ciWith variable point vjAssociated each edge passes to ciLikelihood information and vjBetween priori mutual trust
Breath: IAc(i, j)=IEv(i, j)；
Update checkpoint ciPass to variable point vjLikelihood information and vjBetween external mutual information:σEcThe expression formula of (i, j) are as follows:
More new variables point vjWith checkpoint ciAssociated each edge passes to vjLikelihood information and vjPriori mutual trust before
Cease IAv(i, j): IAv(i, j)=IEc(i, j)
More new variables point vjPosteriority likelihood information and variable point vjBetween posteriority mutual information:
Wherein,
(33) judge whether to meet iteration stopping criterion:
If satisfied, then end step (3), export d, then using d as the distance threshold of the current round of corresponding protographIf not satisfied, then updating d=d+ Δ d, Δ d is preset increment step-length, then return step (32)；
(4) P is obtained to step (2)numIt is total to obtain P after a H-matrix executes step (3) respectivelynumA distance threshold
From this PnumA distance thresholdIn find out the global optimum of epicycle iteration
(5) dual vector updates:
Update particle rapidity:
Wherein,
Wherein, λ, c1、c2、η1、η2For coefficient, wherein λ=1, c1=c2=2, η1、η2It is uniformly to divide in section (0,1)
The random number of cloth； The equally distributed random number between (0,1)；
(6) the number of iterations r=r+1 is updated, judges whether the number of iterations meets r > Imax；If satisfied, then stop iteration, root
Optimal QC-LDPC code matrix, and the global optimum that will be obtained are constructed according to the current corresponding dual vector of particleAssignment
To d；Otherwise, P is regenerated according to step (5) updated particlenumA H-matrix returns to (3) step.
The present invention also proposes a kind of magnetic induction communication system designed using the design method, comprising:
Transmitting terminal and receiving end；Transmitting terminal includes encoder, modulator and transmitting coil；Receiving end include receiving coil,
Demodulator and decoder；The distance between transmitting coil and receiving coil d are obtained by design method described in claim 1；
Transmitting terminal executes following steps: the H-matrix for the quasi-cyclic LDPC code that encoder is constructed using claim 1 is to source number
According to being encoded, after the data after coding are sent to modulators modulate, then transmitting coil is arrived, transmitting coil generates induced current；
Receiving end executes following steps: excitation receiving coil generates induced current, then related data is sent to demodulator solution
It adjusts, serial data is sent to decoder again after demodulation, the data after being rebuild after decoder for decoding.
The utility model has the advantages that compared with prior art, present invention has the advantage that
The magnetic induction communication system that the present invention constructs effectively improves the communication of system by using quasi-cyclic LDPC code
Performance, the quasi-cyclic LDPC code that optimization design obtains can be encoded directly with H-matrix, and encoder complexity is low, and energy consumption is small, non-
Often suitable for the low energy consumption requirement of underground wireless sensing net node.
It is magnetic in addition, can also predict the coil maximum communication distance in direct magnetic induction communication by the algorithm proposed
Induction Communications System design coil deployment provides good guidance.
Detailed description of the invention
Fig. 1 show the flow chart of design method of the present invention；
Fig. 2 is the architecture diagram using the magnetic induction communication system of design method of the present invention design；
Fig. 3 is the basic matrix schematic diagram of three kinds of A codes described in embodiment；
Fig. 4 is the BER performance comparison schematic diagram of A2 code and 802.16e code；
Fig. 5 is the H-matrix schematic diagram of B1 code and B2 code；
Fig. 6 is basic matrix and the H-matrix signal of B3 code；
Fig. 7 is the BER performance comparison schematic diagram of three B codes.
Specific embodiment
The present invention will be further explained in the following with reference to the drawings and specific embodiments.
The present invention proposes a kind of magnetic induction Communication System Design method based on quasi-cyclic LDPC code, and this method is using discrete
Particle swarm optimization algorithm and PEXIT algorithm generate the basic matrix of quasi-cyclic LDPC code and calculate to be emitted in magnetic induction communication system
The setting distance of coil and receiving coil, comprising steps of
(1) setting the number of iterations is r, initializes r=1, initializes local optimum pbest at this timerFor a very little
Value；Maximum number of iterations I is setmaxWith particle number Pnum；
(2) P is generated at randomnumA dimension is the dual vector of M × KEach dual vectorThe letter of corresponding basic matrix B
Cease part, i.e. B (HI), B (HI) size be M × K, B (HI) with size be K × K basic matrix check part B (Hp) group symphysis
The binary basic matrix for being M × N at size: B (H)=[B (HI)|B(Hp)], wherein N=M+K；B(Hp) it is shown below:
Wherein, the position of " 1 " of first row and most next two columns is fixed, and first row column weight is 3, remaining column column weight is 2, upper right
Diagonal line is " 1 ", and in addition to first row and most next two columns, the position of another " 1 " of remaining each column is random, but is to ensure that column weight
It is 2；
PnumA dual vectorP is constructed altogethernumA basic matrix after obtaining each basic matrix, reconstructs its corresponding H square
Battle array, that is, have PnumA H-matrix；
With PnumA dual vector is PnumA particle initializes population；NoteFor the corresponding binary arrow of p-th of particle
AmountT-th of bit；
(3) it to the H-matrix of each quasi-cyclic LDPC code, is calculated in the magnetic induction communication system and is sent out by PEXIT algorithm
It penetrates, the distance between receiving coil, specific steps include:
(31) the distance between initialization transmitting, receiving coil d=d0；Initialize j-th of variable point vjIt is verified with i-th
Point ciAssociated each edge passes to vjLikelihood information and vjPriori mutual information I beforeAv(i, j) is 0；Calculate vjIt is initial
Likelihood information are as follows:
Wherein,
σ*=1.6363
aJ, 1=-0.0421061, bJ, 1=0.209252
cJ, 1=-0.00640081,
aJ, 2=0.00181491, bJ, 2=-0.142675
cJ, 2=-0.0822054, dJ, 2=0.0549608
Wherein,R indicates the code rate of quasi-cyclic LDPC code, f (d)=Pt-LMI-Pn, PtTo emit function
Rate, PnFor noise power, LMIFor the path loss of magnetic induction communication, LMIExpression formula are as follows:
Nt、NrRespectively the number of turns of transmitting coil and receiving coil, at、arRespectively the half of transmitting coil and receiving coil
Diameter；
(32) it updates:
More new variables point vjPass to checkpoint ciLikelihood information and vjBetween related mutual information IEv(i, j):σEvThe expression formula of (i, j) are as follows:
I*=0.3646
aσ, 1=1.09542, bσ, 1=0.214217, cσ, 1=2.33727
aσ, 2=0.706692, bσ, 2=0.386013, cσ, 2=-1.75017
Wherein bI, jFor in the corresponding bipartite graph of H-matrix of quasi-cyclic LDPC code, link variable node vjWith check-node ci
Between side, be 1 when position is the element of (i, j) in H-matrix, show vjAnd ciBetween have Bian Xianglian, then its corresponding bI, j=
1, otherwise bI, j=0,
Update checkpoint ciWith variable point vjAssociated each edge passes to ciLikelihood information and vjBetween priori mutual trust
Breath: IAc(i, j)=IEv(i, j)；
Update checkpoint ciPass to variable point vjLikelihood information and vjBetween external mutual information:σEcThe expression formula of (i, j) are as follows:
More new variables point vjWith checkpoint ciAssociated each edge passes to vjLikelihood information and vjPriori mutual trust before
Cease IAv(i, j): IAv(i, j)=IEc(i, j)
More new variables point vjPosteriority likelihood information and variable point vjBetween posteriority mutual information:
Wherein,
(33) judge whether to meet iteration stopping criterion:
If satisfied, then end step (3), export d, then using d as the distance threshold of the current round of corresponding protographIf not satisfied, then updating d=d+ Δ d, Δ d is preset increment step-length, then return step (32)；
(4) P is obtained to step (2)numIt is total to obtain P after a H-matrix executes step (3) respectivelynumA distance threshold
From this PnumA distance thresholdIn find out the global optimum of epicycle iteration
(5) dual vector updates:
Update particle rapidity:
Wherein,
Wherein, λ, c1、c2、η1、η2For coefficient, wherein λ=1, c1=c2=2, η1、η2It is uniformly to divide in section (0,1)
The random number of cloth； The equally distributed random number between (0,1)；
(6) the number of iterations r=r+1 is updated, judges whether the number of iterations meets r > Imax；If satisfied, then stop iteration, root
Optimal QC-LDPC code matrix, and the global optimum that will be obtained are constructed according to the current corresponding dual vector of particleAssignment
To d；Otherwise, P is regenerated according to step (5) updated particlenumA H-matrix returns to (3) step.
The present invention also proposes a kind of magnetic induction communication system designed using the design method, comprising:
Transmitting terminal and receiving end；Transmitting terminal includes encoder, modulator and transmitting coil；Receiving end include receiving coil,
Demodulator and decoder；The distance between transmitting coil and receiving coil d are obtained by the design method；
Transmitting terminal executes following steps: the H-matrix pair for the quasi-cyclic LDPC code that encoder is constructed using the design method
Source data is encoded, and after the data after coding are sent to modulators modulate, then arrives transmitting coil, transmitting coil generates induced electricity
Stream；
Receiving end executes following steps: excitation receiving coil generates induced current, then related data is sent to demodulator solution
It adjusts, serial data is sent to decoder again after demodulation, the data after being rebuild after decoder for decoding.
Unlike traditional electromagnetic communication, system coding, modulated data are sent out in such a way that coil couples
Receiving coil is given, then carries out demodulation decoding.And the relative parameters setting of the distance between coil and coil, transmission power
It is had an important influence Deng all being transmitted for signal.
Unlike existing magnetic induction communication system, the present invention considers quasi-cyclic LDPC code as channel coding side
Case is just encoded before signal transmission, and there are no magnetic induction communication systems to be followed using standard in other existing researchs
The scheme of ring LDPC code.
Below by the specific embodiment technical effect that the present invention is further explained.
Embodiment: the present embodiment is constructed 2 groups of quasi-cyclic LDPC codes (A code and B code) using the design method, and totally 6
A code design parameter related to B code is set forth to table 3 in code, table 1, and table 2 and table 4 provide method construct proposed by the present invention
The corresponding system coil of quasi-cyclic LDPC code between maximum transmission distance threshold value, i.e. maximum between transmitting coil and receiving coil
Distance, the basic matrix of code and the diagram of H-matrix are as shown in Figure 3, Figure 5 and Figure 6.
The relevant parameter of 1 A code of table
The maximum transmission distance threshold value of 2 A code of table
| Title | A1 code | A2 code | A3 code |
| Maximum transmission distance (rice) | 8.4 | 8.42 | 8.37 |
The relevant parameter of 3 B code of table
The maximum transmission distance threshold value of 4 B code of table
| Title | B1 code | B2 code | B3 code |
| Maximum transmission distance (rice) | 17.16 | 18.48 | 36.91 |
Fig. 4 is the quasi-cyclic LDPC code (parameters such as code rate, code length are identical as A2 code) in A2 code and 802.16 standards in magnetic
The BER performance emulated in Induction Communications System.As can be seen that the standard of method construct proposed by the present invention follows LDPC code in figure
When BER is (- 4) 10^, transmission range ratio 802.16e code improves 4 meters or so.
Fig. 7 is B1, B2 and B3 code, as it can be seen that transmission distance of three quasi-cyclic LDPC codes in magnetic induction communication system in figure
From having larger improvement than A code, illustrate that the parameters such as transmission power, coil radius have larger impact for performance etc..
Fig. 4 and Fig. 7 can also be seen that the maximum transmission distance threshold value predicted by the method for the invention in table 2 and table 4 and
The maximum transmission distance that BER is emulated coincide fine.
It is said that the bright one aspect of the present invention of the figures above sum number can predict very well based on quasi-cyclic LDPC code maximum transmitted away from
From, the quasi-cyclic LDPC code of magnetic induction communication system on the other hand can be suitble to Optimal Construction, be very suitable for as underground without
The design alternatives of line Sensor Network.
The above is only a preferred embodiment of the present invention, it should be pointed out that: for the ordinary skill people of the art
For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also answered
It is considered as protection scope of the present invention.
Claims (2)
1. a kind of magnetic induction Communication System Design method based on quasi-cyclic LDPC code, which is characterized in that this method is using discrete
Particle swarm optimization algorithm and PEXIT algorithm generate the basic matrix of quasi-cyclic LDPC code and calculate to be emitted in magnetic induction communication system
The setting distance of coil and receiving coil, comprising steps of
(1) setting the number of iterations is r, initializes r=1, initializes local optimum pbest at this timerFor the value of a very little；
Maximum number of iterations I is setmaxWith particle number Pnum；
(2) P is generated at randomnumA dimension is the dual vector of M × KEach dual vectorThe information portion of corresponding basic matrix B
Point, i.e. B (HI), B (HI) size be M × K, B (HI) with size be K × K basic matrix check part B (Hp) combination producing is big
The small binary basic matrix for M × N: B (H)=[B (HI)|B(HP)], wherein N=M+K；B(Hp) it is shown below:
Wherein, the position of " 1 " of first row and most next two columns is fixed, and first row column weight is 3, remaining column column weight is 2, and upper right is diagonal
Line is " 1 ", and in addition to first row and most next two columns, the position of another " 1 " of remaining each column is random, but is to ensure that column weight is 2；
PnumA dual vectorP is constructed altogethernumA basic matrix after obtaining each basic matrix, reconstructs its corresponding H-matrix, i.e.,
There is PnumA H-matrix；
With PnumA dual vector is PnumA particle initializes population；NoteFor the corresponding dual vector of p-th of particle
T-th of bit；
(3) it to the H-matrix of each quasi-cyclic LDPC code, is calculated in the magnetic induction communication system by PEXIT algorithm and emits, connects
The distance between take-up circle, specific steps include:
(31) the distance between initialization transmitting, receiving coil d=d0；Initialize j-th of variable point vjWith i-th of checkpoint ci
Associated each edge passes to vjLikelihood information and vjPriori mutual information I beforeAv(i, j) is 0；Calculate vjInitial likelihood
Information are as follows:
Wherein,
σ*=1.6363
aJ, 1=-0.0421061, bJ, 1=0.209252
cJ, 1=-0.00640081,
aJ, 2=0.00181491, bJ, 2=-0.142675
cJ, 2=-0.0822054, dJ, 2=0.0549608
Wherein,R indicates the code rate of quasi-cyclic LDPC code, f (d)=Pt-LMI-Pn, PtFor transmission power, Pn
For noise power, LMIFor the path loss of magnetic induction communication, LMIExpression formula are as follows:
Nt、NrRespectively the number of turns of transmitting coil and receiving coil, at、arThe respectively radius of transmitting coil and receiving coil；
(32) it updates:
More new variables point vjPass to checkpoint ciLikelihood information and vjBetween related mutual information IEv(i, j):σEvThe expression formula of (i, j) are as follows:
I*=0.3646
aσ, 1=1.09542, bσ, 1=0.214217, cσ, 1=2.33727
aσ, 2=0.706692, bσ, 2=0.386013, cσ, 2=-1.75017
Wherein, bI, jFor in the corresponding bipartite graph of H-matrix of quasi-cyclic LDPC code, link variable node vjWith check-node ciBetween
Side, be 1 when position is the element of (i, j) in H-matrix, show vjAnd ciBetween have Bian Xianglian, then its corresponding bI, j=1, it is no
Then bI, j=0,
Update checkpoint ciWith variable point vjAssociated each edge passes to ciLikelihood information and vjBetween priori mutual information:
IAc(i, j)=IEv(i, j)；
Update checkpoint ciPass to variable point vjLikelihood information and vjBetween external mutual information:σEcThe expression formula of (i, j) are as follows:
More new variables point vjWith checkpoint ciAssociated each edge passes to vjLikelihood information and vjPriori mutual information I beforeAv
(i, j): IAv(i, j)=IEc(i, j)
More new variables point vjPosteriority likelihood information and variable point vjBetween posteriority mutual information:
Wherein,
(33) judge whether to meet iteration stopping criterion:
If satisfied, then end step (3), export d, then using d as the distance threshold of the current round of corresponding protographIf
It is unsatisfactory for, then updates d=d+ Δ d, Δ d is preset increment step-length, then return step (32)；
(4) P is obtained to step (2)numIt is total to obtain P after a H-matrix executes step (3) respectivelynumA distance thresholdFrom this
PnumA distance thresholdIn find out the global optimum of epicycle iteration
(5) dual vector updates:
Update particle rapidity:
Wherein,
Wherein, λ, c1、c2、η1、η2For coefficient, wherein λ=1, c1=c2=2, η1、η2It is equally distributed in section (0,1)
Random number； The equally distributed random number between (0,1)；
(6) the number of iterations r=r+1 is updated, judges whether the number of iterations meets r > Imax；If satisfied, then stop iteration, according to working as
The corresponding dual vector of preceding particle constructs optimal QC-LDPC code matrix, and the global optimum that will be obtainedIt is assigned to
d；Otherwise, P is regenerated according to step (5) updated particlenumA H-matrix returns to (3) step.
2. a kind of magnetic induction communication system designed using design method described in claim 1 characterized by comprising
Transmitting terminal and receiving end；Transmitting terminal includes encoder, modulator and transmitting coil；Receiving end includes receiving coil, demodulation
Device and decoder；The distance between transmitting coil and receiving coil d are obtained by design method described in claim 1；
Transmitting terminal execute following steps: encoder using claim 1 construct quasi-cyclic LDPC code H-matrix to source data into
Row encodes, and after the data after coding are sent to modulators modulate, then arrives transmitting coil, transmitting coil generates induced current；
Receiving end executes following steps: excitation receiving coil generates induced current, then related data is sent to demodulator demodulation, solution
Serial data is sent to decoder again after tune, the data after being rebuild after decoder for decoding.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910709725.0A CN110311759B (en) | 2019-08-01 | 2019-08-01 | Magnetic induction communication system design method based on quasi-cyclic LDPC code and communication system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910709725.0A CN110311759B (en) | 2019-08-01 | 2019-08-01 | Magnetic induction communication system design method based on quasi-cyclic LDPC code and communication system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110311759A true CN110311759A (en) | 2019-10-08 |
| CN110311759B CN110311759B (en) | 2021-07-02 |
Family
ID=68082829
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910709725.0A Active CN110311759B (en) | 2019-08-01 | 2019-08-01 | Magnetic induction communication system design method based on quasi-cyclic LDPC code and communication system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110311759B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115801063A (en) * | 2022-10-21 | 2023-03-14 | 华中科技大学 | Underground magnetic induction wireless communication method, device and system |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104104393A (en) * | 2013-04-02 | 2014-10-15 | 盐城师范学院 | Quasi-cycle LDPC code design method with simple iterative code structure |
| US20190181882A1 (en) * | 2016-08-12 | 2019-06-13 | Telefonaktiebolaget Lm Ericsson (Publ) | Determining elements of base matrices for quasi-cyclic ldpc codes having variable code lengths |
-
2019
- 2019-08-01 CN CN201910709725.0A patent/CN110311759B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104104393A (en) * | 2013-04-02 | 2014-10-15 | 盐城师范学院 | Quasi-cycle LDPC code design method with simple iterative code structure |
| US20190181882A1 (en) * | 2016-08-12 | 2019-06-13 | Telefonaktiebolaget Lm Ericsson (Publ) | Determining elements of base matrices for quasi-cyclic ldpc codes having variable code lengths |
Non-Patent Citations (2)
| Title |
|---|
| YL GUAN: ""Finite-length EXIT Analysis and Design of Protograph LDPC Codes for Ultra-High-Density Magnetic Recording Channels"", 《IEEE》 * |
| 徐华: ""一种改进的基于有限域的准循环LDPC码设计方法"", 《微电子学与计算机》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115801063A (en) * | 2022-10-21 | 2023-03-14 | 华中科技大学 | Underground magnetic induction wireless communication method, device and system |
| CN115801063B (en) * | 2022-10-21 | 2024-05-14 | 华中科技大学 | Underground magnetic induction wireless communication method, device and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110311759B (en) | 2021-07-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106230489B (en) | Polarization code encoding modulation method suitable for any high order modulation | |
| CN109921882B (en) | Deep learning-based MIMO decoding method, device and storage medium | |
| CN1983861B (en) | Apparatus and method for transmitting/receiving data in a communication system | |
| CN105554865B (en) | A kind of MIMO-SCMA system down link design method based on STBC | |
| CN101141133B (en) | Method of encoding structured low density check code | |
| CN107026656A (en) | A kind of long Polar codes efficient coding method of short code in CRC auxiliary based on disturbance | |
| CN110504997B (en) | MIMO-based rapid coding LDPC optical communication system | |
| CN103916220A (en) | Network coding cooperative communication method based on polarization codes | |
| CN105897379A (en) | Polarization code cascade space-time code system and cascade polarization code coding method thereof | |
| CN109302267A (en) | Interpretation method, equipment and the storage medium of mimo system based on LDPC | |
| CN109889208A (en) | LDPC channel decoding method based on normalization minimum-sum algorithm in NAVDAT | |
| CN101997652A (en) | Acceptance detection method and device based on LDPC-MIMO (low density parity check-multiple input multiple output) communication system | |
| CN107294885A (en) | Allied signal detection and the method for estimation of channel in a kind of MIMO relay system | |
| CN102723975A (en) | Signal detection method and device of MIMO (multiple input multiple output) system | |
| CN110311759A (en) | A kind of magnetic induction Communication System Design method and communication system based on quasi-cyclic LDPC code | |
| CN106487446B (en) | It is suitable for the light space-time grid code encoding method of intensity modulated/direct-detecting mode | |
| CN106059594B (en) | A kind of combined channel state-detection and interpretation method based on state transition model | |
| CN100574127C (en) | A kind of BLAST receive-transmit system | |
| CN102487309B (en) | Signal detecting method and device under MIMO (Multiple Input Multiple Output) system | |
| CN104836607A (en) | Fast antenna selection method for united receiving and transmitting end of MIMO system | |
| CN100373841C (en) | Multiuser space hour group coding detection method | |
| CN103401647B (en) | The relay transmission method of LDPC code is extended based on bilayer | |
| CN106899388A (en) | Joint-detection and coding/decoding method of the LDPC code under mimo channel | |
| CN105827289A (en) | Limited feedback pre-coding method of generalized spatial modulation system | |
| CN109361491A (en) | Interpretation method, equipment and the storage medium of multi-input multi-output system based on LDPC |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |