CN108092931B - Time-delay multi-carrier modulation and demodulation method based on time-frequency pulse shaping - Google Patents

Time-delay multi-carrier modulation and demodulation method based on time-frequency pulse shaping Download PDF

Info

Publication number
CN108092931B
CN108092931B CN201711278898.9A CN201711278898A CN108092931B CN 108092931 B CN108092931 B CN 108092931B CN 201711278898 A CN201711278898 A CN 201711278898A CN 108092931 B CN108092931 B CN 108092931B
Authority
CN
China
Prior art keywords
wave
information
interval
wavelet
waves
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.)
Active
Application number
CN201711278898.9A
Other languages
Chinese (zh)
Other versions
CN108092931A (en
Inventor
梁德群
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN201711278898.9A priority Critical patent/CN108092931B/en
Publication of CN108092931A publication Critical patent/CN108092931A/en
Application granted granted Critical
Publication of CN108092931B publication Critical patent/CN108092931B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2602Signal structure
    • H04L27/2605Symbol extensions, e.g. Zero Tail, Unique Word [UW]
    • H04L27/2607Cyclic extensions

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

The invention discloses a method based on time frequencyThe pulse-forming time-delay multi-carrier modulation and demodulation method is characterized in that a plurality of periodic symbol waveforms capable of carrying information are constructed at a sending end; the symbol waveform is composed of H wavelets, the starting point of each wavelet is required to be placed in the life cycle of the first wavelet, and the symbol waveform is obtained by linear superposition after the starting point of the symbol waveform sequentially moves for a certain time period called time delay; the wavelet is a fundamental wavelet carrying information, and the fundamental wavelet comprises three intervals of waveforms which are respectively in an interval Tθi1Inner prepulse information wave in interval Tθi2Main wave sum without information in interval Tθi3An internal and a back pulse information wave; at a receiving end, the received symbol waveform is equalized, then amplitude information of information carried by the equalized symbol waveform is solved in each front-rear pulse information wave interval in sequence, and then the amplitude information is converted into binary number to complete a signal demodulation process. The invention has the characteristics of higher transmission efficiency, lower demodulation complexity and the like.

Description

Time-delay multi-carrier modulation and demodulation method based on time-frequency pulse shaping
Technical Field
The invention relates to a Modulation technology in digital communication, in particular to a Time-Delay Multi-carrier Modulation and Demodulation method (TFP-TDMC-A Time-Frequency Pulse-Shaping Based Time-Delay Multi-Carriers Modulation and Demodulation for short) Based on Time-Frequency Pulse Shaping, belonging to the technical field of Multi-carrier Modulation in digital communication.
Background
OFDM is a modulation technique in 4G as well as 5G mobile communications. The inventor has proposed a symbol waveform with a special structure and a corresponding modulation and demodulation method in a patent of time-frequency phase-mixing multi-carrier modulation method (application number: 2008101194121), and is called a time-shift orthogonal multi-carrier modulation technique (NMT or TS-NMT for short) later in published documents. The inventor has proved from both theory and experiment that the transmission efficiency is much higher than that of OFDM, and has the advantages of small peak-to-average power ratio, strong frequency shift resistance and no need of cyclic prefix. However, the demodulation complexity of TS-NMT is high, and when the wavenumber of wavelets is large, the demodulation equation set presents a large ill-conditioned behavior, so the inventor has proposed "a time-delay multi-carrier modulation and demodulation method", abbreviated as TDMC (patent application number: 2015100924095). TDMC keeps the characteristics of TS-NMT time delay multi-carrier, greatly simplifies demodulation method on the basis of inheriting the advantage of high transmission efficiency of TS-NMT, but has the defect of error propagation. In digital communication, the structure of a symbol waveform and a pulse shaping technique of the waveform are two important means for improving the efficiency of modulation. In terms of waveform structure, the multi-carrier structure has higher efficiency than the single carrier structure, for example, a composite wave (called QAM) composed of sine and cosine waves is more efficient than a single carrier of single sine or cosine waves, and OFDM composed of multiple QAMs has higher efficiency. The concept of pulse shaping was originally proposed in baseband transmission by nyquist, and demodulation was completed by shaping a transmitted square pulse into a SINC function waveform that does not cause interference at one point by an ideal filter that equalizes a communication system to a 1/(2T) bandwidth at a receiving end. Various pulse shaping methods in the time domain and the frequency domain appear later to meet different performance requirements in the time domain and the frequency domain. It should be noted that in the present invention, TS-NMT and TDMC are modified by using the idea of pulse shaping, and it is proposed that TFP-TDMC also uses multicarrier and pulse shaping as important means for improving transmission efficiency.
Disclosure of Invention
In view of the defects in the prior art, the present invention aims to provide a time-frequency pulse shaping-based time-delay multi-carrier modulation and demodulation method, so that the modulation technology has the characteristics of higher transmission efficiency, low demodulation complexity, etc.
In order to achieve the purpose, the technical scheme of the invention is as follows:
a time-delay multi-carrier modulation and demodulation method based on time-frequency pulse shaping is characterized in that:
at a sending end, binary information is converted into corresponding quantization amplitude and then sequentially given to a front information wave and a back information wave of a base wavelet to form corresponding wavelets to generate H wavelets, all wavelet starting points are required to be placed in the survival period of the first wavelet, all wavelets are sequentially moved for a time delay from the starting point of a symbol waveform and then linearly superposed to form a synthetic wave, namely, a periodic symbol waveform with T as a period is formed to complete a modulation process, namely, the method of sequentially moving for a time delay and then linearly superposed is taken as a time delay rule, different modulation methods have different time delay rules, and the difference is expressed as the rule of the time delay interval of adjacent waveletsDifferent; at a sending end, converting binary information into corresponding quantization amplitude, and sequentially giving a front information wave and a back information wave of the base wavelet to form corresponding wavelets; wherein the quantization amplitude complies with traditional value-taking rules in digital communications; the fundamental wave is composed of waveforms of three sections in which a front information wave, a main wave and a rear information wave are distributed in sequence, and the waveform structures of the front information wave and the rear information wave are completely the same; the wavelets are waveforms formed by respectively changing the amplitudes of a front information wave and a rear information wave of the fundamental wave to enable the front information wave and the rear information wave to respectively carry different binary information and enable the amplitude of the main wave to be unchanged; the three sections of the fundamental wave are respectively Tθ1、Tθ2And Tθ3The term "life time of the wavelet" means that the synthesis interval, in which the three intervals are sequentially combined into one fundamental wavelet, is called the life time of the fundamental wavelet, and the interval of the fundamental wavelet is also called the life time of the wavelet, and is denoted as Td,Td=Tθ1+Tθ2+Tθ3All intervals in the formula are half-open and half-closed intervals, the waveform of each basic wavelet only exists in the corresponding life cycle, and the basic wavelet is not defined or takes a value of 0 outside the life cycle; correspondingly, the tokens associated with different wavelets or base wavelets are distinguished by the subscript i ═ 1, …, H, i.e., TdiIndicating the lifetime of the ith wavelet or base wavelet such that Tdi=Td(i+1)Or Tdi≠Td(i+1)Generally, take Tdi=Td(i+1)=Td,Tθi1、Tθi2、Tθi3Respectively representing the interval of the ith front information wave, the main wave and the back information wave; and the fundamental wave is divided into two types, namely type 1 fundamental wave and type 2 fundamental wave which are expressed by formula
Figure BDA0001497145760000031
The corresponding wavelets are classified into type 1 wavelets and type 2 wavelets, and are formulated as
Figure BDA0001497145760000032
Superscript 1\2 for distinguishing type 1 baseThe wave or type 2 fundamental wave, the symbol "\\" represents the meaning of "or" in the following description, and the following formulas are expressed according to the above expression forms; wherein ω iskRepresenting the angular frequencies of the wavelet and the fundamental wavelet main information wave and the reference main wave, so that the wavelets with the same angular frequency form a subchannel in the frequency domain, and the center frequency of each subchannel is defined by the angular frequency omegakDetermining that K is 1, …, K, i.e. K subchannels are formed in total; and each sub-channel is not required to be orthogonal in frequency domain, and the number of co-frequency sub-waves contained in each sub-channel can be equal or unequal, that is, if the number of sub-waves with the same frequency contained in each sub-channel is denoted as LkK1, …, K, allowing Lk=Lk+1Or Lk≠Lk+1Angular frequency omegakThe subscript k and the wavelet subscript i have the following relation: wavelet sub-order i ═ L for each k-th sub-channelk-1+1,…,Lk-1+Lk
At the receiving end, firstly, the received waveform is restored to the waveform close to the transmitted waveform or the most suitable demodulation waveform as much as possible through equalization, and then the amplitudes of all wavelets are solved to complete the demodulation process.
Compared with the prior art, the invention has the beneficial effects that:
(1) compared with OFDM and other multi-carriers, TFP-TDMC keeps the characteristics of time-frequency phase mixed multi-carrier modulation method and TDMC high transmission efficiency on transmission efficiency;
(2) compared with the TDMC proposed by the inventor, as the TFP-TDMC information is only carried by the information wave, the amplitude of the main wave of the wavelet is not changed and is known during demodulation, the demodulation of each information wave is independently carried out, and the problem of error propagation is not introduced like the TDMC;
(3) the modulation and demodulation method of the dual overlapping of the packet delay in the TFP-TDMC has the similarity with a 'multi-carrier modulation method of time-frequency phase mixing', but the number of sub-waves contained in each group is reduced, the ill-posed characteristic of an equation set can be effectively reduced, and through the matching of alpha and beta values, the demodulation signal-to-noise ratio can be improved by improving the alpha value, so that the bit number carried by the information wave is increased, but the power of the synthesized wave is not improved.
(4) The information wave in TFP-TDMC occupies a much smaller interval than the wavelet period, and when the pipeline hardware structure is adopted for demodulation, the demodulation delay can be reduced, which is expected by 5G.
Drawings
FIG. 1 is an illustration of a transmitting end type 1 wavelet of the present invention;
FIG. 2 is an illustration of a transmitting end type 2 wavelet of the present invention;
FIG. 3 illustrates an exemplary type 2 wavelet through a 75 th order filter according to the present invention:
FIG. 4 is an illustration of a first exemplary transmission wave according to the present invention;
FIG. 5 is an illustration of a separate representation of sine and cosine waves of a wavelet pair in a second transmitted wave in accordance with the present invention;
FIG. 6 is an illustration of a synthesis of sine and cosine wavelet pairs in a second transmitted wave in accordance with the present invention;
fig. 7 is an illustration of a third transmission wave according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In the present invention, the term "symbol" is a term specific to a symbol waveform of one period, the period of the symbol being denoted by T, and T being used for each of three modulation methods to be described below1、T2And T3The symbol periods of the characters are represented, and other Pinyin or English letters (including the ones added with upper and lower marks) are called signs; in the following description, all the interval marks may be used as a mark of a time segment or a mark of a time point of an interval starting point depending on the application; the invention changes the amplitude of the information wave before and after the base wave to make them each otherThe basic wavelets carrying information are called wavelets, wherein the basic wavelets carry different binary information but do not carry information because the amplitude of the main wave is unchanged; let i be 1, …, H as a subscript to distinguish between different wavelets and the fundamental wavelets, thus, TdiIndicating the lifetime, T, of the ith wavelet or fundamental waveletθi1、Tθi2And Tθi3Respectively representing the interval of the ith front information wave, main wave and back information wave, and using T when different sub-waves are not needed to be distinguishedθ1、Tθ2And Tθ3These three intervals represent all wavelets; in the following formula, f and a are respectively used as subscripts or superscripts in different descriptions to distinguish front and back information waves, the information waves are general names of the front and back information waves, the information waves can be various waveforms such as square waves, raised cosine waves, sine waves or cosine waves, and the like, which are called pulse shaping in digital communication, however, the research of the inventor finds that the square waves and the raised cosine waves are used as the information waves, and a synthetic wave with high efficiency in both frequency domain and time domain is difficult to construct, and in order to simplify the description, in the following specific description, only type 1 and type 2 base waves of the information waves composed of the sine waves and the cosine waves respectively are included, and a designer can select other types of waveforms as the information waves according to different applications according to the following description;
based on the principle, the time-frequency pulse forming-based time-delay multi-carrier modulation and demodulation method is characterized by comprising the following steps of: at the transmitting end, the information is added to each base wavelet to become corresponding wavelet, that is, the quantization amplitude corresponding to the binary number representing the information is given to each base wavelet (as the common technology, the conversion step is omitted below), and then all the wavelets are combined into different synthetic waves to form a symbol waveform, thus completing the modulation process; at a receiving end, firstly, restoring the received waveform to a waveform close to the transmitted waveform or an optimum demodulation waveform as much as possible through equalization, and solving the amplitudes of all wavelets to complete the demodulation process; because the equalization can only obtain a waveform similar to the transmission, in order to distinguish different shapes, subscripts t and r are used for distinguishing a transmitting end from a receiving end, and t \ r is used for representing by an abbreviation of t or r; in the following, x \ y is an abbreviated representation of x or y, and x or y can be replaced by any number and letter;
the basic structure of the fundamental waves of types 1 and 2
Figure BDA0001497145760000051
Wherein u (t) is as defined above1)|v(t2)|w(t3) The form gives the sequential connection of three waveforms separated by vertical lines marked with "|", u (t)1)、v(t2) And w (t)3) The following general terms refer to three different waveforms, and the listed waveform notation may be substituted for any waveform in the same representation in the formula of the present invention, specifically in formula (1):
Figure BDA0001497145760000052
and
Figure BDA0001497145760000053
all the intervals are half-open and half-closed intervals, subscripts 1, 2 and 3 of the intervals give a connection sequence on a time coordinate, superscripts 1\2 are used for distinguishing type 1 or type 2 basewaves,
Figure BDA0001497145760000054
and
Figure BDA0001497145760000055
referred to as front or rear information waves, distinguished by subscript f or a respectively,
Figure BDA0001497145760000056
referred to as the reference dominant wave, i is the index of the different fundamental wavelets, and τ i is the time delay of the ith fundamental wavelet relative to the origin of a symbol (τ when a symbol contains multiple sub-composite waves)iIs rewritten as
Figure BDA0001497145760000057
Is indicated as
Figure BDA0001497145760000058
The fundamental wave is relative to the sub-synthetic waveThe time delay of the point),
Figure BDA0001497145760000059
meaning that the intervals of the front and back information waves are equal, subscript t \ r represents a sending end or a receiving end, and delay expansion is generated when the sending waveform reaches the receiving end, so that a section of empty interval is reserved for the waveforms of the three intervals of the sending end, and delta taui=τi+1iRepresenting the time delay interval between the (i + 1) th wavelet and the ith wavelet, typically by Δ τi=Δτh=Δτ,i&H1, …, H, notation&Expressing the sum, Δ τ generally refers to the time delay interval, and τ is specifiedi+1i、τ1=0、τi∈Td1That is, the origin of the 1 st wavelet coincides with the origin of the symbol, and the origins of the other wavelets are all within the lifetime of the 1 st wavelet, allowing Tdi≠Tdh,Δτi≠ΔτhOr Tdi=Tdh=TdWhen equal, the subscript may be omitted, and when T ═ 2TdOr T ═ 2TdΔ τ, in general, taking Tdi=Tdj,Δτi=Δτj,T=2Td(ii) a The time delay tau is required in the above formulaiThe use of (b) is specifically explained: definition of tauiIs the time delay of a base wavelet from a symbol starting point, the time delay of a front information wave is the time delay of the base wavelet, but the time delay of a main wave should be larger than that of the front information wave, the time delay of a back information wave should be larger than that of the front main wave, so that three-section waveforms in the base wavelet should be expressed as
Figure BDA0001497145760000061
However, such representation is too complicated to read, and therefore, in the present invention, only the time delays between the fundamental waves are distinguished, the time delays of the fundamental waves for all three waveforms in the fundamental waves are not distinguished, and different time delays are distinguished by the sequential connection relationship, that is, the sequential connection relationship given by the symbol "|" of the vertical line in the above expression (1).
For ease of design, the respective regions of type 1 and type 2 waveletsThe same value is used for the interval, so the superscript 1\2, omega of the interval mark can be omittedkThe angular frequency of main information wave representing wavelet and base wavelet and reference main wave, the wavelets with same angular frequency form a sub-channel in frequency domain, each sub-channel can be orthogonal or non-orthogonal in frequency domain, the same-frequency wavenumbers of sub-channels can be equal or not, and the angular frequency omegakIs the center frequency of the K-th 1, …, K subchannels, i.e. K subchannels are formed in total; in the following description, the frequency of the information wave is the same as the frequency of the main wave, and the spectrum of the sub-channel is a composite of the information wave and the main wave spectrum; let L be the number of wavelets containing the same frequency in each subchannelk,Lk=LjOr Lk≠Lj,k&j is 1, …, K, meaning that the co-channel sub-wave numbers contained in each sub-channel may be equal or unequal; angular frequency omegakThe subscript k and the wavelet subscript i have the following relation: the same-frequency wavelet sub-sequence i of k-th sub-channel is Lk-1+1,…,Lk-1+LkIn all the following descriptions, the subchannel indices will be used for the case where there are sub-composite waves
Figure BDA0001497145760000062
Or
Figure BDA0001497145760000063
Denotes k and i,
Figure BDA0001497145760000064
And
Figure BDA0001497145760000065
and
Figure BDA0001497145760000066
and
Figure BDA0001497145760000067
all of the same; the sending end type 1 wavelet shown in FIG. 1 has a pre-information wave interval Tθt11750, the front information wave is the expansion reserved interval T of the receiving endθt1275 (meaning)The order of the filter is 75), the main wave interval Tθt214500, the main wave is the expansion reserved interval T of the receiving endθt2275, the post-information wave interval Tθt31750, the post information wave is the expansion reserved interval T of the receiving end θt3275, 2, 0.1, wavelet lifetime Td1=6225,ω 12 pi/750 x 6 x pi/4500 x 0.08378; the transmitting end type 2 wavelet shown in FIG. 2, the front master information wave interval Tθt112750, the front main information wave rear end moves to the section T of the copy part of the front information wave front endθt11175, the front end of the front main information wave moves to the section T of the copy part of the rear end of the front information waveθt11375, the front information wave is the expansion reserved interval T of the receiving endθt12Main wave interval T75θt24500, the main wave is the expansion reserved interval T of the receiving endθt2275, the post-master information wave interval Tθt312When the post-main information wave front end moves to the section T of the copied part of the post-main information wave front end, 750θt311When the post-main information wave front end moves to the section T of the copy part of the post-information wave rear end, 75θt31375, the rear information wave is the expansion reserved interval T of the receiving endθt32=75,α=2,β=0.5,
Figure BDA0001497145760000071
Angular frequency omega of information wave and main wave 12 pi/750 2 x 8 x pi/4500 0.08378, wavelet lifetime Td16225. In the present document, the specific values of the intervals in the various illustrations are represented by dimensionless sampling points, and the corresponding specific values of the frequencies are also calculated values using the points as time units, so that the calculated results can be easily converted into actual time or frequency values by a designer as long as the actual time units of the sampling intervals are given.
More specifically, 1\2 type waves of the transmitting end and the receiving end are represented as
At the transmitting end
The information wave and the main wave of the type 1 fundamental wave are represented as waveforms of expressions (1-1) and (1-2), respectively:
Figure BDA0001497145760000072
Figure BDA0001497145760000073
Figure BDA0001497145760000074
Figure BDA0001497145760000075
z(t(1\3)2i)=z(t21i)=0
the information wave and the main wave both include two intervals, interval Tθti(1\3)1And Tθti(1\3)2The interval T of the information wave and the main waveθti(1\3)2And Tθti22The subscript (1\3) is used for distinguishing two abbreviated forms of intervals of front or rear information waves, but on a time axis, the two intervals are respectively distributed at the front end and the rear end of a main wave;
the type 2 fundamental wave information wave and the dominant wave are represented by waveforms of the formulae (1-3) and (1-4)
Figure BDA0001497145760000081
Figure BDA0001497145760000082
Figure BDA0001497145760000083
Scale formula (1-3-1) middle interval Tθti(1\3)12The waveform of (a) is a main information wave,
Figure BDA0001497145760000084
is sin omegak(t(1\3)12i) A replica of the waveform of the back end portion,
Figure BDA0001497145760000085
is sin omegak(t(1\3)12i) A replica of the waveform of the front-end portion,
Figure BDA0001497145760000086
is cos omegak(t(1\3)12i) A replica of the waveform of the back end portion,
Figure BDA0001497145760000087
is cos omegak(t(1\3)12i) A replica of the waveform to which the front portion corresponds, where Tθti(1\3)1=Tθti(1\3)11+Tθti(1\3)12+Tθti(1\3)13The waveform of the front and back end parts is similar to the cyclic prefix in OFDM, except that in OFDM, only sin ω is usedk(t(1\3)12i) And cos omegak(t(1\3)12i) The section width of the copied part is determined by the channel length, while in the present invention, the section width of the copied part is determined by the order TP of the equalization filter, T is a bidirectional copy in which the front end is copied to the back end and the back end is copied to the front endθti(1\3)11=Tθti(1\3)13=TP/2;
Figure BDA0001497145760000088
z(t(1\3)2i)=z(t22i) 0 represents Tθti(1\3)2And Tθti22The value of the waveform in the receiving end is 0, and a reserved interval is expanded for the receiving end;
at a transmitting end, a modulation process of a symbol waveform comprises: firstly, converting binary information into corresponding quantization amplitudes according to a general method of digital communication, sequentially endowing the quantization amplitudes to front and rear information waves of each base wavelet to form a type 1 wavelet or a type 2 wavelet of a transmitting end as described in the following formula (2),
Figure BDA0001497145760000089
in the formula (2), the reaction mixture is,
Figure BDA0001497145760000091
respectively the quantized amplitude of the preceding and following information waves,
Figure BDA0001497145760000092
y belongs to integer field, that is, the amplitude of front and back information wave can select one quantization amplitude from Y grades, which is a universal method in digital communication; note that the main wave in equation (2) is not given quantization amplitude, which indicates that only front and back information waves carry information in one wavelet, and the main wave does not carry information, so that the TFP-TDMC avoids the problem of error propagation in TDMC; alpha and beta are amplitude adjustment parameters which are fixed and unchangeable under the condition that the system operation environment is unchanged after the design is finished, and the designer performs certain adjustment only when needed; in order to improve the signal-to-noise ratio of the information wave, alpha is required to be more than or equal to 1, and only beta is required to be taken so that the amplitude of the main wave does not exceed the maximum value allowed by a system, under different conditions, the different values of alpha and beta are matched to bring larger difference of the power of the synthesized wave, a proper value is selected, alpha is taken as a larger value under the condition of not increasing the power, so that the signal-to-noise ratio of the information wave is increased, and an optimized value can be obtained through computer simulation; t is t1\3∈Tθti(1\3)Is an abbreviated form, which indicates that different time variables belong to different intervals to indicate different positions of the front and back information waves on the time axis, and is used for indicating the position of the front and back information waves on the time axis
Figure BDA0001497145760000093
Indicating different quantized amplitudes of the front or back information waves;
further, the composite wave of the transmitting end is composed by selecting one of the following three methods to complete modulation: for convenience of description, in the following description of the three methods, similar symbols to those described above are used, and only symbols given special meanings are additionally described herein, and no special description follows the meaning described above;
the method comprises the following steps: direct addition modulation method based on type 1 wavelet
Specifically, the 1-type wavelets with different time delays are superposed to form a composite wave, as shown in the following formula
Figure BDA0001497145760000094
τi∈Td110 means that all the H wavelets have their origins in the lifetime of the first wavelet and the time delay interval is equal to the information wave interval, i.e. Δ τ Tθti1=Tθti3That means that all the information waves occupy the corresponding delay interval sections without overlapping each other, but there is a partial overlap between the information waves and the main wave and between the main wave and the main wave;
Figure BDA0001497145760000101
the expression of (2) is shown in the formula, and the corresponding basic wavelets are type 1 basic wavelets; as in the first transmit wave illustration of fig. 4: this is a superposition of 7 type 1 wavelets, only the three wavelets 1, 2 and 7 are shown in the figure, which are respectively indicated by solid lines, broken lines and dotted lines, and the information waves are not overlapped; for clarity of illustration, the lifetimes T of the 1 st wavelets are omittedd16225, all wavelets α -2, β -0.1, ai1 is ═ 1; the first transmission wave is referred to as equation (2-1).
The method 2 comprises the following steps: modulation method based on type 2 wavelet pair
The method is specifically divided into three steps:
the first step is as follows: constructing same-frequency sine-cosine wavelet pairs on the basis of formulas (1-3), (1-4) (1-3-1) and (1-4-1)
Figure BDA0001497145760000102
Figure BDA0001497145760000103
Figure BDA0001497145760000104
Scale formula (2-3-1) middle interval Tθti(1\3)12The waveform of (a) is a master information wave pair and orthogonal sine and cosine pairs,
Figure BDA0001497145760000105
is sin omegak(t(1\3)12i) A replica of the waveform of the back end portion,
Figure BDA0001497145760000106
is sin omegak(t(1\3)12i) A replica of the front-end portion waveform,
Figure BDA0001497145760000107
is cos omegak(t(1\3)12i) A replica of the waveform of the back end portion,
Figure BDA0001497145760000108
is cos omegak(t(1\3)12i) A replica of the front-end portion waveform; wherein T isθti(1\3)1=Tθti(1\3)11+Tθti(1\3)12+Tθti(1\3)13The partial waveform of the front end and the back end is similar to the cyclic prefix in OFDM, except that in OFDM, only sin omega is usedk(t(1\3)12i) And cos omegak(t(1\3)12i) The section width of the copied part is determined by the channel length, while in the present invention, the section width of the copied part is determined by the order TP of the equalization filter, T is a bidirectional copy in which the front end is copied to the back end and the back end is copied to the front endθti(1\3)11=Tθti(1\3)13=TP/2
The second step is that: constructing a type 2 sub-synthetic wave consisting of the same-frequency type 2 sub-wave pairs according to a time delay rule:
Figure BDA0001497145760000111
Figure BDA0001497145760000112
representing the lifetime of the first wavelet in the sub-composite wave, with a delay interval equal to the information wave pair interval, i.e.
Figure BDA0001497145760000113
The sub-synthetic waves all comprise wavelet angular frequencies omegak
Let K be 1, …, K then obtain K sub-composite waves of the formula (2-3-3), and all pairs of main information waves with the same time delay in the respective sub-composite waves are defined to be orthogonal in the frequency domain, i.e., for K, j be 1, …, K,
Figure BDA0001497145760000114
and
Figure BDA0001497145760000115
interval Tθti(1\3)12Are orthogonal.
The third step: composition type 2 synthetic wave:
Figure BDA0001497145760000116
the formula (2-3-4) is called as a second transmitting wave; fig. 5 and 6, where fig. 5 is a graphical illustration of a sine and cosine wave separated (not added together) representation of a wavelet pair in a second transmitted wave: the solid line in the figure represents sine wavelets, the dotted line represents cosine wavelets, and the parameter values of each interval are the same as the type 2 wavelets in FIG. 2; FIG. 6 is an illustration of a wavelet pair obtained by adding sine and cosine wavelets in a second transmitted wave, each interval having the same value of parameter as the type 2 wavelet in FIG. 2. The arrangement of wavelets given in the table below may help to understand method 2 more clearly:
Figure BDA0001497145760000117
the table merely gives an example of the positional relationship among wavelets for pairs of information waves of different frequency bands (note that this is only an example of the positional relationship and is not a description of the complete waveform), where
Figure BDA0001497145760000118
Representing a sine-cosine pair of information waves, the main pairs of information waves in all columns being orthogonal in frequency domain and superposed without delay, the pairs of information waves in a row in the transverse direction having the same frequency and being arranged in such a way that
Figure BDA0001497145760000121
Given the superposition of delays (which may actually allow different frequencies to be taken and do not require orthogonality, but are cumbersome to design).
The method 3 comprises the following steps: modulation method of packet delay double overlapping based on 1-type wavelet
The method is to compose a synthetic wave from wavelets of type 1, and to divide all wavelets into U groups, each group of wavelets being
Figure BDA0001497145760000122
The starting points of the wavelets in the u-th group are sequentially shifted by a time delay in the first previous information wave interval of the group
Figure BDA0001497145760000123
Referred to as intra-group delay, and further having an intra-group delay interval
Figure BDA0001497145760000124
Thereby making
Figure BDA0001497145760000125
The origin of the type 1 wavelet being in the front information wave interval T of the first wavelet of the grouptθi1With internal sequential time delay of one
Figure BDA0001497145760000126
Post-superposition, the individual information waves within a group being partially overlapping, whereupon
Figure BDA0001497145760000127
The front information wave of the 1-type wavelet occupies an interval of
Figure BDA0001497145760000128
This is also called information band, and the main wave and the post-information wave are overlapped and added in the same way, and the post-information wave occupies the interval
Figure BDA0001497145760000129
This constitutes the u-th composite wave, called the sub-composite wave, which is the first wave; defining interclass delay intervals
Figure BDA00014971457600001210
Then, each sub-synthesis is sequentially moved by an inter-group time delay interval and then added to form a symbol synthesis wave, which is the second time, the information wave bands of each sub-synthesis wave are not overlapped, but the main wave band is still partially overlapped; in methods 1 and 2, the time delay interval between the wavelets is greater than or equal to the information wave interval, i.e. τi+1i=Δτi≥Ttθi(1\3)And for method 3, the sum U is determined at H>Under the condition of 1. the method comprises the following steps,
Figure BDA00014971457600001211
the following formula is a sub-composite wave formed by overlapping and adding information waves in the u group:
Figure BDA00014971457600001212
the corresponding wavelets have the same basic structure as equation (2), except that they follow the subscript
Figure BDA00014971457600001213
With variations, subscripts
Figure BDA00014971457600001214
Indicates the u in the group
Figure BDA00014971457600001215
The sign of the correlation of the individual wavelets,
Figure BDA00014971457600001216
is the number of sub-waves within a set,
Figure BDA00014971457600001217
is each wavelet lifetime of the u-th group, but the subscript u may distinguish the positional variation of the wavelet lifetimes of the different groups, T in widthd(u,1)=TdBut the position varies with U-1, …, U,
Figure BDA0001497145760000131
is at the position of
Figure BDA0001497145760000132
And U is 1, …, varying from U,
Figure BDA0001497145760000133
indicating that the range of action of the time variable is extended by a period of time longer than the wavelet lifetime
Figure BDA0001497145760000134
Figure BDA0001497145760000135
Time-shifting intervals for the information waves within the group and defining
Figure BDA0001497145760000136
Allow for
Figure BDA0001497145760000137
Or
Figure BDA0001497145760000138
In order to reduce the design complexity, the design is typically done in an equal manner,
Figure BDA0001497145760000139
an interval representing a front information wave of the first wavelet in the u-th group,
Figure BDA00014971457600001310
meaning that the start points of all wavelets belonging to the u-th group are in the interval
Figure BDA00014971457600001311
In the above, the superscript f corresponds to the meaning of the front information wave, and it is noted that here, in order not to make the subscript number too much, the subscript f is changed to the superscript f, Δ τu=τu+1uThe time delay interval between groups is also the interval occupied by all information waves of the u-th group; the interval distribution position and the symmetry characteristic of the information waves before and after the wavelet structure are determined, and after the relative position between the current information waves is determined, the starting points of the information waves after all the wavelets in the group are naturally in the interval of the information waves after the first wavelet in the group; all the previous information waves in the group are in the first previous information wave interval and sequentially move by a time delay, and then the whole wavelet corresponding to the previous information wave complies with the time delay rule; the sub-composite wave thus composed is divided into three intervals: the first interval is
Figure BDA00014971457600001312
A partial overlap interval of the front information wave, which also includes a partial main wave component and is called as the overlap interval of the front information wave and marked as Tu1The second interval is the main wave overlapping region, denoted as Tu2And the third interval is
Figure BDA00014971457600001313
The post information wave part overlaps the interval, also includes part of main wave component, marked as Tu3
Then, the sub-synthesized waveforms of all groups are sequentially shifted by a time delay tauu(U-1, …, U) are added to form the final composite wave, as shown in the following equation
Figure BDA00014971457600001314
Based on the symmetry of the wave structure of the front and back information, the grouping of the wavelets will naturally result in the front and back informationThe same grouping of waves, so that the expressions (2-4-1) and (2-4-2) do not particularly use a representation for distinguishing between preceding and following information waves, H being the total wavenumber, where τuIs a group u delay, so that the information waves of the wavelets in each group have a delay in turn
Figure BDA0001497145760000141
Making the information waves in the u-th group partially overlapped, but the information waves between the groups are not overlapped, and only the main waves still keep partially overlapped; the formula (2-4-2) is called as a third transmitting wave; third transmission wave illustration as shown in fig. 7: two sets of overlapping information waves are shown, each set comprising 5 partially overlapping information waves, two adjacent sets of information waves do not overlap, only the information waves in the sets overlap, the time delay interval of the information waves in the sets is 150 points, the two sets of sub-waves have the same angular frequency omegau1In practice, each group may have a different frequency, which is not clearly visible in the figure, the wavelets are not added up in order to make the wavelet distribution clearly visible in the figure, and only the first two groups are shown, and two groups of sub-composite waves may be arranged in the remaining interval, which are omitted for clarity of the drawing. For the above three methods, the symbol periods are T respectively1、T2、T3
At a receiving end, firstly, restoring the waveform subjected to channel distortion to be close to a transmitting waveform as much as possible or balancing the waveform to be a waveform suitable for demodulation through balancing; in general digital communication, equalization needs to follow the principle of obtaining the maximum signal-to-noise ratio of demodulation or the minimum distortion of a waveform, the invention allows the adoption of various existing equalization methods, but follows the principle of receiving the minimum distortion of the waveform or the most suitable demodulation of the waveform, the expanded width of an information wave is ensured not to exceed the wavelet interval, and the specific expanded width is determined by training; the specially proposed band-limited inverse filter equalization method of the present invention is an equalization method for achieving the above-mentioned purpose by controlling the parameters of the equalization filter.
Specifically, the mathematical expression of the band-limited inverse filtering equalization method is expressed as
Figure BDA0001497145760000142
HP(u)=BP(u)/H(u) (3-1)
Figure BDA0001497145760000143
hp(t)=IFFT[BP(u)/H(u)] (3-2)
Equations (3-1) and (3-2) are the frequency and time domain representations, respectively, of band-limited inverse filtering, where Sr(u)/H (u) represents inverse filtering, which has a certain ill-conditioned property, and the addition of the filter BP (u) effectively reduces the ill-conditioned property, u and t are frequency domain and time domain variables respectively,
Figure BDA0001497145760000144
and
Figure BDA0001497145760000145
representing equalized frequency-domain and time-domain waves, St(u) and st(t) symbol waveforms, S, representing frequency and time domain transmissions, respectivelyr(u) and sr(t) represents the received symbol waveform in frequency domain and time domain, respectively, H (u), and N (u) are the channel transfer function or channel model, low pass or band pass filter, and noise in frequency domain, respectively, h (t) and bp (t) are the channel impulse response and impulse response of band pass or low pass filter, respectively, in time domain, are the convolved symbols,
Figure BDA0001497145760000151
is a deconvolution symbol, [ N (u)/H (u)]BP (u) and
Figure BDA0001497145760000152
is a frequency domain time domain representation of the equalized noise term. The present invention requires the filter bp (u) to have symmetry, and requires the bandwidth of the filter to be smaller than the system-limited bandwidth to reduce the ill-conditioned nature of the inverse filtering, and does not limit other specific forms, such as the FIR or FIR filter, or other specially designed pulse or spectrum shaping filters, and also requires the filter to have symmetry, and the shaping filters are also called equalization filtering filters in the followingA machine; due to the band-limiting effect of the filter, so that sr(t) relative to st(t) in the time domain, the spreading widths of the waveforms caused by them are collectively called shaping spreading or equalization spreading, for FIR filters, the shaping spreading is equal to the order of the filter, for specially designed filters the shaping width is determined by the design parameters, for convenience of description, this parameter is also called the order, denoted TP, as illustrated in fig. 3 for the wavelet of type 2 passing through a 75 th order filter: front end extended wave interval T of front information waveθr1175, the former master information wave interval Tθr12750, the rear end of the front information wave is spread the wave interval Tθr1375, main interval T without expansionθr214500, its extension interval Tθr2275, front end spreading wave section T of rear information waveθr3175, the post-master information wave interval Tθr32750, rear end spreading wave section T of rear information waveθr33Three spreading waveforms can be seen in three spreading intervals, called tail waveforms, α is 2, β is 0.5, ai1, due to the trailing waveform, some high-frequency components are added in the side lobe of the wavelet spectrum, but the central angular frequency is not changed by omegak=2×6×π/4500=0.08378;
After equalization, the corresponding base wavelet waveform and base wavelet of the receiving end
Figure BDA0001497145760000153
The equation (4-1) corresponds to the equation (1-1), and the difference is only that the reserved empty interval of each segment of waveform is occupied by the expanded trailing waveform, so that the waveform can only be an approximate sine/cosine wave, and correspondingly, the subscript t representing the transmitting end in the equation is changed into the subscript r representing the receiving end.
The information wave and the main wave of the type 1 fundamental wave are expressed by the formulas (4-1-1) and (4-1-2)
Figure BDA0001497145760000154
Figure BDA0001497145760000155
The information wave and the main wave of the type 2 fundamental wave are expressed by the formulas (4-1-3) and (4-1-4)
Figure BDA0001497145760000156
Wherein,
Figure BDA0001497145760000157
Figure BDA0001497145760000161
and
Figure BDA0001497145760000162
are no longer partial replicas of sine and cosine waves, but instead are deformation waves with time delay spread of these replicas;
Figure BDA0001497145760000163
wherein,
Figure BDA0001497145760000164
the corresponding wavelet is
Figure BDA0001497145760000165
A in the formula (4-2)fiAnd aaiIs the actual amplitude; the composite wave of a symbol at the receiving end is composed of wavelets shown in (4-2); specifically, for different transmitted symbol waveforms shown by the formulas (2-1), (2-3-4) and (2-4-2), corresponding different symbol composite wave waveform structures exist after equalization is performed on a receiving end, so that the following three demodulation methods are correspondingly adopted; the receiving waveforms of different structures are respectively as follows: compared with the transmitting wave, the subscript t is changed into r, and other marks have the same meaning;
the first method comprises the following steps: received waves corresponding to the first type of transmitted waves:
Figure BDA0001497145760000166
the corresponding wavelets and the fundamental waves are shown as formulas (4-2) and (4-1);
and the second method comprises the following steps: received wave corresponding to second transmitted wave
Figure BDA0001497145760000167
Figure BDA0001497145760000168
Wherein,
Figure BDA0001497145760000169
Figure BDA00014971457600001610
Figure BDA00014971457600001611
and
Figure BDA00014971457600001612
no longer are partial replicas of sine and cosine waves, but instead are deformation waves with extended time delay of those replicas, with only intervals
Figure BDA0001497145760000171
The waveforms in (a) are orthogonal,
Figure BDA0001497145760000172
and
Figure BDA0001497145760000173
is the actual amplitude of the wavelet at the receiving end, sinceA certain difference is caused between the quantization amplitude and the noise and the balance error, and finally the quantization amplitude is classified into the corresponding quantization amplitude according to a certain principle (such as 4 SHE 5 input);
Figure BDA0001497145760000174
the main wave pair is changed into approximate sine wave and cosine wave with time delay expansion from sine wave and cosine wave
And the third is that: received wave corresponding to third transmitted wave
Figure BDA0001497145760000175
Figure BDA0001497145760000176
For the three received waveforms, there are three demodulation methods, which are respectively:
the method comprises the following steps: for coherent demodulation of the first received wave, specifically,
in the interval Tθri1And Tθri3Respectively performing operations in the following steps:
Figure BDA0001497145760000177
Figure BDA0001497145760000178
Figure BDA0001497145760000179
the expression (6-1-3) means that the main wave part which does not carry information is subtracted in the information wave interval, and only the information wave which carries information is left, wherein the mark x in the expressions (6-1-1) and (6-1-2) is a vector dot product;
when i is 1, …, H, the amplitude of the information wave before and after each wavelet corresponding to the symbol waveform can be obtained by repeating the equations (6-1-1) and (6-1-2), that is, the demodulation of the first type of received wave is completed;
the method 2 comprises the following steps: for the demodulation of the second type of received wave, the specific steps are as follows:
the first step is as follows: dividing the waveform of formula (4-3-1) into
Figure BDA00014971457600001710
Section, front
Figure BDA00014971457600001711
Each section interval in a section is Trθi1After, after
Figure BDA0001497145760000181
Each section interval in a section is Trθi3Wherein, the first
Figure BDA0001497145760000182
The waveform within a segment may be represented as
Figure BDA0001497145760000183
Figure BDA0001497145760000184
Is the corresponding dominant wave component;
the second step is that: subtracting the main wave component in (6-2-1)
Figure BDA0001497145760000185
Thus, in the formula (6-2-2)
Figure BDA0001497145760000186
Intermediate interval T of segmentrθi(1\3)2Is a combination of K orthogonal pairs of sine-cosine waves,
the third step: fast Fourier operation is performed on the waveform shown in the formula (6-2-2)
Figure BDA0001497145760000187
K front or back information wave amplitudes can be obtained
Figure BDA0001497145760000188
Superscripts f and a indicate the front or back information wave, respectively; order to
Figure BDA0001497145760000189
Repeating the formulas (6-2-1) to (6-2-3) to obtain all information wave amplitudes;
the method 3 comprises the following steps: for the packet demodulation of the third type of received wave, specifically, the following four steps are performed:
the first step is as follows: the received wave represented by the formula (5-3-1)
Figure BDA00014971457600001810
Carrying out segmentation interception to obtain waveforms of 2U groups of sections, wherein the front U group of sections comprises all front information waves and main wave components in corresponding intervals, and the rear U group of sections corresponds to all rear information waves and main wave components in corresponding intervals; each group of segment waveforms is expressed as follows
Figure BDA00014971457600001811
Wherein the superscript f \ a indicates the front or back information wave, TruThe interval representing the waveform of the u-th group of segments
Figure BDA00014971457600001812
The information wave interval before and after the wavelet is
Figure BDA00014971457600001813
The second step is that: in the interval
Figure BDA00014971457600001814
All main wave components are removed:
Figure BDA00014971457600001815
the third step: obtaining a demodulation equation set
In the intervalPerforming coherent operation as shown in the following formula:
Figure BDA0001497145760000192
order to
Figure BDA0001497145760000193
Repeating the formula (6-3-3) to obtain the coherent operation column vector
Figure BDA0001497145760000194
Further form a linear equation system
Figure BDA0001497145760000195
This is called the demodulation equation set, where
Figure BDA0001497145760000196
For the magnitude column vector to be solved,
Figure BDA0001497145760000197
in order to coherently operate on the column vectors,
Figure BDA0001497145760000198
Figure BDA0001497145760000199
in order to demodulate the matrix, the demodulation matrix,
Figure BDA00014971457600001910
the subscripts in the above formulas are all increased by u as group differences;
of a matrixThe element is obtained by coherent operation of the fundamental wave obtained by training and a part of waveforms corresponding to other fundamental waves in the interval; where, x is the vector point-by-sign, ωkAnd ωjIndicating that different information waves may have different angular frequencies, k&j=1,…,K;。
The fourth step: solving the equation set shown in equation (6-3-4)
Figure BDA00014971457600001911
Obtaining amplitude vectors of the U-th group of front/back information waves, enabling U to be 1, … and U, and repeatedly performing the formulas (6-3-3) and (6-3-4) to obtain the amplitude vectors of the front/back information waves of all groups;
further, in order to implement demodulation, when the amplitude of information carried by a symbol waveform is solved, firstly, the correlation parameters of a fundamental wave required for demodulation and a demodulation matrix required for demodulation of a third received wave need to be obtained through training; specifically, the obtaining of the relevant parameters of the training base wavelet is obtained through band-limited inverse filtering training; meanwhile, because the same-frequency wavelets have the same training result, each subchannel only takes one base wavelet, and because the training is carried out on a single base wavelet, the time delay tau of each wavelet does not need to be considerediThe subscripts i and k may be modified
Figure BDA00014971457600001912
And
Figure BDA00014971457600001913
or
Figure BDA00014971457600001914
And
Figure BDA00014971457600001915
the same applies to the following training steps;
the band-limited inverse filtering training is characterized in that:
(1) training for 1\2 type base wavelets: for K1, …, K, the following is done:
selecting a 1\2 type base wavelet from the k-th sub-channel
Figure BDA0001497145760000201
Let 1\2 type base wave
Figure BDA0001497145760000202
Is formulated as
Figure BDA0001497145760000203
Then in the interval
Figure BDA0001497145760000204
From
Figure BDA0001497145760000205
Is separated out
Figure BDA0001497145760000206
Rear information wave
Figure BDA0001497145760000207
And a reference main wave
Figure BDA0001497145760000208
t2∈Tθrk2Three sections of waveforms are used for completing the training of a fundamental wave;
making K equal to 1, … and K, and repeating the first step to obtain the trained waveforms of the 1\2 type base wavelets of all the sub-channels;
(2) training of demodulation matrix required for demodulation of third received wave
The demodulation matrix in the rewritten equation (6-3-4) is as follows
Figure BDA0001497145760000209
Figure BDA00014971457600002010
And
Figure BDA00014971457600002011
the angular frequencies of the wavelets in the u-th group, i.e. different indices, to indicate the crossing between the information waves
Figure BDA00014971457600002012
As long as it is obtained
Figure BDA00014971457600002013
Can form a demodulation matrix to obtain
Figure BDA00014971457600002014
Is mainly to do
Figure BDA00014971457600002015
A about
Figure BDA00014971457600002016
The operation of (1); for this purpose, the following matrix is constructed
Figure BDA00014971457600002017
Figure BDA00014971457600002018
Elements in a matrix
Figure BDA00014971457600002019
Is that
Figure BDA00014971457600002020
To (1) a
Figure BDA00014971457600002021
And (4) sampling points, wherein 0 in the matrix is a time delay interval in the front of the wavelet, and in order to form the matrix, the tail part of the right side of each row needs to be correspondingly supplemented with 0.
Further, the band-limited filteringThe device BP (u) is designed according to the following principle: i. filter bandwidth W of the filterbpW ≦ W, W is the symbol bandwidth, and W isbp>>Wsc,WscThe order TP of the filter determines the trailing length of an information wave, if the TP is large, the trailing is lengthened, the time delay interval is enlarged, the number of sub-waves in a symbol is reduced, and the transmission rate is influenced, so that the problem of engineering compromise is solved, and the optimal value can be selected through computer simulation.
The specific implementation method of the technical scheme of the invention comprises the following steps:
the technical principle and some key technologies of implementation of the present invention are further described below by taking an implementation of the TFP-TDMC for wireless broadband access as an example; the scheme mainly describes special technical problems in the implementation of the TFP-TDMC, and conventional technologies in digital communication are not explained in detail.
The following embodiments are all referenced to the 802.11a international standard. The following environmental parameters are specifically taken: when the channel bandwidth of 5MHz, the number of subcarriers is 52, the data carrier is 48, the OFDM symbol duration is 16 mus, the guard interval is 3.2 mus, the occupied bandwidth is 4.15MHz, the modulation method is 64QAM, and the highest transmission rate is 13.5 Mbps.
The first example is an explanation of an embodiment example of a second transmitting and receiving composite wave composed of type 2 wavelets:
1. the system gives the conditions:
system sampling rate fs120 MHz; if the order of the equalization filter is 20, the total expansion point number of the waveform is 20 points, the corresponding time value of the total expansion tailing is 0.167us, the expansion point numbers at two ends are respectively 10 points, and the corresponding time values are respectively 0.083 us;
2. arrangement of wavelet intervals:
let the symbol period of OFDM be TO=TO1+ CP is (16+3.2) ═ 19.2us, CP is cyclic prefix (or time isolation zone), the wavelet interval of TFP-TDMC:
Td=Tθ1+Tθ2+Tθ3=TO=(16+3.2)=19.2us
Tθ1=Tθ11+Tθ12+Tθ13=TO/10=1.92us
Tθ11+Tθ13=0.167us,Tθ12=1.92-0.167=1.75us
wavelet lifetime Td=TO,Tθ1、Tθ2、Tθ3The front information wave, the main wave and the back information wave are respectively intervals, the information wave interval is 1/10 of the OFDM symbol interval, Tθ1=Tθ11+Tθ12+Tθ13=TO1.92us (/ 10 ═ Tθ11+Tθ130.167us is the extension interval),
3. the relationship between the coherent demodulation signal-to-noise ratio and the demodulation interval,
according to the literature, the relation between the signal-to-noise ratio of coherent demodulation and the demodulation interval and amplitude under the condition of narrow-band white noise can be formulated
Figure BDA0001497145760000221
Meaning that the signal-to-noise ratio SNR is proportional to the product of the interval and the amplitude. Thus, a wavelet (i.e., a sine or cosine wave) of OFDM is related to the coherent demodulation signal-to-noise ratio of an information wave of TFP-TDMC as follows:
Figure BDA0001497145760000222
there is the ratio of the coherent demodulation signal-to-noise ratio of OFDM and TPF-TDMC single wavelet
γ=SNRTPF/SNRO=Tθ12α2/TO1=0.109α2 (7‐2)
Description of the drawings: t isO1、Tθ12Respectively, OFDM and TPF-TDMC coherent demodulation interval, aOAnd aTFPRespectively OFDM wavelet and TPF-TDMC information wave amplitude, n0Is the power spectral density, SNR, of the noiseOAnd SNRTFPRespectively, the signal-to-noise ratio of OFDM and TPF-TDMC coherent demodulation; it can be seen that increasing T is required to improve the signal-to-noise ratio of TPF-TDMC coherent demodulationθ12α2Or, when T is reducedθ12The signal-to-noise ratio can be kept unchanged by increasing alpha.
4. Effect of alpha, beta on TFP-TDMC wavelet power and relation to OFDM wavelet power
Since the information wave and the main wave are both sine and cosine waves, and the amplitude of the wave is normalized to 1, the power of one wavelet is
Figure BDA0001497145760000223
Wherein, Tθ1、Tθ2、Tθ3、TdThe intervals of the pre-information wave, the main wave, the post-information wave and the wavelet lifetime are respectively. The amplitude normalization value makes the OFDM wavelet power be
Figure BDA0001497145760000224
TO=TdIs the symbol period and is also the interval occupied by the wavelet. Taking the sum of normalized values in amplitudeO=TdUnder the condition of (1), 1/10 indicating that the information wave interval is the OFDM symbol interval and the power of the two is equal, there are
Figure BDA0001497145760000225
Let Tθ1=Tθ3=TO/10,Tθ2=Td-2Tθ1=(4/5)TOThen there is
0.2α2+0.8β2=1 (7-4)
Since the main wave does not carry information and a high signal-to-noise ratio of the information wave is desired, the value of β is reduced while the value of α is increased. If γ of the formula (7-2) is 1, α is guaranteed29.173, it can ensure that the information wave carries 3 bits as one sub-channel of OFDM; however, according to the formula (7-4), β2(1-1.83)/0.8, which results in β being an imaginary number; for this purpose, 0.2 α is used2<1,α2<5, e.g. by taking alpha24.9, then beta2This corresponds to an amplitude reduction of about one time by 0.025, which reduces the number of carrying bits by 1, so that each information wave carries 2 bits and an information wave pair carries 4 bits.
4. Forming a second received wave
According to the above analysis, the following parameters were used:
(1) each section of the wavelet
Td=Tθ1+Tθ2+Tθ3=TO(16+3.2)=19.2us
Tθ1=Tθ11+Tθ12+Tθ13=TO/10=1.92us
Tθ11+Tθ13=0.167us,Tθ12=1.92-0.167=1.75us
Wherein, the main information wave interval is 1.6us, and the two ends of the main information wave interval are respectively added with a 0.083us interval and a 0.083us expansion interval for storing the copy waveform of sine and cosine, so that the information wave interval is Tθri1=Tθri32.0332us, main wave interval Tθri2=12.766us
(2) Amplitude parameter: alpha is alpha24.9, β 0.05, information wave quantization width 22The number of the grades is 4,
(3) construction of symbols
i. Number of wavelets N in a frequency bandw=Td/Tθ1-1-10-1-9; according to equation (7-3), the power of one wavelet of TFP-TDMC is equal to the power of one wavelet of OFDM (1/2 corresponding to the power of one subchannel), and thus the power of one band is equal to the power of 5 subchannels of OFDM;
number of bands: the condition is that the power of the TFP-TDMC and OFDM wavelet pair is equal; the condition that the OFDM comprises 52 sub-channels is regulated according to the design and the standard of the interval and the values of alpha and beta, corresponding to TFP-TDMC, the frequency range number is 52/5 and 10.4, and 10 frequency ranges are taken;
frequency arrangement of the center of each band:
the first frequency band has a central angular frequency of
Figure BDA0001497145760000241
The angular frequencies of the other frequency bands being 2 omega respectively1,3ω1,…,10ω1
The configuration of the waveform: rewriting outputWave receiving formula
Figure BDA0001497145760000242
Wherein,
Figure BDA0001497145760000243
5. demodulation (omit training and equalization)
The first step is as follows: subtracting the dominant wave component in the received wave:
Figure BDA0001497145760000244
the second step is that: performing fast Fourier operation
Figure BDA0001497145760000245
K front or back information wave amplitudes can be obtained
Figure BDA0001497145760000246
Superscripts f and a indicate the front or back information wave, respectively; order to
Figure BDA0001497145760000247
Repeating the first step and the second step to obtain all information wave amplitudes;
6. analysis of transmission rate and spectral efficiency
(1) TFP-TDMC: since the information wave pairs are densely populated in the composite wave, the transmission rate of one information wave pair interval is the transmission rate of the entire symbol. According to the analysis of wavelet power and signal-to-noise ratio of information wave, one information wave pair carries 4 bits, 40 bits in total are carried in 10 frequency bands, and the transmission rate R is (10 x 4 bits)/Tθ1=40/1.92us=20.83Mbps;
(2) OFDM: transmission rate of 13.5Mbps
(3) And (3) comparison: the transmission rate of the second received wave of the TFP-TDMC is 1.5 times that of OFDM.
(4) Analysis of spectral efficiency: OFDM occupies 52 sub-channels, while TFP-TDMC occupies only 10 sub-channels, so the spectral efficiency of TFP-TDMC is 7.8 times that of OFDM.
The second example is an implementation case of a third sending and receiving synthetic wave for the wavelet type 1 configuration:
the system setting conditions are the same as the first example described above.
1. Composition of synthetic wave
(1) Wavelet region:
Figure BDA0001497145760000251
(2) each interval in the composite wave
Time delay interval of information wave in group
Figure BDA0001497145760000252
Meaning that a group contains 100 wavelets, the inter-group delay interval Δ τu=2×T θri12 × 4 ═ 8us, the number of groups in one composite wave is 2;
(3) selection of values of α and β: for a single wavelet, the smaller beta is, the smaller the wavelet power is, but for different synthesized waves, because each wavelet has time delay, the larger alpha can be obtained under the condition of ensuring that the power of the synthesized wave is unchanged due to the fact that the values of alpha and beta and the positive and negative changes of the alpha and beta are different, so that the improvement of the signal-to-noise ratio of the information wave is facilitated, and the optimization in a certain degree can be obtained through computer simulation; to reduce space, the following gives a case: α -16 and β -12, the composite wave power is 51.854, the power of OFDM is 52;
then there are
Figure BDA0001497145760000253
Corresponding to 18 decibels of the audio signal,
in digital communication, for each 6 db improvement of the noise ratio, one bit is added, that is, the information wave of the TFP-TDMC can carry 2 more bits than one sine wave of the OFDM, and carries 5 bits in total, and then the transmission rate of 100 information waves is (5 × 100)/8 × 62.5Mbps, which is 4.6 times of the OFDM standard.
(4) With the interval, the values of alpha and beta and the number of bits carried by the information wave, a composite wave can be created according to the formulas (2-4-1) and (2-4-2);
2. analysis of transmission rate and spectral efficiency
The transmission rate of the third waveform of TFP-TDMC is (5 × 100)/8 ═ 62.5Mbps, which is 4.6 times that of the OFDM standard.
Spectral efficiency: if all wavelets are at the same frequency, the spectral zero bandwidth is 0.042MHz, the spectral efficiency is 62.5/0.042 ═ 1.5Mbps/Hz, and if 52 wavelets are at the same frequency, the spectral efficiency is 62.5/(0.042 × 52) ═ 28.8bps/Hz
The above description is only for the preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art should be considered to be within the technical scope of the present invention, and the technical solutions and the inventive concepts thereof according to the present invention should be equivalent or changed within the scope of the present invention.

Claims (3)

1. A time-delay multi-carrier modulation and demodulation method based on time-frequency pulse shaping is characterized in that:
at a sending end, binary information is converted into corresponding quantization amplitude and then sequentially given to a front information wave and a back information wave of a base wavelet to form corresponding wavelets, H wavelets are generated, all wavelet starting points are required to be placed in the survival period of the first wavelet, all wavelets are sequentially moved by different time delays from the starting point of a symbol waveform and then linearly superposed to form a synthetic wave, namely, a periodic symbol waveform with T as a period is formed to complete a modulation process, namely, the method of sequentially moving by one time delay and then linearly superposed is called as a time delay rule, different modulation methods have different time delay rules, and the difference is mainly expressed as different time delay intervals of adjacent wavelets; the fundamental wave is composed of waveforms corresponding to respective corresponding sections of a front information wave, a main wave and a rear information wave which are distributed in sequence, and the waveform structures of the front information wave and the rear information wave are consistent; the wavelets are waveform structures formed by respectively changing the amplitudes of the front information wave and the rear information wave of the fundamental wave so that the front information wave and the rear information wave respectively carry different binary information and the amplitude of the main wave is kept unchanged; front information wave of the fundamental waveThe sections corresponding to the main wave and the post information wave are respectively Tθ1、Tθ2And Tθ3Expressing that the three intervals are combined into a fundamental wavelet or a wavelet interval in turn, namely the life cycle of the fundamental wavelet/wavelet, and is marked as Td,Td=Tθ1+Tθ2+Tθ3All intervals in the formula are half-open and half-closed intervals, the waveform corresponding to each base wavelet only exists in the life cycle corresponding to each base wavelet, and the base wavelets outside the life cycle are not defined or take the value of 0; the tokens associated with different fundamental wavelets are distinguished by the subscript i ═ 1, …, H, then TdiRepresents the lifetime of the ith base wavelet, such that Tdi=Td(i+1)Or Tdi≠Td(i+1)Taking Tdi=Td(i+1)=TdAnd T isθi1、Tθi2、Tθi3Respectively representing the sections corresponding to the ith front information wave, the main wave and the back information wave; the symbol "\\" means "or";
more specifically, the fundamental waves are divided into two types, i.e., type 1 fundamental wave and type 2 fundamental wave, which are formulated as
Figure FDA0002586387160000011
The corresponding wavelets are classified into type 1 wavelets and type 2 wavelets, and are formulated as
Figure FDA0002586387160000012
Wherein the superscript 1\2 is used to distinguish type 1 basewaves or type 2 basewaves, where ω iskRepresenting the angular frequency of the main wave of the wavelet and the fundamental wave, so that the wavelets with the same angular frequency form a sub-channel on the frequency domain, and the center frequency of each sub-channel is formed by the angular frequency omega of the main wavekDetermining that K is 1, …, K, i.e. K subchannels are formed in total; and each sub-channel is not required to be orthogonal in frequency domain, so that the same-frequency wavelets contained in each sub-channel are equal or unequal, that is, if the wavelets with the same frequency contained in each sub-channel are denoted by LkK1, …, K, allowing Lk=Lk+1Or Lk≠Lk+1Angular frequency omegakThe subscript k and the wavelet subscript i have the following relation: i.e. the wavelet index i of each k-th sub-channel is Lk-1+1,…,Lk-1+Lk
The fundamental wave can be uniformly expressed by formula (1)
Figure FDA0002586387160000021
In the formula, the vertical line symbol "|" represents
Figure FDA0002586387160000022
And
Figure FDA0002586387160000023
three waveforms are connected in sequence, the interval corresponding to each waveform is a half-open half-close interval, subscripts 1, 2 and 3 of the three waveforms give the connection sequence on a time coordinate, namely the three waveforms are expressed by a formula
Figure FDA0002586387160000024
Representation of affiliation of time variables
Figure FDA0002586387160000025
So that the front and rear information waves are distinguished by subscripts f, a, respectively, i.e.
Figure FDA0002586387160000026
In the form of a front information wave,
Figure FDA0002586387160000027
in order to be the rear information wave,
Figure FDA0002586387160000028
for the reference dominant wave, i denotes the subscript, τ, of the different fundamental wavesiThe delay of the ith base wavelet relative to a symbol starting point is used for distinguishing various waves of a transmitting end and a receiving end by subscripts t and rThe shape and the time parameters of the shape and the time,
Figure FDA0002586387160000029
the interval representing the front information wave and the rear information wave is equal, and the width of the interval is taken according to the following principle: is provided with
Figure FDA00025863871600000210
c≤1,
Figure FDA00025863871600000211
Determining through computer simulation; delta taui=τi+1iRepresenting the time delay interval between the (i + 1) th wavelet and the (i) th wavelet, and defining taui+1i,τ10 and τi∈Td1That is, the starting point of the 1 st wavelet is coincident with the starting point of the symbol, and the starting points of other wavelets are all in the life cycle of the 1 st wavelet, allowing Tdi≠Tdj,Δτi≠ΔτjOr Tdi=Tdj,Δτi=Δτj,i&j-1, …, H, notation&To indicate the sum, when taken equal, the subscripts are omitted and when T is 2TdOr T ═ 2Td- Δ τ, taking Tdi=Tdj,Δτi=Δτj,T=2Td
The performance of the base wavelet at the sending end and the receiving end is different, specifically:
at the transmitting end, the information wave and the main wave of the type 1 wavelet are embodied as waveforms of the following formulas (1-1) and (1-2): namely, it is
Figure FDA00025863871600000212
Figure FDA00025863871600000213
Figure FDA00025863871600000214
Figure FDA00025863871600000215
z(t(1\3)2i)=z(t22i) 0 represents Tθti(1\3)2And Tθti22The value of the waveform in the receiving end is 0, and a reserved interval is expanded for the receiving end;
the information wave and the main wave both include two intervals, interval Tθti(1\3)1And Tθti(1\3)2The interval T of the information wave and the main waveθti(1\3)2And Tθti22The subscript (1\3) is used for distinguishing two kinds of interval abbreviation forms of front or back information waves;
the type 2 fundamental wave information wave and the dominant wave are represented by waveforms of the following formulae (1-3) and (1-4):
Figure FDA0002586387160000031
Figure FDA0002586387160000032
Figure FDA0002586387160000033
scale formula (1-3-1) middle interval Tθti(1\3)12The waveform of (a) is a main information wave,
Figure FDA0002586387160000034
is sin omegak(t(1\3)12i) A replica of the waveform to which the trailing portion corresponds,
Figure FDA0002586387160000035
is sin omegak(t(1\3)12i) A replica of the waveform to which the front end portion corresponds,
Figure FDA0002586387160000036
is cos omegak(t(1\3)12i) A replica of the waveform to which the trailing portion corresponds,
Figure FDA0002586387160000037
is cos omegak(t(1\3)12i) A replica of the waveform to which the front portion corresponds, where Tθti(1\3)1=Tθti(1\3)11+Tθti(1\3)12+Tθti(1\3)13,Tθti(1\3)11=Tθti(1\3)13=TP/2;
Figure FDA0002586387160000038
z(t(1\3)2i)=z(t22i) 0 represents Tθti(1\3)2And Tθti22The value of the waveform in the equalization filter is 0, the reserved interval is expanded for a receiving end, and TP is the order of the equalization filter;
ω in the above formulakThe wavelets with the same angular frequency form a sub-channel on a frequency domain, each sub-channel is orthogonal or non-orthogonal on the frequency domain, so that the same frequency sub-wave number contained in each sub-channel is equal or not, wherein the central frequency of the kth sub-channel is formed by the angular frequency omegakDeciding that K subchannels are formed in total;
at the receiving end, the received base wavelet is represented as
Figure FDA0002586387160000039
Wherein the information wave and the main wave of the type 1 fundamental wave are expressed as the following formulas (1-5-1) and (1-5-2), respectively:
Figure FDA0002586387160000041
Figure FDA0002586387160000042
wherein, the information wave and the main wave of the type 2 fundamental wave are expressed as the formulas (1-5-3) and (1-5-4) respectively
Figure FDA0002586387160000043
Figure FDA0002586387160000044
Wherein,
Figure FDA0002586387160000045
scale formula (1-5-5) middle interval Tθri(1\3)12The waveform of (a) is a main information wave,
Figure FDA0002586387160000046
is that
Figure FDA0002586387160000047
The spread wave of (a) is generated,
Figure FDA0002586387160000048
is that
Figure FDA0002586387160000049
The spread wave of (a) is generated,
Figure FDA00025863871600000410
is that
Figure FDA00025863871600000411
The spread wave of (a) is generated,
Figure FDA00025863871600000412
is that
Figure FDA00025863871600000413
The spreading wave of (1).
2. The time-frequency pulse shaping-based time-delay multi-carrier modulation and demodulation method according to claim 1, characterized in that:
in the modulation process of the transmitting end, firstly, binary information is converted into corresponding quantization amplitudes according to a general method of digital communication, each quantization amplitude is sequentially given to a front pulse information wave and a rear pulse information wave of each base wavelet, so as to form 1 type wavelets or 2 type wavelets of the transmitting end as shown in the following formula (2),
Figure FDA00025863871600000414
in the formula (2), the reaction mixture is,
Figure FDA00025863871600000415
the information pulse amplitude quantization method comprises the steps that quantization amplitudes corresponding to front or rear information waves are respectively, Y belongs to an integer field, and Z is a quantization amplitude level number, namely the amplitude of front and rear information pulses can select one quantization amplitude from Y levels; alpha and beta are amplitude adjustment parameters, and alpha is taken>1, the value of beta only needs to ensure that the dominant wave amplitude does not exceed the allowable value of the system, the specific values of alpha and beta are relatively unchanged, and the values are changed only when the working state of the system needs to be adjusted, and are generally determined by computer simulation;
specifically, the modulation process at the transmitting end adopts any one of the following three methods:
the method comprises the following steps: modulation method based on type 1 wavelet direct addition
Specifically, the 1-type wavelets are superimposed according to a time-delay rule to form a synthetic wave, as shown in the following formula
Figure FDA0002586387160000051
τi∈Td110 means that all the H wavelets have their origins in the lifetime of the first wavelet and the time delay interval is equal to the information wave interval, i.e. Δ τ Tθti1=Tθti3,H=(Tdi-Tθti3)/Tθti1(ii) a The formula (2-1) is called a first transmission wave;
Figure FDA0002586387160000052
the expression of (2) is shown in the formula, and the corresponding basic wavelets are type 1 basic wavelets;
the method 2 comprises the following steps: modulation method based on type 2 wavelet pair
The method is specifically divided into three steps:
the first step is as follows: constructing same-frequency sine-cosine wavelet pairs on the basis of formulas (1-3), (1-4), (1-3-1) and (1-4-1)
Figure FDA0002586387160000053
Figure FDA0002586387160000054
Figure FDA0002586387160000055
Scale formula (2-3-1) middle interval Tθti(1\3)12The inner waveform is the master information wave pair,
Figure FDA0002586387160000056
is sin omegak(t(1\3)12i) A replica of the waveform of the back end portion,
Figure FDA0002586387160000057
is sin omegak(t(1\3)12i) A replica of the front-end portion waveform,
Figure FDA0002586387160000058
is cos omegak(t(1\3)12i) A replica of the waveform of the back end portion,
Figure FDA0002586387160000059
is cos omegak(t(1\3)12i) A replica of the front-end portion waveform; wherein, in the interval Tθti(1\3)12The main information wave pair in the system is a same-frequency sine and cosine orthogonal wavelet pair;
the second step is that: constructing a type 2 sub-synthetic wave consisting of the same-frequency type 2 sub-wave pairs according to a time delay rule:
Figure FDA00025863871600000510
Figure FDA00025863871600000511
representing the lifetime of the first wavelet in the sub-composite wave, with a delay interval equal to the information wave pair interval, i.e.
Figure FDA00025863871600000512
The sub-synthetic waves all comprise wavelet angular frequencies omegak(ii) a Let K be 1, …, K then obtain K sub-composite waves of the formula (2-3-3), which defines that all pairs of sub-information waves with the same time delay in different sub-composite waves are orthogonal in the frequency domain, i.e., for K, j be 1, …, K,
Figure FDA0002586387160000061
and
Figure FDA0002586387160000062
the pair of master information waves contained in (1) is orthogonal;
the third step: composition type 2 synthetic wave:
Figure FDA0002586387160000063
the formula (2-3-4) is called as a second transmitting wave;
the method 3 comprises the following steps: modulation method of packet delay double overlapping based on 1-type wavelet
The method is to form a synthetic wave from wavelets of type 1, and specifically, to divide all wavelets into groups of U1, …, U, where each group of wavelets is
Figure FDA0002586387160000064
The starting points of the wavelets in the u-th group are sequentially shifted by a time delay in the first previous information wave interval of the group
Figure FDA0002586387160000065
Figure FDA0002586387160000066
Balance
Figure FDA0002586387160000067
The time delay in the group of the ith wavelet in the u group is called
Figure FDA0002586387160000068
Is an inter-group delay interval, thereby
Figure FDA0002586387160000069
The origin of the wavelets of type 1 being in the front information wave interval T of the first wavelet of the grouptθi1Therein is thus
Figure FDA00025863871600000610
The information wave of the 1-type wavelet occupies an interval of
Figure FDA00025863871600000611
The wave is called an information wave band, and meanwhile, the main wave and the post information wave are overlapped and added according to the same time delay rule, so that a u-th combined wave is formed, which is called a sub-combined wave and is the first wave; provision for
Figure FDA00025863871600000612
The sub-combinations are sequentially separated for the time delay interval between groups
Figure FDA00025863871600000613
Moving a time delay and adding them to form a composite wave of symbols, this being the second term
Figure FDA00025863871600000614
Group delay, and make the information bands of the sub-composite waves not overlap, but make the main band still partially overlap, and the inter-group delay interval
Figure FDA00025863871600000615
The following formula is a sub-composite wave formed by overlapping and adding the time delays of the wavelets in the u group:
Figure FDA00025863871600000616
the corresponding wavelet has the same basic structure as the wavelet in equation (2-1) except that it follows the subscript
Figure FDA00025863871600000617
With variations, subscripts
Figure FDA00025863871600000618
Indicates the u in the group
Figure FDA00025863871600000619
The sign of the correlation of the individual wavelets,
Figure FDA00025863871600000620
the wavenumbers within a group are equal in width for different groups of wavelet lifetimes, but the position of the wavenumbers in different groups varies with the subscript U1, …, U, i.e., U
Figure FDA00025863871600000621
Is at the position of
Figure FDA00025863871600000622
And U is 1, …, varying from U,
Figure FDA00025863871600000623
is the u-th group
Figure FDA00025863871600000624
The time delay of the individual wavelets is,
Figure FDA0002586387160000071
an interval representing a front information wave of the first wavelet in the u-th group,
Figure FDA0002586387160000072
means that all the information wave origins of the wavelets belonging to the u-th group are in the interval
Figure FDA0002586387160000073
In the interior of said container body,
Figure FDA0002586387160000074
will represent subscripts of the front information wavefIs changed into a superscriptfAnd specify
Figure FDA0002586387160000075
Recording the time shift interval of information wave in the u-th group
Figure FDA0002586387160000076
And stipulate
Figure FDA0002586387160000077
Indicating equal spacing of information waves within a group, then
Figure FDA0002586387160000078
The section of the u-th group of the sub-synthetic waves; the symmetry characteristics of the information wave before and after the base wave structure are determined, and the description is simultaneously suitable for the information wave after the base wave structure; the start of the information waves after all wavelets in the set also being in the interval of the information waves after the first wavelet in the set
Figure FDA0002586387160000079
And
Figure FDA00025863871600000710
respectively representing the intervals of the first sub-wavefront information wave and the rear information wave in the group, and the superscripts f and a respectively correspond to the meanings of the front information wave and the rear information wave; each sub-composite wave is divided into three intervals: the first interval is
Figure FDA00025863871600000711
A front information wave time delay overlapping interval, which also contains part of main wave components and is called front information wave overlapping interval and marked as Tu1The second interval is the main wave overlapping region, denoted as Tu2And the third interval is
Figure FDA00025863871600000712
The post-information wave time delay overlapping interval also contains partial main wave component, which is marked as Tu3
Then, the sub-synthesized waveforms of all groups are separated in sequence
Figure FDA00025863871600000713
Moving one delay by one delay and adding them together to form the final composite wave, as shown in the following formula
Figure FDA00025863871600000714
The formula (2-4-2) is called as a third transmitting wave; the periods of the above three transmission waves are respectively T1、T2And T3Represents;
after the transmission wave passes through the channel, the waveform will change, and after proper equalization, the waveform is specifically represented as:
corresponding to three kinds of transmission waves transmitted by the transmitting end, there are three kinds of reception waves after the receiving end is equalized:
the first method comprises the following steps: received waves corresponding to the first type of transmitted waves:
Figure FDA00025863871600000715
the corresponding wavelet:
Figure FDA0002586387160000081
wherein, afi\aaiα is the actual amplitude of the preceding or following information wave;
and the second method comprises the following steps: received wave corresponding to second transmitted wave
Figure FDA0002586387160000082
Figure FDA0002586387160000083
Figure FDA0002586387160000084
Figure FDA0002586387160000085
Figure FDA0002586387160000086
And
Figure FDA0002586387160000087
all are deformation waves with time delay expansion,
Figure FDA0002586387160000088
and
Figure FDA0002586387160000089
is the actual amplitude of the wavelet at the receiving end;
Figure FDA00025863871600000810
Figure FDA00025863871600000811
and the third is that: received wave corresponding to third transmitted wave
Figure FDA00025863871600000812
Wherein, the sub-synthetic wave:
Figure FDA00025863871600000813
Figure FDA00025863871600000814
the corresponding wavelet has the same basic structure as equation (4-2), but the range of the time quantum varies with the change of the subscript.
3. The time-frequency pulse shaping-based time-delay multi-carrier modulation and demodulation method according to claim 2, characterized in that:
there are corresponding respective demodulation methods for the three received waveforms: namely, it is
The method comprises the following steps: for coherent demodulation of the first received wave, specifically,
in the interval Tθri1And Tθri3Respectively performing operations in the following steps:
Figure FDA0002586387160000091
Figure FDA0002586387160000092
wherein,
Figure FDA0002586387160000093
in the formulae (5-1-1) and (5-1-2), the notation x is a vector dot product; wherein
Figure FDA0002586387160000094
The main wave component in the interval is obtained;
repeating the equations (5-1-1) and (5-1-2) for i, j ∈ {1, …, H }, and obtaining the amplitude of the information wave before and after each wavelet corresponding to the symbol waveform, thereby completing demodulation of the first type of received wave;
the method 2 comprises the following steps: for the demodulation of the second type of received wave, the specific steps are as follows:
the first step is as follows: dividing the waveform of formula (4-3-1) into
Figure FDA0002586387160000095
Section, front
Figure FDA0002586387160000096
Each section interval in a section is Trθi1After, after
Figure FDA0002586387160000097
Each section interval in a section is Trθi3Wherein, the first
Figure FDA0002586387160000098
The waveform within a segment may be represented as
Figure FDA0002586387160000099
Figure FDA00025863871600000910
Is the corresponding dominant wave component;
the second step is that: subtracting the main wave component in (5-2-1)
Figure FDA00025863871600000911
Thus, in the formula (5-2-2)
Figure FDA00025863871600000912
Intermediate interval T of segmentrθi(1\3)2Is a combination of K orthogonal pairs of sine-cosine waves,
the third step: fast Fourier operation is performed on the waveform shown in the formula (5-2-2)
Figure FDA00025863871600000913
Repeating the formula (5-2-3) for K equal to 1, … to obtain K front or back information wave amplitudes
Figure FDA00025863871600000914
Superscripts f and a indicate the front or back information wave, respectively; order to
Figure FDA0002586387160000101
Repeating the formulas (5-2-1) to (5-2-3) to obtain all information wave amplitudes;
the method 3 comprises the following steps: for the packet demodulation of the third type of received wave, the following four steps are specifically performed:
the first step is as follows: the received wave represented by the formula (4-4-1)
Figure FDA0002586387160000102
Carrying out segmentation interception to obtain waveforms of 2U groups, namely waveforms of front U groups and rear U groups, wherein the front U groups comprise all front information waves and main waves in corresponding intervals, and the rear U groups correspond to all rear information waves and main waves in corresponding intervals; each group of segment waveforms is expressed as follows
Figure FDA0002586387160000103
Wherein the superscript f \ a represents the front or rear information wave, TruThe interval representing the waveform of the u-th group of segments
Figure FDA0002586387160000104
The information wave interval before and after the wavelet is
Figure FDA0002586387160000105
The second step is that: in the interval
Figure FDA0002586387160000106
All main wave components are removed:
Figure FDA0002586387160000107
Figure FDA0002586387160000108
the action interval representing the time variable is in the front information wave interval of the u group section;
the third step: obtaining a demodulation equation set
In the interval
Figure FDA0002586387160000109
Performing coherent operation as shown in the following formula:
Figure FDA00025863871600001010
order to
Figure FDA00025863871600001011
Repeating the formula (5-3-3) to obtain the column vector related to the coherent operation of the front information wave
Figure FDA00025863871600001012
Further forming a linear system of equations for the front information wave
Figure FDA00025863871600001013
This is called the demodulation equation set for the front information wave, where,
Figure FDA00025863871600001014
for the magnitude column vector to be solved,
Figure FDA00025863871600001015
in order to coherently operate on the column vectors,
Figure FDA00025863871600001016
Figure FDA0002586387160000111
in order to demodulate the matrix, the demodulation matrix,
Figure FDA0002586387160000112
according to the symmetry of the front and rear information waves, the formula (5-3-3) is changed into:
Figure FDA0002586387160000113
order to
Figure FDA0002586387160000114
Repeating the formula (5-3-5) to obtain the column vector related to the post-information wave coherent operation
Figure FDA0002586387160000115
Further form a linear equation system
Figure FDA0002586387160000116
Wherein,
Figure FDA0002586387160000117
for the magnitude column vector to be solved,
Figure FDA0002586387160000118
in order to coherently operate on the column vectors,
Figure FDA0002586387160000119
in order to demodulate the matrix, the demodulation matrix,
Figure FDA00025863871600001110
obtaining a demodulation equation set related to the post information wave; the combined writing of formulae (5-3-4) and (5-3-6) is:
Figure FDA00025863871600001111
the subscripts in the above formulas are each increased by u as a group difference,
Figure FDA00025863871600001112
and ωjIndicating that different information waves may have different angular frequencies, k&j=1,…,K;
The elements of the matrix are obtained by training and are obtained by carrying out coherent operation on the fundamental wave and a part of waveforms corresponding to all fundamental waves including the fundamental wave in the interval; wherein, x is a vector point-by-symbol;
the fourth step: solving the system of equations shown in equation (5-3-7)
Figure FDA00025863871600001113
Obtaining the amplitude vectors of all the front and back information waves of the u-th group;
let U be 1, …, U, repeat the magnitude vector of the preceding and following information waves of all groups.
CN201711278898.9A 2017-12-06 2017-12-06 Time-delay multi-carrier modulation and demodulation method based on time-frequency pulse shaping Active CN108092931B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711278898.9A CN108092931B (en) 2017-12-06 2017-12-06 Time-delay multi-carrier modulation and demodulation method based on time-frequency pulse shaping

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711278898.9A CN108092931B (en) 2017-12-06 2017-12-06 Time-delay multi-carrier modulation and demodulation method based on time-frequency pulse shaping

Publications (2)

Publication Number Publication Date
CN108092931A CN108092931A (en) 2018-05-29
CN108092931B true CN108092931B (en) 2020-11-10

Family

ID=62174549

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711278898.9A Active CN108092931B (en) 2017-12-06 2017-12-06 Time-delay multi-carrier modulation and demodulation method based on time-frequency pulse shaping

Country Status (1)

Country Link
CN (1) CN108092931B (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1943197A (en) * 2004-04-16 2007-04-04 数据流技术公司 Single and multiple sinewave modulation and demodulation techniques, apparatus, and communications systems
CN101662437A (en) * 2008-08-29 2010-03-03 北京良桑通信技术有限责任公司 Time-frequency-phase mixed multicarrier modulation method
CN101316158B (en) * 2007-05-29 2010-11-17 北京良桑通信技术有限责任公司 Additive waveshape pretreatment method in digital communication modulation
US8687734B1 (en) * 2013-01-03 2014-04-01 Harris Corporation Non-causal orthogonal frequency division multiplexing tapered peak suppression
CN104468454A (en) * 2014-12-29 2015-03-25 大连海事大学 Multi-orthogonal frequency division multiplexing modulation and demodulation method
CN104601517A (en) * 2015-02-28 2015-05-06 大连海事大学 Time-delay multi-carriers modulation and demodulation method
CN106961405A (en) * 2016-01-11 2017-07-18 中兴通讯股份有限公司 Data modulation, demodulation method, data transmission method and the node of multicarrier system

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10411843B2 (en) * 2012-06-25 2019-09-10 Cohere Technologies, Inc. Orthogonal time frequency space communication system compatible with OFDM

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1943197A (en) * 2004-04-16 2007-04-04 数据流技术公司 Single and multiple sinewave modulation and demodulation techniques, apparatus, and communications systems
CN101316158B (en) * 2007-05-29 2010-11-17 北京良桑通信技术有限责任公司 Additive waveshape pretreatment method in digital communication modulation
CN101662437A (en) * 2008-08-29 2010-03-03 北京良桑通信技术有限责任公司 Time-frequency-phase mixed multicarrier modulation method
US8687734B1 (en) * 2013-01-03 2014-04-01 Harris Corporation Non-causal orthogonal frequency division multiplexing tapered peak suppression
CN104468454A (en) * 2014-12-29 2015-03-25 大连海事大学 Multi-orthogonal frequency division multiplexing modulation and demodulation method
CN104601517A (en) * 2015-02-28 2015-05-06 大连海事大学 Time-delay multi-carriers modulation and demodulation method
CN106961405A (en) * 2016-01-11 2017-07-18 中兴通讯股份有限公司 Data modulation, demodulation method, data transmission method and the node of multicarrier system

Also Published As

Publication number Publication date
CN108092931A (en) 2018-05-29

Similar Documents

Publication Publication Date Title
Michailow et al. Generalized frequency division multiplexing: Analysis of an alternative multi-carrier technique for next generation cellular systems
JP4959791B2 (en) Orthogonal frequency division multiplexing using subsymbol processing
CN104601517B (en) A kind of time delay multi-carrier modulation demodulation method
CN108283025B (en) Multiple access in orthogonal time-frequency space communication systems
Kattoush et al. The performance of multiwavelets based OFDM system under different channel conditions
CN104519004B (en) The forming method of the preserved sub-carrier position pattern of NGB W systems
CN105049386B (en) A kind of active interference clearance method in UFMC systems
CN112532351A (en) Interleaving transmission method for weighted fractional Fourier transform frequency domain two-component signal
KR20160039219A (en) Method and apparatus for generating a phy header field
CN110430152A (en) Time-frequency compresses multi-carrier transmission method, method of reseptance, transmitter and receiver
CN105874734A (en) Transmission device, reception device, and communication system
CN101867548A (en) Blind frequency tracking algorithm based on multi-carrier of filter bank
CN106789809A (en) A kind of non-orthogonal multi-carrier transmission method
CN110213191A (en) FBMC-OQAM timing and channel estimation training sequence design method
CN107070836A (en) The design method of FBMC receive-transmit systems based on FRM technologies in a kind of 5G systems
CN112787790B (en) Frequency hopping multi-carrier communication method of power communication system
CN106464273A (en) Signal processing method, transmitter and compressive sampling receiver
CN112702298A (en) OFDM signal transmission method for expanding mixed carrier wave precoding
CN102710404A (en) Transmission method for low transmitting power and single carrier-frequency division multiplexing access system
CN101807954A (en) Time domain synchronous frequency division multiple access method for uplink multi-users
CN108092931B (en) Time-delay multi-carrier modulation and demodulation method based on time-frequency pulse shaping
CN101958866B (en) Pilot frequency insertion method and module
CN101662437B (en) Time-frequency-phase mixed multicarrier modulation method
CN106101047B (en) A kind of anonymous communication Waveform generation method based on OFDM
CN111385234B (en) OFDM radar communication integrated waveform processing method, device, terminal and storage medium

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