KR102120021B1 - 수신 방법, 수신 장치, 송신 방법, 송신 장치, 송수신 시스템 - Google Patents

수신 방법, 수신 장치, 송신 방법, 송신 장치, 송수신 시스템 Download PDF

Info

Publication number
KR102120021B1
KR102120021B1 KR1020197009104A KR20197009104A KR102120021B1 KR 102120021 B1 KR102120021 B1 KR 102120021B1 KR 1020197009104 A KR1020197009104 A KR 1020197009104A KR 20197009104 A KR20197009104 A KR 20197009104A KR 102120021 B1 KR102120021 B1 KR 102120021B1
Authority
KR
South Korea
Prior art keywords
shift
equation
frequency
time
signal
Prior art date
Application number
KR1020197009104A
Other languages
English (en)
Other versions
KR20190089850A (ko
Inventor
토루 코다
Original Assignee
라디우스 가부시키가이샤
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 라디우스 가부시키가이샤 filed Critical 라디우스 가부시키가이샤
Publication of KR20190089850A publication Critical patent/KR20190089850A/ko
Application granted granted Critical
Publication of KR102120021B1 publication Critical patent/KR102120021B1/ko

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
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2668Details of algorithms
    • H04L27/2669Details of algorithms characterised by the domain of operation
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/18Phase-modulated carrier systems, i.e. using phase-shift keying
    • H04L27/22Demodulator circuits; Receiver circuits
    • H04L27/233Demodulator circuits; Receiver circuits using non-coherent demodulation
    • H04L27/2334Demodulator circuits; Receiver circuits using non-coherent demodulation using filters
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/69Spread spectrum techniques
    • H04B1/707Spread spectrum techniques using direct sequence modulation
    • H04B1/7073Synchronisation aspects
    • H04B1/7075Synchronisation aspects with code phase acquisition
    • H04B1/70757Synchronisation aspects with code phase acquisition with increased resolution, i.e. higher than half a chip
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/69Spread spectrum techniques
    • H04B1/707Spread spectrum techniques using direct sequence modulation
    • H04B1/709Correlator structure
    • H04B1/7093Matched filter type
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00Arrangements for detecting or preventing errors in the information received
    • H04L1/004Arrangements for detecting or preventing errors in the information received by using forward error control
    • H04L1/0045Arrangements at the receiver end
    • H04L1/0054Maximum-likelihood or sequential decoding, e.g. Viterbi, Fano, ZJ algorithms
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/18Phase-modulated carrier systems, i.e. using phase-shift keying
    • H04L27/20Modulator circuits; Transmitter circuits
    • 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
    • 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/2626Arrangements specific to the transmitter only
    • H04L27/2627Modulators
    • H04L27/2639Modulators using other transforms, e.g. discrete cosine transforms, Orthogonal Time Frequency and Space [OTFS] or hermetic transforms
    • 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/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2657Carrier synchronisation
    • 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/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2662Symbol synchronisation
    • 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/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2668Details of algorithms
    • H04L27/2673Details of algorithms characterised by synchronisation parameters
    • H04L27/2676Blind, i.e. without using known symbols
    • H04L27/2679Decision-aided
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/32Carrier systems characterised by combinations of two or more of the types covered by groups H04L27/02, H04L27/10, H04L27/18 or H04L27/26

Landscapes

  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Artificial Intelligence (AREA)
  • General Physics & Mathematics (AREA)
  • Discrete Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)
  • Image Processing (AREA)
  • Superheterodyne Receivers (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

본 발명의 일 형태에 따른 수신 방법은, 신호를 수신하는 수신 방법으로서, 여차원 2를 갖는 비가환 시프트의 파라미터 공간을 참조해, 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정하는 추정 스텝을 포함하고 있다.

Description

수신 방법, 수신 장치, 송신 방법, 송신 장치, 송수신 시스템
본 발명은 수신 방법, 수신 장치, 송신 방법, 송신 장치 및 송수신 시스템에 관한 것이다.
통신 시스템 기술 및 이와 관련된 기술에 대해 이미 많은 연구나 개발이 이루어지고 있는데, 발명자들은 시간 분할 방식 및 주파수 분할 방식 등을 이용한 송수신 시스템을 연구 및 개발해 왔다(예를 들면, 특허 문헌 1∼6, 비특허 문헌 1∼32를 참조).
특허 문헌 1: 일본 특허공개 2016-189500호 공보 특허 문헌 2: 국제 특허공개 WO2012/153732 A1 공보 특허 문헌 3: 국제 특허공개 WO2014/034664 A1 공보 특허 문헌 4: 국제 특허공개 WO2013/183722 A1 공보 특허 문헌 5: 일본 특허공개 2016-189501호 공보 특허 문헌 6: 일본 특허공개 2016-189502호 공보 특허 문헌 7: 일본 특허공개 2013-251902호 공보 특허 문헌 8: 일본 특허공개 2012-170083호 공보
비특허 문헌 1: D. Gabor, 「Theory of Communication」, in Proc. Inst. Elect. Engr., pt. Ⅲ, 93, 429-41, 1946 비특허 문헌 2: D. Gabor, 「Lectures on communication theory」, in Technical Report., No., 238, Research Laboratoty of Electronics, Massachusettes Institute of technology Cambribge, 1-48, 1952. 비특허 문헌 3: L. W. Couch, Ⅱ, 「Digital and Analog Communication Systems」, 8th Ed., Pearson, 2013. 비특허 문헌 4: G. B. Folland, Harmonic Analysis in Phase Space, Annals of Mathematics Studies Princeton Univ. Press, 1989. 비특허 문헌 5: Howe, R. On the role of the Heisenberg groups in harmonic analysis, Bulltein of the American Mathematical Society, 3-2, 821-840(1980) 비특허 문헌 6: Rosenberg, J. A selective history of the Stone-von Neumann theorem, Contemporary Mathematics, 365. American Mathematical Society(2004) 비특허 문헌 7: I. Daubechies, Time-frequency localization operators: a geometric phase space approach, IEEE Trans. Information Theory, 34-4, 605-612, 1988. 비특허 문헌 8: I. Daubechies, The wavelet transform, time-frequency localization and signal analysis, IEEE Trans. Information Theory, 36-5, 961-1005, 1990. 비특허 문헌 9: B. Le. Floch, M. Alard & C. Berrou, 「coded OFDM: Coded Orthgonal Frequency Division Multiplex」, Proc. of IEEE, 83-6, 982-996, 1995. 비특허 문헌 10: P. Siohan, C. Sciclet & N. Lacaille, 「OFDM/OQAM: Analysis and Design of OFDM/ OQAM Systems, based on Filterbank Theory」, IEEE, Trans. Sig. 50-5, 1170-1183, May, 2002. 비특허 문헌 11: Sciclet, C., Siohan, P. & Pinchon, D. Perfect reconstruction conditions and design of oversampled DFT-based transmultiplexers, EURASIP J. on Applied Signal Processing, 2006, Article ID 15756, 1-14, 2006. 비특허 문헌 12: B. F. Boroujeny, 「OFDM Versus Filter Bank Multicarrier, IEEE Signal Processing Magazine, 28-3, 92-112, 2011. 비특허 문헌 13: P. P. Vaidyanathan, 「Multirate Systems and Filter Banks」, Prentice-Hall, 1993. 비특허 문헌 14: J. Ville, 「Th'eorie et application de la notion de signal analytique」, C'ables et transmission, no. 2, pp. 61-74. 1948. (J. Ville:「Theory and Applications of the notion of complex signal」, translated from the French by I. Stein, T-92, 8-1-58, U.S.Air Force Project Rand, 1958) 비특허 문헌 15: P. M. Woodward, Probability and Information Theory, with Applications to Radar, Pergamon Press, New York, 1953. 비특허 문헌 16: C. H. Wilcox, 「The synthesis problem for radar ambiguity function」, MRC Technical Report, No. 157, pp. 1-46. Mathematics Reaearch Center, U.S.Army, Univ. Wisconsin, Madison, 1960. 비특허 문헌 17: L. Auslander and R. Tolimieri, 「Radar Ambiguity Functions and Group Theory, SIAM J. Math. Anal., 16-3, 577-601, 1985. 비특허 문헌 18: C. W. Helstrom, Elements of Signal Detection and Estimation, PTR Prentice-Hall, 1995. 비특허 문헌 19: N. Levanson and E. Mozeson, 「Radar Signals」, Wiley interscience, 2004 비특허 문헌 20: Sakurai, J. J. Modern quantum mechanics, S. F. Tuan editor, Rev. ed., Addison-Wesley Pub. Comp. 1994. 비특허 문헌 21: J. von Neumann, The Geometry of Operators, vol. Ⅱ(Ann. Math. Studies, no. 22), 1950. 비특허 문헌 22: Youla, D. C Generalized image restoration by the method of alternating orthogonal projections, IEEE Trans. Circuits and Systems, CAS-25-9, 694-702, 1978. 비특허 문헌 23: Stark, H., Cahana, D. & Webb, H. Restoration of arbitrary finiteenergy optical objects from limited spatial and spectral information, J. Opt. Soc. Amer., 71-6, 635-642, 1981. 비특허 문헌 24: Kohda, T., Jitsumatsu, Y. & Aihara, K. Separability of time-frequency synchronization, Proc. Int. Radar Symp., 964-969, 2013. 비특허 문헌 25: T. Kohda, Y. Jitsumatsu and K. Aihara, 「Gabor division/spread spectrum system is separable in time and frequency synchronization」, Proc. VTC 2013 Fall, 1-5, 2013. 비특허 문헌 26: Y. Jitsumatsu, T. Kohda and K. Aihara, 「Spread Spectrum-based Cooperative and individual time-frequency synchronization」, Proc. (ISWCS), 1-5, 2013. 비특허 문헌 27: Jitsumatsu, Y., Kohda, T. & Aihara, K. Delay-Doppler space divisionbased multiple-access solves multiple-target detection, Jonnsson, M., et al, (eds.) MACOM2013, LNCS8310, Springer, 39-53, 2013. 비특허 문헌 28: T. Kohda, Y. Jitsumatsu and K. Aihara, 「Recovering noncoherent MPSK signal with unknown delay and Doppler using its ambiguity function」, 4th International workshop on recent Advanced in BroadbandAccess Network, (RABAN2013), 251-256, 2013. 비특허 문헌 29: T. Kohda, Y. Jitsumatsu and K. Aihara, 「Phase-tuned layers with multiple 2D SS codes realize 16PSK communication」, 2014 2014 IEEE Wireless Commun. Networking Conference, WCNC 2014, 469-474 (2014). 비특허 문헌 30: Jitsumatsu, Y. & Kohda, T. Digital phase updating loop and delay-Doppler space division multiplexing for higher order MPSK, Jonnsson, M., et al, (eds.) MACOM2014, LNCS8715, Springer, 1-15, 2014. 비특허 문헌 31: T. Kohda, Y. Jitsumatsu and K. Aihara, 「Frequency-division spread spectrum makes frequency synchronisation easy」, Proc. IEEE Globecom 2012, 3952-3958, 2012. 비특허 문헌 32: T. Kohda, Y. Jitsumatsu and K. Aihara, 「Frequency synchronisation using SS technique」, Proc. The ninth Int. Sympo. on Wireless Communication Systems, 855-859, 2012. 비특허 문헌 33: J. F. Daughman, 「Two-dimensional analysis of cortical receptive field profiles」, Vision Research, 20, 846-856, 1980. 비특허 문헌 34: J. F. Daughman, 「Image analysis by local 2-D spectral signatures」, J. Opt. Soc. Amer. (A), 2, p. P74, 1985. 비특허 문헌 35: J. F. Daughman, 「Complete Discrete 2-D Gabor Transform by Neural Networks for Image Analysis and Compression」, IEEE Trans. Acoustics, Speech and Signal Processing, 36-7, 1169-1179, 1988. 비특허 문헌 36: Movella, Javier R., 「Tutorial on 가보 필터s」, Archived from on 2009-04-19, Retrieved 2008-05-14. 비특허 문헌 37: Hans G. Feichtinger and Thomas Strohmer: Gabor Analysis and Algorithms, Birkhauser, 1998. 비특허 문헌 38: Jones, J. P. and Palmer, L. A., 「An evaluation of the two-dimensional 가보 필터 model of simple receptive fields in cat striate cortex」, J. Neurophysiol. 58(6): 1233-1258. 1987 비특허 문헌 39: Tai Sing Lee, 「Image representation using 2d Gabor wavelets」, IEEE Trans. on pattern analysis and machine intelligence, 18-10, 1-13, 1996. 비특허 문헌 40: W. D. Montgomery, 「Optical applications of von Neumann's alternating projection theorem」, Optics Letters, 7-1, 1-3, 1982. 비특허 문헌 41: W. D. Montgomery, 「Restoration of images processing a finite Fourier series」, Optics Letters, 7-2, 54-56, 1982. 비특허 문헌 42: A. V. Oppenheim and J. S. Lim, 「Importance of phase in Signals」, 1980, Proc of the IEEE, 96-5, 529-541, 1981
본 발명자는 통신 기술의 이론적 측면에 관한 연구를 진행시켜, 시간 주파수 위상 공간(Time Frequency Plane, 단순히 'TFP'라고도 기재한다)상에서의 시프트 연산자의 본질적인 비가환성(非可換性)에 착안해야 한다는 지견에 이르렀다.
본 발명의 일 형태는 상기 지견에 기초해 이루어진 것으로, 고효율의 수신 방법, 수신 장치, 송신 방법, 송신 장치 및 송수신 시스템을 실현하는 것을 목적으로 한다.
상기 과제를 해결하기 위해, 본 발명의 일 형태는, 신호를 수신하는 수신 방법으로서, 여차원 2를 갖는 비가환 시프트의 파라미터 공간을 참조해, 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정하는 추정 스텝을 포함한다.
또한, 상기 과제를 해결하기 위해, 본 발명의 일 형태는, 신호를 수신하는 수신 장치로서, 여차원 2를 갖는 비가환 시프트의 파라미터 공간을 참조해, 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정하는 추정부를 포함한다.
또한, 상기 과제를 해결하기 위해, 본 발명의 일 형태는, 신호를 송신하는 송신 방법으로서, 여차원 2를 갖는 비가환 시프트의 파라미터 공간을 참조해, 송신 대상 신호의 시간 및 주파수를 시프트하는 시프트 스텝을 포함한다.
또한, 상기 과제를 해결하기 위해, 본 발명의 일 형태는, 신호를 송신하는 송신 장치로서, 여차원 2를 갖는 비가환 시프트의 파라미터 공간을 참조해, 송신 대상 신호의 시간 및 주파수를 시프트하는 시프트부를 포함한다.
또한, 상기 과제를 해결하기 위해, 본 발명의 일 형태는, 화상 신호를 수신하는 수신 방법으로서, 수신한 화상 신호가 나타내는 공간 시프트 및 공간 주파수 시프트를 파라미터 공간을 참조해 추정하는 추정 스텝을 포함하고, 상기 공간 시프트 및 공간 주파수 시프트의 각각은 2 이상의 차원을 갖고 있다.
또한, 상기 과제를 해결하기 위해, 본 발명의 일 형태는, 화상 신호를 송신하는 송신 방법으로서, 파라미터 공간을 참조해 송신 대상의 화상 신호의 공간 및 주파수를 시프트하는 시프트 스텝을 포함하고, 상기 공간의 시프트 및 상기 주파수의 시프트의 각각은 2 이상의 차원을 갖고 있다.
본 발명의 구성에 의하면, 고효율의 수신 방법, 수신 장치, 송신 방법, 송신 장치 및 송수신 시스템을 실현할 수 있다.
도 1은 실시 형태에 따른 TFP에서의 3종의 분할법을 나타낸다. (a)는 시간 분할(Time Division(TD))을 나타내고, (b)는 주파수 분할(Frequency Division(FD))을 나타내고, (c)는 가보 분할(Gabor Division(GD))[비특허 문헌 1]을 나타낸다. (a)에서의 실선은 정보 데이터의 시간폭 T의 구분을 나타내고, 파선은 TD-Phase code(PC)에 의한 세분할을 나타낸다. (b)에서의 굵은 파선은 정보 데이터의 대역폭 F의 구분을 나타내고, 가는 파선은 FD-Phase code(PC)에 의한 세분할을 나타낸다.
도 2는 실시 형태에 따른 시간 주파수 시프트(time-frequency shift)의 비가환성을 모식적으로 나타낸 도면이다. 비가환성은, 시프트 연산자의 곱(積)에 관한 관계식 Tτ,0·T0,ν=e-i2πτνT0,νTτ,0(좌변은 도면에서의 삼각형, 우변은 도면에서의 사각형에 대응)에서 위상 왜곡(Phase Distortion(PD)) e-i2πτν로서 나타난다. 도면에서 '○'는, 대칭 시간 주파수 연산자(Symmetrical Time-Frequency Shift Operator: Symmetrical TFSO)
[수학식 1]
Figure 112019032179095-pct00001

을 나타낸다[비특허 문헌 26].
도 3의 a)는 실시 형태에 따른 TFP상에 배치한 chip level의 가보 함수 및 그에 관련된 정보를 나타낸다. 특히, a-0)은 TFP상에 배치한 chip level의 가보 함수 gmm'(t)와 그 푸리에 변환(FT)인 Gmm'(f)를 나타내고, a-1)은 gmm'(t)를 (FD-Phase code(FD-PC))X'm'로 가중 가산한 TD-template의 실부, 허부를 나타내고, a-2)는 Gmm'(f)를 (TD-Phase Code(TD-PC))Xm으로 가중 가산한 FD-template의 실부, 허부를 나타낸다. 도 3의 b)는 TFP상의 NN'개의 상호 상관 함수(Cross Correlation Function(CCF)) 및 N'개의 열의 합의 TD-CCF값과 N개의 행의 합의 FD-CCF값을 나타낸다. 도 3의 c)는 Alternating Projection Theorem: APT에 기초하는 시간 제한 공간(TL-TD: Time Limited Time Domain), 대역 제한 공간 BL-FD(Band Limited Frequency Domain)상으로의 직교 투영에 의한 추정값
[수학식 2]
Figure 112019032179095-pct00002

의 교호 갱신 과정과 수렴값 td, fD를 나타낸다.
도 4는 실시 형태에 따른 TD-, FD-PCs Xm, X'm' 및 m'-번째의 TD template
[수학식 3]
Figure 112019032179095-pct00003

을 갖는 식 (25), (67)의 TD signature v[k]를 생성하는 Synthesis Filter Bank(SFB)를 나타낸다.
도 5는 TD-, FD-PCs Xm, X'm' 및 m'-번째의 TD template
[수학식 4]
Figure 112019032179095-pct00004

를 갖는 식 (25), (67)의 FD signature
[수학식 5]
Figure 112019032179095-pct00005

를 생성하는 SFB를 나타낸다.
도 6은 실시 형태에 따른 복소수값 데이터
[수학식 6]
Figure 112019032179095-pct00006

을 입력 신호로 하는 식 (27), (71)의 TD-Complex Envelope(TD-CE)
[수학식 7]
Figure 112019032179095-pct00007

을 생성하는 SFB를 나타낸다.
도 7은 실시 형태에 따른 복소수값 데이터
[수학식 8]
Figure 112019032179095-pct00008

을 입력 신호로 하는 식 (27), (71)의 FD-Complex Envelope(FD-CE)
[수학식 9]
Figure 112019032179095-pct00009

를 생성하는 SFB를 나타낸다.
도 8은 실시 형태에 따른 복소수값 데이터
[수학식 10]
Figure 112019032179095-pct00010

복호를 위한 TD 상관기 어레이를 구비한 Analysis Filter Bank(AFB)를 나타낸다.
도 9는 실시 형태에 따른 복소수값 데이터
[수학식 11]
Figure 112019032179095-pct00011

복호를 위한 FD 상관기 어레이를 구비한 AFB를 나타낸다.
도 10의 (a)는 실시 형태에 따른 TD 상관기 어레이를 나타내고, 도 10의 (c)는 FD 상관기 어레이의 Analysis Filter Bank(AFB)를 나타내고, 도 10의 (b)는 von Neumann의 APT에 의한 두 어레이의 최대 우도 추정값의 교호 갱신을 모식적으로 나타낸다.
도 11은 실시 형태에 따른 von Neumann의 APT(Alternating Projection Theorem)를 모식적으로 나타낸 도면이다. 도면에서의 TL-TD(Time Limited Time Domain), BL-FD(Band Limited Frequency Domain)는, 각각 힐베르트 공간(Hilbert space)의 2개의 부분 공간을 나타낸다. 또한, 부분 공간 TL-TD는 LΔt-(또는 Ts)-시간 제한(Time-limited: TL) 공간을 의미하고, 부분 공간 BL-FD는 LΔf-(또는 Fs)-대역 제한(Band-limited: BL) 공간을 의미한다. 도면에서 화살표는 각각의 부분 공간 상으로의 직교 투영을 나타내고, 최대 우도 추정값과 최대 우도값의 최대값을 제공하는 CCF 번호가 검출된다.
도 12는 실시 형태에 따른
[수학식 12]
Figure 112019032179095-pct00012

-PSK 통신 가능하고, 고속·고정도 측거(測距) 가능한 송신기(부호화기)의 구성을 나타내는 블록도이다.
도 13은 실시 형태에 따른
[수학식 13]
Figure 112019032179095-pct00013

-PSK 통신 가능하고, 고속·고정도 측거 가능한 수신 동기기(복호화기)의 구성을 나타내는 블록도이다.
도 14는 실시 형태에 따른 Main channel(MC)의 delay τ-Doppler shift ν 공간상에서 CCF 실부의 크기의 분포예를 모식적으로 나타내는 도면이다.
도 15는 실시 형태에 따른 Artificial Channel(AC)이 Main channel(MC)에 중첩되었을 경우의 τ-ν 공간상에서 CCF 실부의 크기의 분포예를 모식적으로 나타내는 도면이다.
도 16은 실시 형태에 따른 여차원 2의, 비가환성을 갖는 AC 시프트의 파라미터 공간을 이용한 신호의 TFP 분할을 나타내는 도면이다. 신호(시간폭 Ts, 대역폭 Fs)의 시간·주파수 평면(TFP) S의 분할(가보 분할)(S(0), S(1), S(2), S(3)) 평면과 직교하는 축에 AC의 비가환 시프트의 양
[수학식 14]
Figure 112019032179095-pct00014

를 나타내는 눈금, 및 TFP 분할과 대응하는 2차원 PC 부호
[수학식 15]
Figure 112019032179095-pct00015

를 나타낸다.
도 17은 실시 형태에 따른 여차원 2의, 비가환성을 갖는 AC 시프트의 파라미터 공간을 이용한 신호의 TFP 분할을 나타내는 도면이다. AC의 비가환 시프트의 양에 따라 TFP는 각각 AC0의 TFP, AC1의 TFP, AC2의 TFP, AC3의 TFP와 같이 시프트된다.
도 18은 제1 실시 형태에 따른 송수신 시스템의 구성예를 나타내는 블록도이다.
도 19는 제1 실시 형태에 따른 송수신 시스템을 이용한 데이터 송수신 처리의 흐름을 나타내는 플로우차트이다.
도 20은 제1 실시 형태에 따른 송수신 시스템의 구성예를 나타내는 블록도이다.
도 21은 제1 실시 형태에 따른 송수신 시스템을 이용한 데이터 송수신 처리의 흐름을 나타내는 플로우차트이다.
〔제1 실시 형태〕
이하, 본 실시 형태에 따른 송수신 시스템에 대해 도면을 참조해 설명한다. 우선, 본 실시 형태에 따른 송수신 시스템의 이론적 측면 및 구체적 구성예에 대해 설명하고, 그 후, 특허 청구의 범위에 기재된 내용과의 대응 설명을 행한다.
한편, 이하의 설명에서 인용되는 특허 문헌 및 비특허 문헌은 참조에 의해 본원에 원용된다.
또한, 이하의 설명에서 인용되는 특허 문헌 및 비특허 문헌은 어디까지나 본 명세서에 기재되는 내용에 다소의 관련성을 갖는다는 관점에서 인용하는 것에 지나지 않는다. 따라서, 이들 문헌의 인용은 본 명세서에 기재된 발명의 특허성에 영향이 미치지 않는다.
《이론적 측면의 개요: 시간 지연과 주파수 시프트의 2개 연산자의 비가환성을 이용한 통신(부제: 고정도 측거 가능한 송신 신호 설계와 그 수신기 설계)》
시간폭 Ts, 대역폭 Fs 데이터의 고효율 전송을 위해 시간(TD)·주파수 영역(FD)을 분할하는 다중화 통신 방식의 동기는 어려운 문제이다. 발사 전파의 에코 신호로부터 delay td와 Doppler fD를 추정하는 레이더 문제도 현재 미해결인 상태이다. 이들은 양자역학의 위치, 운동량 작용소[비특허 문헌 4]와 마찬가지로, 시간·주파수 시프트 작용소(Time-Frequency Shift Operator: TFSO) 유래의 위상 왜곡(phase distortion: PD)
[수학식 16]
Figure 112019032179095-pct00016
이나
[수학식 17]
Figure 112019032179095-pct00017
에 원인이 있다. 측거 문제는 비가환 TFSO의 PD에 포함되는 2개의 미지 파라미터 추정 문제이므로, Weyl-Heisenberg Group: WHG에 기초하는 신호 검출·파라미터 추정론의 원용 없이는 고정도(高精度)·고속화를 도모할 수 없다. 본 명세서에서 설명하는 내용의 개요는 이하의 5점으로 정리할 수 있다.
(개요 1)
첫째, 송신 신호 설계에서는 TD 신호와 FD 신호를 동등하게 대상으로 하는 것 외에, 시간·주파수 대칭성(time-frequency symmetrical property: TFSP)을 만족하는 대칭 시프트 작용소(symmetrical TFSO)가, 시간 지연과 주파수 시프트의 비가환성에 의해 다중화 신호의 어드레스 정보를 현재화(顯在化)시키는 연산자라는 점이 제시된다(식 (39), (44), (51), (56) 참조).
(개요 2)
둘째, 시간·주파수 시프트된 가우스 칩 함수(가보 함수)를 주기 N, N'의 TD-, FD-phase code: PC를 이용한 변조, 즉 TD-, FD-Binary phase-shift Keying: 2차원 BPSK가, 파라미터 최대 우도 추정의 M값 검정법에 유용하다는 점이 제시된다.
종래 2종의 PC는 '2차원 확산 부호'라고 호칭되어, 용어의 오용도 많고 확대 해석되었다. 그러나 본 명세서에서는, BPSK에는 지금까지 깨닫지 못했던 다음과 같은 득실의 양면의 기능이 있다는 것을 분명히 한다.
변조로 얻어진 TD-, FD-광대역 신호('signature'라고 부른다)의 CE(식 (27))에는 PC에서 발생한 PD가,
N개의 type-3(또는 type-1)의 TD-template CE(식 (29))(또는 식 (49))(각각의 support는 Ts×LΔf, LΔt×Fs),
N'개의 type-4(또는 type-2)의 FD-template CE(식 (33))(또는 식 (54))(각각의 support는 LΔt×Fs, Ts×LΔf)
로서 자동적으로 매장되므로, template 검출의 가설 검정이 TD, FD로 가능해진다(식 (30), 명제 4 및 식 (35), 명제 5).
(개요 3)
셋째, 위상 정보가 유효하게 이용되지 않는 통상의 최대 우도법과는 달리, type-3(또는 type-1)의 TD-CE template, type-4(또는 type-2)의 FD-CE template matching의 4종의 TD-, FD-cross-correlation function: CCF의 완전식은, ambiguity function(AF) 외에, TD-, FD-PC 변조에 수반해 시간폭 Δt, 주파수폭 Δf로 이산화(離散化)한 L=(ΔtΔf)-1점의 twiddle factor
[수학식 18]
Figure 112019032179095-pct00018
을 포함한다. 또한, 상기 완전식에서는, 상기 W의 DFT, IDFT형의 합에서 비가환성 유래의 PD가 W의 멱수로서 등장하고, 그 결과, 상기 완전식은 3 요소의 곱(積)으로서 표현된다(Lemma 2 식 (41), Lemma 4 식 (45) 및 식 (52), (57)). 이 때문에, 본 실시 형태에 기재의 방법에 의하면, DSP에 의한 고속 계산이 보장된다.
(개요 4)
넷째, 특허 문헌 1에서 정의·도입한 Phase-Updating Loop: PUL의 수렴에 관해, Youla[비특허 문헌 22]의 신호 복원법에 기초하는 증명을 제공한다.
먼저, 힐베르트 공간인 신호 공간의 부분 공간으로서, Ts-(또는 LΔt-) 시간 제한 공간(TL-TD) E3(또는 E1), Fs-(또는 LΔf-) 대역 제한 공간(BL-FD) E4(또는 E2)를 도입 및 정의했다. 계속해서, N, N'개의 TD-, FD-CCF 어레이에 기초해, 이들 공간상으로의 직교 투영 연산자(Projection operator: PO), Ts-(또는 LΔt-) 시간 제한 연산자 P3(또는 P1) 및 Fs-(또는 LΔf-) 대역 제한 연산자 P4(또는 P2)의 4개의 PO를 식 (62)∼(64)로 정의했다.
계속해서, Alternating Projection Theorem: APT[비특허 문헌 21]의 교호 투영 연산자(APTO)
[수학식 19]
Figure 112019032179095-pct00019
를 본 발명의 관점에서 정의했다. 다음으로, 추정값
[수학식 20]
Figure 112019032179095-pct00020
의 함수인 통신로의 감쇠 상수 Ae의 MLE의 갱신식 (59)와
[수학식 21]
Figure 112019032179095-pct00021
의 갱신식 (60), (61)을 이용해, APTO의 수렴 영역(TD-CE template의 support와 FD-CE template의 support의 곱집합, 즉 칩 어드레스(ρ, ρ')의 시간폭 LΔt, 대역폭 LΔf의 직사각형 영역) 내에 td, fD가 있는 것을 최대 우도 추정해, 추정값
[수학식 22]
Figure 112019032179095-pct00022
가 LΔt×LΔf의 정도(精度)이고, 계산량은
[수학식 23]
Figure 112019032179095-pct00023
을 대신해
[수학식 24]
Figure 112019032179095-pct00024
가 되는 것을 증명했다.
즉. 이 APTO는, 시간 주파수 공간(time-frequency plane: TFP)의 어느 영역을 추출하고, 다른 부분을 필터링하는(filter-out) 국소 선택 연산자(localization operator: LO)이며, DSP에서 많이 쓰이는 급준한 필터의 대체물이다.
종래의 통신에서는 나이퀴스트(nyquist) 조건을 만족하지 않는다는 이유로 이용되지 않았던 가우스 함수가, 다른 뛰어난 성질에 의해 본질적 역할을 하는 것이 밝혀졌다.
PUL은, 시간 PTs와 대역 P'Fs의 필요 자원과 데이터 레벨의 어드레스
[수학식 25]
Figure 112019032179095-pct00025
를 이용하면, td, fD의 존재 범위를 한정하지 않는 탐색법이다. 이 때문에, TD-, FD-가우스 함수의 2차원 PC 변조 신호의 송신기와 TD-, FD-CCF 어레이에 PUL을 실장한 수신기의 세트는 고정도·고속의 신호 복원이 가능하여, 고정도이면서 고속으로 파라미터 추정이 가능한 통신 시스템을 실현할 수 있다는 점이 제시된다. 환언하면, 상기의 구성을 이용함으로써, 비가환성을 이용한 통신 시스템의 패러다임의 변화를 제시한다.
(개요 5)
다섯째, 동기를 실시하면서, 초고(超高)의
[수학식 26]
Figure 112019032179095-pct00026
-PSK의 부호화·복호화 시스템을 제공한다.
[수학식 27]
Figure 112019032179095-pct00027
-PSK 통신은, 측거와 데이터 통신이 동시에 가능한 차량 탑재 레이더에 응용된다. 고차(高次)의
[수학식 28]
Figure 112019032179095-pct00028
-PSK 신호
[수학식 29]
Figure 112019032179095-pct00029
의 전송은, 각종 위상 잡음의 영향으로 위상
[수학식 30]
Figure 112019032179095-pct00030
의 변별이 곤란하므로, 주파수 자원 유효 이용(
[수학식 31]
Figure 112019032179095-pct00031
전송)의 면에서 중요하지만, 그 실현은 어렵다고 여겨지고 있다.
이를 해결하기 위해, 도 12의 하단 부분(Switch 1-1과 Switch 1-2 사이의 블록)과 같이
1) 「정보 k」를
[수학식 32]
Figure 112019032179095-pct00032
으로 분해한다. 그리고
2) 시간 지연 및 도플러 시프트의 파라미터 평면('target space'라고 부른다)을
[수학식 33]
Figure 112019032179095-pct00033
개로 균등 배반 분할해, 각각에 2차원 PC(TD-PC, FD-PC)를 할당하고(여기에서,
[수학식 34]
Figure 112019032179095-pct00034
이다), 칩 펄스를
[수학식 35]
Figure 112019032179095-pct00035
번째의 2차원 PC로 변조해, 이들을
[수학식 36]
Figure 112019032179095-pct00036
-code 다중화한다. 그리고,
3) 얻어진 signature 신호를 j번째 target space의 부분 평면의 중심점 대응의 시간 지연
[수학식 37]
Figure 112019032179095-pct00037
및 도플러 시프트
[수학식 38]
Figure 112019032179095-pct00038
의 양(
[수학식 39]
Figure 112019032179095-pct00039
를 인공적 통신로: Artificial Channel(AC)의 시프트라고 부른다)만큼 시프트하고,
4) 시프트된 signature로
[수학식 40]
Figure 112019032179095-pct00040
-PSK 신호
[수학식 41]
Figure 112019032179095-pct00041
을 변조한 신호를 전송한다. 그 결과
[수학식 42]
Figure 112019032179095-pct00042
-시프트된 signature는
[수학식 43]
Figure 112019032179095-pct00043
(Main Channel: MC의 시프트라고 부른다)의 전파로에서 이중의 시프트를 받는다. 수신 CE 신호와의 상관 함수(CCF)에 필요한 추정 수신 CE를 이하의 스텝으로 구한다. 도 13의 중단 부분(Switch 2-1에 접속된 블록)과 같이,
1') 추정값
[수학식 44]
Figure 112019032179095-pct00044
[수학식 45]
Figure 112019032179095-pct00045
로 분해하고,
2') 칩 펄스를
[수학식 46]
Figure 112019032179095-pct00046
번째의 2차원 PC로 변조하고,
3') 그 신호를 추정
[수학식 47]
Figure 112019032179095-pct00047
번째 AC의
[수학식 48]
Figure 112019032179095-pct00048
의 양만큼 시프트해,
4') 얻어진 신호로 추정
[수학식 49]
Figure 112019032179095-pct00049
번째의
[수학식 50]
Figure 112019032179095-pct00050
-PSK 신호
[수학식 51]
Figure 112019032179095-pct00051
을 변조한 신호가 추정 수신 CE가 된다.
N, N'개의 어레이형 TD-, FD-CCF의 실부 최대화 변수는 칩 어드레스(ρ', ρ), 데이터 어드레스
[수학식 52]
Figure 112019032179095-pct00052
로 정해지는 CCF 번호 및 추정 부호
[수학식 53]
Figure 112019032179095-pct00053
이다. 우선
[수학식 54]
Figure 112019032179095-pct00054
종의 2차원 PC 번호
[수학식 55]
Figure 112019032179095-pct00055
를 정하고, 데이터 어드레스
[수학식 56]
Figure 112019032179095-pct00056
의, 위상 보상항
[수학식 57]
Figure 112019032179095-pct00057
이 부가된, 2종의 CCF 실부 최대화를 PUL에 기초하여 행해, 도 13의 하단 부분(Switch 2-2에 접속된 블록)과 같이
[수학식 58]
Figure 112019032179095-pct00058
및 변수
[수학식 59]
Figure 112019032179095-pct00059
의 최대 우도 추정에 의한
[수학식 60]
Figure 112019032179095-pct00060
으로 k를 최대 우도 복호한다. 이는, 저차(低次)
[수학식 61]
Figure 112019032179095-pct00061
-ary PSK-based 부호·복호 시스템을
[수학식 62]
Figure 112019032179095-pct00062
종의 비가환 AC-shift-다중화에 의한, 초고(超高)의
[수학식 63]
Figure 112019032179095-pct00063
-ary PSK 시스템이다. 즉,
[수학식 64]
Figure 112019032179095-pct00064
종의 비가환 시프트의 AC 중에서, 「정보 k」로 정해지는 AC에 종속 접속된,
[수학식 65]
Figure 112019032179095-pct00065
의 MC의 출력으로부터 파라미터 추정을 행하는 동기기(측거기) 및 k의 복호기 겸용 시스템이다. 따라서, 비가환 AC-shift-부호화·복호화에 기초하는 다중화 방식은 통신에서의 PD 활용의 패러다임의 변화이다. 한편,
[수학식 66]
Figure 112019032179095-pct00066
-PSK의 복호의 계산 복잡도는 동기·측거의 계산 복잡도(N+N')의
[수학식 67]
Figure 112019032179095-pct00067
배이다.
《이론적 측면의 상세 및 통신 시스템의 구체적 구성예》
이하에서는, 본 실시 형태에 따른 통신 시스템의 이론적 측면의 상세 및 통신 시스템의 구체적 구성예에 대해 설명한다.
<1. 배경>
통신의 중요 과제의 하나는, 주파수 자원을 효과적으로 활용한 통신 방식의 설계이다[비특허 문헌 1]. 시간 영역(TD) 및 주파수 영역(FD)(도 1 참조)을 분할해 시간폭 Ts와 주파수 대역 Fs의 데이터를 다중화하는 통신 방식(도 1 참조)인 Orthogonal Frequency Division Multiplex(OFDM)는, delay td 나 Doppler fD에 의해 직교성이 무너지는 결점을 갖는다.
한편, 도 1은 TFP에서의 3종의 분할법을 나타낸다. (a)는 시간 분할(Time Division(TD))을 나타내고, (b)는 주파수 분할(Frequency Division(FD))을 나타내고, (c)는 가보 분할(Gabor Division(GD))[비특허 문헌 1]을 나타낸다. (a)에서의 실선은 정보 데이터의 시간폭 T의 구분을 나타내고, 파선은 TD-Phase code(PC)에 의한 세분할을 나타낸다. (b)에서의 굵은 파선은 정보 데이터의 대역폭 F의 구분을 나타내고, 가는 파선은 FD-Phase code(PC)에 의한 세분할을 나타낸다.
한편, td, fD를 갖는 통신로를 개재해 통신을 행하기 전에 최초로 실시하는 절차가 동기이다. 그러나, 다중화에 필요한 2개의 시프트 연산인, 시간, 주파수 시프트 연산(time-frequency shift operator: TFSO)에 수반되는 위상 왜곡(Phase distortion: PD)
[수학식 68]
Figure 112019032179095-pct00068
이, 통신로의 PD
[수학식 69]
Figure 112019032179095-pct00069
에 계속해서 발생하므로, 동기가 쉽지 않다. 에코로부터 td, fD를 추정하는 레이더 문제[비특허 문헌 15]에서조차 유효한 해결책은 찾아내지 못했다.
본 발명자는 우선, td, fD에 관해 TFSP(도 2 참조)를 만족하는[비특허 문헌 24, 26] TD 신호와 그 푸리에 변환(FT)의 FD 신호를 TD와 FD의 위상 부호(phase code: PC)로 변조한, TD-, FD-signature를 생성했다.
한편, 도 2는 시간 주파수 시프트(time-frequency shift)의 비가환성을 모식적으로 나타낸 도면이다. 비가환성은, 시프트 연산자의 곱(積)에 관한 관계식 Tτ,0·T0,ν=e-i2πτνT0,νTτ,0(좌변은 도면에서의 삼각형, 우변은 도면에서의 사각형에 대응)에서 위상 왜곡(Phase Distortion(PD)) e-i2πτν로서 나타난다. 도면에서 '○'는, 대칭 시간 주파수 연산자(Symmetrical Time-Frequency Shift Operator: Symmetrical TFSO)
[수학식 70]
Figure 112019032179095-pct00070
을 나타낸다[비특허 문헌 26].
다음으로, signature에는 PC에 의한 PD가 template로서 매장되므로, TD-, FD-template 검출을 위한, ambiguity function(AF)형[비특허 문헌 15]의 TD-, FD-cross-corrrelation function: CCF 어레이(도 3의 a) 참조)를 정의했다.
그리고, CCF 실부를 최대로 하는 파라미터값과 최대값의 최대를 달성하는 CCF 번호를 검출했다. 그 다음, 파라미터값을 교호 갱신(도 3의 b), c) 참조)하는 td, fD의 무정보 추정법이 레이더 문제의 해결책인 것을 분명히 했다(비특허 문헌 27∼비특허 문헌 32 및 특허 문헌 1∼6).
한편, 도 3의 a)는 TFP상에 배치한 chip level의 가보 함수 및 그에 관련된 정보를 나타낸다. 특히, a-0)은 TFP상에 배치한 chip level의 가보 함수 gmm'(t)와 그 푸리에 변환(FT)인 Gmm'(f)를 나타내고, a-1)은 gmm'(t)를 (FD-Phase code(FD-PC))X'm'로 가중 가산한 TD-template의 실부, 허부를 나타내고, a-2)는 Gmm'(f)를 (TD-Phase Code(TD-PC))Xm으로 가중 가산한 FD-template의 실부, 허부를 나타낸다. 도 3의 b)는 TFP상의 NN'개의 상호 상관 함수(Cross Correlation Function(CCF)) 및 N'개의 열의 합의 TD-CCF값과 N개의 행의 합의 FD-CCF값을 나타낸다. 도 3의 c)는 Alternating Projection Theorem: APT에 기초하는 시간 제한 공간(TL-TD), 대역 제한 공간(BL-FD)상으로의 직교 투영에 의한 추정값
[수학식 71]
Figure 112019032179095-pct00071
의 교호 갱신 과정과 수렴값 td, fD를 나타낸다.
에코 신호로부터 delay td와 Doppler fD를 추정하는 측거 문제는, 시간, 주파수 시프트 작용소의 비가환성 유래의 PD에 포함되는 2개의 미지 파라미터 추정 문제이다. 따라서, 그 범주는 Weyl-Heisenberg Group: WHG에 기초하는 신호 검출·파라미터 추정론이다. 그러나, 일부 예외[비특허 문헌 17]를 제외하고, 많은 레이더 연구자는 비가환성을 무시한 추정론에 의거하므로, 측거의 고정도화에 성공하지 못하였다. 한편, 본 명세서에 기재된 발명은, 레이더 문제를 포함하는 무선 통신에 있어서, 비가환성이야말로 효율을 높이기 위해 중요하다는 본 발명자의 지견에 기초하고 있다. 또한, Wavelet론[비특허 문헌 8]에서는, 함수 f(t)에 대한 시간·주파수 시프트한 가보 함수[비특허 문헌 1]
[수학식 72]
Figure 112019032179095-pct00072
에 의한 전개
[수학식 73]
Figure 112019032179095-pct00073
의 계수 am,m'가 중심 과제이다. 또한, 5G, after 5G의 OFDM/OQAM이나 FBMC[비특허 문헌 9, 10, 12]에서는, am,m'를 전송 정보로 하는 다중화 신호 f(t)의 intersymbol interference(ISI) 및 interchannel interference(ICI)가 0(零)이 되는 함수 g(t)의 설계와 gmm'(t)의 비직교성의 해결이 주요 과제이다.
한편, 무선 통신에서는, 시간, 주파수 오프셋 내성의 동기법이 필요하지만, td, fD 추정의 시도는 거의 찾아볼 수 없다. 또한, 많은 통신 연구자는 신호 다중화를 위한 시간·주파수 시프트 mτ0, m'ν0에 수반하는 PD
[수학식 74]
Figure 112019032179095-pct00074
는 무시할 수 있다고 생각하고 있다. 그러나, WHG의 시간·주파수 시프트의 군론적 성질에 의해, 통신로의 PD
[수학식 75]
Figure 112019032179095-pct00075
에 이어 시프트의 PD
[수학식 76]
Figure 112019032179095-pct00076
외에 멀티 캐리어화에 수반하는 PD
[수학식 77]
Figure 112019032179095-pct00077
도 동시 발생하므로, 그 메카니즘은 단순하지 않다.
먼저, 레이더 문제에 잠재하는 이하의 3개의 과제에 임해야 할 것이다. 한편, 종래 레이더의 발사 신호로서 처프(chirp) 신호의 LFM-CW(Linear FM Continuous Wave)를 이용하거나, 혹은 짧은 시간의 펄스로 위상 변조한 압축 레이더나 그 멀티 캐리어 버전이 최근 제안되고 있다[비특허 문헌 19].
(과제 1)
첫번째 과제로서, 레이더는 원래 td, fD 2개의 미지 변수 문제이지만, 많은 수신기는 시간 shift τ와 주파수 shift ν의 2변수 복소수값 상관 함수, Ambiguity Function(AF)의 표시를 위해, AF의 절대값의 피크치 탐색이나 혹은 처프 신호의 AF 특성 이용에 기초하고 있는 것을 들 수 있다. 2개의 미지수 문제에는 2개 이상의 함수를 이용해야만 하는 것이다.
(과제 2)
두번째 과제로서, 레이더 문제에는 양자역학의 위치, 운동량 작용소와 마찬가지로, 비가환의 시간, 주파수 시프트에 의한 PD
[수학식 78]
Figure 112019032179095-pct00078
나 시간폭 Tp, 주파수 전이폭 Fp의 처프 펄스열에 의한 PD
[수학식 79]
Figure 112019032179095-pct00079
가 발생하는 것을 들 수 있다. 또한, td나 fD를 갖는 통신로를 개재해 시간폭 Ts와 주파수 대역 Fs의 데이터를 TFP상에서 무중첩 다중화 전송하는 방식(도 1 참조)에서도,
[수학식 80]
Figure 112019032179095-pct00080
에 이어 PD
[수학식 81]
Figure 112019032179095-pct00081
이 존재한다. 대부분의 연구자는 이 PD를 간과하고 있지만, PD의 존재는 정보 데이터를 담당하는 실수값 수신 신호의 공통의 문제이다.
(과제 3)
세번째 과제는, 두번째 과제와 관련되는데, PD 발생의 메카니즘을 알기 힘들다는 점에 있다. 즉, 통신이나 레이더 분야에서는 통상적으로 시간 시프트 연산자나 주파수 시프트 연산자는 각각
[수학식 82]
Figure 112019032179095-pct00082
로 정의되고, M(v)S(u)=e-i2πuvS(u)M(v)로부터 M(v)와 S(u)의 비가환성은 PD e-i2πuv의 지수의 어깨에 두 시프트량의 곱으로서 여실히 나타나 있다. 고정도 측거나 동기법의 해결의 단서는 여기에 집약되어 있다. 이 PD e-i2πuv의 취급(이후, TFSP(식 4, 24, 도 2 참조)라고 부른다)이 본 명세서에서 중요한 고찰 대상이다.
<2. 시간·주파수 대칭인 시간·주파수 시프트 연산자>
전형적인 에코 신호는
[수학식 83]
Figure 112019032179095-pct00083
으로 주어진다. 단,ψ(t)는 펄스의 complex envelope: CE라고 불리는 복소수값 신호이고,
[수학식 84]
Figure 112019032179095-pct00084
는, 각각 진폭, 도래 시간, 캐리어 주파수, 캐리어의 위상, 캐리어 주파수의 변위, 즉 참조 주파수 Ωr=2πfr로부터의 시프트(도플러 시프트)이다. 간단히 하기 위해 당분간 Ωr=0으로 둔다.
[수학식 85]
Figure 112019032179095-pct00085
로 푸리에 변환(Fourier Transform(FT))을 표기하고,
[수학식 86]
Figure 112019032179095-pct00086
을 ψ(t)의 FT라고 하면, re(t; td, fD)의 FT는
[수학식 87]
Figure 112019032179095-pct00087
이 된다. 여기에서, re(t; td, fD)와 Re(f; td, fD)의 쌍은 td, fD에 관해 대칭이 아니다. 미지수 td, fD의 곱이 TD 함수의 PD에만 나타나고 있다.
그러나, 약간 수정한 식[비특허 문헌 24, 26]
[수학식 88]
Figure 112019032179095-pct00088
은, TD 함수 x(t)와 그 FD 함수
[수학식 89]
Figure 112019032179095-pct00089
에 대해, td, fD에 관해, time-frequency symmetry property: TFSP를 만족하는 대칭인 시간 주파수 시프트 연산자(symmetrical time-frequency shift operators(symmetrical TFSO))의 정의식(도 2 참조)
[수학식 90]
Figure 112019032179095-pct00090
및, 2개의 시프트 연산자간의 항등식
[symmetrical TFSO의 성질 1]
[수학식 91]
Figure 112019032179095-pct00091
을 제공한다. 통상 이용되는, 시간 시프트 연산자 S(-td)x(t)=x(t-td)나 주파수 시프트 연산자
[수학식 92]
Figure 112019032179095-pct00092
와 달리, 식 (4)의 TD 신호 x(t)의 하프 시프트
[수학식 93]
Figure 112019032179095-pct00093
이나 FD 신호 X(f)의 하프 시프트
[수학식 94]
Figure 112019032179095-pct00094
는, 신호의 시간 주파수 표현론에서 이용되는
[수학식 95]
Figure 112019032179095-pct00095
또는
[수학식 96]
Figure 112019032179095-pct00096
과 비교해 약간의 수정처럼 보이지만, 이것은 식 (39), (44)와 같이 레이더 신호나 수신 신호의 위상 정보가 완전하게 TD, FD로 추적 가능한 표현법이다.
식 (4)의 TFSO는 양자역학의 von Neumann의 Canonical Commutative Relations: CCR[비특허 문헌 4, 6]과 동일하므로, 이후 von Neumann의 TFSO라고 부른다. 유명한 Heisenberg의 불확정성 원리라고 불리는 관계식[비특허 문헌 4, 5]
[수학식 97]
Figure 112019032179095-pct00097
을 통신의 시간 주파수 시프트 연산자[Tτ,0, T0,ν]에 대응시키면
[수학식 98]
Figure 112019032179095-pct00098
이 얻어지므로, 고전역학 극한[비특허 문헌 20]
[수학식 99]
Figure 112019032179095-pct00099
는, 통신의 무왜곡 조건
[수학식 100]
Figure 112019032179095-pct00100
과 대응한다. 단,
[수학식 101]
Figure 112019032179095-pct00101
은, 각각, 플랑크 상수(Planck constant)
[수학식 102]
Figure 112019032179095-pct00102
, 및 양자역학에서의 교환자를 나타낸다.
한편, TFSO의 연쇄법칙:
[수학식 103]
Figure 112019032179095-pct00103
으로부터 유도되는
[symmetrical TFSO의 성질 2]:
[수학식 104]
Figure 112019032179095-pct00104
에 있어서, 식 (7)은 식 (9)의 제1식의 예제이다.
symmetrical TFSO의 보다 중요한 성질은, TD, FD 신호의 대칭성으로부터 미지수의 곱 td·fD가 TD, FD 함수의 PD에 대칭적으로 나타나, 후술하는 바와 같이, 다중화 신호의 어드레스 정보가 PD로서 현재화되어 파라미터 추정에 중요한 역할을 한다는 점에 있다. 예를 들면, 무선 통신에서는[비특허 문헌 9, 10], OFDM 신호
[수학식 105]
Figure 112019032179095-pct00105
가 중심 과제이다. 단, 계수 am,n은 전송할 복소 데이터를 나타내고, x(t)는 시간 파형 함수이다. OFDM 신호는
[수학식 106]
Figure 112019032179095-pct00106
이라고 표현할 수 있다.
따라서, 무왜곡 조건 τ0ν0=1(식 (8) 참조)이 만족되고 있으면, PDeiπnτ0mν0은 amn의 부호를 별도로 하고 영향을 주지 않는다. 그러나, 오프셋 τ'=τ0τ, ν'=ν0ν이 있으면
[수학식 107]
Figure 112019032179095-pct00107
이 td, fD의 통신로의 PD
[수학식 108]
Figure 112019032179095-pct00108
에 이어 발생하므로, 위상 왜곡 PD를 포함하는 신호
[수학식 109]
Figure 112019032179095-pct00109
를 처리해야만 한다. 식 (10)과 같이 관용의 TFP상의 신호의 무중첩 중합법에서는[비특허 문헌 1, 2, 9, 10, 12], 식 (12)와 같이, 군론적 성질로부터, PD
[수학식 110]
Figure 112019032179095-pct00110
의 집적을 초래하므로, 각종 PD가 수신기 출력을 약하게 함으로써, 디지털 통신[비특허 문헌 3]에서 중요한 동기의 열화로 이어진다. 그럼에도 불구하고, 심볼간 간섭(inter-symbol interference: ISI)이나 채널간 간섭(inter-channel interference: ICI)을 삭감하는 시간 함수 x(t)의 설계는 중심적 과제이다. 이 관찰이 본 명세서의 출발점이다.
<3. 우도 함수와 CCF>
레이더 이론은, 레이더의 최적 시스템의 해석과 설계의 통계론을 논한 Woodward[비특허 문헌 15]나, 신호 검출 및 추정론의 포괄적인 연구를 실시한 Helstrom[비특허 문헌 18]의 교과서가 기초이다.
Wilcox의 연구[비특허 문헌 16]에 대한 Auslander와 Tolimieri에 의한 주의[비특허 문헌 17]와 같이, 이들은 비가환성을 고려하고 있지 않다는 것에 유의할 필요가 있다.
그러나, 레이더 이론의 기초는 이하의 신호 검출과 파라미터 추정론에 있다.
<3.1 신호 검출>
수신기에 에코 신호가 도달하면, 당해 에코 신호는 잡음에 섞여 있으므로, 당해 에코 신호의 결정은 불확정이 될 수 밖에 없다.
어느 정해진 형식의 신호 s(t)가 가우스 잡음 n(t) 중에 정해진 시간에 도달했는지 여부를 판정하는 문제인 신호 검출 문제를 생각한다. 관측 시간
[수학식 111]
Figure 112019032179095-pct00111
중에 측정한 수신기 입력 w(t)를 기초로 해, 수신기는 2개의 가설 검정:
H0, 「신호 없음」, 즉, w(t)=n(t)와,
H1, 「신호 있음」, 즉, w(t)=s(t)+n(t)
를 행한다. 관측 시간 중의 시간 t=tk에 측정치 wk=w(tk)가 얻어졌다고 하면, n개의 표본 wk는 가설 Hi, i=0, 1하에서의 결합 밀도 함수 joint probability density function(p.d.f.) pi(w)를 갖는 랜덤 변수이다. 수신기에서, 우도비
Λ(w)=p1(w)/p0(w), w=(w1, … , wn)
에 기초해 관측자의 최적 결정을 행한다.
우선, 임의의 결정 레벨 Λ0에 대해 관측자는, Λ(w)<Λ0이면 H0를, Λ(w)>Λ0이면 H1을 선택한다.
레이더 신호는
[수학식 112]
Figure 112019032179095-pct00112
로 표현된다[비특허 문헌 15, 18]. 단, ψ(t)는 CE(complex envelope)라고 불리는 복소수값 신호이고, Ω=2πfc는 캐리어 주파수이다. s(t)의 스펙트럼
[수학식 113]
Figure 112019032179095-pct00113
은 fc 근방과 -fc 근방에 좁은 피크를 갖고, 그 밴드폭이 Ω에 비해 작을 때, 신호는 narrowband(NB) 또는 quasi-harmonic이라고 불린다.
수신기 입력
[수학식 114]
Figure 112019032179095-pct00114
는 NB라고 가정하고, 그 CE
[수학식 115]
Figure 112019032179095-pct00115
는 변조기로 측정 가능하다고 가정한다.
공분산 함수
[수학식 116]
Figure 112019032179095-pct00116
을 갖는 정상 NB 백색 가우스 잡음하의 NB 신호
[수학식 117]
Figure 112019032179095-pct00117
의 최적 검출기는, 대수 우도 함수(logarithm of likelihood functional: LLF[비특허 문헌 18, p. 106]
[수학식 118]
Figure 112019032179095-pct00118
을 갖는다. 단, N0은 백색 잡음의 편측 스펙트럼 밀도를 나타내고, g 및 d2은 우도 함수(likelihood functional: LF) Λ[ψw(t)]의 통계량(statistic) 및 signal-to-noise ratio(SNR)를 나타낸다.
<3.2 파라미터 추정>
가설 검정의 원리는 다중 가설 검정에도 응용할 수 있다. 송신기로부터 M개의 신호 중 1개가 송신되었다고 하고, 수신기는 관측 시간(0, T)중에 M개의 신호 중 어느 것인지를 결정한다. 즉, 가설 Hk, 「신호 sk(t)가 송신되었다.」는 가설하에서 수신기 입력은
[수학식 119]
Figure 112019032179095-pct00119
이다. 단, ψk(t)는 NB CE를 나타내고, fk는 캐리어를 나타내고,
[수학식 120]
Figure 112019032179095-pct00120
은 sk(t)의 위상을 나타내고, n(t)는 랜덤 잡음을 나타낸다.
수신기에서는, 입력 w(t)의 측정치에 기초해 M개의 가설 중 1개를 선택한다. N개의 측정치 w1, … , wn에 대해 pk(w)를 가설 Hk 하에서의 p.d.f. ζk를, 가설 Hk 의 사전 확률이라고 한다. 간단히 하기 위해, ζk=M-1라고 하면, 신호 sk(t)의 직교성
[수학식 121]
Figure 112019032179095-pct00121
을 전제(Ei는 i번째 신호 에너지)로 한 다음에, 모든 k≠j에 대해 Λk(w)>Λj(w)일 때, 수신기는 Hk를 선택한다.
신호의 미지 파라미터를 θ1, … , θm이라고 하고, 이들을 m-차원 파라미터 공간 Θ중의 벡터 θ=(θ1, … , θm)로 나타낸다. 레이더 에코 신호는
[수학식 122]
Figure 112019032179095-pct00122
로 표기된다. 단, Ae는 감쇠 상수이고, A, td, fc, κ, fD는 진폭, 도착 시간, 캐리어 주파수, 위상, 도플러 시프트이다. 에코 신호의 미지 파라미터는 θ=(A, κ, td, fD)가 된다.
[수학식 123]
Figure 112019032179095-pct00123
이고, 잡음이 스펙트럼 밀도 N0를 갖는 백색일 때, LLF[비특허 문헌 18, p. 251]는
[수학식 124]
Figure 112019032179095-pct00124
이다. 변수 변환
[수학식 125]
Figure 112019032179095-pct00125
(여기에서,
[수학식 126]
Figure 112019032179095-pct00126
은 허부를 나타낸다)에 의해, θ'=(td, fD)를 갖는 A, κ의 최대 우도 추정값(maximum likelihood estimate(MLE))
[수학식 127]
Figure 112019032179095-pct00127
이 얻어진다. 따라서 나머지 미지 파라미터 θ'의 MLE는,
[수학식 128]
Figure 112019032179095-pct00128
을 최대화하는 파라미터값을 구하면 된다[비특허 문헌 18, p. 251]. 따라서, 수신기는 θ'의 추정에 전념할 수 있다. MLE θ'는 도플러 시프트
[수학식 129]
Figure 112019032179095-pct00129
에 근접하는 값의 집합의 1개에 대한 신호
[수학식 130]
Figure 112019032179095-pct00130
에 매칭한 필터 뱅크를 구성함으로써 얻어질 수도 있다. 단, WD는 최대의 fD이다. 그러나, 병렬의 NB 필터의 구성에서 통계량을 조사하는 것은 사실상 불가능하다. 이것이 미지 2변수 문제를 2개의 미지 1변수 문제로 분해한 최대의 이유이다.
식 (15)와 그 FT를 고쳐 쓰면
[수학식 131]
Figure 112019032179095-pct00131
단,
[수학식 132]
Figure 112019032179095-pct00132
는 td, fD를 갖는 von Neumann의 TFSO를 나타내고,
[수학식 133]
Figure 112019032179095-pct00133
은 설계할 CE, 위상(상세는 생략)이고,
[수학식 134]
Figure 112019032179095-pct00134
는 기저 대역 복소수값 신호 CE를 참조 주파수 fr로 통과역 신호로 시프트(변조)한 신호이다.
(td, fD)의 TFP상의 정확한 위치를 파악하기 위해, 주기 N의 TD-phase code(PC)와 주기 N'의 FD-PC를 이용해, TFP상의 2-차원 격자
[수학식 135]
Figure 112019032179095-pct00135
상의 위치 특정을 하자. 단, Tc=Ts/N, Fc=Fs/N', Ts 및 Fs는 칩 펄스 간격, 칩(서브) 캐리어 주파수 간격, 데이터 신호 시간폭 및 데이터 신호 대역폭이다. (td, fD)의 파라미터 공간 Θ'를 데이터 어드레스
[수학식 136]
Figure 112019032179095-pct00136
을 갖는 NN'개의 직사각형 영역
[수학식 137]
Figure 112019032179095-pct00137
로 분할해, 가설 「θ'는 영역
[수학식 138]
Figure 112019032179095-pct00138
중에 있다」를 Hm, m'로 나타낸다.
그러나, NN'개의 가설은, 후술하는 바와 같이, TD 신호 secho(t; A, κ, td, fD)의 fD를 추정하는 N'개의 가설, 및 FD 신호 Secho(f; A, κ, td, fD)의 td를 추정하는 N개의 가설로 분해할 수 있다.
신호 sk(t; θ')가 서로 직교하는 것을 전제로, 식 (14)의 가설 Hk 중의 CE ψk(t)와 위상
[수학식 139]
Figure 112019032179095-pct00139
를 갖는 k번째의 NB 에코 신호
[수학식 140]
Figure 112019032179095-pct00140
을 생각하자. 잡음이 스펙트럼 밀도 N0를 갖는 백색 가우스라고 하면, 그 LLF는[비특허 문헌 18, p. 129, p. 251]
[수학식 141]
Figure 112019032179095-pct00141
이 된다.
임의의 결정 레벨 r0에 대해
[수학식 142]
Figure 112019032179095-pct00142
를 만족하는 정수를 k=k0이라고 하면, 수신기는 k번째의 신호가 도달했다고 판정하고, 모든 gk가 r0 이하이면, 수신기는 신호 없음으로 판정한다. 이것은 ML 수신기이다. 따라서, 직교 신호의 설계가 중요하다. 식 (18)은 |gk(θ')|를 최대화하는 방법(하나는 피적분 함수의 최대화, 그 외는 캐리어의 위상 e 및 신호의 위상
[수학식 143]
Figure 112019032179095-pct00143
을 보상한다)을 나타내고 있다. 그러나, 통상, 위상 인자는
[수학식 144]
Figure 112019032179095-pct00144
의 재정의에 의해 흡수된다. 또한, 통상적으로, gk(θ')가 아니라 |gk(θ')|를 단순히 평가한다. 이와 같은 방법은, 위상 정보를 소거하게 된다. Woodward[비특허 문헌 15]는 Ville[비특허 문헌 14]의 ambiguity function(AF)이라 불리는 2차원의 CCF
[수학식 145]
Figure 112019032179095-pct00145
를 레이더 해석에 이용했다. 시간·주파수 시프트한 함수의 WHG-based의 비가환적이고 군론적 성질은 식 (9)의 위상 함수에 나타나고 있다. 또한,
[symmetrical TFSO의 성질 3]:
TD 신호 z(t)와 그
[수학식 146]
Figure 112019032179095-pct00146
에 대해, TD와 FD의 시간·주파수 시프트한 함수간의 내적(inner product(IP))은,
[수학식 147]
Figure 112019032179095-pct00147
이다. 단,
[수학식 148]
Figure 112019032179095-pct00148
은, r(t)와 s(t) 간의 TD-IP를 나타내고,
[수학식 149]
Figure 112019032179095-pct00149
[수학식 150]
Figure 112019032179095-pct00150
[수학식 151]
Figure 112019032179095-pct00151
사이의 FD-IP를 나타낸다.
식 (19)는, i) TD-, FD-IP의 실부는 t2=t1과 f2=f1일 때 최대가 되고, 그 때, AF의 최대값도 달성된다. ⅱ) IP의 좌우의 항을 송신·수신 신호라고 보면, 송신 신호의 비가환성 유래의 PD는 잘 설계된 수신기의 PD로 보상 가능한 것을 의미하고 있다. 식 (19)는, 변조 신호의 위상이 전기 공학의 교류의 전압, 전류의 「phasor」[비특허 문헌 18]와 같이 중요한 것과, 2개의 양 td, fD가 항상 PD에 나타나는 것을 시사하고 있다. 이는 TD의 matched filter 의존의 통상법과는 달리, WHG를 기초로 하는 TD와 FD에서의 (td, fD)-추정법을 실현하기 위한 큰 스텝이며, 본 명세서의 중요한 근간 부분이다. TD, FD에서 추적 용이한 위상항을 갖는 신호를 설계함으로써, 송신·수신기의 모두에서 위상을 효과적으로 이용 가능해진다.
Gabor[비특허 문헌 1]는, TFP상의 신호 해석과 시간과 주파수의 불확정성 관계 하한을 달성하는 가우스 함수의 중요성을 지적해, 이하의 함수 f의 시간·주파수 표현을 제시했다.
[수학식 152]
Figure 112019032179095-pct00152
Gabor의 가우스 함수 gm, n(t)의 집합은, TFP상에서 신호를 국재화하는 성질을 갖는다. 그러나, 이 기저는 직교하지 않고, 프레임[비특허 문헌 7]도 아니다. 또한, 많은 통신 연구자는 나이퀴스트 조건을 만족하지 않는 가우스 함수를 이용하지 않는다. 그러나, 본 명세서의 (td, fD)-추정법에서는, 가우스 함수의 많은 양호한 성질이 중요한 역할을 한다.
<4. TD-, FD-signature와 template>
TD-PC(phase code), 즉 스펙트럼 확산(spreading spectrum)[비특허 문헌 3]은, code-division multiple access(CDMA)를 실현한다. CE
[수학식 153]
Figure 112019032179095-pct00153
이 TFSP를 만족하도록, 그 FTΨ(f)도 위상 변조한다. 한편, 송신 신호 s(t)는 독립한 펄스 부호 c(t)로 전송 신호 m(t)를 변조해 s(t)=m(t)c(t)로 한다. 각 사용자는 신호가 직교하도록 부호가 할당됨으로써 넓은 주파수대 대역의 동시 사용을 가능하게 한다.
이하에서는,
[수학식 154]
Figure 112019032179095-pct00154
에 의해 표현되는 연속 시간 신호 대신
[수학식 155]
Figure 112019032179095-pct00155
에 의해 표현되는 이산 시간 신호를 고찰한다. TD 신호 s(t)를 Δt의 간격으로 표본화하고, 이산 FD 신호는 L-점의 이산 푸리에 변환(DFT)으로 정의한다. 즉, FD의 인접하는 bin의 주파수 간격은 Δf=1/(LΔt)이다.
[수학식 156]
Figure 112019032179095-pct00156
(분수
[수학식 157]
Figure 112019032179095-pct00157
의 절상)
[수학식 158]
Figure 112019032179095-pct00158
(분수
[수학식 159]
Figure 112019032179095-pct00159
의 절상)
의 이산 변수, 및, 칩 펄스의 직교성을 위해
[수학식 160]
Figure 112019032179095-pct00160
으로 하고
[수학식 161]
Figure 112019032179095-pct00161
이라고 하면, L-점의 twiddle factor
[수학식 162]
Figure 112019032179095-pct00162
가 정해진다.
이하, 7 종류의 이산 시간, 이산 주파수 신호를 정의한다.
[수학식 163]
Figure 112019032179095-pct00163
단, X=(X0, … , XN-1)∈{-1, 1}N은, 주기 N의 TD-PC이고, X'=(X'0, … , X'N'-1)∈{-1, 1}N'는 주기 N'의 FD-PC이고, χ=(X, X')이다.
support[-LΔt/2, LΔt/2](즉, 시간폭 LΔt)를 갖는 연속 시간 칩 펄스 g(t)에 대해, 지연 (D/2)Δt, D=L-1, L=(ΔtΔf)-1=MM'를 갖는[비특허 문헌 10], 인과적 이산 시간의 LΔt-시간 제한(time-limited: TL) 칩 펄스 g[k]
[수학식 164]
Figure 112019032179095-pct00164
가 얻어진다.
또한, support[-LΔf/2, LΔf/2], 즉, 대역폭 LΔf를 갖는, 이산 주파수의 LΔf-대역 제한(band-limited: BL) 칩 펄스 G[l]은, g[k]의 DFT
[수학식 165]
Figure 112019032179095-pct00165
에 의해 얻어진다. 한편, 칩 펄스 g[k](또는 G[l])의 펄스 간격 Tc=MΔt(또는 Fc=M'Δf)에 대해, 그 support 시간폭 LΔt(또는 대역폭 LΔf) 중에 좌우 양측으로부터 M'/2(또는 M/2)개의 칩 펄스의 간섭을 받는다. guard interval을 이용하지 않기 때문에, 이것이 관용법과 가장 다른 점이다.
다음 식으로 정의되는 이산 시간 TD-signature v[k;χ]와 그 FD-signature
[수학식 166]
Figure 112019032179095-pct00166
에는, 다음 식의 type-3의 TD-template
[수학식 167]
Figure 112019032179095-pct00167
과 type-4의 FD-template
[수학식 168]
Figure 112019032179095-pct00168
이 매장(embedded)되어 있다. 단,
[수학식 169]
Figure 112019032179095-pct00169
[수학식 170]
Figure 112019032179095-pct00170
은, 식 (4)의 von Neumann의 TD와 FD의 TFSO의 이산 버전으로,
[수학식 171]
Figure 112019032179095-pct00171
이다.
TD-signature v[k;χ]는 N'개의 TD-templates
[수학식 172]
Figure 112019032179095-pct00172
를 포함하고, 한편, FD-signature
[수학식 173]
Figure 112019032179095-pct00173
은 N개의 FD-templates
[수학식 174]
Figure 112019032179095-pct00174
를 포함하므로, PC의 이용에 의해, signature와 매장된 template 사이의 CCF는 큰 값을 갖는다. TD-template
[수학식 175]
Figure 112019032179095-pct00175
는 TFP상에서 Ts×LΔf의 직사각형 support를 갖고, FD-template
[수학식 176]
Figure 112019032179095-pct00176
은, LΔt×Fs의 직사각형 support를 갖는다. TFSO의 연쇄
[수학식 177]
Figure 112019032179095-pct00177
(또는
[수학식 178]
Figure 112019032179095-pct00178
)을 식 (22)에 적용하면, TD-signature와 FD-signature
[수학식 179]
Figure 112019032179095-pct00179
가 완전 대칭인 것을 알 수 있다.
[수학식 180]
Figure 112019032179095-pct00180
과 캐리어 주파수
[수학식 181]
Figure 112019032179095-pct00181
을 갖는 레이더 TD-신호 s[k;χ]와 그 FT의 FD-신호 S[l;χ]는,
[수학식 182]
Figure 112019032179095-pct00182
로 구해진다.
레이더 TD-신호의 CE와 그 DFT
[수학식 183]
Figure 112019032179095-pct00183
에 의해, TD-, FD-신호
[수학식 184]
Figure 112019032179095-pct00184
를 설계할 수 있다. 이는, 시간폭 Ts=NMΔt, 캐리어 간격 Fs=N'M'Δf의 signatures
[수학식 185]
Figure 112019032179095-pct00185
(또는,
[수학식 186]
Figure 112019032179095-pct00186
)를 무중첩으로, PP'개 중합한 2-차원 계열이다. 단,
[수학식 187]
Figure 112019032179095-pct00187
은 TFP상의 격자
[수학식 188]
Figure 112019032179095-pct00188
상의 데이터 어드레스
[수학식 189]
Figure 112019032179095-pct00189
를 갖는 데이터 심볼이다. 즉, 레이더 시스템에서는, 사전에는 그 범위를 알 수 없는 지연 td∈(0, PTs)와 Doppler
[수학식 190]
Figure 112019032179095-pct00190
을 탐색하기 위해, PTs×P'Fs의 시간폭·대역폭을 준비한다. 한편, 데이터 통신에서는 P·P'개의
[수학식 191]
Figure 112019032179095-pct00191
, 즉,
[수학식 192]
Figure 112019032179095-pct00192
를 송신한다.
여기에서,
[수학식 193]
Figure 112019032179095-pct00193
을 갖는 통신로를 개재해 신호 s[k;χ]를 송신했다고 가정한다. 이 때,
[수학식 194]
Figure 112019032179095-pct00194
에서 얻어지는 수신 TD-신호는
[수학식 195]
Figure 112019032179095-pct00195
로 표현된다. 단,
[수학식 196]
Figure 112019032179095-pct00196
은 수신 신호의 신호 성분의 CE이고, η[k] 및 ξ[k]는 간섭 성분 및 가우스 잡음이다. FD-신호 및 DFT
[수학식 197]
Figure 112019032179095-pct00197
의 상세는 생략했다. 비가환의 변·복조에 수반하는 PD
[수학식 198]
Figure 112019032179095-pct00198
은, e의 재정의에 의해 흡수시킬 수도 있지만, 후술하는 바와 같이, 상관 수신기에서 보상되어야 한다. 이와 같은 TD, FD 수신 신호는 수신기에서의 관측값 w와 그 DFT W를 제공한다.
독립 동일 분포(i.i.d.)의 TD-, FD-PC는, M-종 검정법에서 필요한 독립적인 N'개의 TD-templates
[수학식 199]
Figure 112019032179095-pct00199
및 N개의 FD-templates
[수학식 200]
Figure 112019032179095-pct00200
을 생성한다. PC에는 2개의 기능(신호의 난잡화 및 TFSO 유래의 PD 발생)이 있다. 다행히, PD 자체는 송신 신호에 추적 용이성을 부여한다는 의미에서, 파라미터 추정을 위한 좋은 표적이 된다. 이것은 PC 도입의 공과를 나타내고 있다. 실제로, 위상 변조 시스템의 대역폭은 고전적 레이더 시스템의 대역폭보다 N배 커지고, 또한 멀티 캐리어화, 즉 FD-PC는 sub-carriers 수 N'배의 대역폭을 더 필요로 한다.
<5. M종 가설 검정에 의한 TD-, FD-신호 검출과 추정>
드디어, 식 (27)의 CE
[수학식 201]
Figure 112019032179095-pct00201
(또는,
[수학식 202]
Figure 112019032179095-pct00202
)를 갖는 식 (26)의 레이더 신호
[수학식 203]
Figure 112019032179095-pct00203
(또는,
[수학식 204]
Figure 112019032179095-pct00204
)를 검출하는, M-종 가설 검정을 적용할 준비가 갖추어졌다. 수신기가 식 (17)에 관련지워진 NN'개의 가설 Hmm'를 선택하는 전략을 생각하자. 그것은, TD(또는 FD)의 LLF(또는 그 관련의 CCF)를 최대화하는 파라미터 θ',d를 찾아내면 충분하다. 먼저, TFP의 격자
[수학식 205]
Figure 112019032179095-pct00205
상의 어드레스
[수학식 206]
Figure 112019032179095-pct00206
을 갖는, 수신 TD-template CE
[수학식 207]
Figure 112019032179095-pct00207
의 검출 문제부터 생각한다. 단,
[수학식 208]
Figure 112019032179095-pct00208
이고,
[수학식 209]
Figure 112019032179095-pct00209
는 칩 어드레스 ρ'의 FD-PC TD-template((22) 참조)이고,
[수학식 210]
Figure 112019032179095-pct00210
은 kd의 추정 정수값이고,
[수학식 211]
Figure 112019032179095-pct00211
[수학식 212]
Figure 112019032179095-pct00212
를 추정하는 제어용 정수값이다. 이 CE는 식 (27)(식 (28)의
[수학식 213]
Figure 112019032179095-pct00213
참조)의 추정 수신 CE
[수학식 214]
Figure 112019032179095-pct00214
중에 매장되어 있다(도 3의 a-1) 참조). 단, 관계식
[수학식 215]
Figure 112019032179095-pct00215
를 이용했다. 식 (29)는, 다음 장에 기재하는 바와 같이, CE가 여러 가지 PD 유래의 의미 있는 위상을 포함하는 것을 나타낸다.
[수학식 216]
Figure 112019032179095-pct00216
과 상보적 CE
[수학식 217]
Figure 112019032179095-pct00217
을,
[수학식 218]
Figure 112019032179095-pct00218
로 나타낸다. 식 (22)와 식 (27)로부터, N'개의 TD-template
[수학식 219]
Figure 112019032179095-pct00219
가 이용 가능하므로, 수신기는 N'개의 CE
[수학식 220]
Figure 112019032179095-pct00220
을 이용해 N'개의 LLF 중에서 어느 것이 최대인지를 결정한다(도 3의 b) 참조).
i) TD의 이산 시간 신호 검출 및 Doppler-shift-ML 추정 문제:
NT개의 확률 변수 w=(w[0], … , w[NT-1])에 기초해 수신기는 2개의 가설
[수학식 221]
Figure 112019032179095-pct00221
을 선택한다. 단,
[수학식 222]
Figure 112019032179095-pct00222
는 시각 k의 관측값
[수학식 223]
Figure 112019032179095-pct00223
의 NB CE이고,
[수학식 224]
Figure 112019032179095-pct00224
는 가우스 잡음
[수학식 225]
Figure 112019032179095-pct00225
의 NB CE이고,
[수학식 226]
Figure 112019032179095-pct00226
은 관측 시간(0, T)의 샘플수이다. 가설 H1의 신호 성분은
[수학식 227]
Figure 112019032179095-pct00227
과 같다. N'개의 CE
[수학식 228]
Figure 112019032179095-pct00228
은 등(等)에너지이며,
[수학식 229]
Figure 112019032179095-pct00229
의 의미에서 의사(擬似)적 직교인 것을 전제로 하여, 다음의 결과를 얻는다.
명제 4:
스펙트럼 밀도 N0의 가우스 잡음하의 관측값 w=(w[0], … , w[NT-1])에 기초하는
[수학식 230]
Figure 112019032179095-pct00230
의 검출·추정의
[수학식 231]
Figure 112019032179095-pct00231
번째의 대수 우도는[비특허 문헌 18]
[수학식 232]
Figure 112019032179095-pct00232
이다. 단,
[수학식 233]
Figure 112019032179095-pct00233
이다.
임의의 결정 레벨 r0에 대해, 다음 식을 만족하는 정수값을
[수학식 234]
Figure 112019032179095-pct00234
라고 한다.
[수학식 235]
Figure 112019032179095-pct00235
단,
[수학식 236]
Figure 112019032179095-pct00236
은, 식 (31)에서의
[수학식 237]
Figure 112019032179095-pct00237
의, 식 (16)에서의 Ae-iκ의 MLE
[수학식 238]
Figure 112019032179095-pct00238
에 의한 보상 버전이다. 이 때, 수신기는 TFP상의 격자
[수학식 239]
Figure 112019032179095-pct00239
상의 어드레스
[수학식 240]
Figure 112019032179095-pct00240
의 CE가 도달했다고 판정하고, 만약, 모든
[수학식 241]
Figure 112019032179095-pct00241
이 r0 이하일 때, 수신기는 신호 없음으로 판정한다. 따라서,
[수학식 242]
Figure 112019032179095-pct00242
는 주어진
[수학식 243]
Figure 112019032179095-pct00243
하에서의
[수학식 244]
Figure 112019032179095-pct00244
의 MLE 이다. 한편,
[수학식 245]
Figure 112019032179095-pct00245
는, TFSO
[수학식 246]
Figure 112019032179095-pct00246
에 수반하는 위상 함수이다.
다음으로, 관측값 w의 DFT
[수학식 247]
Figure 112019032179095-pct00247
에 기초해, FD에서의 신호 검출·지연 추정 문제로 이동한다. 격자
[수학식 248]
Figure 112019032179095-pct00248
상의 어드레스
[수학식 249]
Figure 112019032179095-pct00249
를 갖는 수신 FD-template CE
[수학식 250]
Figure 112019032179095-pct00250
의 검출 문제를 생각한다. 단
[수학식 251]
Figure 112019032179095-pct00251
이고,
[수학식 252]
Figure 112019032179095-pct00252
는, 칩 어드레스 ρ의 TD phase-coded FD-template((22) 참조)이다. 이 CE는 식 (27)의 DFT
[수학식 253]
Figure 112019032179095-pct00253
의 수신 FD-CE
[수학식 254]
Figure 112019032179095-pct00254
중에 매장되어 있다(도 3의 a-2) 참조). 단
[수학식 255]
Figure 112019032179095-pct00255
[수학식 256]
Figure 112019032179095-pct00256
의 추정 정수값이고, kσ는 kd의 추정 제어용 정수값이고, 관계식
[수학식 257]
Figure 112019032179095-pct00257
을 이용했다.
[수학식 258]
Figure 112019032179095-pct00258
과 상보적 집합
[수학식 259]
Figure 112019032179095-pct00259
를,
[수학식 260]
Figure 112019032179095-pct00260
으로 나타낸다.
식 (22)와 식 (27)로부터 N개의 FD-templates
[수학식 261]
Figure 112019032179095-pct00261
이 이용 가능하므로, 수신기는 N개의 FD-CE
[수학식 262]
Figure 112019032179095-pct00262
를 이용해, N개의 LLF 중 어느 것이 최대인지를 결정하면 된다(도 3의 b) 참조).
ⅱ) FD에서의 신호 검출 및 지연-ML 추정 문제:
관측값 W=(W[0], … , W[NT-1])에 기초해, 수신기는 FD에서 이하의 2개의 가설
[수학식 263]
Figure 112019032179095-pct00263
을 선택한다. 단
[수학식 264]
Figure 112019032179095-pct00264
는, NB CE
[수학식 265]
Figure 112019032179095-pct00265
를 갖는 관측값 w[k]의 DFT이고,
[수학식 266]
Figure 112019032179095-pct00266
은 NB CE
[수학식 267]
Figure 112019032179095-pct00267
을 갖는 잡음 n[k]의 DFT이고,
[수학식 268]
Figure 112019032179095-pct00268
은 대역폭 B 중의 샘플수이며, 간단히 하기 위해 NB=NT로 한다.
[수학식 269]
Figure 112019032179095-pct00269
는 식 (33)의 CE
[수학식 270]
Figure 112019032179095-pct00270
을 갖는
[수학식 271]
Figure 112019032179095-pct00271
번째 template의 신호 스펙트럼이고,
[수학식 272]
Figure 112019032179095-pct00272
는 식 (34)의 CE
[수학식 273]
Figure 112019032179095-pct00273
을 갖는 상보 신호 스펙트럼이다. 가설 H'1의 신호 성분은
[수학식 274]
Figure 112019032179095-pct00274
와 같다. N개의 CE
[수학식 275]
Figure 112019032179095-pct00275
가 등에너지이며,
[수학식 276]
Figure 112019032179095-pct00276
의 의미에서 의사적 직교인 것을 전제로 하여, 다음의 결과를 얻는다.
명제 5:
스펙트럼 밀도 N0의 백색 가우스 잡음중의 관측값
[수학식 277]
Figure 112019032179095-pct00277
에 기초해,
[수학식 278]
Figure 112019032179095-pct00278
을 검출 및 추정하는
[수학식 279]
Figure 112019032179095-pct00279
번째 대수 우도 함수(LLF)는
[수학식 280]
Figure 112019032179095-pct00280
로 표시된다.
단,
[수학식 281]
Figure 112019032179095-pct00281
이고, 임의의 결정 레벨 r'0에 대해, ρ=ρ0
[수학식 282]
Figure 112019032179095-pct00282
를, 다음 식을 만족하는 정수로 한다.
[수학식 283]
Figure 112019032179095-pct00283
이 때, 수신기는 TFP의 격자
[수학식 284]
Figure 112019032179095-pct00284
상의 어드레스
[수학식 285]
Figure 112019032179095-pct00285
의 신호가 도착했다고 판정하고, 모든
[수학식 286]
Figure 112019032179095-pct00286
이 r'0 이하이면, 수신기는 신호 없음으로 판정한다.
[수학식 287]
Figure 112019032179095-pct00287
은 주어진
[수학식 288]
Figure 112019032179095-pct00288
하에서의 kd의 MLE이다. 한편,
[수학식 289]
Figure 112019032179095-pct00289
는 TFSO
[수학식 290]
Figure 112019032179095-pct00290
에 의한,
[수학식 291]
Figure 112019032179095-pct00291
의 FD 버전의 위상 함수이다.
<6. 파라미터 추정용 TD-, FD-CCF>
<6.1 TD-, FD-CCF>
관측값
[수학식 292]
Figure 112019032179095-pct00292
는 NB이고, 그 CE
[수학식 293]
Figure 112019032179095-pct00293
의 실부 및 허부는 각각 측정할 수 있다고 가정한다[비특허 문헌 18, 3]. 신호
[수학식 294]
Figure 112019032179095-pct00294
(또는,
[수학식 295]
Figure 112019032179095-pct00295
)를 검출하기 위해, 그리고, 식 (31)의
[수학식 296]
Figure 112019032179095-pct00296
을 갖는
[수학식 297]
Figure 112019032179095-pct00297
(또는, 식 (36)의
[수학식 298]
Figure 112019032179095-pct00298
을 갖는
[수학식 299]
Figure 112019032179095-pct00299
) 대신에 시간·주파수 대칭인 통계량을 얻기 위해, 이하에서 정의되는, 각각 관련지워진 CCF를 이용한다(도 3의 b) 참조).
Lemma 1: 식 (30)의 2개의 가설 H0, H1의 CE
[수학식 300]
Figure 112019032179095-pct00300
은 가우스 잡음이라고 가정하면, CCF
[수학식 301]
Figure 112019032179095-pct00301
이 성립한다. 단 <·,·>d,k는 이산 시간 함수의 공간
[수학식 302]
Figure 112019032179095-pct00302
에서의 내적(inner product(IP))이다. 따라서, (수신기 입력 CE
[수학식 303]
Figure 112019032179095-pct00303
대신) 감쇠 상수
[수학식 304]
Figure 112019032179095-pct00304
를 갖는 수신 CE
[수학식 305]
Figure 112019032179095-pct00305
(즉 식 (28)의 수신 CE
[수학식 306]
Figure 112019032179095-pct00306
의 신호 성분)와, 식 (29)의 어드레스
[수학식 307]
Figure 112019032179095-pct00307
에서의 추정 TD-template CE
[수학식 308]
Figure 112019032179095-pct00308
의 복소 임펄스 응답의 NB matched filter 사이의 CCF(type-3의 상관기라고 부른다)를,
[수학식 309]
Figure 112019032179095-pct00309
로 정의한다. 여기에서, 식 (29)의 부호 X'ρ', X 대신에, Y'ρ', Y를 갖는
[수학식 310]
Figure 112019032179095-pct00310
을 이용했다.
조금 계산하면, 이 CCF는
[수학식 311]
Figure 112019032179095-pct00311
로 표현되는 것을 알 수 있다.
유감스럽지만, AF θgg(τ,ν), ΘGG(ν,-τ)는, 일반적으로 다수의 사이드 로브(sidelobe)를 갖는다. 그러나, 가우스 칩 펄스 g(t)는, 그 AF의 τ, ν에 관한 변수 분리성과 지수 함수적 감쇠 특성
[수학식 312]
Figure 112019032179095-pct00312
에 의해, 추정 문제의 혁신적인 해답을 제공한다. 여기에서,
[수학식 313]
Figure 112019032179095-pct00313
이다.
[수학식 314]
Figure 112019032179095-pct00314
가 큰 값을 갖기 위해서는, N, N'≫1의 경우, θgg[·,·]의 제1 및 제2 변수는 모두 작아야만 한다. 즉, 식 (39)의
[수학식 315]
Figure 112019032179095-pct00315
의 항은 전부 무시 가능하다. 가우스 함수의 이 성질이
[수학식 316]
Figure 112019032179095-pct00316
을 최대화하는 lμ의 결정에 중심적 역할을 한다.
twiddle factor W의 멱수의 PD와 관련되는 PC의 3중합을 평가하기 위해, IDFT-형 합
[수학식 317]
Figure 112019032179095-pct00317
을 각괄호
[수학식 318]
Figure 112019032179095-pct00318
과 둥근 괄호(a)m'(편리한 기술법 (a)m'=W-am'로 정의한다)의 쌍으로 기호적으로 표시하고, DFT-형 합
[수학식 319]
Figure 112019032179095-pct00319
[수학식 320]
Figure 112019032179095-pct00320
의 쌍으로 기호적으로 표시하고, 기호
[수학식 321]
Figure 112019032179095-pct00321
을 이용하면,
Lemma 2: type-3의 수신기의 어드레스가,
[수학식 322]
Figure 112019032179095-pct00322
이고, Y=X, Y'=X'라면, CCF는
[수학식 323]
Figure 112019032179095-pct00323
으로 표현된다.
따라서,
[수학식 324]
Figure 112019032179095-pct00324
가 큰 값을 갖기 위해서는
[수학식 325]
Figure 112019032179095-pct00325
가 성립할 필요가 있다. 단, 식 (39)의 2번째 twiddle factor의 6개 항 중의 5개를 다시 나열하고
[수학식 326]
Figure 112019032179095-pct00326
나머지 1개를 1번째의 twiddle factor로 이동했다. 한편, FD에서는,
Lemma 3: 식 (35)의 2개의 가설 H'0, H'1
[수학식 327]
Figure 112019032179095-pct00327
은 가우스 잡음이라고 하면, CCF
[수학식 328]
Figure 112019032179095-pct00328
이 성립한다. 단, <·,·>d,l은 이산 주파수의 FD-함수의 공간
[수학식 329]
Figure 112019032179095-pct00329
의 IP이다. (FD-CE
[수학식 330]
Figure 112019032179095-pct00330
대신에) 감쇠 상수
[수학식 331]
Figure 112019032179095-pct00331
을 갖는 수신 FD-CE
[수학식 332]
Figure 112019032179095-pct00332
, 즉 수신 CE
[수학식 333]
Figure 112019032179095-pct00333
의 신호 성분의 FT와, 식 (33)의 어드레스
[수학식 334]
Figure 112019032179095-pct00334
에서의 추정 FD-template CE
[수학식 335]
Figure 112019032179095-pct00335
의 NB matched filter의 복소 임펄스 응답의 CCF(type-4의 상관기라고 부른다)를,
[수학식 336]
Figure 112019032179095-pct00336
으로 정의한다. 이 CCF는
[수학식 337]
Figure 112019032179095-pct00337
로 표현된다. 마찬가지로
[수학식 338]
Figure 112019032179095-pct00338
이 큰 값을 갖기 위해
[수학식 339]
Figure 112019032179095-pct00339
로 가정할 수 있으므로, N, N'≫1일 때, 식 (44)의
[수학식 340]
Figure 112019032179095-pct00340
의 모든 항은 무시할 수 있다. twiddle factor의 3개의 항은 합쳐져, PC의 3중합을 평가하면,
Lemma 4: type-4의 상관 수신기가 어드레스
[수학식 341]
Figure 112019032179095-pct00341
을 갖고, Y=X, Y'=X'이면, CCF는
[수학식 342]
Figure 112019032179095-pct00342
로 표현된다. 따라서
[수학식 343]
Figure 112019032179095-pct00343
이 큰 값을 갖기 위해서는,
[수학식 344]
Figure 112019032179095-pct00344
이어야만 한다. 단, 식 (44)의 2번째 twiddle factor의 6개 항 중의 5개를 다시 나열하고
[수학식 345]
Figure 112019032179095-pct00345
나머지 1개를 식 (44)의 1번째의 twiddle factor로 이동시키고, 가우스 함수의 변수 분리성을 이용하면, 식 (45)가 얻어진다. 식 (41) 및 (45)는
[수학식 346]
Figure 112019032179095-pct00346
[수학식 347]
Figure 112019032179095-pct00347
[수학식 348]
Figure 112019032179095-pct00348
에 관해 완전 대칭인 것을 나타내고 있다.
[수학식 349]
Figure 112019032179095-pct00349
및,
[수학식 350]
Figure 112019032179095-pct00350
의 쌍(complementary pair(CP)[비특허 문헌 27]라고 부른다)의 TD-PC와 FD-PC의 역할을 바꾸면
[수학식 351]
Figure 112019032179095-pct00351
[수학식 352]
Figure 112019032179095-pct00352
의 쌍(original pair(OP)[비특허 문헌 24, 30]라고 부른다)이 다음와 같이 얻어진다.
[OP의 TD-, FD-CCF:]
v[k;χ] 및
[수학식 353]
Figure 112019032179095-pct00353
은 다른 분해
[수학식 354]
Figure 112019032179095-pct00354
가 가능하다. 단, TD-template
[수학식 355]
Figure 112019032179095-pct00355
및 FD-template
[수학식 356]
Figure 112019032179095-pct00356
[수학식 357]
Figure 112019032179095-pct00357
이다.
[수학식 358]
Figure 112019032179095-pct00358
을 다음 식으로 정의되는, 격자
[수학식 359]
Figure 112019032179095-pct00359
상의 어드레스
[수학식 360]
Figure 112019032179095-pct00360
을 갖는 TD-template
[수학식 361]
Figure 112019032179095-pct00361
을 갖는 추정 수신 TD-CE
[수학식 362]
Figure 112019032179095-pct00362
로 하면, type-1 상관기라고 불리는 CCF
[수학식 363]
Figure 112019032179095-pct00363
이 얻어지고,
[수학식 364]
Figure 112019032179095-pct00364
로 표현된다. Y=X, Y'=X'으로 두고,
[수학식 365]
Figure 112019032179095-pct00365
의 항을 제외한 모든
[수학식 366]
Figure 112019032179095-pct00366
의 항은 무시할 수 있으므로
[수학식 367]
Figure 112019032179095-pct00367
이 얻어진다.
단, 식 (51)의 제2 twiddle factor의 6개 항 중의 5개를 다시 나열하여
[수학식 368]
Figure 112019032179095-pct00368
로 하고, 나머지 1개를 제1 twiddle factor로 이동시키고, 가우스 함수의 AF의 변수 분리성을 이용하면, 식 (52)가 얻어진다.
다음으로,
[수학식 369]
Figure 112019032179095-pct00369
를 다음 식으로 정의되는, 격자
[수학식 370]
Figure 112019032179095-pct00370
상의 어드레스
[수학식 371]
Figure 112019032179095-pct00371
의 FD-template
[수학식 372]
Figure 112019032179095-pct00372
를 갖는 추정 수신 FD-CE를,
[수학식 373]
Figure 112019032179095-pct00373
이라고 하면, type-2 상관기라고 부르는 관련 FD 상관기
[수학식 374]
Figure 112019032179095-pct00374
가 얻어지고,
[수학식 375]
Figure 112019032179095-pct00375
가 얻어진다. Y=X, Y'=X'로 하고
[수학식 376]
Figure 112019032179095-pct00376
을 만족하는
[수학식 377]
Figure 112019032179095-pct00377
을 제외하고 모든 항을 무시하면,
[수학식 378]
Figure 112019032179095-pct00378
이 얻어진다. 단, 식 (56)의 제2 twiddle factor의 6개 항 중의 5개를 다시 나열하고
[수학식 379]
Figure 112019032179095-pct00379
나머지 1개를 제1 twiddle factor로 이동시키면, 식 (57)이 얻어진다.
<6.2 PUL와 von Neumann의 APT>
만일, Lemma 2의
[수학식 380]
Figure 112019032179095-pct00380
이나 Lemma 4의
[수학식 381]
Figure 112019032179095-pct00381
에 관해, 칩 펄스(chip pulse)의 시간폭 LΔt나 대역폭 LΔf의 정도(精度)로 2개의 정확한 추정값
[수학식 382]
Figure 112019032179095-pct00382
가 얻어진다면, 2개의 CCF
[수학식 383]
Figure 112019032179095-pct00383
[수학식 384]
Figure 112019032179095-pct00384
에 각각 포함되어 있는 간섭 성분
[수학식 385]
Figure 112019032179095-pct00385
를 필터 제거해, 관용되는 첨예한 필터를 이용하지 않고
[수학식 386]
Figure 112019032179095-pct00386
을 복원할 수 있다. 이는 종래 통신에서 이용되고 있는 디지털 신호 처리와는 근본적으로 다르다. 두 쌍의 상관기의 추정값
[수학식 387]
Figure 112019032179095-pct00387
을 갱신하는 간단한 방법이 Phase-Updating Loop(PUL)라고 부르는 이하의 절차이며, 종래 통신의 동기에서 관용되는 「Phase-Locked Loop」와는 전혀 다르다.
[감쇠 상수의 MLE 갱신 부가 PUL의 알고리즘]:
[수학식 388]
Figure 112019032179095-pct00388
[수학식 389]
Figure 112019032179095-pct00389
를 각각 type-3 및 type-4의 상관기 어레이의 Complementary Pair(CP)라고 하고,
[수학식 390]
Figure 112019032179095-pct00390
[수학식 391]
Figure 112019032179095-pct00391
을, 각각 type-1 및 type-2의 상관기 어레이의 Original Pair(OP)라고 한다. 간단히 하기 위해 PUL의 알고리즘이 수렴할 때까지,
[수학식 392]
Figure 112019032179095-pct00392
로 둔다. 파라미터 θ'=(td, fD)의 s-step의 이산화 추정값
[수학식 393]
Figure 112019032179095-pct00393
및 식 (16)의
[수학식 394]
Figure 112019032179095-pct00394
대신에, 각각
[수학식 395]
Figure 112019032179095-pct00395
를 이용한 감쇠 상수
[수학식 396]
Figure 112019032179095-pct00396
의 s-step의 MLE
[수학식 397]
Figure 112019032179095-pct00397
[수학식 398]
Figure 112019032179095-pct00398
이라고 정의한다. 정수값의 쌍
[수학식 399]
Figure 112019032179095-pct00399
의 갱신은
[수학식 400]
Figure 112019032179095-pct00400
으로 정해진다.
[수학식 401]
Figure 112019032179095-pct00401
을 각각
[수학식 402]
Figure 112019032179095-pct00402
로 한다. 단, 초기값
[수학식 403]
Figure 112019032179095-pct00403
은 자유롭게 선택 가능하다. 예를 들면,
[수학식 404]
Figure 112019032179095-pct00404
로 설정할 수 있다. 칩 펄스 g[k], 및 G[l]의 시간폭 LΔt 및 대역폭 LΔf에 대해, 만일
[수학식 405]
Figure 112019032179095-pct00405
가 성립했을 때, (s+1)-스텝은 종료한다. 얻어지는 추정값은 MLE가 되고, 두 상관기는 ML 수신기가 된다.
왜곡이 있는 신호의 복원 문제는, 신호 처리의 중요한 한 분야이다. Youla[비특허 문헌 22]는 복원 문제의 해답을 주었다. Youla의 방법이나 기술법을 이용함으로써, PUL 알고리즘의 수렴이 어떻게 von Neumann의 APT[비특허 문헌 21]에 의존하고 있는지를 알 수 있다.
내적(inner product(IP))
[수학식 406]
Figure 112019032179095-pct00406
(또는
[수학식 407]
Figure 112019032179095-pct00407
), 및 놈(norm)
[수학식 408]
Figure 112019032179095-pct00408
(또는
[수학식 409]
Figure 112019032179095-pct00409
)을 갖는, 2승 가산가능의 연속적인 이산 시간 함수(또는, 이산 주파수 함수)로 이루어지는 힐베르트 공간
[수학식 410]
Figure 112019032179095-pct00410
을 생각하자. ε을
[수학식 411]
Figure 112019032179095-pct00411
에서의 임의의 폐다양체(closed linear manifold(CLM))로 한다. 투영 정리[비특허 문헌 22]에 의해, ε' 및 ε"를 서로 직교하는
[수학식 412]
Figure 112019032179095-pct00412
의 부분 공간으로 하면, 임의의 f∈ε은 일의적(一意的)인 분해 f=g+h, g∈ε', h∈ε"를 갖는다. 단, g, h는 f의 ε', ε" 상으로의 투영을 가리키고, g=Pf, h=Qf라고 표기한다. P는 ε' 상으로의 투영 연산자(projection operator(PO))이고, Q=I-P는 ε" 상으로의 PO이고, I는 항등 연산자이다.
ε1(또는 ε3)을 LΔt-(또는 Ts-) 시간 제한(Time limited(TL)) 신호인
[수학식 413]
Figure 112019032179095-pct00413
의 전부로 이루어지는 집합으로 한다. 한편, LΔf-(또는 Fs-) 대역 제한(band limited(BL)) 신호인
[수학식 414]
Figure 112019032179095-pct00414
의 전부로 이루어지는 집합을 ε2(또는 ε4)로 한다.
[수학식 415]
Figure 112019032179095-pct00415
는 CLM[비특허 문헌 22]이 된다.
Pi
[수학식 416]
Figure 112019032179095-pct00416
상으로 투영하는 투영 연산자로 하고, Qi=I-Pi를 εi의 직교보공간
[수학식 417]
Figure 112019032179095-pct00417
상으로 투영하는 투영 연산자로 한다. CCF는 이하의 의미에서 PO의 역할을 한다. 즉, 임의의 신호 r 및 s에 대해 신호 r∈ε는 일의적 분해
[수학식 418]
Figure 112019032179095-pct00418
을 갖는다. 단,
[수학식 419]
Figure 112019032179095-pct00419
는 s의 직교 공간이다. 상관계수
[수학식 420]
Figure 112019032179095-pct00420
은 r을 ε' 상으로 투영하는 투영 연산자라고 볼 수 있다.
type-3, type-4의 상관기의 쌍 CP(Complementary pair)와, type-1, type-2의 상관기의 쌍 OP(Original Pair)는, 다음 식과 같이 직교 투영 연산자(PO)가 된다.
[수학식 421]
Figure 112019032179095-pct00421
단,
[수학식 422]
Figure 112019032179095-pct00422
는 각각
[수학식 423]
Figure 112019032179095-pct00423
의 직교보공간이다.
Alternating Projection Theorem(APT)(도 11 참조)을 이용하면 이하의 결과가 얻어진다. 한편, APT에 의하면[비특허 문헌 21, p. 55, theorem 13.7], E, F를 각각 힐베르트 공간의 CLM ε,
[수학식 424]
Figure 112019032179095-pct00424
상으로의 투영 연산자로 할 때, 연산자의 계열 E, FE, EFE, FEFE, … 은 극한 G를 갖는다. 또한, 계열 F, EF, FEF, …도 같은 극한 G를 갖는다. 또한, G는
[수학식 425]
Figure 112019032179095-pct00425
상으로의 투영 연산자가 된다.(대칭 조건 EF=FE는 불요.)
[Theorem: PUL 알고리즘의 수렴 정리](도 3의 c) 참조):
식 (60)의 s-스텝의 추정값
[수학식 426]
Figure 112019032179095-pct00426
(또는 (61))을 갖는
[수학식 427]
Figure 112019032179095-pct00427
및, 마찬가지로 식 (60)의 s-스텝의 추정값
[수학식 428]
Figure 112019032179095-pct00428
(또는 (61))을 갖는 (ρ, kσ)(또는 (ρ', kσ))에 관해 행하는 argmax-연산을 생각한다. 당해 argmax-연산에서 정해지는 (s+1)-스텝의 최대 우도 추정값을 포함하는 TD 및 FD의 4개의 PO를, 간략화한 기호로
[수학식 429]
Figure 112019032179095-pct00429
[수학식 430]
Figure 112019032179095-pct00430
으로 표기한다. PUL 알고리즘은 수렴한다.
[증명]:
CP의 PO 쌍(Ts-TL-PO P3, Fs-BL-PO P4)(OP의 PO 쌍(LΔt-TL-PO P1, LΔf-BL-PO P2))의 적용순으로 2개의 상이한 재귀식이 얻어진다. 한편, TD-PC X와 FD-PC X', 즉 첨자(3, 4)와 (1, 2)를 바꿔 넣으면 OP는 CP와 같으므로, CP의 증명만을 제공한다.
먼저,
[수학식 431]
Figure 112019032179095-pct00431
의 복원 알고리즘을 제공한다.
[수학식 432]
Figure 112019032179095-pct00432
라고 하면,
[수학식 433]
Figure 112019032179095-pct00433
이므로,
[수학식 434]
Figure 112019032179095-pct00434
가 얻어지고,
[수학식 435]
Figure 112019032179095-pct00435
는 연산자 방정식
[수학식 436]
Figure 112019032179095-pct00436
을 만족하므로, 방정식
[수학식 437]
Figure 112019032179095-pct00437
은, 이하의 TD의 재귀식[비특허 문헌 22, 23]
[수학식 438]
Figure 112019032179095-pct00438
을 제공한다. APT에 의해
[수학식 439]
Figure 112019032179095-pct00439
[수학식 440]
Figure 112019032179095-pct00440
의 CLM εc=⊥ ε3∩ε4 상으로의 PO가 된다. εc는 항등적으로 0(零)인 함수만 포함하므로[비특허 문헌 22, p. 699][비특허 문헌 23, p. 637],
[수학식 441]
Figure 112019032179095-pct00441
이 된다. 이는, Youla[비특허 문헌 22]의 주요한 결과의 하나이다. 따라서, 결합 PO
[수학식 442]
Figure 112019032179095-pct00442
는, 격자 공간
[수학식 443]
Figure 112019032179095-pct00443
상의 chip과 데이터 어드레스
[수학식 444]
Figure 112019032179095-pct00444
의 LΔt×LΔf의 직사각형 영역(CE
[수학식 445]
Figure 112019032179095-pct00445
의 support와
[수학식 446]
Figure 112019032179095-pct00446
의 support의 곱집합:
[수학식 447]
Figure 112019032179095-pct00447
)을 뽑아내고, TFP의 나머지를 필터 제거한다. 단, Δa, Δb, Δa', Δb'는, Δa+Δb=M, Δa'+Δb'=M'를 만족하는 정수이다. 이와 같은 PO 쌍은 위상 공간(또는 시간·주파수 공간(TFP))의 국소 선택 연산자(localization operator)[비특허 문헌 7]라고 불린다. CE ψ[k]가 신호 복원 가능하므로, 파라미터 kd
[수학식 448]
Figure 112019032179095-pct00448
은 LΔt 및 LΔf의 정도(精度)로 추정 가능해진다. 통상의 엄밀한 TL(또는 BL) 연산자[비특허 문헌 23]와 달리, TL-PO P3(또는 BL-PO P4)은 N개의 위상 변조한 TD-가우스 칩 펄스 g[k](또는 N'개의 위상 변조한 FD-가우스 칩 펄스
[수학식 449]
Figure 112019032179095-pct00449
)를 시간 주파수 시프트해, guard interval을 마련하지 않고 무중첩으로 중합한 신호를 template로서 정의하고 있다.
반대로, FD에서는, Ψ∈ε3이면,
[수학식 450]
Figure 112019032179095-pct00450
[수학식 451]
Figure 112019032179095-pct00451
이 성립한다. 즉, Ψ는 연산자 방정식[비특허 문헌 23]
[수학식 452]
Figure 112019032179095-pct00452
를 만족하므로, FD의 재귀식
[수학식 453]
Figure 112019032179095-pct00453
이 얻어진다. APT에 의해,
[수학식 454]
Figure 112019032179095-pct00454
는 (Ψ0-Ψ)의 CLM
[수학식 455]
Figure 112019032179095-pct00455
상으로의 PO가 된다. 이 CLM는 영함수만이므로, limi→∞ Ψi=Ψ이다.
FD CE
[수학식 456]
Figure 112019032179095-pct00456
이 복원되므로, 파라미터
[수학식 457]
Figure 112019032179095-pct00457
은 LΔt, LΔf의 정도로 추정 가능해진다.
[수학식 458]
Figure 112019032179095-pct00458
은 다른 하나의 국소 선택 연산자이며, 격자 공간
[수학식 459]
Figure 112019032179095-pct00459
상의 chip과 데이터 어드레스
[수학식 460]
Figure 112019032179095-pct00460
의 LΔt×LΔf의 직사각형 영역(TD-CE
[수학식 461]
Figure 112019032179095-pct00461
의 직사각형 support와 FD-CE
[수학식 462]
Figure 112019032179095-pct00462
의 직사각형 support의 곱집합:
[수학식 463]
Figure 112019032179095-pct00463
을 뽑아내고, 나머지를 필터 제거한다.(증명 끝).
<7. Twinned FBMC>
<7.1 SFB: signature 송신기와 레이더 신호 송신기>
식 (25)의 TD-signature v[k;χ] 및 FD-signature
[수학식 464]
Figure 112019032179095-pct00464
를 생성하는 multicarrier filter bank(FBMC)를 제공한다. 주기 N의 TD-PC Xm을 반복해 이용함으로써 시간 t의 표본점 -∞, … , -Tc, 0, Tc, … , ∞로 무한 계열을 생성하고, 한편, 주기 N'의 FD-PC X'm'를 반복해 이용함으로써 주파수 f의 표본점 -∞, … , -Fc, 0, Fc, … ,∞로 무한 계열을 생성하면,
명제 6:
식 (25)의 v[k;χ] 및
[수학식 465]
Figure 112019032179095-pct00465
는 각각 다음 식으로 고쳐 쓸 수 있다.
[수학식 466]
Figure 112019032179095-pct00466
단,
[수학식 467]
Figure 112019032179095-pct00467
은 다음 식으로 정의되는 변조 필터이다.
[수학식 468]
Figure 112019032179095-pct00468
한편, Vaidyanathan[비특허 문헌 13]은, 입력 x[k] 및 출력 y[n]에 대해 필터 계수 h[·]를 갖는 3종류의 multi-rate filter: 축소율 Mf의 데시메이션(decimation) 필터:
[수학식 469]
Figure 112019032179095-pct00469
, 확대율 Lf의 보간(interpolation) 필터:
[수학식 470]
Figure 112019032179095-pct00470
, 및, 축소율 Mf/Lf의 데시메이션 필터:
[수학식 471]
Figure 112019032179095-pct00471
을 정의했다.
[증명] 식 (25)의 TD-signature v[k;χ](또는 FD
[수학식 472]
Figure 112019032179095-pct00472
)는, N'(또는 N)개의 부대역(sub-band) 및 확대율 M(또는 M')을 갖고, 각 m'th 부대역(또는 mth 부대역)이 FD-PC X'm'(또는 TD-PC Xm)로 위상 변조된 synthesis filter bank(SFB)[10, 12]의 출력에서 얻어지는 신호라고 볼 수 있다. 실제로,
[수학식 473]
Figure 112019032179095-pct00473
이 입력 신호,
[수학식 474]
Figure 112019032179095-pct00474
가 해당 SFB의 (m, m')th 부대역 필터라면, 출력
[수학식 475]
Figure 112019032179095-pct00475
는 식 (67)이 된다.(증명 끝).
한편, 기호 lcm[M, N']=M0N'=MN'0, lcm[M', N]=M'0N=M'N0을 이용하면(비특허 문헌[13], [10], [11]), M0N', M'0N개의 polyphase component
[수학식 476]
Figure 112019032179095-pct00476
의 도입에 의해, 각각 polyphase filter(Vaidyanathan [13]의 Type 1 polyphase)
[수학식 477]
Figure 112019032179095-pct00477
이 정의된다. 한편, lm=lcm[M, N'], gd=gcd[M, N']에 대해, lm=MN'0=N'M0을 만족하는 각각 N'0, M0이 존재하고, 한편, 항등식 lm·gd=M·N'로부터 M=gd·M0, N'=gdN'0이 성립한다.
따라서, 2차원의 TD-, FD-PC로 위상 변조된 도 4 및 도 5에 나타낸 signature v[k], V [l]의 SFB가 얻어진다.
한편, 도 4는, TD-, FD-PCs Xm, X'm' 및 m'-번째의 TD template
[수학식 478]
Figure 112019032179095-pct00478
을 갖는 식 (25), (67)의 TD signature v[k]를 생성하는 Synthesys Filter Bank(SFB)를 나타낸다.
또한, 도 5는, TD-, FD-PCs Xm, X'm' 및 m'-번째의 TD template
[수학식 479]
Figure 112019032179095-pct00479
를 갖는 식 (25), (67)의 FD signature
[수학식 480]
Figure 112019032179095-pct00480
을 생성하는 SFB를 나타낸다.
한편, 식 (27)은
[수학식 481]
Figure 112019032179095-pct00481
[수학식 482]
Figure 112019032179095-pct00482
를 생성하는 SFB를 제공한다.
명제 7:
CE ψ[k;χ] 및
[수학식 483]
Figure 112019032179095-pct00483
은 다음 식과 같이 기술할 수 있다.
[수학식 484]
Figure 112019032179095-pct00484
단,
[수학식 485]
Figure 112019032179095-pct00485
는 다음 식으로 정의되는 변조 필터이다.
[수학식 486]
Figure 112019032179095-pct00486
[증명] 식 (71)은 이하를 의미한다.
[수학식 487]
Figure 112019032179095-pct00487
(또는,
[수학식 488]
Figure 112019032179095-pct00488
)이 입력 신호이고,
[수학식 489]
Figure 112019032179095-pct00489
(또는,
[수학식 490]
Figure 112019032179095-pct00490
)가 해당 SFB의 (q, q')th 부대역의 필터이면,
[수학식 491]
Figure 112019032179095-pct00491
(또는 FD-CE
[수학식 492]
Figure 112019032179095-pct00492
)은 PP'개의 부대역을 갖고, 한편, 각 부대역에서 확대율 NM(또는 N'M')의 SFB의 출력이다.
따라서, 도 6, 7의 정보 데이터가 매장된 송신 신호의 SFB가 얻어진다. 통상의 SFB에서는 도 4, 5의 signature 생성 과정이 없다. 도 6, 7에서 N=N'=1의 경우가, 통상의 SFB에 대응한다. 한편, 데이터 전송의 modulated filter(MF) 특성
[수학식 493]
Figure 112019032179095-pct00493
을 이용하면, 각각 polyphase filter(Vaidyanathan[비특허 문헌 13]의 Type 1 polyphase)
[수학식 494]
Figure 112019032179095-pct00494
가 정의된다. 단, P'0, P0은, P'0MN=lcm[P', MN], P0MN=lcm[P, MN]를 만족하는 정수이다.(증명 끝).
한편, 도 6은, 복소수값 데이터
[수학식 495]
Figure 112019032179095-pct00495
를 입력 신호로 하는 식 (27)의 TD-Complex Envelope(TD-CE)
[수학식 496]
Figure 112019032179095-pct00496
을 생성하는 SFB를 나타낸다.
또한, 도 7은, 복소수값 데이터
[수학식 497]
Figure 112019032179095-pct00497
을 입력 신호로 하는 식 (27)의 FD-Complex Envelope(FD-CE)
[수학식 498]
Figure 112019032179095-pct00498
을 생성하는 SFB를 나타낸다.
<7.2 AFB: 수신기와 복호기>
감쇠 상수 Ae 외에 비가환인 변복조의 PD
[수학식 499]
Figure 112019032179095-pct00499
를 포함하는, 식 (28)의 수신 신호
[수학식 500]
Figure 112019032179095-pct00500
(또는, 그 FT
[수학식 501]
Figure 112019032179095-pct00501
)(
[수학식 502]
Figure 112019032179095-pct00502
로 간략히 기재한다)와, 식 (29)의 추정 template TD-CE
[수학식 503]
Figure 112019032179095-pct00503
(또는 식 (33)의 추정 template FD-CE
[수학식 504]
Figure 112019032179095-pct00504
) 사이의 type-3(또는 type-4) CCF를 각각
[수학식 505]
Figure 112019032179095-pct00505
로 한다. 단, PD
[수학식 506]
Figure 112019032179095-pct00506
(또는,
[수학식 507]
Figure 112019032179095-pct00507
)은, r[k](식 (28), 식 (38)의
[수학식 508]
Figure 112019032179095-pct00508
참조)(또는,
[수학식 509]
Figure 112019032179095-pct00509
(식 (28), 식 (43)의
[수학식 510]
Figure 112019032179095-pct00510
참조))의 신호 성분의 변복조 PD
[수학식 511]
Figure 112019032179095-pct00511
의 보상항이다.
TD-(또는 FD-) analysis filter bank(AFB)[비특허 문헌 10, 12]가 다음과 같이 얻어진다.
명제 8:
type-3과 type-4의 CCF는 각각
[수학식 512]
Figure 112019032179095-pct00512
이다. 단, D=L-1이고,
[수학식 513]
Figure 112019032179095-pct00513
은 각각 TD-, FD-AFB의 변조 필터
[수학식 514]
Figure 112019032179095-pct00514
이다. 한편, N=N'=1의 경우, 통상의 AFB[비특허 문헌 10, 11]에 상당하고, P×P'개의
[수학식 515]
Figure 112019032179095-pct00515
이 전송되면,
[수학식 516]
Figure 112019032179095-pct00516
을 추정하기 위해 상관기 어레이의 어드레스
[수학식 517]
Figure 112019032179095-pct00517
(p-번째 대역, p'-번째 시간)의 출력값을 평가한다.
[증명]: 해당 AFB가 각 부대역에서 축소율 NM(또는 N'M')의 P개의 부대역(또는 P'개의 부대역)을 갖고,
[수학식 518]
Figure 112019032179095-pct00518
번째(또는,
[수학식 519]
Figure 112019032179095-pct00519
번째) 부대역의 입력 신호가, 식 (75) (또는 식 (76))의 위상 변조 수신 TD-신호
[수학식 520]
Figure 112019032179095-pct00520
(또는 FD-신호
[수학식 521]
Figure 112019032179095-pct00521
)이라고 하면, 다음의 2개의 사실:
(i) TD 신호의 대칭성[비특허 문헌 10]:
[수학식 522]
Figure 112019032179095-pct00522
(식 (21) 참조);
(ⅱ) TD 신호의 대칭성으로부터 계승된 FD 신호의 성질
[수학식 523]
Figure 112019032179095-pct00523
으로부터, 필터 계수
[수학식 524]
Figure 112019032179095-pct00524
(또는,
[수학식 525]
Figure 112019032179095-pct00525
)는 식 (77)이라고 할 수 있다.(증명 끝).
또한, 식 (77)의 AFB의 필터 특성과 P'0MN, P0M'N'개의 polyphase components
[수학식 526]
Figure 112019032179095-pct00526
을 이용하면, 각각 N', N개의 polyphase filter(Vaidyanathan[비특허 문헌 13]의 Type 2 polyphase)
[수학식 527]
Figure 112019032179095-pct00527
이 정의되고, 도 8, 9의 AFB가 얻어진다. N=N'=1, N0=N'0=1의 경우가 통상의 AFB에 대응한다.
TD-AFB(또는 FD-AFB)의 2값 심볼은,
[수학식 528]
Figure 112019032179095-pct00528
(또는,
[수학식 529]
Figure 112019032179095-pct00529
)의 부호로 정해진다. 단
[수학식 530]
Figure 112019032179095-pct00530
은, MLE값
[수학식 531]
Figure 112019032179095-pct00531
(식 (59) 참조)을 갖는 감쇠 상수 Ae의 MLE이다. TD 및 FD에 관해 대칭인 도 4∼도 7의 SFB, 및, 도 8, 9의 TD-AFB, FD-AFB의 쌍을 「twinned-FBMC」라고 부른다.
한편, 도 8은, 복소수값 데이터
[수학식 532]
Figure 112019032179095-pct00532
복호를 위한 TD 상관기 어레이를 구비한 Analysis Filter Bank(AFB)를 나타낸다.
또한, 도 9는, 복소수값 데이터
[수학식 533]
Figure 112019032179095-pct00533
복호를 위한 FD 상관기 어레이를 구비한 AFB를 나타낸다.
한편, 도 8, 도 9에서의
[수학식 534]
Figure 112019032179095-pct00534
[수학식 535]
Figure 112019032179095-pct00535
는, 각각, 전술한
[수학식 536]
Figure 112019032179095-pct00536
[수학식 537]
Figure 112019032179095-pct00537
에 대응하고 있다.
PUL 알고리즘을, 도 10과 같이 TD-AFB 및 FD-AFB 간의 인터페이스로서 실장하면, 얻어진 시변·적응형 AFB는 레이더에 대해서는 파라미터 추정기가 되고, 다른 통신 시스템에 대해서는 PUL 알고리즘의 수렴 후 동기기가 된다. 또한 이것은, 데이터 통신 시스템에 대해서는 복소 CCF값
[수학식 538]
Figure 112019032179095-pct00538
[수학식 539]
Figure 112019032179095-pct00539
에 기초하는
[수학식 540]
Figure 112019032179095-pct00540
의 복호기가 된다. 한편, 이 FBMC의 식 (68), 식 (72) 및 식 (77)에서 정의되는 필터는, 모두 다상 필터(Vaidyanathan의 type1- 또는 type2-다상형[비특허 문헌 13]이라고 불린다)로 실현될 수 있는데, 상세는 생략한다.
한편, 도 10의 (a)는 TD 상관기 어레이를 나타내고, (c)는 FD 상관기 어레이의 Analysis Filter Bank(AFB)를 나타내고, (b)는 von Neumann의 APT에 의한 두 어레이의 최대 우도 추정값의 교호 갱신을 모식적으로 나타낸다.
<8. 시간, 주파수 시프트의 비가환성을 이용한 통신의 다른 예>
비가환성을 이용한 통신의 전형예의 레이더 문제에서는, 통신로의 delay td, Doppler shift fD를 파라미터로 하는 TFSO
[수학식 541]
Figure 112019032179095-pct00541
및, 수신기에서 필요한 지연(delay), 도플러 시프트(Doppler shift)의 추정값
[수학식 542]
Figure 112019032179095-pct00542
를 갖는 TFSO
[수학식 543]
Figure 112019032179095-pct00543
의 쌍에 관해, 각각의 송신·수신 TD, FD 신호로의 사전의 하프 시프트가 td, fD의 파라미터 추정에 유용하다는 것을 분명히 했다. 본 절에서는 비가환성을 이용한 통신의 구체예를 두세 가지 설명한다.
지금까지는 간략화를 위해 단일 타깃을 논의했다. 복수 개 타깃으로의 확장은 어렵지 않다. 예를 들면 식 (32), 식 (37)의 결정 레벨 r0, r'0을 이용함으로써 타깃 탐색 공간 중의 복수 개 타깃이 검출 가능해진다. 그 밖에도,
<8.1 CDMT에 의한 복수 개 타깃 검출>
예를 들면, 타깃 탐색 공간 Θ'=[0, T)×[0, F)를 4분할:
[수학식 544]
Figure 112019032179095-pct00544
하고, 각각의 공간에 TD-, FD-PC:
[수학식 545]
Figure 112019032179095-pct00545
를 할당하고, 2차원 PC
[수학식 546]
Figure 112019032179095-pct00546
의 칩 어드레스를
[수학식 547]
Figure 112019032179095-pct00547
로 나타낸다. 또한, signature를
[수학식 548]
Figure 112019032179095-pct00548
로 한다. 단, χ(1), χ(2), χ(3), χ(4)
[수학식 549]
Figure 112019032179095-pct00549
에 대응한다. 레이더 송신 신호의 CE는
[수학식 550]
Figure 112019032179095-pct00550
이고, 다중 타깃 검출용이다.
Npath개의 doubly dispersive chanel의 delay, Doppler, 감쇠 상수
[수학식 551]
Figure 112019032179095-pct00551
의 통신로를 거친 수신 신호의 신호 성분은
[수학식 552]
Figure 112019032179095-pct00552
가 된다. 한편, 수신기에서는, 식 (29), 식 (33)의 TD-, FD-CE
[수학식 553]
Figure 112019032179095-pct00553
에, 2차원 PC χ(또는
[수학식 554]
Figure 112019032179095-pct00554
)를 대입해, 식 (38), 식 (43)의 type-3, type-4의 CCF
[수학식 555]
Figure 112019032179095-pct00555
[수학식 556]
Figure 112019032179095-pct00556
의 가동 범위를, 각각
[수학식 557]
Figure 112019032179095-pct00557
(또는
[수학식 558]
Figure 112019032179095-pct00558
)로 하는 최대 우도 추정으로 타깃 검색을 행한다. 다른 부분 영역에서의 타깃 검색도 마찬가지이다. 이는, code division multiple access의 고찰을 따라 2차원 PC로 다중 타깃 검색을 행하므로, Code Division Multiple Target(CDMT)이라고 부른다.
N=N'=64, Npath=4의 수치 시뮬레이션 결과에서는, PUL 알고리즘을 적용하면 SNR 5 dB 이상에서 80%의 확률로 3, 4개의 타깃 검출에 성공하고 있다.
<8.2 인공적인 비가환 시프트를 매장하는 delay-Doppler space Division multiplexing에 기초하는 초고 MPSK>
PUL의 수렴 증명에서 중요한 역할을 한, 힐베르트 공간의 부분 공간, Ts-TL TD 공간, 및 Fs-BL FD 공간 각각의 위의 직교 투영 연산자 P3, P4에서 이용하는 식 (38), 식 (43)의 TD-, FD-CCF
[수학식 559]
Figure 112019032179095-pct00559
는, 비가환 통신 유래의 각종 TFSO 중에서 송·수신 신호 생성의 TFSO 쌍
[수학식 560]
Figure 112019032179095-pct00560
(또는,
[수학식 561]
Figure 112019032179095-pct00561
)이 본질인 것을 나타내고 있다. 왜냐하면, 이들에는 미지수
[수학식 562]
Figure 112019032179095-pct00562
와, 최대 우도값
[수학식 563]
Figure 112019032179095-pct00563
의 탐색용 제어 파라미터
[수학식 564]
Figure 112019032179095-pct00564
의 어느 것인가가 포함되어 있기 때문이다. 변·복조의 TFSO
[수학식 565]
Figure 112019032179095-pct00565
나 데이터, 칩 어드레스
[수학식 566]
Figure 112019032179095-pct00566
을 파라미터로 하는 TFSO
[수학식 567]
Figure 112019032179095-pct00567
은, 이들의 어느 것과도 관계가 없다.
한편,
[수학식 568]
Figure 112019032179095-pct00568
로 지정되는, TD-, FD-CCF 쌍
[수학식 569]
Figure 112019032179095-pct00569
의 출력(식 (41), (45))에는, 각종 PD를 포함하는 형태로 정보 데이터
[수학식 570]
Figure 112019032179095-pct00570
이 양(陽)으로 나타난다. 이것은,
[수학식 571]
Figure 112019032179095-pct00571
의 위상항의 인공적인 파라미터로서 파라미터값
[수학식 572]
Figure 112019032179095-pct00572
(간단히 하기 위해 감쇠 상수는
[수학식 573]
Figure 112019032179095-pct00573
이라고 한다)가, 상기 복수 개 타깃 검출법을 이용하면 이용 가능한 것을 나타낸다. 이를 감안하여, 본 실시 형태에 따른 송수신 시스템은, 구체적으로, 송·수신 신호의 Npath개의 TFSO 쌍 및 TD-, FD-CCF 쌍
[수학식 574]
Figure 112019032179095-pct00574
를 이용하는 비가환 통신이며, 파라미터값
[수학식 575]
Figure 112019032179095-pct00575
를 매장한 송신 신호의 Npath 쌍의 TD-, FD-template을 이용한다. 이 경우의 PUL을 송신 신호에 인공적으로 매장한 비가환 시프트량
[수학식 576]
Figure 112019032179095-pct00576
을 수신기에서 복원한다는 의미에서, 통상의 PUL와 조금 다르다. 한편, 비특허 문헌 [25, 30]에서는, transmitter, receiver간의 파라미터 갱신을 「active PUL」라고 불렀지만, PUL은 수렴 증명의 문제상 transmitter에서의 갱신에는 적용할 수 없기 때문에, 본 명세서에서는 transmitter에 기존의 비가환 시프트량
[수학식 577]
Figure 112019032179095-pct00577
을 매장하고, 수신기(reciever)에서 그 시프트량을 이용해 최대 우도 추정을 행한다. 이는 통상적인 PUL의 응용예이다.
복수 개의 서로 독립적인 2차원 PC를 이용하는 CDMT는,
[수학식 578]
Figure 112019032179095-pct00578
phase shift keying(PSK) 데이터
[수학식 579]
Figure 112019032179095-pct00579
의 데이터 통신에 효과적이다. N', N개의 어레이형 TD-, FD-CCF의 입력 부분에 각각 위상 보상항
[수학식 580]
Figure 112019032179095-pct00580
을 전치한, phase tuned layer(PTL)는 TD-, FD-CCF 쌍
[수학식 581]
Figure 112019032179095-pct00581
의 출력(식 (41), (45))의
[수학식 582]
Figure 112019032179095-pct00582
를 각각
[수학식 583]
Figure 112019032179095-pct00583
으로 치환하는 위상 보상과, 식 (60)의 PUL의 최대 우도 추정의 연산
[수학식 584]
Figure 112019032179095-pct00584
의 2변수를 대신해, PTL용 변수
[수학식 585]
Figure 112019032179095-pct00585
를 추가한 3변수의
[수학식 586]
Figure 112019032179095-pct00586
으로 수정하는 것에 상당한다. 이에 따라, 최대 우도 추정의 복소출력의 위상 요동은 평균 0(零)의 가우스 분포를 닮은 분포로 바뀐다.
본 실시 형태에 따른 송수신 시스템은, N=N'=16,
[수학식 587]
Figure 112019032179095-pct00587
SNR 30 dB에서 양호한 수치 시뮬레이션 결과를 나타내고 있다. 그러나,
[수학식 588]
Figure 112019032179095-pct00588
의 경우, 분포의 메인로브(main lobe)의 좌우 양측에 사이드로브가 생기므로 복호 에러의 원인이 된다. 이 때문에, 단순한 위상 보상을 이용하는 한, 복호 에러의 관점에서는,
[수학식 589]
Figure 112019032179095-pct00589
까지로 하는 것이 바람직하다.
계속해서, 본 실시 형태에 따른 송수신 시스템으로서,
[수학식 590]
Figure 112019032179095-pct00590
의 데이터 통신을 행하기 위해 delay-Doppler의 파라미터 공간 Θ'를
[수학식 591]
Figure 112019032179095-pct00591
개로 균등 분할해, 각각의 부분 공간에 2차원 PC
[수학식 592]
Figure 112019032179095-pct00592
를 할당하고, 또한 인공적인 시프트량
[수학식 593]
Figure 112019032179095-pct00593
이 그 i번째 부분 공간의 중심점 근방에 있다고 간주함으로써, 8-PSK를 기본 단위로 하는 TD-, FD-CCF의 PTL의 위상 보상항을 각각
[수학식 594]
Figure 112019032179095-pct00594
로 하는, delay-Doppler space division muliplex(dDSDM)를 실현하는 비가환 통신을 예시한다. 이 경우,
[수학식 595]
Figure 112019032179095-pct00595
에 대한 PTL 및 그 최대 우도 추정의 4변수 연산은 각각
[수학식 596]
Figure 112019032179095-pct00596
이 된다. 따라서, 본 실시 형태에 따른 송수신 시스템은, 합계
[수학식 597]
Figure 112019032179095-pct00597
개의 CCF를 구비하고 있다. 16, 32개의 2차원 PC를 이용한 128 PSK, 256 PSK의 복호 결과가 얻어지고 있다. 그러나, 이것은 본질적으로 위상 WM의 변별 정도를 강요하는 것으로, 초고(超高)의
[수학식 598]
Figure 112019032179095-pct00598
-PSK의 부호화·복호화의 개선으로 연결되지 않는다.
이하, 동기를 하면서, 초고의
[수학식 599]
Figure 112019032179095-pct00599
-PSK의 부호화·복호화하는 시스템, 및, 동기·측거 겸용의
[수학식 600]
Figure 112019032179095-pct00600
-ary PSK-부호화·복호화기에 대해, 도 12∼도 17을 참조해 설명한다.
우선, 도 12는,
[수학식 601]
Figure 112019032179095-pct00601
-PSK 통신 가능/고속·고정도 측거 가능한 송신기(부호화기)의 구성을 나타내는 블록도이다.
동기·측거의 경우(도 12의 Switch 1-1, 1-2가 상측에 접속되어 있는 상태)와
[수학식 602]
Figure 112019032179095-pct00602
-PSK의 데이터 통신을 행하는 경우(도 12의 Switch 1-1, 1-2가 하측에 접속되어 있는 상태)에서, 송신기 및 수신부에서 각각 절환을 실시한다. 입력 데이터
[수학식 603]
Figure 112019032179095-pct00603
[수학식 604]
Figure 112019032179095-pct00604
-PSK의 경우(환언하면, '데이터 통신유'의 경우), 각 도면에 나타내는 스위치가 아래 쪽의 경로로 절환된다.
이하에서는, 우선 도 12의 송신기의 블록도에서 각 블록에 대해 설명한다. 도 12에서, 좌단의 입력은,
[수학식 605]
Figure 112019032179095-pct00605
이고, 도 12에 나타내는 Switch 1-1 및 Switch 1-2가 상측에 접속된 상태의 동기·측거 모드에서는, 송신기는
[수학식 606]
Figure 112019032179095-pct00606
을 행한다. 한편, 도 12에 나타내는 Switch 1-1 및 Switch 1-2가 하측에 접속된 상태(환언하면, '데이터 통신유'의 상태)에서는, 송신기는 chip 파형을 입력으로 하여
[수학식 607]
Figure 112019032179095-pct00607
을 행한다. 다음으로, 「k의 부호화기」 블록에서는 0th-AC부터,
[수학식 608]
Figure 112019032179095-pct00608
th-AC가 병렬로 배치되고, 그 중간에는 복수의 검은 원이
[수학식 609]
Figure 112019032179095-pct00609
와 같이 세로로 배치되어 있다. 여기에서, 각 검은 원은 각 AC를 생략적으로 표현한 것에 불과하다.
본 실시 형태에 따른 부호화기는, 「k의 부호화기」 블록에서, j의 값에 따라, 「k의 부호화기」 블록에서의 좌우 양측의 Switch 1-3 및 Switch 1-4를 절환하고, 그 후, 「k의 부호화기」 블록에서의 하류측에서
[수학식 610]
Figure 112019032179095-pct00610
의 변조를 행해,
[수학식 611]
Figure 112019032179095-pct00611
로 부호화한다.
그 후, 송신기가 구비하는 Switch 1-2의 상하에는 각각 송신
[수학식 612]
Figure 112019032179095-pct00612
, 각각 송신
[수학식 613]
Figure 112019032179095-pct00613
이 얻어진다. 본 실시 형태에 따른 부호화기에서는, 또한
[수학식 614]
Figure 112019032179095-pct00614
변조 후,
[수학식 615]
Figure 112019032179095-pct00615
-MC를 거쳐 노이즈가 더해지고,
[수학식 616]
Figure 112019032179095-pct00616
복조를 거쳐, 수신 CE가 얻어진다.
이상이 도 12에 나타낸 각 블록의 설명이다.
고차(高次)의
[수학식 617]
Figure 112019032179095-pct00617
-PSK 데이터
[수학식 618]
Figure 112019032179095-pct00618
전송은 위상 잡음의 영향으로 위상
[수학식 619]
Figure 112019032179095-pct00619
의 변별이 곤란하므로, 우선 k의 부호화기(도 12에서의 하단 오른쪽의 실선으로 둘러싸인 부분)의 중간 블록에서
[수학식 620]
Figure 112019032179095-pct00620
으로 부호화하고, 저차(低次)의
[수학식 621]
Figure 112019032179095-pct00621
-PSK 데이터
[수학식 622]
Figure 112019032179095-pct00622
전송을 실시하는 것을 생각한다. 한편,
[수학식 623]
Figure 112019032179095-pct00623
의 전송은, 시간 시프트, 주파수 시프트가 정확하게 복원 가능한 본 특허 제안의 동기·측거법을 원용한다. 이를 위해 시간 지연·도플러 시프트의 파라미터 평면 Θ를
[수학식 624]
Figure 112019032179095-pct00624
개로 균등 배반 분할한 부분 평면 Θ(i)에 2차원 PC χ(i)를 할당한다(도 17 참조). chip 파형을
[수학식 625]
Figure 112019032179095-pct00625
에 의한 2차원 BPSK 변조로 이들을
[수학식 626]
Figure 112019032179095-pct00626
다중화하고, 또한, 이 신호가 Θ(i)의 중심점에 상당하는 시간 시프트, 주파수 시프트의 인공적 통신로(Artificial Channel: AC)를 통과했다고 상정한다. 이를 위해, 0th AC부터
[수학식 627]
Figure 112019032179095-pct00627
th AC 중에서 j에 부합하는 AC를 선택하고, 그 출력 신호에서
[수학식 628]
Figure 112019032179095-pct00628
[수학식 629]
Figure 112019032179095-pct00629
-PSK 변조한다. 이것이 송신 CE가 되어,
[수학식 630]
Figure 112019032179095-pct00630
변조를 거쳐 전파로(Main Channel: MC)의
[수학식 631]
Figure 112019032179095-pct00631
-MC로 송출된다.
[수학식 632]
Figure 112019032179095-pct00632
복조를 거쳐 노이즈가 인가되고 수신 CE가 된다.
계속해서, 도 13에 나타내는 각 블록을 설명한다.
도 13은,
[수학식 633]
Figure 112019032179095-pct00633
-PSK 통신 가능하고, 고속·고정도 측거 가능한 수신 동기기(복호화기)의 구성을 나타내는 블록도이다.
도 13에 나타낸 Switch 2-1이 상측에 접속된 상태에서는, 수신기는
[수학식 634]
Figure 112019032179095-pct00634
를 행한다. 한편, Switch 2-1이 하측에 접속된 상태(환언하면, '데이터 통신유' 상태)에서는, 우선, 「k의 복호화기」 블록에서
[수학식 635]
Figure 112019032179095-pct00635
으로 복호화해,
[수학식 636]
Figure 112019032179095-pct00636
된 신호를 좌측의 블록에 입력한다. 그 블록은, 0th-AC,
[수학식 637]
Figure 112019032179095-pct00637
th-AC AC가 병렬로 배치되고, 그 중간에는 복수의 검은 원이
[수학식 638]
Figure 112019032179095-pct00638
과 같이 세로로 배치되어 있다. 여기에서, 각 검은 원은 각 AC를 생략적으로 표현한 것에 불과하다.
본 실시 형태에 따른 부호화기는, 「k의 복호화기」 블록에서
[수학식 639]
Figure 112019032179095-pct00639
의 값에 따라 좌우 양측의 Switch 2-3 및 Switch 2-4를 절환한 후, 좌측에서
[수학식 640]
Figure 112019032179095-pct00640
복조를 행해, 하단의 블록 COR에 좌단의 수신 CE와 함께 입력된다.
단, 도 13에서의 「COR」은 상관기(correlator)를 가리킨다.
도 13에서의 「COR」 우측의 큰 블록은 (ρ', ρ) 선택에 의한
[수학식 641]
Figure 112019032179095-pct00641
, Ae-최대 우도 추정 및
[수학식 642]
Figure 112019032179095-pct00642
복원, k의 복호를 행하는 블록이다. Switch 2-2가 상측에 접속된 상태에서 수신기는
[수학식 643]
Figure 112019032179095-pct00643
의 처리를 행하고, Switch 2-2가 하측에 접속된 상태(환언하면, '데이터 통신유' 상태)에서 수신기는
[수학식 644]
Figure 112019032179095-pct00644
의 처리를 행한다.
이상이 도 13에 나타낸 각 블록의 설명이다.
수신 CE를, chip 파형의 2차원 PCχ에 의한 2차원 BPSK 복조를 실시하는 경우를 생각한다. 수신기에서 Switch 2-1이 하측에 접속된 상태의 '데이터 통신유'의 모드에서는, chip 파형을 j의 추정값
[수학식 645]
Figure 112019032179095-pct00645
의 2차원 PC
[수학식 646]
Figure 112019032179095-pct00646
에 의한 2차원 BPSK 복조해 얻어지는 신호를, 0th-AC부터
[수학식 647]
Figure 112019032179095-pct00647
th-AC 중에서 선택한
[수학식 648]
Figure 112019032179095-pct00648
th-AC에 입력하고, 그 출력 신호로
[수학식 649]
Figure 112019032179095-pct00649
[수학식 650]
Figure 112019032179095-pct00650
-PSK 복조한다. 마지막 단의 파라미터 추정 블록(즉 「COR」의 출력이 입력되는 블록)에서는
[수학식 651]
Figure 112019032179095-pct00651
에 대해 각각 4변수
[수학식 652]
Figure 112019032179095-pct00652
에 관해 연산 argmax를 실시해, PUL에 의한
[수학식 653]
Figure 112019032179095-pct00653
-최대 우도 추정(PUL의 수렴값은
[수학식 654]
Figure 112019032179095-pct00654
)을 행한다.
여기에서, 최대 우도 추정은 Switch 2-2가 상측에 접속되었을 경우와 하측에 접속되었을 경우에서 공통이지만, Switch 2-2가 상측에 접속된 상태에서는, 복호기는
[수학식 655]
Figure 112019032179095-pct00655
의 복원을 행하고, Switch 2-2가 하측에 접속된 상태인 '데이터 통신유' 모드에서는
[수학식 656]
Figure 112019032179095-pct00656
-최대 우도 추정값에 의한 k의 복호
[수학식 657]
Figure 112019032179095-pct00657
을 행한다. 한편, 통상의 동기·측거 모드(상측의 switch 상태)에서도 데이터
[수학식 658]
Figure 112019032179095-pct00658
은 저차의
[수학식 659]
Figure 112019032179095-pct00659
-PSK 변·복조는 가능(예를 들면,
[수학식 660]
Figure 112019032179095-pct00660
)하다.
계속해서, 도 14는, Main channel(MC)의 delay τ-Doppler shift ν 공간상에서 CCF 실부의 크기의 분포예를 모식적으로 나타내는 도면이다.
동기·측거 모드에서 메인로브는
[수학식 661]
Figure 112019032179095-pct00661
에 위치해, 다른 사이드로브와 구별할 수 있다. 한편, 물결선은 백그라운드 노이즈를 나타낸다.
도 15는, AC가 Main channel(MC)에 중첩된 경우의 τ-ν 공간상에서 CCF 실부의 크기의 분포예를 모식적으로 나타내는 도면이다.
MC의 앞에
[수학식 662]
Figure 112019032179095-pct00662
-PSK 신호의 정보 k의 부호화로 도입된 3종의 Artificial channel(AC)이 종속 접속된 경우의 CCF 실부 크기의 분포를 나타낸다. MC는
[수학식 663]
Figure 112019032179095-pct00663
에 위치한 채로 있지만, AC가 중첩된 경우 MC+AC0, MC+AC1, MC+AC2는 각각
[수학식 664]
Figure 112019032179095-pct00664
에 위치한다. 시간 시프트, 주파수 시프트의 양은 가법적(加法的)이 되지만, PD는 군론적 성질을 수반해 전파한다. 도 15와 같은 분포도를 얻기 위해서는, 본 명세서에서 분명히 한 바와 같이, PD의 엄밀 평가와 함께 CCF 실부의 극대화를 행하는 최대 우도 추정으로 CCF 실부의 극대치를 제공하는 변수의 수를 증가시켜 정확하게 추측할 필요가 있다. 한편, 추가한 변수 j, j'의 종류수는 각각
[수학식 665]
Figure 112019032179095-pct00665
이다.
도 16은, 여차원 2의, 비가환성을 갖는 AC 시프트의 파라미터 공간을 이용한 신호의 TFP 분할을 나타내는 도면이다(파라미터 공간이므로 「여차원」이라는 표현을 이용하고 있다).
도 16에서는, 신호(시간폭 Ts, 대역폭 Fs)의 시간·주파수 평면(TFP) S의 분할(가보 분할)(S(0), S(1), S(2), S(3)) 평면과 직교하는 축에 AC의 비가환 시프트의 양
[수학식 666]
Figure 112019032179095-pct00666
을 나타내는 눈금을 나타내고 있다.
도 17도, 여차원 2의, 비가환성을 갖는 AC 시프트의 파라미터 공간을 이용한 신호의 TFP 분할을 나타내는 도면이다.
도 17과 같이, AC의 비가환 시프트의 양에 따라 TFP는 각각 AC0의 TFP, AC1의 TFP, AC2의 TFP, AC3의 TFP와 같이 시프트된다. 한편, 신호의 S는
[수학식 667]
Figure 112019032179095-pct00667
의 MC의 target space의 단위 평면
[수학식 668]
Figure 112019032179095-pct00668
과 동일시할 수 있다. 즉, 타깃이 데이터 어드레스
[수학식 669]
Figure 112019032179095-pct00669
근방의 target space의 부분 평면[(p-1)Ts, pTs)×[(p'-1)Fs, p'Fs)에 존재하는 경우, 어드레스
[수학식 670]
Figure 112019032179095-pct00670
의 특정이 필요하므로, 도 13의 2종의 CCF 실부의 최대화 변수로서 모든 경우에서(데이터 통신 유무에 관계없이)
[수학식 671]
Figure 112019032179095-pct00671
을 채택하고 있다.
전술한 바와 같이, 도 12 및 도 13의 송수신기에서 상측 및 하측의 시스템이 절환된다.
도 12 및 도 13에서 Switch 1-1, 1-2, 2-1, 2-2가 상측에 접속된 시스템이 특허 문헌 6에 개시된 시스템이며, 이들 스위치가 하측에 접속된 시스템이 비가환 AC-shift-부호화·복호화에 기초하는 다중화 방식이다.
이와 같이, 도 12 및 도 13은, 종래의 시스템과의 차이, 신규성이 분명한 도면이다.
본 실시 형태에 따른 송수신 시스템에서는, 다치 위상 변조(
[수학식 672]
Figure 112019032179095-pct00672
-ary PSK)
[수학식 673]
Figure 112019032179095-pct00673
의 데이터 신호 전송에서는, 상기 송신 신호에 「정보 k」를 효율적으로 매장하기 위해, 정수
[수학식 674]
Figure 112019032179095-pct00674
에 대해, 2차원 위상 편이 변조(BPSK)용의, 주기 N의 TD 위상 변조 부호(TD-PC)
[수학식 675]
Figure 112019032179095-pct00675
와 주기 N'의 FD 위상 변조 부호(FD-PC)
[수학식 676]
Figure 112019032179095-pct00676
의 쌍을
[수학식 677]
Figure 112019032179095-pct00677
세트 준비한다. 시간 지연, 도플러 시프트의 파라미터 평면
[수학식 678]
Figure 112019032179095-pct00678
및, 신호의 TFP S=[0, Ts)×[0, Fs)를
[수학식 679]
Figure 112019032179095-pct00679
개 균등 분할한 부분 평면을 각각
[수학식 680]
Figure 112019032179095-pct00680
으로 한다.(Ts=NMΔt, Fs= N'M'Δf는 데이터 심볼
[수학식 681]
Figure 112019032179095-pct00681
의 시간, 대역폭, Tmax, Fmax는 검출할 시간 지연, 도플러 시프트의 최대값, Ts, Fs의 정수배 P, P'로 한다. Δt, Δf는 양자화폭.) 다음으로 정보 k를
[수학식 682]
Figure 112019032179095-pct00682
로 분해해,
[수학식 683]
Figure 112019032179095-pct00683
과 그 나머지
[수학식 684]
Figure 112019032179095-pct00684
로 부호화한다. 부호(j, j') 전송을 위해, 첫번째로 시간 Tc= MΔt, 대역폭 Fc= M'Δf의 칩 펄스를 2차원 PC X(i)로 변조한 TD-signature, FD-signature를 (Ts, Fs) 레벨의 시프트를 행해
[수학식 685]
Figure 112019032179095-pct00685
-signature-다중화한다. 두번째로 그 다중화 신호가 부분 파라미터 평면 Θ(j)의 중심:
[수학식 686]
Figure 112019032179095-pct00686
의 시프트를 갖는 인공적 전파로(Artificial Channel: AC)를 통과했다고 상정한다. 이것을 시뮬레이션하기 위해, TD-, FD-signature 신호에 제4 시프트 연산자
[수학식 687]
Figure 112019032179095-pct00687
을 작용시킨다. 세번째로, 얻어진 신호로
[수학식 688]
Figure 112019032179095-pct00688
-ary 데이터 신호
[수학식 689]
Figure 112019032179095-pct00689
를 변조한다.
[수학식 690]
Figure 112019032179095-pct00690
이 신호를 TFP상에서 심볼
[수학식 691]
Figure 112019032179095-pct00691
의 (P/P')의 무중첩 중합 신호를 송신 TD-, FD-CE로 한다(도 12 참조).
[수학식 692]
Figure 112019032179095-pct00692
-ary 데이터 통신의 TD-신호의 CE와 그 DFT는, 식 (27) 대신에
[수학식 693]
Figure 112019032179095-pct00693
으로 정해진다. 단,
[수학식 694]
Figure 112019032179095-pct00694
는 TD-signature와 FD-signature
[수학식 695]
Figure 112019032179095-pct00695
이다(식 (25) 참조). 여기에서,
[수학식 696]
Figure 112019032179095-pct00696
은, 도 16 및 도 17에 나타내는 부분 평면 S(i)에 부수된 칩 어드레스(m, m')의 가동 범위(S의 어드레스의 대응 부분)를 나타낸다.
[수학식 697]
Figure 112019032179095-pct00697
이고, 동기법이나 측거시의 N, N'에 비해, 데이터 통신에서는 신호의 시간·주파수 평면 및 타깃 평면을 PC 부호 분할하므로, 정도(精度)를 얻기 위해서
[수학식 698]
Figure 112019032179095-pct00698
배 필요하다(도 12 및 도 13을 참조('데이터 통신유'의 스위치 상하는, 동기법이나 측거시, 데이터 통신에 대응한다)).
식 (87)에서는 i번째의 2차원 PC X(i)로 변조한 signatureν[k;χ(i)](또는
[수학식 699]
Figure 112019032179095-pct00699
)를
[수학식 700]
Figure 112019032179095-pct00700
-(χ(i)&signature)-다중화한 신호가
[수학식 701]
Figure 112019032179095-pct00701
-ary 데이터
[수학식 702]
Figure 112019032179095-pct00702
의 k의 부호(j, j') 대응의 시프트
[수학식 703]
Figure 112019032179095-pct00703
의 AC를 통과하고, 이것에
[수학식 704]
Figure 112019032179095-pct00704
-ary 데이터
[수학식 705]
Figure 112019032179095-pct00705
의 정보를 부과한 CE
[수학식 706]
Figure 112019032179095-pct00706
(또는 DFT
[수학식 707]
Figure 112019032179095-pct00707
)이 새로운 송신 TD-CE(또는 FD-CE)이다. 단,
[수학식 708]
Figure 112019032179095-pct00708
은 크기 1의 실수이다(도 12의 송신기 참조).
수신 시스템에서는, k의 복호를 위해 추정값
[수학식 709]
Figure 112019032179095-pct00709
대응의 추정 부호
[수학식 710]
Figure 112019032179095-pct00710
과 대응의 시프트
[수학식 711]
Figure 112019032179095-pct00711
의 AC를 시뮬레이션한, 시프트 연산자
[수학식 712]
Figure 112019032179095-pct00712
를 작용시켜, 다시 추정
[수학식 713]
Figure 112019032179095-pct00713
의 PD 보상을 행한다(도 13의 수신기 참조). 구체적으로는 AC의 시프트를 이용한, type-3, type-4의 상관기는, 식 (38), (43) 대신에, 2변수 j, j'와 그에 따른 최대화 변수
[수학식 714]
Figure 112019032179095-pct00714
가 추가되고,
[수학식 715]
Figure 112019032179095-pct00715
상단 식 2종의 상관 함수에서는, 식 (87)의 송신 CE 신호에서 연산자
[수학식 716]
Figure 112019032179095-pct00716
이 파라미터 추정용 2차원 PC 연산자
[수학식 717]
Figure 112019032179095-pct00717
[수학식 718]
Figure 112019032179095-pct00718
의 전파로(Main Channel:MC)의 연산자
[수학식 719]
Figure 112019032179095-pct00719
의 사이에 삽입된다(도 12 참조). 한편, 수신측에서는 식 (89), (90)의
[수학식 720]
Figure 112019032179095-pct00720
-ary PSK를 위한
[수학식 721]
Figure 112019032179095-pct00721
종의
[수학식 722]
Figure 112019032179095-pct00722
의 AC의 추정 AC 연산자
[수학식 723]
Figure 112019032179095-pct00723
(또는
[수학식 724]
Figure 112019032179095-pct00724
)이
[수학식 725]
Figure 112019032179095-pct00725
의 추정 MC 연산자
[수학식 726]
Figure 112019032179095-pct00726
(또는
[수학식 727]
Figure 112019032179095-pct00727
)과 추정 2차원 PC 연산자
[수학식 728]
Figure 112019032179095-pct00728
의 사이에 삽입된다(도 13 참조).
도 12 및 도 13의
[수학식 729]
Figure 112019032179095-pct00729
-PSK 통신 가능/고속·고정도 측거 가능한 송수신기와 부호 복호화기는, 시간 지연·도플러 시프트의 여차원 2의 파라미터 공간 Θ을 배반 분할한 delay-Doppler Space Division Multiplexing(dD-SDM)(비특허 문헌 30, 특허 문헌 6)의 실현예이다. 이는
[수학식 730]
Figure 112019032179095-pct00730
-PSK의 데이터 k 대신에 부호
[수학식 731]
Figure 112019032179095-pct00731
을 전송하는 것으로 하여, 신호가 Θ(j)의 중심의 시간 시프트
[수학식 732]
Figure 112019032179095-pct00732
, 주파수 시프트
[수학식 733]
Figure 112019032179095-pct00733
의 AC를 통과함으로써 생기는 PD 검출(도 14, 15 참조)로 j, j'를 복호한다. 수신측에서 j를 검출하기 위해서는 송신측에서 상이한 i의 PC χ(i)로 2차원 BPSK하고, 이들을
[수학식 734]
Figure 112019032179095-pct00734
다중화해, 이것을
[수학식 735]
Figure 112019032179095-pct00735
th-AC에 입력하고, 그 출력 신호로
[수학식 736]
Figure 112019032179095-pct00736
[수학식 737]
Figure 112019032179095-pct00737
-PSK 복조를 행한다.
[수학식 738]
Figure 112019032179095-pct00738
[수학식 739]
Figure 112019032179095-pct00739
종 있으므로
[수학식 740]
Figure 112019032179095-pct00740
개의 AC와 2차원 BPSK의 PC
[수학식 741]
Figure 112019032179095-pct00741
이 필요하다. 수신측에서는 chip 파형의 TFP상의 4중 다중화로 매장된 template의 검출을 행한다. 추정
[수학식 742]
Figure 112019032179095-pct00742
th-AC의 선택으로,
[수학식 743]
Figure 112019032179095-pct00743
의 최대 우도 추정을 행한다.
[수학식 744]
Figure 112019032179095-pct00744
-PSK 통신을
[수학식 745]
Figure 112019032179095-pct00745
종의 비가환 시프트의 AC 삽입에 의한 다중화로 고차(高次)의
[수학식 746]
Figure 112019032179095-pct00746
-PSK 통신을 실현하고 있다.
이 dD-SDM은, 여차원 2의, 비가환 시프트의 파라미터 공간을 이용한 다중화 통신 방식이다. 이것은, [0, Ts)×[0, Ts)의 심볼 평면의 TFP(타깃 평면의 단위 부분 평면과 동일시)와의 직교축이 AC의 비가환 시프트축인, 3차원화(AC의 시프트로 계층화)된 시간·주파수 공간(Time-Frequency Space: TFS)을 가지므로(도 16, 도 17 참조), 종래의 신호 TFP 분할의 TDMA, FDMA와는 다른 방식이다. 파라미터의 고속·고정도 추정을 보증하는 APT-based의 PUL은 근간 기술이다. 한편, 비가환성 시프트를 갖는 AC는 각종 전파 선로로 대체 가능하므로, 초고차(超高次)의 PSK 통신이 요구되는 광통신에의 응용이 기대된다.
비특허 문헌 30이나 특허 문헌 6에서는,
[수학식 747]
Figure 112019032179095-pct00747
의 다중 타깃 검출 문제를 CDMT(Code-Division Multiple Targets: CDMT)(비특허 문헌 27)로서 파악해, Npath개의 2차원 PC와
[수학식 748]
Figure 112019032179095-pct00748
의 시프트 연산과 2개의 위상 회전 인자의 곱
[수학식 749]
Figure 112019032179095-pct00749
에 의한 위상 보상으로 타깃 검출을 행하고 있다. 이것은 실질적으로 미소한 위상량
[수학식 750]
Figure 112019032179095-pct00750
의 변별을 강요하는 방법으로, 큰
[수학식 751]
Figure 112019032179095-pct00751
에 대해 위상 잡음에 의한 복호 오차가 생겨 한계가 있었다. 또한,
[수학식 752]
Figure 112019032179095-pct00752
-MC의 연산자
[수학식 753]
Figure 112019032179095-pct00753
이나 그 보상의 시프트
[수학식 754]
Figure 112019032179095-pct00754
를 이용하지 않고,
[수학식 755]
Figure 112019032179095-pct00755
-ary의 복호를 행하고 있으므로, MC를 통과한 2중의 시프트된 신호의 다중 타깃 검출 문제가 아니라 파라미터 추정의 수렴의 보장도 없었다. 이것을 해결하기 위해, 본 명세서에 따른 발명에서는 새롭게 2종의 상관 함수
[수학식 756]
Figure 112019032179095-pct00756
을 정의하고, 각종 시프트 연산자에 기초한 정확한 PD 평가와 그 보상을 행하였다.
식 (89), (90)의 2종의 상관 함수의 식 (87)의
[수학식 757]
Figure 112019032179095-pct00757
(또는
[수학식 758]
Figure 112019032179095-pct00758
)에 포함되는
[수학식 759]
Figure 112019032179095-pct00759
를 복호하기 위한 수신부(도 12 참조)에 매장한
[수학식 760]
Figure 112019032179095-pct00760
에 의한, PD 보상의 계산은 매우 번잡하다. 그러나, 상세한 계산 과정은 생략하지만, 고차의
[수학식 761]
Figure 112019032179095-pct00761
-ary PSK의 복호를 위한 상관 함수는,
[수학식 762]
Figure 112019032179095-pct00762
종의 2차원 PC
[수학식 763]
Figure 112019032179095-pct00763
을 이용함으로써 이하의 변수의 치환:
[수학식 764]
Figure 112019032179095-pct00764
(수신측에서의 치환:
[수학식 765]
Figure 112019032179095-pct00765
에 수반해, 다수의 위상항이 발생하고, 또한 상관 함수 실부를 최대화하는 최대 우도 추정 파라미터가 2변수
[수학식 766]
Figure 112019032179095-pct00766
(또는
[수학식 767]
Figure 112019032179095-pct00767
)으로부터 4변수
[수학식 768]
Figure 112019032179095-pct00768
(또는
[수학식 769]
Figure 112019032179095-pct00769
)로 확장되고(도 13 참조. 단,
[수학식 770]
Figure 112019032179095-pct00770
(또는
[수학식 771]
Figure 112019032179095-pct00771
), 식 (39), (44)를 일반화한 형식으로 각각
[수학식 772]
Figure 112019032179095-pct00772
단,
[수학식 773]
Figure 112019032179095-pct00773
이다. 또한
[수학식 774]
Figure 112019032179095-pct00774
는 신호 부분 평면
[수학식 775]
Figure 112019032179095-pct00775
에 부수된 칩 어드레스 n, n'의 가동 범위,
[수학식 776]
Figure 112019032179095-pct00776
이다. 식 (91), (92)의 ambiguity function(AF) θgg[·],ΘGG[·]는, g, G가 가우스 함수의 경우 지수 함수적 감소 함수인 것, 및, 서로 배반인 칩 어드레스 분할
[수학식 777]
Figure 112019032179095-pct00777
이나
[수학식 778]
Figure 112019032179095-pct00778
이 Θ(j)의 중심인 것으로부터, AF의 제1 변수, 제2 변수가 큰 데이터 어드레스 (p-q)MN1, (p'-q')M'N1'의 p'≠q', p≠q의 항이나 칩 어드레스
[수학식 779]
Figure 112019032179095-pct00779
의 항은 무시할 수 있어, p'=q', p=q나
[수학식 780]
Figure 112019032179095-pct00780
의 항만이 남아,
[수학식 781]
Figure 112019032179095-pct00781
중의
[수학식 782]
Figure 112019032179095-pct00782
가 추출되고
[수학식 783]
Figure 112019032179095-pct00783
이 특정된다. 이것들이 Θ를 배반 분할해 2차원 PC χ(i)로 변조한 신호를
[수학식 784]
Figure 112019032179095-pct00784
-다중화해, 굳이 PD 발생원의 AC를
[수학식 785]
Figure 112019032179095-pct00785
개 병렬화하고, 신호를 그 하나의 AC에 통과시키는 최대의 이유이다. 한편,
[수학식 786]
Figure 112019032179095-pct00786
의 특정은
[수학식 787]
Figure 112019032179095-pct00787
-ary의 위상 보상
[수학식 788]
Figure 112019032179095-pct00788
으로 달성된다(도 12 및 도 13 참조).
상기 2종의 상관 함수는 서로 독립적인 2차원 PC
[수학식 789]
Figure 112019032179095-pct00789
를 이용해 k의 부호
[수학식 790]
Figure 112019032179095-pct00790
과 그 추정 부호
[수학식 791]
Figure 112019032179095-pct00791
-의존의 비가환 시프트의
[수학식 792]
Figure 112019032179095-pct00792
-AC의 PD를 이용한, PUL에 기초하는 최대 우도 추정으로, 위상의 최소 변별을
[수학식 793]
Figure 112019032179095-pct00793
으로 하는, 위상 잡음에 내성을 가진 복호법이다.
[수학식 794]
Figure 112019032179095-pct00794
-ary PSK의 데이터 통신을 행하기 전에 동기(동기 포착, 동기 유지)나 측거를 확립해 두는 것은 효과적이다. 그러나, 현실의 통신에서는
[수학식 795]
Figure 112019032179095-pct00795
의 MC를 개재한 통신을 전제로 하는, 상기 식의 2개의 상관 함수가 필요하다. 두 함수는 도 12 및 도 13의 상측 스위치 상태에서 동기 포착·동기 유지를 행하고, 도 12 및 도 13의 하측 스위치 상태에서
[수학식 796]
Figure 112019032179095-pct00796
-ary PSK의 부호화·복호화를 병행적으로 가능하게 한다. 즉, 상기 2종의 상관 함수는, 레이더에 대해서는 시간 지연·도플러 시프트의 측거기로서 동작하고,
[수학식 797]
Figure 112019032179095-pct00797
-ary 데이터 무선 통신(도 12 및 도 13의 하측 스위치 상태)에 대해서는 시간 지연·도플러 시프트의 동기기 및 부호·복호기로서 동작하는, 진폭 변조 불요 다치 위상 변조 복조 시스템이다. 이것은 2차원 PC
[수학식 798]
Figure 112019032179095-pct00798
을 비밀키로 하는 무선의 비밀 통신 시스템의 실현예, 혹은, 비밀 데이터 통신 및 측거가 가능한 차량 탑재 레이더·데이터 통신 시스템을 제공한다.
종래법과 달리, 진폭에 「정보」를 부과하지 않는 이유는, (1) 통신로의 감쇠 상수 Ae의 수시 추정, (2) 복호기의 단순화 등에 있다. 정확한 PD 평가와 그 보상 및 PD 활용은, 주파수 자원 유효 이용의 통신 방식으로서의 초고차의
[수학식 799]
Figure 112019032179095-pct00799
-ary 통신의 실현을 용이하게 한다. 한편,
[수학식 800]
Figure 112019032179095-pct00800
개의 2차원 PC를 이용하므로
[수학식 801]
Figure 112019032179095-pct00801
다중 통신을 행하는 것도 가능하다.
<8.3 다차원 비가환성을 이용한 신호 처리>
Daughman은 가보 함수
[수학식 802]
Figure 112019032179095-pct00802
를 2차원으로 확장하였다[비특허 문헌 33, 34]. 사전 준비로서 Howe[비특허 문헌 5]를 따라,
[수학식 803]
Figure 112019032179095-pct00803
에 대해 2개의 시프트 연산자
[수학식 804]
Figure 112019032179095-pct00804
를 도입한다. 단, ν·t는 내적을 나타낸다. 계속해서,
[수학식 805]
Figure 112019032179095-pct00805
에 대해 scalar product
[수학식 806]
Figure 112019032179095-pct00806
을 도입하면, 집합
[수학식 807]
Figure 112019032179095-pct00807
[수학식 808]
Figure 112019032179095-pct00808
상의 유니터리 연산의 군이 되어,
[수학식 809]
Figure 112019032179095-pct00809
(H를 차수 n의 Heisenberg group이라고 부른다)상에서의 결합법칙
[수학식 810]
Figure 112019032179095-pct00810
이 성립한다. 또한,
[수학식 811]
Figure 112019032179095-pct00811
을 도입하고 푸리에 변환을 이용하면, 2개의 H의 유니터리 표현
[수학식 812]
Figure 112019032179095-pct00812
사이에
[수학식 813]
Figure 112019032179095-pct00813
의 관계가 성립한다(Howe('80))
Daughman[비특허 문헌 35]은, 다음 식의 서로 동일 형식의 2차원 공간 영역(space domain: SD)의 공간 변수
[수학식 814]
Figure 112019032179095-pct00814
의 complex Gabor elementary function g(x)와 2차원 공간 주파수 영역(frequency domain: FD)의 공간 주파수
[수학식 815]
Figure 112019032179095-pct00815
의 2-D 푸리에 변환(FT)
[수학식 816]
Figure 112019032179095-pct00816
[비특허 문헌 35, 36]
[수학식 817]
Figure 112019032179095-pct00817
의 2-D 가보 필터족을 정의했다. 한편, P는 직류분이 0(零)이 되는 초기 위상이다.
[수학식 818]
Figure 112019032179095-pct00818
은 타원형 Gaussian을 시계 방향으로 θ 회전시킨
[수학식 819]
Figure 112019032179095-pct00819
SD 함수 g(x), FD 함수
[수학식 820]
Figure 112019032179095-pct00820
은 완전 대칭으로 각각의 시프트량
[수학식 821]
Figure 112019032179095-pct00821
의 위치에서 피크가 되는 2-D Gaussian envelope를 포함해, 2차 모멘트
[수학식 822]
Figure 112019032179095-pct00822
사이의 불확정성 관계[비특허 문헌 35, 36]
[수학식 823]
Figure 112019032179095-pct00823
을 제공한다(등식은 식 (97)의 2-D 가보 필터족에 의해 달성된다.). 한편,
[수학식 824]
Figure 112019032179095-pct00824
는 변조, 스케일 파라미터라고 불린다.
식 (99)의 조건하에서 최적인 해상도를 갖는 공간·공간 주파수 동시 표현(화상 표현은 음성의 스펙트로그램의 4차원 확장판)이 화상 처리에서 중요하다.
임의의 주어진 2-D 신호 I[x](예를 들면, 256×256 픽셀상의 화상, i.e., (256)2=65536차원 벡터 공간의 벡터)의 표현을 위해, I[x]를 임의의 선택된 2-D 벡터의 집합
[수학식 825]
Figure 112019032179095-pct00825
상으로의 투영에 의한 최적 표현을 생각한다[비특허 문헌 35]. 이는 선형합
[수학식 826]
Figure 112019032179095-pct00826
과의 오차 에너지 ||I[x]-H[x]||2을 최소화하는 전개 계수
[수학식 827]
Figure 112019032179095-pct00827
의 결정에 귀착된다. Daughman은 비직교계
[수학식 828]
Figure 112019032179095-pct00828
에 대한, ai를 구하는 구배법이 n×n(상기 예에서 n=65536)의 희소행렬 문제가 되어 그 실행은 비현실적이므로, 신경망(neural network)에 기초하는 방법을 제안하고 있다. 그 외에, 화상의 에지 검출이나 특징 추출 등의 화상 해석에 관한 여러 가지 방법이 제안되고 있다[비특허 문헌 37, 38, 39].
본 특허에서는, 식 (100)이 공간·공간 주파수 시프트(x0i, u0i)를 갖는 2D 가우스 함수
[수학식 829]
Figure 112019032179095-pct00829
에 의한 2D 가보 전개이므로, 이전 절의 1차원 비가환 연산에 기초하는 신호 처리법을 모방해, SD 표현(100)의 FD 표현을 대칭적으로 고찰하기 위해, I[x], H[x]의 2-D 푸리에 변환
[수학식 830]
Figure 112019032179095-pct00830
을 이용해 화상의 SD-FD 표현 문제
[수학식 831]
Figure 112019032179095-pct00831
을 생각한다. 화상을 힐베르트 공간의 부분 공간 ε,ε'의 요소
[수학식 832]
Figure 112019032179095-pct00832
라고 생각하고, 오차
[수학식 833]
Figure 112019032179095-pct00833
의 최소화를 각각 n쌍의 SD, FD의 직교 투영 연산자
[수학식 834]
Figure 112019032179095-pct00834
에 의한
[수학식 835]
Figure 112019032179095-pct00835
의 직합 분해(Direct Sum Decomposition)
[수학식 836]
Figure 112019032179095-pct00836
으로 행한다. template 집합
[수학식 837]
Figure 112019032179095-pct00837
의 선택이 중요하다. 예를 들면,
[수학식 838]
Figure 112019032179095-pct00838
을 각각 (L1Δx×L2Δy)-공간 제한의 가우스 함수나 (L1Δu×L2Δv)-공간 주파수 대역 제한의 가우스 함수로 하면,
[수학식 839]
Figure 112019032179095-pct00839
는, 각각 힐베르트 공간의 부분 공간 상으로의 직교 투영(orthogonal projection: OP)이다. 두 OP에 2차원의 von Neumann의 APT를 적용할 수 있다. 두 OP로부터 APT의 교호 투영인, 이들의 결합 연산자(localization operator)
[수학식 840]
Figure 112019032179095-pct00840
이 유도되므로
[수학식 841]
Figure 112019032179095-pct00841
이들은, 2×2차원의 SD-FD 공간중의 피크-어드레스(x0i, u0i)의
[수학식 842]
Figure 112019032179095-pct00842
의 부분 영역을 추출해 다른 영역의 신호를 필터링한다(filter-out). 이 신호 재생법은, 종래법과 완전히 다른, 2차원 비가환성을 이용한 신호 처리이다. 단,
[수학식 843]
Figure 112019032179095-pct00843
이다.
[수학식 844]
Figure 112019032179095-pct00844
[수학식 845]
Figure 112019032179095-pct00845
의 직교 공간이다. 한편, 이 경우, Howe의 시프트 연산자(93)와 달리, SD, FD 신호
[수학식 846]
Figure 112019032179095-pct00846
의 대칭성 때문에, 하프 시프트
[수학식 847]
Figure 112019032179095-pct00847
을 갖는 2차원의 von Neumann의 SD, FD의 symmetrical space-space frequency shift operator: SFSO
[수학식 848]
Figure 112019032179095-pct00848
을 이용하면, 2×2차원의 SD-FD 공간중의 피크-어드레스(x0i, u0i)를 갖는 2-D 가보 함수
[수학식 849]
Figure 112019032179095-pct00849
및 그
[수학식 850]
Figure 112019032179095-pct00850
[수학식 851]
Figure 112019032179095-pct00851
DFT를 생각한다. 공간 좌표
[수학식 852]
Figure 112019032179095-pct00852
의 샘플링 간격을 (M, N), 공간 변수
[수학식 853]
Figure 112019032179095-pct00853
의 SD 함수
[수학식 854]
Figure 112019032179095-pct00854
에 대해, 공간 주파수
[수학식 855]
Figure 112019032179095-pct00855
의 샘플링 간격
[수학식 856]
Figure 112019032179095-pct00856
의 이산화 공간 주파수
[수학식 857]
Figure 112019032179095-pct00857
의 FD 함수
[수학식 858]
Figure 112019032179095-pct00858
은, SD 함수
[수학식 859]
Figure 112019032179095-pct00859
와의 사이의 twiddle factor
[수학식 860]
Figure 112019032179095-pct00860
을 갖는 2차원 DFT
[수학식 861]
Figure 112019032179095-pct00861
의 관계
[수학식 862]
Figure 112019032179095-pct00862
로 정해진다. 따라서, SD 함수
[수학식 863]
Figure 112019032179095-pct00863
(또는, FD 함수
[수학식 864]
Figure 112019032179095-pct00864
)의 support는
[수학식 865]
Figure 112019032179095-pct00865
(또는,
[수학식 866]
Figure 112019032179095-pct00866
)이다. SD 공간의 피크-어드레스 x0i 간격을 (MxΔx, MyΔy), FD 공간의 피크-어드레스 u0i의 피크 간격을
[수학식 867]
Figure 112019032179095-pct00867
로 하고, 정규화 조건
[수학식 868]
Figure 112019032179095-pct00868
을 요청하면,
[수학식 869]
Figure 112019032179095-pct00869
를 얻는다. 또한,
[수학식 870]
Figure 112019032179095-pct00870
의 공간 시프트, 주파수 시프트도
[수학식 871]
Figure 112019032179095-pct00871
로 이산화하면, von Neumann의 이산 TFSO의 2차원 버전은, 다음 식의 2차원의 대칭인 공간 시프트·공간 주파수 시프트 연산자 symmetrical SFSO
[수학식 872]
Figure 112019032179095-pct00872
로서 정의된다. 화상 처리에는 이것에 기초하는 계산법이 효과적일 것이다. 왜냐하면, 식 (108)의 두 시프트 연산자에 의해 하프 시프트
[수학식 873]
Figure 112019032179095-pct00873
의 위상항이 각각 SD, FD 신호
[수학식 874]
Figure 112019032179095-pct00874
에 매장되어 있으므로, 식 (102)의 내적에서
[수학식 875]
Figure 112019032179095-pct00875
가 유효하게 추출되기 때문이다. 한편, 시간 변화가 있는 화상 해석은 3차원 비가환성을 이용한 신호 처리의 대상이라고 생각할 수 있다.
<9. 본 명세서에 기재된 발명의 특징>
본 명세서에서 구축한, 시간 지연, 도플러 시프트의 고정도 추정법은. 비가환성을 이용한 다원 통신 시스템의 일례에 지나지 않는다. 이것은, 통신 이론의 창시자 가보(Gabor)가 1946년의 논문[비특허 문헌 1]에 이어지는 논문[비특허 문헌 2]에서, 양자역학의 2개의 비가환 작용소에 주목해 제창한 「비가환성을 이용한 다원 통신」의 하나의 구체적 예제라고 볼 수 있다.
TFP상에서 실시하는 신호의 중첩에 수반하는, 신호의 시간 및 주파수 시프트의 비가환성의 이해가 진행되지 않은 현재의 상태에서, 비가환성을 이용한 통신 시스템의 실현이 요구된다.
본 명세서에서 제안한 시간, 주파수 시프트의 비가환성을 이용한 통신 시스템은 적어도 다음의 관점을 갖는다.
(관점 1)
특허 문헌 1 및 참고 문헌[비특허 문헌 26, 27, 30]에서 정의·도입한, 시간 주파수 대칭인 시프트 연산자 symmetrical TFSO(식 4, 24)는,
(i) 그 하프 시프트에 의해 파라미터 추정에서 중요한 시프트량을 TD, FD의 PD로서 현재화하고;
(ⅱ) 비가환인 변복조에 수반하는 PD
[수학식 876]
Figure 112019032179095-pct00876
의 발생이 자명하고;
(ⅲ) 비가환성에 수반하는 PD 평가는 twiddle factor
[수학식 877]
Figure 112019032179095-pct00877
의 거듭제곱 연산에 귀착될 수 있으므로, 그 PD 보상법이 용이하게 되는 등의 특징을 갖는다.
(관점 2)
TD-, FD-PC 변조(2차원 BPSK)는 주지의 통신 기술이지만, 시간, 주파수 시프트 연산의 일종이라고 이해함으로써, 각각 TD-, FD-PC 변조한 광대역 송신 신호에는 PC에 의한 PD가 TD-, FD-template으로 매장된다.
(관점 3)
주파수 자원의 효과적 이용의 상투적 수단인, TFP상의 신호의 무중첩 중합에서 이용하는 시간, 주파수 시프트 연산은, 데이터 레벨의 PD를 유발한다. 이 PD는 파라미터 추정에 중요한 요소이다.
(관점 4)
파라미터 추정을 위한 M종 검정법에서 필요한 TD-, FD-template 검출의 TD-, FD-우도 함수는 각각 TD-, FD-CCF 어레이를 유도한다. TD-, FD-CCF는 힐베르트 공간의 부분 공간인 TL-TD 공간, BL-FD 공간 상으로의 PO P3, P4(또는 P1, P2)를 정의하고, APT의 적용의 프레임을 제공한다.
(관점 5)
APT의 교호 투영을 나타내는 PO쌍의 결합 연산자
[수학식 878]
Figure 112019032179095-pct00878
(또는,
[수학식 879]
Figure 112019032179095-pct00879
)은, 나이퀴스트 조건을 만족하지 않는다는 이유로 통신에서는 이용되지 않는 가우스 함수의 AF의 변수 분리성이나 지수 함수적 감쇠 특성에 의해, TFP의 특정 부분을 국소 선택하는 localization operator가 된다. 그 결과. TD-, FD-CCF 어레이는 양질인 수신기가 된다.
전술한 각 관점은, 특허 문헌 및 비특허 문헌의 공표 이후에 분명히 한 새로운 관점이며, 본 명세서에 기재된 발명의 포인트는 이하와 같이 표현할 수도 있다.
(포인트 1)
시간·주파수 대칭인 시간 주파수 시프트 연산자 symmetrical TFSO(식 4, 24)에 기초해, 시간 지연과 주파수 시프트의 2개의 비가환성을 이용한 다중 통신을 위한 TFP상에서 실시하는 비가환 조작에 수반하는 PD 평가법 및 PD 보상법의 알고리즘과 그 이론적 근거를 제공한다.
(포인트 2)
TD, FD의 우도 함수를 구성함으로써 최대 우도 상관기에 의한 신호 복원이나 파라미터 최대 우도 추정을 실현한다.
(포인트 3)
고전적 위상 편이 변조(BPSK)는 그 비가환성에 의한 PD를 유발하는데, 이 PD는 파라미터 추정이나 신호 복원에서의 중요한 단서가 된다.
(포인트 4)
통신의 상투적 수단인, TFP상에서의 신호의 무중첩 중합이나 변복조에 수반해 발생하는 비가환의 PD는 데이터 레벨의 신호 검출에 유효하다.
(포인트 5)
이들 PD를 모두 보상하는 TD의 CCF(또는 FD의 CCF)는, 각각 von Neumann의 APT의 힐베르트 공간의 2개의 불가결한 부분 공간(시간 제한 공간, 대역 제한 공간)상으로의 직교 투영 연산자를 정의한다.
(포인트 6)
TD, FD의 직교 투영 연산자의 결합 연산자가 TFP상에서의 국소 선택 연산자가 되어, 관용의 DSP의 첨예한 필터의 역할을 한다.
(포인트 7)
화상 신호 표현에서는 2×2의 공간·공간 주파수 평면(Space-Frequency plane: SFP)상의 신호의 분리성이 뛰어난 성질을 갖는 2차원 가보 함수(97)에 의한 가보 전개(100)가 많이 쓰이고 있다. 그러나, 종래의 기술에서는, 가보 함수에 내포되는 2차원의 공간·주파수 시프트 연산자가 1차원 신호의 비가환 시간·주파수 시프트 연산자와 마찬가지로, 비가환의 시프트 유래의 위상 왜곡의 발생을 알아차리지 못한 듯하다. SD, FD 신호를 대칭적으로 취급하는 화상의 공간·공간 주파수의 동시 표현 문제(101)를 힐베르트 공간의 SD-영역 제한 함수나 FD-대역 제한 함수상으로의 직교 투영 연산자
[수학식 880]
Figure 112019032179095-pct00880
에 의한 직합 분해(102)로 풀어, SD, FD에 관해 대칭인 공간·주파수 시프트 연산자 SFSO(105)에 의한 2차원 가보 함수(106)나 이산 버전 symmetrical SFSO(108)를 정의하고, localization operator(103)를 유도해, von Neumann의 APT에 기초하는 화상 표현이나 화상 복원 등의 비가환성을 이용한 화상 처리법의 이론적 프레임을 제공한다.
《송수신 시스템의 구성예 1》
이하에서는, 전술한 이론적 측면에 준거한 송수신 시스템의 제1 구성예에 대해 도면을 참조해 설명한다.
도 18은, 본 구성예에 따른 송수신 시스템(1)의 구성예를 나타내는 블록도이다. 도 18에 나타낸 바와 같이, 통신 시스템(1)은 송신 장치(100) 및 수신 장치(200)를 구비하고 있다. 송수신 시스템(1)은 레이더로서 이용되어도 되고, 레이더 이외의 통신 시스템(예를 들면, 주로 음성 데이터나 화상 데이터 등을 송수신하는 데이터 송수신 시스템)으로 이용되어도 된다.
<송신 장치(100)>
송신 장치(100)는, 전술한 설명에서 「송신기」의 동작으로서 기재된 사항을 실제로 실행하는 장치이다.
일례로서 도 18에 나타낸 바와 같이, 송신 장치(100)는 송신용 데이터 취득부(101), 송신 신호 생성부(102) 및 송신부(103)를 구비하고 있다.
(송신용 데이터 취득부(101))
송신용 데이터 취득부(101)는 송신 대상의 데이터를 취득한다. 송수신 시스템(1)을 데이터 송수신 시스템으로서 이용하는 경우에는, 송신 대상의 데이터는, 예를 들면 음성 데이터, 화상 데이터 및 텍스트 데이터 중 적어도 어느 하나를 디지털 데이터화한 것이라도 되고, 그 외의 데이터라도 무방하다.
또한, 송수신 시스템(1)을 레이더로서 이용하는 경우에는, 송신 대상의 데이터는 레이더용 펄스파여도 된다.
(송신 신호 생성부(102))
송신 신호 생성부(102)는, 송신용 데이터 취득부(101)가 취득한 송신 대상 데이터에 대해, 송신 신호 생성 처리를 행함으로써 송신 신호를 생성한다.
송신 신호 생성부(102)는 일례로서, 도 4∼도 7에 나타낸 SFB의 적어도 어느 하나를 구비해 구성된다.
송신 신호 생성 처리의 예로는, 전술한 주기 N, N'의 TD-, FD-phase code: PC를 이용한 변조 처리, 즉 TD-, FD-Binary phase-shift Keying: 2차원 BPSK 등을 들 수 있다.
송신 신호 생성부(102)가 생성하는 송신 신호로는, 예를 들면, 전술한
[수학식 881]
Figure 112019032179095-pct00881
등을 들 수 있다. 또한, 송신 신호 생성부(102)에 의한 구체적인 처리로서 <4. TD-, FD-signature와 template>에서 설명한 각 처리를 들 수 있다.
또한, 송신 신호 생성부(102)는 전술한 <7.1 SFB: signature 송신기와 레이더 신호 송신기>에서 설명한 각 처리를 실행하는 구성으로 해도 된다. 또한, 송신 신호 생성부(102)는 전술한 <8.2 인공적인 비가환 시프트를 매장하는 delay-Doppler space Division multiplexing에 기초하는 초고 MPSK>에서 설명한 각 처리를 실행하는 구성으로 해도 된다.
(송신부(103))
송신부(103)는 송신 신호 생성부(102)가 생성한 송신 신호를 송신한다.
<수신 장치(200)>
수신 장치(200)는, 전술한 설명에서 「수신기」의 동작으로서 기재된 사항을 실제로 실행하는 장치이다.
일례로서 수신 장치(200)는, 도 18에 나타낸 바와 같이, 수신부(201), 시프트 추정 및 수신 데이터 추출부(202), 수신 데이터 출력부(203)를 구비하고 있다.
(수신부(201))
수신부(201)는 송신 장치(100)가 송신한 신호를 수신한다. 수신부(201)가 수신하는 신호는, 예를 들면 전술한
[수학식 882]
Figure 112019032179095-pct00882
를 들 수 있다.
(시프트 추정 및 수신 데이터 추출부(202))
시프트 추정 및 수신 데이터 추출부(202)(단순히, 추정부라고도 한다)는, 수신부(201)가 수신한 신호에 대해 시프트 추정 처리를 행함과 함께, 수신 데이터를 추출한다.
시프트 추정 및 수신 데이터 추출부(202)는, 일례로서 도 8∼도 10에 나타낸 AFB의 적어도 어느 하나를 구비해 구성된다.
시프트 추정 및 수신 데이터 추출부(202)는, 예를 들면 전술한 <5. M종 가설 검정에 의한 TD-, FD-신호 검출과 추정> 및 <6. 파라미터 추정용 TD-, FD-CCF>에서 설명한 각 처리를 행한다.
또한, 시프트 추정 및 수신 데이터 추출부(202)는, 전술한 <7.2 AFB: 수신기와 복호기>에서 설명한 각 처리를 실행하는 구성으로 해도 된다. 또한, 송신 신호 생성부(102)는 전술한 <8.2 인공적인 비가환 시프트를 매장하는 delay-Doppler space Division multiplexing에 기초하는 초고 MPSK>에서 설명한 각 처리를 실행하는 구성으로 해도 된다.
따라서, 시프트 추정 및 수신 데이터 추출부(202)는, 신호를 수신하는 수신 방법(도 13 참조)으로서, 여차원 2를 갖는 비가환 시프트의 파라미터 공간(도 17 참조)을 참조해, 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정하는 추정 스텝(도 13 및 도 13의 Switch 2-1,2-2의 상하측 접속에 따라 각각 식 (38), (43) 및 식 (89), (90) 참조)을 실행한다.
또한, 송신 신호 생성부(102)는, 여차원 2를 갖는 비가환 시프트의 파라미터 공간을 참조해, 송신 대상 신호의 시간 및 주파수를 시프트하는 시프트 스텝(도 12 및 도 12의 Switch 1-1,1-2의 상하측 접속에 따라 각각 식 (25), (27) 및 식 (88), (87) 참조)을 실행한다.
또한, 상기 여차원 2를 갖는 비가환 시프트의 파라미터 공간은, 시간을 나타내는 제1 축 및 주파수를 나타내는 제2 축에 의해 펼쳐지는 평면을, 시간 및 주파수 시프트를 나타내는 제3 축에 의해 3차원화해 얻어지는 공간(도 17 참조)이다.
한편, 전술한
[수학식 883]
Figure 112019032179095-pct00883
-MC 자신도 AC와 마찬가지로, 여차원 2의 파라미터 평면이며, 시간·주파수의 신호 평면 TFP상에 실시되는 시프트 연산에 수반하는 파라미터 평면이다. 「shift 평면」과 시간 t와 주파수 f의 TFP(도 16 참조)와 구별할 필요가 있다. 신호의 변수 t, f에 대해서 행해지는 4칙연산, 미분, 적분 등과는 별개로, 위상에 관련되는 비가환의 시프트 연산(도 17의 제3축 참조)을 특별히 다뤘다.
또한, 특허 문헌 6∼8에 관련해 다음과 같이 부언한다.
a) 변조·복조시의
[수학식 884]
Figure 112019032179095-pct00884
(또는 그 이산 버전
[수학식 885]
Figure 112019032179095-pct00885
) 및 감쇠 상수 Ae의 PD e를 모두 포함시켜 연립시킬 필요가 있다. 왜냐하면, 시프트의 군론적 성질로부터 각 PD는 독립적이 아니고 전파하므로 정확한 PD 보상이 되지 않는다. 시간 시프트, 및 주파수 시프트 추정값의 갱신 과정에서는, 이들 PD의 보상과의 연동(동시 병행 처리: 식 59)에 유의할 필요가 있다.
b) 본 특허의 시간 시프트, 주파수 시프트는 각각 시간, 주파수 변수에 대응하는, 여차원 2의 파라미터 공간의 변수의 호칭이며, 이들의 비가환 시프트 연산은 위상 함수(위상 왜곡)를 유발한다. 한편, 통신에서의 시간 오프셋 및 주파수 오프셋은 주기나 반송 주파수의 「어긋남」이다.
(1) 특허 문헌 7, 8의 추정된 시간 오프셋 및 주파수 오프셋이나, 특허 문헌 7의 OFDM 방식의 통신 방식에서의 시간, 주파수 오프셋은, 각각 펄스 파형의 시간폭이나 반송 주파수의 요동, 본 특허의 시간 시프트 td, 주파수 시프트 fD에 대응한다. 두 특허 문헌은 td, fD 사이의 비가환성을 무시하고 있으므로, 통상의 동기 포착·동기 유지법으로서 이해된다.
(2)「의도적인」 시프트와 오프셋은 혼용되기 쉽기 때문에, 본 명세서에서는 「오프셋」이 아니라, 「지수 함수 어깨의 하프 시프트」(종래의 비가환 연산의 시프트의 반)라는 표현을 이용하고 있다. 또한, 본 명세서에 기재된 발명의 일 특징은 송신 신호에 「사전의 하프 시프트된 위상 함수」를 포함하는 것이나, 비가환성을 갖는 시간, 주파수 시프트의 군론적 성질을 고려한 보상법이라고 표현할 수 있다.
c) 특허 문헌 8의 채널 추정을 위한 수신된 프리엠블(preamble)은 시간 영역 template 신호이며, 주파수 영역 신호는 언급하고 있지 않고, 비가환성에 근거하지 않는, 통상의 채널 추정법이다.
(1) 특허 문헌 6의 PUL 알고리즘에 대해, 1) 힐베르트 공간인 신호 공간의 부분 공간으로서, 각종 template이 속하는 힐베르트 공간 εi, 1≤i≤4를 정의하고, 2) 상기 제1의 N'개의 상관 함수(type3- 또는 type1-CCF)가 NMΔt(또는 LΔt)-TL-TD 함수의 힐베르트 공간 ε3(또는 ε1)상으로의 직교 투영 연산자 P3(또는 P1)인 것 및 상기 제2의 N개의 상관 함수(type4- 또는 type2-CCF)가 N'M'Δf(또는 LΔf)-BL-FD 함수의 힐베르트 공간 ε4(또는 ε2)상으로의 직교 투영 연산자 P4(또는 P2)인 것, 3) TFP에서 2차원 template-matching을 행하는, 국소 선택 연산자 LO의 교대 투영 정리: Alternating Projection Theorem(APT) Operator: APTO P3F-1,dP4Fd(또는 P4FdP3F-1,d)에 의한 갱신 과정은, 두 힐베르트 공간의 합병 집합내의 함수에 수렴한다(von Neumann의 APT: F-1d, Fd는 IDFT, DFT), 4) 그 TL-TD 함수 및 BL-FD 함수의 합병 집합은 공집합, 즉 영함수이다(Youla의 정리) 등을 분명히 하고, 또한, PUL는 상관 함수의 종류수 N, N'의 곱 N·N'가 아니라, 합 N+N'에 비례하는 계산 복잡도와 가우스 칩 펄스의 시간폭 LΔt, 대역폭 LΔf의 정도의 시간 시프트, 주파수 시프트의 추정값을 제공하는 것을 증명하였다(L=(ΔtΔf)-1= MM').
(2) 특허 문헌 8의 채널 추정의 MMSE 알고리즘은, Bayes 법칙에 기초하는 2 파라미터 추정이므로, 계산 복잡도는 곱(N·N')이 된다. 한편, 본 명세서에서의 우도 함수는 Bayes 법칙과 관련이 있지만, 시간 시프트 추정의 N'개의 우도 함수와 주파수 시프트 추정의 N개의 우도 함수로 이루어져, 각각 TL-TD 공간, BL-FD 공간에서 교대로 최대 우도 추정을 행하므로 계산 복잡도는 합(N+N')이 되어, 문헌 3의 방법과 전혀 다른 파라미터 추정법이다.
또한, 시프트 추정 및 수신 데이터 추출부(202)는, 식 (84) 및 관련 기재를 이용해 설명한 바와 같이, 일례로서 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정하는 추정 스텝을 실행하고, 상기 추정 스텝은, 시간 시프트의 추정값 및 주파수 시프트를 표현하는 제1 시프트 연산자
[수학식 886]
Figure 112019032179095-pct00886
과, 시간 시프트 및 주파수 시프트의 추정값을 표현하는 제2 시프트 연산자
[수학식 887]
Figure 112019032179095-pct00887
을 이용해, 상기 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정한다.
상기 추정 스텝은, 시간 시프트의 추정값, 주파수 시프트 및 추정 시간 시프트의 절반인 하프 시프트의 위상항을 표현하는 제1 시프트 연산자
[수학식 888]
Figure 112019032179095-pct00888
과, 시간 시프트, 주파수 시프트의 추정값 및 추정 주파수 시프트의 절반인 하프 시프트의 위상항을 표현하는 제2 시프트 연산자
[수학식 889]
Figure 112019032179095-pct00889
와, 추정 대상의 시간 시프트 및 추정 대상의 주파수 시프트의 관측값, 그리고, 추정 대상의 시간 시프트의 절반인 하프 시프트의 위상항을 표현하는 제3 시프트 연산자(
[수학식 890]
Figure 112019032179095-pct00890
또는 그 주파수쌍의 짝
[수학식 891]
Figure 112019032179095-pct00891
)와, 시간 시프트 및 주파수 시프트의 추정값, 및, 추정 시간 시프트의 절반인 하프 시프트의 위상항을 표현하는 제4 시프트 연산자(
[수학식 892]
Figure 112019032179095-pct00892
또는 그 주파수쌍의 짝
[수학식 893]
Figure 112019032179095-pct00893
)를 이용해, 상기 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정한다.
송신 시간 신호의 위상 함수에는, 시간에 관한 1 또는 복수의 시프트 파라미터의 시프트량의 절반인 하프 시프트의 위상량, 및, 송신 주파수 신호의 위상 함수에는 주파수에 관한 1 또는 복수의 시프트 파라미터의 시프트량의 절반인 하프 시프트의 위상량을 포함하고 있다.
상기와 같이 구성된 수신 장치(200)에 의하면, 고효율의 수신 장치를 실현할 수 있다. 일례로서, 수신 장치(200)를 레이더 수신 장치로서 구성했을 경우, 고속의 수신 장치를 실현할 수 있다.
또한, 시프트 추정 및 수신 데이터 추출부(202)가 실행하는 상기 추정 스텝에서는, 식 (84) 및 관련 기재를 이용해 설명한 바와 같이, 일례로서 상기 제1 시프트 연산자를 이용해 표시되는 제1 상관 함수
[수학식 894]
Figure 112019032179095-pct00894
, 및, 상기 제2 시프트 연산자를 이용해 표시되는 제2 상관 함수
[수학식 895]
Figure 112019032179095-pct00895
를 참조해 상기 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정한다.
또한, 전술한 바와 같이, 상기 제1 상관 함수는
[수학식 896]
Figure 112019032179095-pct00896
으로 표현되고, 상기 제2 상관 함수는
[수학식 897]
Figure 112019032179095-pct00897
으로 표현된다.
상기 추정 스텝에서는, 상기 제1 시프트 연산자, 제3 시프트 연산자, 및 제4 시프트 연산자를 이용해 표시되는 제1 상관 함수,
[수학식 898]
Figure 112019032179095-pct00898
및, 상기 제2 시프트 연산자, 제3 시프트 연산자 및 제4 시프트 연산자를 이용해 표시되는 제2 상관 함수
[수학식 899]
Figure 112019032179095-pct00899
를 참조해 상기 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정한다.
이와 같이, 제1 상관 함수(식 (89) 참조) 및 제2 상관 함수(식 (90) 참조)는, 변조 주파수 fc의 변조 및 복조에서의, 시간 시프트 td에 수반하는 비가환 연산에 의한 위상 왜곡
[수학식 900]
Figure 112019032179095-pct00900
또는 그 이산 버전
[수학식 901]
Figure 112019032179095-pct00901
및, 전송시의 위상 왜곡 e를 포함하고 있다.
또한, 시프트 추정 및 수신 데이터 추출부(202)는, 식 (60), 식 (61) 및 관련 기재를 이용해 설명한 바와 같이, 시간 시프트의 추정값 및 주파수 시프트의 추정값
[수학식 902]
Figure 112019032179095-pct00902
또는,
[수학식 903]
Figure 112019032179095-pct00903
의 갱신 처리를,
[수학식 904]
Figure 112019032179095-pct00904
로 정해지는
[수학식 905]
Figure 112019032179095-pct00905
를 각각
[수학식 906]
Figure 112019032179095-pct00906
으로 설정함으로써 행한다.
이와 같이, 시프트 추정 및 수신 데이터 추출부(202)가 행하는 추정 처리에는, 제1 상관 함수
[수학식 907]
Figure 112019032179095-pct00907
(또는 식 (89))을 참조해
[수학식 908]
Figure 112019032179095-pct00908
을 구하고, 그 값을
[수학식 909]
Figure 112019032179095-pct00909
로 설정함으로써 주파수 시프트의 추정값을 갱신하는 공정과, 제2 상관 함수
[수학식 910]
Figure 112019032179095-pct00910
(또는 식 (90))을 참조해
[수학식 911]
Figure 112019032179095-pct00911
을 구하고, 그 값을
[수학식 912]
Figure 112019032179095-pct00912
로 설정함으로써 시간 시프트의 추정값을 갱신하는 공정을 교대로 반복하는 교호 갱신 스텝이 포함되어 있다.
또한, 상기 추정 스텝에서는, 전술한 바와 같이, N'개의 TD-template 신호 검출의 우도 함수에 의한 주파수 시프트 최대 우도 추정과, N개의 FD-template 신호 검출의 우도 함수에 의한 시간 시프트 최대 우도 추정에 기초해, 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정한다.
또한, 상기 최대 우도 추정에 있어서, 식 (89)의 상관 함수
[수학식 913]
Figure 112019032179095-pct00913
및 식 (90)의 상관 함수
[수학식 914]
Figure 112019032179095-pct00914
를 참조하는 경우에는, 식 (60)의 제1식 우변
[수학식 915]
Figure 112019032179095-pct00915
및 제2식
[수학식 916]
Figure 112019032179095-pct00916
대신에 각각 적용하고, 또한, argmax 연산에서는,
[수학식 917]
Figure 112019032179095-pct00917
에 대해서도 최대 우도 추정의 대상이 된다.
또한, 상기 최대 우도 추정에 있어서, 수신 장치는
[수학식 918]
Figure 112019032179095-pct00918
통신에서 k의 부호화기의 일부의
[수학식 919]
Figure 112019032179095-pct00919
를 통과한 수신
[수학식 920]
Figure 112019032179095-pct00920
의 수신부이다.
이와 같이, 본 실시 형태에 따른 수신 장치는, 신호를 수신하는 수신 장치로서, 여차원 2를 갖는 비가환 시프트의 파라미터 공간을 참조해, 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정하는 추정부를 포함하고 있다.
또한, 상기 최대 우도 추정에 있어서, 식 (87)(또는 (27))의 송신 TD-CE, FD-CE에서는(여기에서, 도 12의 송신기의 Switch 1-1, 1-2의 상하측에 따라 다중화의 레벨이 상이하므로 대응하는 식을 인용한다), 우변의 식 (88)(또는 (25))의 TD-signature, FD-signature가 다중화되고, 두 signature 자신도 가우스 함수나 그 FD 함수를 2차원 PC에 의한 변조를 이용한 다중화로 얻어진다. 한편,
[수학식 921]
Figure 112019032179095-pct00921
통신의 송신 TD-CE, FD-CE는 식 (87)로 나타낸 바와 같이, TD-signature, FD-signature에 서로 독립인 2차원 PC의 집합
[수학식 922]
Figure 112019032179095-pct00922
에 의한 다중화 및 k의 부호
[수학식 923]
Figure 112019032179095-pct00923
에 따른
[수학식 924]
Figure 112019032179095-pct00924
의 시프트 대응의 시프트 연산
[수학식 925]
Figure 112019032179095-pct00925
가 작용하고 있다.
이와 같이, 본 실시 형태에 따른 송신 방법은, 신호를 송신하는 송신 방법으로서, 여차원 2를 갖는 비가환 시프트의 파라미터 공간을 참조해, 송신 대상 신호의 시간 및 주파수를 시프트하는 시프트 스텝을 포함하고 있다.
또한, 상기 최대 우도 추정에 있어서,
[수학식 926]
Figure 112019032179095-pct00926
통신에서의 k의 부호화기의 일부의
[수학식 927]
Figure 112019032179095-pct00927
을 통과한 수신
[수학식 928]
Figure 112019032179095-pct00928
과 추정 수신 template의 상관 함수(식 (89))(우변 제1, 제2항이 각각 송신 신호, 추정 수신 template이다)를 이용하는 경우, 식 (87)~(90)중의 합 Σ가 의미하는 「다중화=무중첩의 중합」은 각종의 시프트 연산자로 시뮬레이션할 수 있다. 도 12의 송신 신호의 생성 과정에서 각종 다중화를 행한다. 식 (88), (87)에서는 각각 TD-, FD-칩 펄스를 주기
[수학식 929]
Figure 112019032179095-pct00929
의 2차원 PC 변조
[수학식 930]
Figure 112019032179095-pct00930
에 의한 다중화로 signature를 생성하고, 얻어진 signature를 데이터 레벨의 시프트를 이용한 데이터 다중화로 송신 TD-CE, FD-CE를 정의하고 있다. 한편, 시프트가 없는 단순 가산에서는 template 검출이 불가능하므로, 시프트 연산자가 불가결한 것이거나 시프트 연산에 수반하는 PD는 template 검출을 위해 치뤄야 할 대가임과 동시에, 정확한 PD 왜곡 평가의 단서가 되는 것도 자명할 것이다. 한편, 식 (89)의 우변에 변복조에 수반하는 PD
[수학식 931]
Figure 112019032179095-pct00931
, k의 부호
[수학식 932]
Figure 112019032179095-pct00932
[수학식 933]
Figure 112019032179095-pct00933
의 PD 보상이 실시되고 있다.
이와 같이, 상기 시프트 스텝에서는, 시간 펄스 파형을 주기 N의 시간 위상 부호 변조하고, 상기 시간 위상 부호 변조한 시간 펄스 파형을 다시 주기 N'의 주파수 영역 위상 부호 변조에 의해 멀티 캐리어화함으로써 송신 신호를 생성한다.
또한, 본 실시 형태에 따른 송신 장치에서는, 도 12의 Swtich1-1, 1-2에서 하측에 접속된 상태에서 시간 펄스 파형을 주기 N의 시간 위상 부호 변조하고, 상기 시간 위상 부호 변조한 시간 펄스 파형을 다시 주기 N'의 주파수 영역 위상 부호 변조에 의해 멀티 캐리어화하는 등의 각종의 다중화를 행한, 데이터 다중화 된 TD-signature, FD-signature를 k의 부호
[수학식 934]
Figure 112019032179095-pct00934
대응의
[수학식 935]
Figure 112019032179095-pct00935
에 통과시킴으로써
[수학식 936]
Figure 112019032179095-pct00936
이 얻어진다.
또한, 도 13의 Switch 2-1, 2-2에서 하측에 접속된 상태에서는, 식 (89), (90)의 상관 함수의 실부의 최대화는 4변수
[수학식 937]
Figure 112019032179095-pct00937
에 관한 argmax 연산으로 행하고, PUL에 기초하는
[수학식 938]
Figure 112019032179095-pct00938
의 교호 갱신에 의한 수렴값으로
[수학식 939]
Figure 112019032179095-pct00939
의 복원이나 k의 복호를 행한다.
또한, 이하의 각종 구성에 의해 규정되는 화상 신호로의 확장은, 8.3절에 기재한 바와 같이, 공간, 공간 주파수 시프트 사이의 비가환성에 수반하는 하프 시프트를 이용한, 1차원 신호로부터 2차원 신호로의 자연스러운 확장이다(식 (102)∼(108) 참조).
(화상 신호로의 확장 구성 1)
화상 신호를 수신하는 수신 방법으로서, 수신한 화상 신호가 나타내는 공간 시프트 및 공간 주파수 시프트를 파라미터 공간을 참조해 추정하는 추정 스텝을 포함하고, 상기 공간 시프트 및 공간 주파수 시프트의 각각은 2 이상의 차원을 갖고 있다.
(화상 신호로의 확장 구성 2)
상기 추정 스텝에서는, 2차원의 대칭적인 공간 시프트 및 공간 주파수 시프트 연산자
[수학식 940]
Figure 112019032179095-pct00940
[수학식 941]
Figure 112019032179095-pct00941
을 참조해, 공간 시프트 및 공간 주파수 시프트를 추정한다.
(화상 신호로의 확장 구성 3)
화상 신호를 송신하는 송신 방법으로서, 파라미터 공간을 참조해, 송신 대상의 화상 신호의 공간 및 주파수를 시프트하는 시프트 스텝을 포함하고, 상기 공간의 시프트 및 상기 주파수의 시프트의 각각은 2 이상의 차원을 갖고 있다.
(화상 신호로의 확장 구성 4)
상기 시프트 스텝에서는, 2차원의 대칭인 공간 시프트 및 공간 주파수 시프트 연산자
[수학식 942]
Figure 112019032179095-pct00942
[수학식 943]
Figure 112019032179095-pct00943
을 참조해, 공간 및 주파수를 시프트한다.
(수신 데이터 출력부(203))
수신 데이터 출력부(203)는 시프트 추정 및 수신 데이터 추출부(202)가 추출한 수신 데이터를 출력한다.
<송수신의 흐름>
도 19는 송수신 시스템(1)을 이용한 데이터 송수신 처리의 흐름을 나타내는 플로우차트이다.
(S101)
우선, 스텝 S101에서, 송신용 데이터 취득부(101)는 송신용 데이터를 취득한다. 송신용 데이터 취득부(101)에 의한 구체적인 처리는 전술한 바와 같다.
(S102)
계속해서, 스텝 S102에서, 송신 신호 생성부(102)는 송신 신호를 생성한다. 송신 신호 생성부(102)에 의한 구체적인 처리는 전술한 바와 같다.
(S103)
계속해서, 스텝 S103에서, 송신부(103)는 송신 신호를 송신한다. 송신부(103)에 의한 구체적인 처리는 전술한 바와 같다.
(S201)
계속해서, 스텝 S201에서, 수신부(201)는 송신부(103)가 송신한 신호를 수신한다. 수신부(201)에 의한 구체적인 처리는 전술한 바와 같다.
(S202)
계속해서, 스텝 S202에서, 시프트 추정 및 수신 데이터 추출부(202)는 시프트 추정 처리를 행함과 함께, 수신 데이터를 추출한다. 시프트 추정 및 수신 데이터 추출부(202)에 의한 구체적인 처리는 전술한 바와 같다.
(S203)
계속해서, 스텝 S203에서, 수신 데이터 출력부(203)는 시프트 추정 및 수신 데이터 추출부(202)가 추출한 수신 데이터를 출력한다. 수신 데이터 출력부(203)에 의한 구체적인 처리는 전술한 바와 같다.
《통신 시스템의 구성예 2》
이하에서는, 전술한 이론적 측면에 준거한 통신 시스템의 제1 구성예에 대해 도면을 참조해 설명한다.
도 20은 본 구성예에 따른 송수신 시스템(1a)의 구성예를 나타내는 블록도이다. 도 20에 나타낸 바와 같이, 통신 시스템(1a)은 송신 장치(100a) 및 수신 장치(200a)를 구비하고 있다. 송수신 시스템(1a)은, 전술한 송수신 시스템(1)과 마찬가지로, 레이더로서 이용되어도 되고, 레이더 이외의 통신 시스템(예를 들면, 주로 음성 데이터나 화상 데이터 등을 송수신 하는 데이터 송수신 시스템)으로서 이용되어도 된다. 한편, 이미 설명한 부재에는 동일한 부호를 부여하고, 그 설명은 생략한다.
<송신 장치(100a)>
송신 장치(100a)는, 전술한 설명에서 「송신기」의 동작으로서 기재된 사항을 실제로 실행하는 장치이다.
일례로서 도 20에 나타낸 바와 같이, 송신 장치(100a)는 송신용 데이터 취득부(101), 송신 신호 생성부(102a) 및 송신부(103)를 구비하고 있다.
(송신 신호 생성부(102a))
송신 신호 생성부(102a)는, 구성예 1에 따른 송신 신호 생성부(102)가 구비하는 구성 외에, 시프트 매장부(110)를 구비하고 있다.
시프트 매장부(110)는, 식 (84)의 직후에 언급한 바와 같이, 파라미터값
[수학식 944]
Figure 112019032179095-pct00944
를 매장한 송신 신호를 생성한다.
따라서, 송신 신호 생성부(102a)는, 시간에 관한 1 또는 복수의 시프트 파라미터, 및, 주파수에 관한 1 또는 복수의 시프트 파라미터
[수학식 945]
Figure 112019032179095-pct00945
를 이용해, 송신 대상 신호의 시간 및 주파수를 시프트한다.
<수신 장치(200a)>
수신 장치(200a)는, 전술한 설명에서 「수신기」의 동작으로서 기재된 사항을 실제로 실행하는 장치이다.
일례로서 수신 장치(200a)는, 도 20에 나타낸 바와 같이, 수신부(201), 시프트 추정 및 수신 데이터 추출부(202a), 수신 데이터 출력부(203)를 구비하고 있다.
(시프트 추정 및 수신 데이터 추출부(202a))
시프트 추정 및 수신 데이터 추출부(202a)는, 일례로서 구성예 1에 따른 시프트 추정 및 수신 데이터 추출부(202)가 갖는 구성과 같은 구성을 구비한다.
<송수신의 흐름>
도 21은 송수신 시스템(1a)을 이용한 데이터 송수신 처리의 흐름을 나타내는 플로우차트이다. 스텝 S101, S103, S201, S203에 대해서는, 도 19를 이용해 설명한 처리와 동일하기 때문에 여기에서는 설명을 생략한다.
(S102a)
스텝 S102a에서, 송신 신호 생성부(102a)는 송신 신호를 생성한다. 송신 신호 생성부(102a)에 의한 구체적인 처리는 전술한 바와 같다.
(S202a)
스텝 S202에서는, 시프트 추정 및 수신 데이터 추출부(202a)가 시프트 추정 처리를 행함과 함께, 수신 데이터를 추출한다. 시프트 추정 및 수신 데이터 추출부(202a)에 의한 구체적인 처리는 전술한 바와 같다.
《통신 시스템의 구성예 3》
도 18∼도 21을 참조해 설명한 송수신 시스템은, <8.3 다차원 비가환성을 이용한 신호 처리>에서 설명한 각 스텝을 실행하는 구성으로 할 수 있다.
일례로서, <8.3 다차원 비가환성을 이용한 신호 처리>에서 설명한 바와 같이, 수신 장치(200)는, 수신한 화상 신호가 나타내는 공간 시프트 및 공간 주파수 시프트를 파라미터 공간을 참조해 추정하는 추정 스텝을 실행하고, 상기 공간 시프트 및 공간 주파수 시프트의 각각은 2 이상의 차원을 갖고 있다.
또한, <8.3 다차원 비가환성을 이용한 신호 처리>에서 설명한 바와 같이, 일례로서 상기 추정 스텝에서는, 2차원의 대칭인 공간 시프트 및 공간 주파수 시프트 연산자, 및, 공간 주파수 시프트 사이의 비가환성에 수반하는 하프 시프트의 위상항을 표현하는 연산자
[수학식 946]
Figure 112019032179095-pct00946
[수학식 947]
Figure 112019032179095-pct00947
을 참조해, 공간 시프트 및 공간 주파수 시프트를 추정한다.
마찬가지로, 본 구성예에 따른 송신 방법은, 화상 신호를 송신하는 송신 방법으로서, 파라미터 공간을 참조해, 송신 대상의 화상 신호의 공간 및 주파수를 시프트하고, 공간 시프트의 절반인 하프 시프트의 위상항(또는, 공간, 주파수를 시프트하고, 공간 주파수 시프트의 절반인 하프 시프트의 위상항)을 표현하는 2개의 연산자를 이용한 시프트 스텝을 포함하고, 상기 공간의 시프트 및 상기 주파수의 시프트의 각각은 2 이상의 차원을 갖고 있다.
상기 시프트 스텝에서는, 2차원의 대칭인 공간 시프트 및 공간 주파수 시프트 연산자
[수학식 948]
Figure 112019032179095-pct00948
[수학식 949]
Figure 112019032179095-pct00949
를 참조해, 공간 및 주파수를 시프트한다.
〔소프트웨어에 의한 실현예〕
송신 장치(100, 100a) 및 수신 장치(200, 200a)의 제어 블록(특히, 송신 신호 생성부(102, 102a), 시프트 추정 및 수신 데이터 추출부(202, 202a))은, 집적회로(IC칩) 등에 형성된 논리 회로(하드웨어)에 의해 실현되어도 되고, 소프트웨어에 의해 실현되어도 무방하다.
후자의 경우, 송신 장치(100, 100a) 및 수신 장치(200, 200a)는, 각 기능을 실현하는 소프트웨어인 프로그램의 명령을 실행하는 컴퓨터를 구비하고 있다. 이 컴퓨터는, 예를 들면 1개 이상의 프로세서를 구비함과 함께, 상기 프로그램을 기억한 컴퓨터 판독 가능한 기록 매체를 구비하고 있다. 그리고, 상기 컴퓨터에 있어서, 상기 프로세서가 상기 프로그램을 상기 기록 매체로부터 읽어내 실행함으로써, 본 발명의 목적이 달성된다. 상기 프로세서로는, 예를 들면 CPU(Central Processing Unit)를 이용할 수 있다. 상기 기록 매체로는 「일시적이 아닌 유형의 매체」, 예를 들면 ROM(Read Only Memory) 등 외에, 테이프, 디스크, 카드, 반도체 메모리, 프로그래머블 논리 회로 등을 이용할 수 있다. 또한, 상기 프로그램을 전개하는 RAM(Random Access Memory) 등을 더 구비하고 있어도 된다. 또한, 상기 프로그램은, 그 프로그램을 전송 가능한 임의의 전송 매체(통신 네트워크나 방송파 등)를 개재해 상기 컴퓨터에 공급되어도 된다. 한편, 본 발명의 일 형태는, 상기 프로그램이 전자적인 전송에 의해 구현화된, 반송파에 내재된 데이터 신호의 형태로도 실현될 수 있다.
또한, 송신 장치(100, 100a) 및 수신 장치(200, 200a)의 각 기능을 실현하는 컴퓨터 프로그램 제품도 본 발명의 범주에 포함된다. 상기 컴퓨터 프로그램 제품은, 임의의 전송 매체를 개재해 제공되는 프로그램을 적어도 1개의 컴퓨터를 개재해 로딩하고, 그 컴퓨터에 적어도 1개의 프로그램 명령을 실행시킨다. 이에 따라, 상기 적어도 1개의 컴퓨터가 구비하는 프로세서가, 상기 프로그램 명령에 부합하는 처리를 실행함으로써, 송신 장치(100, 100a) 및 수신 장치(200, 200a)의 각 기능이 실현된다. 이 컴퓨터 프로그램 제품은, 본 실시 형태에 따른 송신 처리(송신 방법) 및 수신 처리(수신 방법)의 각 스텝을, 프로그램을 로딩한 적어도 1개의 컴퓨터에 실행시키는 것이다.
본 발명은 전술한 각 실시 형태로 한정되는 것이 아니라, 청구항에 기재된 범위에서 여러 가지의 변경이 가능하고, 다른 실시 형태에 각각 개시된 기술적 수단을 적절하게 조합해 얻어지는 실시 형태에 대해서도 본 발명의 기술적 범위에 포함된다. 또한, 각 실시 형태에 각각 개시된 기술적 수단을 조합함으로써, 새로운 기술적 특징을 형성할 수 있다.
《산업상 이용 가능성》
본 발명은 무선 송수신 시스템 및 레이더 시스템 등으로서 바람직하게 이용할 수 있다.
1, 1a 송수신 시스템
100, 100a 송신 장치
101 송신용 데이터 취득부
102, 102a 송신 신호 생성부
103 송신부
200, 200a 수신 장치
201 수신부
202, 202a 시프트 추정 및 수신 데이터 추출부(추정부)
203 수신 데이터 출력부

Claims (19)

  1. 신호를 수신하는 수신 방법으로서,
    여차원 2를 갖는 비가환 시프트의 파라미터 공간을 참조해, 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정하는 추정 스텝을 포함하고,
    상기 추정 스텝에서는, 추정 시간 시프트의 절반인 하프 시프트를 포함하는 연산자 및 추정 주파수 시프트의 절반인 하프 시프트를 포함하는 연산자를 이용해 상기 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정하는 것을 특징으로 하는 수신 방법.
  2. 제1항에 있어서,
    상기 여차원 2를 갖는 비가환 시프트의 파라미터 공간은, 시간을 나타내는 제1 축 및 주파수를 나타내는 제2 축에 의해 펼쳐지는 평면을, 시간 및 주파수 시프트를 나타내는 제3 축에 의해 3차원화해 얻어지는 공간인 것을 특징으로 하는 수신 방법.
  3. 제1항 또는 제2항에 있어서,
    상기 추정 스텝은,
    시간 시프트의 추정값, 주파수 시프트 및 추정 시간 시프트의 절반인 하프 시프트의 위상항을 표현하는 제1 시프트 연산자
    [수학식 1]
    Figure 112019032179095-pct00950

    과,
    시간 시프트, 주파수 시프트의 추정값 및 추정 주파수 시프트의 절반인 하프 시프트의 위상항을 표현하는 제2 시프트 연산자
    [수학식 2]
    Figure 112019032179095-pct00951

    와,
    추정 대상의 시간 시프트 및 추정 대상의 주파수 시프트의 관측값, 및, 추정 시간 시프트의 절반인 하프 시프트의 위상항을 표현하는 제3 시프트 연산자
    [수학식 3]
    Figure 112019032179095-pct00952

    또는 그 주파수쌍의 짝
    [수학식 4]
    Figure 112019032179095-pct00953

    와,
    시간 시프트, 주파수 시프트의 추정값 및 추정 시간 시프트의 절반인 하프 시프트의 위상항을 표현하는 제4 시프트 연산자
    [수학식 5]
    Figure 112019032179095-pct00954

    또는 그 주파수쌍의 짝
    [수학식 6]
    Figure 112019032179095-pct00955

    를 이용해, 상기 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정하는 것을 특징으로 하는 수신 방법.
  4. 제3항에 있어서,
    상기 추정 스텝에서는, 상기 제1 시프트 연산자, 제3 시프트 연산자 및 제4 시프트 연산자를 이용해 표시되는 제1 상관 함수, 그리고, 상기 제2 시프트 연산자, 제3 시프트 연산자 및 제4 시프트 연산자를 이용해 표시되는 제2 상관 함수를 참조해, 상기 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정하는 것을 특징으로 하는 수신 방법.
  5. 제4항에 있어서,
    상기 제1 상관 함수 및 제2 상관 함수는, 변조 주파수 fc의 변조 및 복조에서의, 시간 시프트 td에 수반하는 비가환 연산에 의한 위상 왜곡
    [수학식 7]
    Figure 112019032179095-pct00956

    또는 그 이산 버전
    [수학식 8]
    Figure 112019032179095-pct00957

    및, 전송시의 위상 왜곡 e를 포함하는 것을 특징으로 하는 수신 방법.
  6. 제4항에 있어서,
    상기 추정 스텝에는, 상기 제1 상관 함수를 참조해 주파수 시프트의 추정값을 갱신하는 공정과, 상기 제2 상관 함수를 참조해 시간 시프트의 추정값을 갱신하는 공정을 교대로 반복하는 교호 갱신 스텝이 포함되어 있는 것을 특징으로 하는 수신 방법.
  7. 제1항 또는 제2항에 있어서,
    상기 추정 스텝에서는, N'개의 TD-template 신호 검출의 우도 함수에 의한 주파수 시프트 최대 우도 추정과, N개의 FD-template 신호 검출의 우도 함수에 의한 시간 시프트 최대 우도 추정에 기초해, 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정하는 것을 특징으로 하는 수신 방법.
  8. 신호를 수신하는 수신 장치로서,
    여차원 2를 갖는 비가환 시프트의 파라미터 공간을 참조해, 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정하는 추정부를 포함하고,
    상기 추정부에서는, 추정 시간 시프트의 절반인 하프 시프트를 포함하는 연산자 및 추정 주파수 시프트의 절반인 하프 시프트를 포함하는 연산자를 이용해 상기 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정하는 것을 특징으로 하는 수신 장치.
  9. 신호를 송신하는 송신 방법으로서,
    여차원 2를 갖는 비가환 시프트의 파라미터 공간을 참조해, 송신 대상 신호의 시간 및 주파수를 시프트하는 시프트 스텝을 포함하고,
    상기 시프트 스텝에서는, 시간 및 주파수에 관해 대칭적 시프트인 하프 시프트만큼 송신 대상 신호의 시간 및 주파수를 시프트하는 것을 특징으로 하는 송신 방법.
  10. 제9항에 있어서,
    상기 시프트 스텝에서는, 시간 펄스 파형을 주기 N의 시간 위상 부호 변조하고, 상기 시간 위상 부호 변조한 시간 펄스 파형을 다시 주기 N'의 주파수 영역 위상 부호 변조에 의해 멀티 캐리어화함으로써 송신 신호를 생성하는 것을 특징으로 하는 송신 방법.
  11. 신호를 송신하는 송신 장치로서,
    여차원 2를 갖는 비가환 시프트의 파라미터 공간을 참조해, 송신 대상 신호의 시간 및 주파수를 시프트하는 시프트부를 포함하고,
    상기 시프트부에서는, 시간 및 주파수에 관해 대칭적 시프트인 하프 시프트만큼 송신 대상 신호의 시간 및 주파수를 시프트하는 것을 특징으로 하는 송신 장치.
  12. 송신 장치와 수신 장치를 포함하는 송수신 시스템으로서,
    상기 송신 장치는, 여차원 2를 갖는 비가환 시프트의 파라미터 공간을 참조해, 송신 대상 신호의 시간 및 주파수를 시프트하는 시프트부를 포함하고,
    상기 시프트부에서는, 시간 및 주파수에 관해 대칭적 시프트인 하프 시프트만큼 송신 대상 신호의 시간 및 주파수를 시프트하고,
    상기 수신 장치는, 여차원 2를 갖는 비가환 시프트의 파라미터 공간을 참조해, 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정하는 추정부를 포함하고,
    상기 추정부에서는, 추정 시간 시프트의 절반인 하프 시프트를 포함하는 연산자 및 추정 주파수 시프트의 절반인 하프 시프트를 포함하는 연산자를 이용해 상기 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정하는 것을 특징으로 하는 송수신 시스템.
  13. 화상 신호를 수신하는 수신 방법으로서,
    수신한 화상 신호가 나타내는 공간 시프트 및 공간 주파수 시프트를 여차원 2를 갖는 비가환 시프트의 복수의 파라미터 공간을 참조해 추정하는 추정 스텝을 포함하고,
    상기 공간 시프트 및 공간 주파수 시프트의 각각은 2 이상의 차원을 갖고,
    상기 추정 스텝에서는, 공간 및 공간 주파수에 관해 대칭적 시프트인 하프 시프트를 포함하는 연산자를 이용해 상기 화상 신호가 나타내는 공간 시프트 및 공간 주파수 시프트를 추정하는 것을 특징으로 하는 수신 방법.
  14. 제13항에 있어서,
    상기 추정 스텝에서는, 2차원의 대칭인 공간 시프트 및 공간 주파수 시프트 연산자
    [수학식 9]
    Figure 112019032179095-pct00958


    [수학식 10]
    Figure 112019032179095-pct00959

    을 참조해, 공간 시프트 및 공간 주파수 시프트를 추정하는 것을 특징으로 하는 수신 방법.
  15. 화상 신호를 송신하는 송신 방법으로서,
    여차원 2를 갖는 비가환 시프트의 복수의 파라미터 공간을 참조해, 송신 대상의 화상 신호의 공간 및 주파수를 시프트하는 시프트 스텝을 포함하고,
    상기 공간의 시프트 및 상기 주파수의 시프트의 각각은 2 이상의 차원을 갖고,
    상기 시프트 스텝에서는, 공간 및 공간 주파수에 관해 대칭적 시프트인 하프 시프트만큼, 송신 대상의 신호의 화상의 공간 및 공간 주파수를 시프트하는 것을 특징으로 하는 송신 방법.
  16. 제15항에 있어서,
    상기 시프트 스텝에서는, 2차원의 대칭인 공간 시프트 및 공간 주파수 시프트 연산자
    [수학식 11]
    Figure 112019077716259-pct00960

    및, 그 주파수쌍의 짝
    [수학식 12]
    Figure 112019077716259-pct00961

    를 참조해, 공간 및 주파수를 시프트하는 것을 특징으로 하는 송신 방법.
  17. 신호를 수신하는 수신 방법으로서,
    여차원 2를 갖는 비가환 시프트의 파라미터 공간을 참조해, 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정하는 추정 스텝을 포함하고,
    상기 추정 스텝은,
    시간 시프트의 추정값, 주파수 시프트 및 추정 시간 시프트의 절반인 하프 시프트의 위상항을 표현하는 제1 시프트 연산자
    [수학식 13]
    Figure 112019136169279-pct00983

    과,
    시간 시프트, 주파수 시프트의 추정값 및 추정 주파수 시프트의 절반인 하프 시프트의 위상항을 표현하는 제2 시프트 연산자
    [수학식 14]
    Figure 112019136169279-pct00984

    와,
    추정 대상의 시간 시프트 및 추정 대상의 주파수 시프트의 관측값, 및, 추정 시간 시프트의 절반인 하프 시프트의 위상항을 표현하는 제3 시프트 연산자
    [수학식 15]
    Figure 112019136169279-pct00985

    또는 그 주파수쌍의 짝
    [수학식 16]
    Figure 112019136169279-pct00986

    와,
    시간 시프트, 주파수 시프트의 추정값 및 추정 시간 시프트의 절반인 하프 시프트의 위상항을 표현하는 제4 시프트 연산자(
    [수학식 17]
    Figure 112019136169279-pct00987

    또는 그 주파수쌍의 짝
    [수학식 18]
    Figure 112019136169279-pct00988

    )를 이용해, 상기 수신 신호가 나타내는 시간 시프트 및 주파수 시프트를 추정하는 것을 특징으로 하는 수신 방법.
  18. 화상 신호를 수신하는 수신 방법으로서,
    수신한 화상 신호가 나타내는 공간 시프트 및 공간 주파수 시프트를 파라미터 공간을 참조해 추정하는 추정 스텝을 포함하고,
    상기 공간 시프트 및 공간 주파수 시프트의 각각은 2 이상의 차원을 갖고,
    상기 추정 스텝에서는,
    2차원의 대칭인 공간 시프트 및 공간 주파수 시프트 연산자
    [수학식 19]
    Figure 112019136169279-pct00989


    [수학식 20]
    Figure 112019136169279-pct00990

    을 참조해, 공간 시프트 및 공간 주파수 시프트를 추정하는 것을 특징으로 하는 수신 방법.
  19. 화상 신호를 송신하는 송신 방법으로서,
    파라미터 공간을 참조해, 송신 대상의 화상 신호의 공간 및 주파수를 시프트하는 시프트 스텝을 포함하고,
    상기 공간의 시프트 및 상기 주파수의 시프트의 각각은 2 이상의 차원을 갖고,
    상기 시프트 스텝에서는,
    2차원의 대칭인 공간 시프트 및 공간 주파수 시프트 연산자
    [수학식 21]
    Figure 112019136169279-pct00991

    및, 그 주파수쌍의 짝
    [수학식 22]
    Figure 112019136169279-pct00992

    를 참조해, 공간 및 주파수를 시프트하는 것을 특징으로 하는 송신 방법.
KR1020197009104A 2018-01-22 2018-06-28 수신 방법, 수신 장치, 송신 방법, 송신 장치, 송수신 시스템 KR102120021B1 (ko)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JPPCT/JP2018/001735 2018-01-22
JP2018001735 2018-01-22
PCT/JP2018/024592 WO2019142372A1 (ja) 2018-01-22 2018-06-28 受信方法、受信装置、送信方法、送信装置、送受信システム

Publications (2)

Publication Number Publication Date
KR20190089850A KR20190089850A (ko) 2019-07-31
KR102120021B1 true KR102120021B1 (ko) 2020-06-05

Family

ID=67301880

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020197009104A KR102120021B1 (ko) 2018-01-22 2018-06-28 수신 방법, 수신 장치, 송신 방법, 송신 장치, 송수신 시스템

Country Status (11)

Country Link
US (2) US10666486B2 (ko)
EP (1) EP3544249A4 (ko)
JP (2) JP6602506B1 (ko)
KR (1) KR102120021B1 (ko)
CN (1) CN110301121B (ko)
AU (1) AU2018338597B2 (ko)
CA (3) CA3052035C (ko)
IL (1) IL268326A (ko)
RU (1) RU2710926C1 (ko)
SG (2) SG11201902729VA (ko)
WO (1) WO2019142372A1 (ko)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019142372A1 (ja) * 2018-01-22 2019-07-25 ラディウス株式会社 受信方法、受信装置、送信方法、送信装置、送受信システム
US11764940B2 (en) 2019-01-10 2023-09-19 Duality Technologies, Inc. Secure search of secret data in a semi-trusted environment using homomorphic encryption
RU2734699C1 (ru) * 2020-05-26 2020-10-22 Федеральное государственное автономное образовательное учреждение высшего образования "Дальневосточный федеральный университет" (ДВФУ) Способ передачи информации с помощью широкополосных сигналов
CN114338297B (zh) * 2021-11-26 2023-06-30 河南工程学院 一种非相干LoRa***下的联合定时同步与频偏估计方法
CN114928520B (zh) * 2022-04-22 2023-12-12 南京邮电大学 一种用于补偿频偏的广义频分复用时频同步方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013126508A (ja) 2011-12-19 2013-06-27 Toshiba Corp 磁気共鳴イメージング装置
JP2016189502A (ja) 2013-10-04 2016-11-04 国立大学法人九州大学 受信器、通信方法及びプログラム
US20170207817A1 (en) 2010-05-28 2017-07-20 Cohere Technologies, Inc. Modulation and equalization in an orthonormal time-frequency shifting communications system

Family Cites Families (47)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09321721A (ja) 1996-05-27 1997-12-12 Sony Corp 送信方法,送信装置,受信方法及び受信装置
US7243064B2 (en) 2002-11-14 2007-07-10 Verizon Business Global Llc Signal processing of multi-channel data
CN1317863C (zh) * 2003-07-02 2007-05-23 矽统科技股份有限公司 无线局域网络的频率位移补偿估计***及其方法
KR100719112B1 (ko) * 2005-07-19 2007-05-17 삼성전자주식회사 Ofdm 시스템에 적용되는 샘플링 주파수 오차 추정장치및 그 방법
JP4991871B2 (ja) * 2006-10-26 2012-08-01 クゥアルコム・インコーポレイテッド 無線通信システムにおけるパケット検出のための方法および装置
CN102160310B (zh) 2008-09-22 2014-06-18 松下电器产业株式会社 无线通信装置及信号分割方法
EP2222011A3 (en) * 2009-02-19 2014-10-29 Samsung Electronics Co., Ltd. Method for performing hybrid automatic repeat request operation in a wireless mobile communication system
US8385222B2 (en) * 2009-10-26 2013-02-26 Telefonaktiebolaget Lm Ericsson (Publ) Method and apparatus for channel quality derivation
US9668148B2 (en) 2010-05-28 2017-05-30 Cohere Technologies, Inc. OTFS methods of data channel characterization and uses thereof
US8547988B2 (en) 2010-05-28 2013-10-01 Ronny Hadani Communications method employing orthonormal time-frequency shifting and spectral shaping
US9071286B2 (en) 2011-05-26 2015-06-30 Cohere Technologies, Inc. Modulation and equalization in an orthonormal time-frequency shifting communications system
US9071285B2 (en) 2011-05-26 2015-06-30 Cohere Technologies, Inc. Modulation and equalization in an orthonormal time-frequency shifting communications system
US9130638B2 (en) 2011-05-26 2015-09-08 Cohere Technologies, Inc. Modulation and equalization in an orthonormal time-frequency shifting communications system
US9083595B2 (en) 2010-05-28 2015-07-14 Cohere Technologies, Inc. Signal modulation method resistant to echo reflections and frequency offsets
US9444514B2 (en) 2010-05-28 2016-09-13 Cohere Technologies, Inc. OTFS methods of data channel characterization and uses thereof
US8902978B2 (en) 2010-05-30 2014-12-02 Lg Electronics Inc. Enhanced intra prediction mode signaling
KR101585447B1 (ko) 2011-02-15 2016-01-18 삼성전자주식회사 무선통신 시스템에서 옵셋 보상을 고려한 엠엠에스이 채널추정 장치 및 방법
JP5751553B2 (ja) 2011-05-10 2015-07-22 国立大学法人九州大学 拡散装置、通信装置、送信装置、通信方法及びプログラム
US9590779B2 (en) 2011-05-26 2017-03-07 Cohere Technologies, Inc. Modulation and equalization in an orthonormal time-frequency shifting communications system
US9294315B2 (en) 2011-05-26 2016-03-22 Cohere Technologies, Inc. Modulation and equalization in an orthonormal time-frequency shifting communications system
US9031141B2 (en) 2011-05-26 2015-05-12 Cohere Technologies, Inc. Modulation and equalization in an orthonormal time-frequency shifting communications system
CN102236268B (zh) 2011-07-20 2012-10-10 中国科学院上海光学精密机械研究所 基于空间像频谱的光刻投影物镜波像差检测方法
EP2670102B1 (en) 2012-05-31 2016-03-30 Mitsubishi Electric R&D Centre Europe B.V. Method and device for synchronizing a receiver on received preamble symbol
EP2860622B1 (en) 2012-06-08 2020-04-08 NEC Corporation Electronic device and controlling method and program therefor
US9912507B2 (en) 2012-06-25 2018-03-06 Cohere Technologies, Inc. Orthogonal time frequency space communication system compatible with OFDM
US10090972B2 (en) 2012-06-25 2018-10-02 Cohere Technologies, Inc. System and method for two-dimensional equalization in an orthogonal time frequency space communication system
US9967758B2 (en) 2012-06-25 2018-05-08 Cohere Technologies, Inc. Multiple access in an orthogonal time frequency space communication system
US10003487B2 (en) 2013-03-15 2018-06-19 Cohere Technologies, Inc. Symplectic orthogonal time frequency space modulation system
CN104662855B (zh) 2012-06-25 2018-10-23 科希尔技术股份有限公司 正交时频移动通信***中的调制和均衡
US9929783B2 (en) * 2012-06-25 2018-03-27 Cohere Technologies, Inc. Orthogonal time frequency space modulation system
US10469215B2 (en) 2012-06-25 2019-11-05 Cohere Technologies, Inc. Orthogonal time frequency space modulation system for the Internet of Things
US10411843B2 (en) 2012-06-25 2019-09-10 Cohere Technologies, Inc. Orthogonal time frequency space communication system compatible with OFDM
JP5725458B2 (ja) 2012-08-27 2015-05-27 国立大学法人九州大学 受信器、受信方法及びプログラム
US9294949B2 (en) 2012-08-28 2016-03-22 Intel Deutschland Gmbh Interference and noise estimation of a communications channel
JP2016189501A (ja) 2013-08-28 2016-11-04 国立大学法人九州大学 通信方法及び受信器
JP2016189500A (ja) 2013-08-27 2016-11-04 国立大学法人九州大学 同期装置、同期方法及びプログラム
US10116476B2 (en) 2014-01-22 2018-10-30 European Space Agency Receiving method and receiver for satellite-based automatic identification systems
JP6508539B2 (ja) * 2014-03-12 2019-05-08 ソニー株式会社 音場収音装置および方法、音場再生装置および方法、並びにプログラム
US10090973B2 (en) 2015-05-11 2018-10-02 Cohere Technologies, Inc. Multiple access in an orthogonal time frequency space communication system
WO2016183230A1 (en) 2015-05-11 2016-11-17 Cohere Technologies Systems and methods for symplectic orthogonal time frequency shifting modulation and transmission of data
CN112532558A (zh) 2015-06-27 2021-03-19 凝聚技术股份有限公司 与ofdm兼容的正交时频空间通信***
US10693581B2 (en) 2015-07-12 2020-06-23 Cohere Technologies, Inc. Orthogonal time frequency space modulation over a plurality of narrow band subcarriers
KR20190008827A (ko) 2015-09-07 2019-01-25 코히어 테크널러지스, 아이엔씨. 직교 시간 주파수 공간 변조를 이용한 다중액세스
EP3437197B1 (en) * 2016-04-01 2022-03-09 Cohere Technologies, Inc. Tomlinson-harashima precoding in an otfs communication system
MY194217A (en) * 2016-05-13 2022-11-22 Ericsson Telefon Ab L M Wireless communication device, transmitter and methods therein
CN108075996A (zh) * 2016-11-11 2018-05-25 索尼公司 控制装置、设备及方法、信号处理装置和方法及移动终端
WO2019142372A1 (ja) * 2018-01-22 2019-07-25 ラディウス株式会社 受信方法、受信装置、送信方法、送信装置、送受信システム

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170207817A1 (en) 2010-05-28 2017-07-20 Cohere Technologies, Inc. Modulation and equalization in an orthonormal time-frequency shifting communications system
JP2013126508A (ja) 2011-12-19 2013-06-27 Toshiba Corp 磁気共鳴イメージング装置
JP2016189502A (ja) 2013-10-04 2016-11-04 国立大学法人九州大学 受信器、通信方法及びプログラム

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Tai Sing Lee, Image representation using 2D Gabor wavelets, IEEE Transactions on Pattern Analysis and Machine Intelligence Volume 18, Issue 10, 1996.10.*
Yutaka Jitsumatsu et al, Spread spectrum-based cooperative and individual time-frequency synchronization, ISWCS 2013 The Tenth International Symposium on Wireless Communication Systems, 2013.10.15.*

Also Published As

Publication number Publication date
JP6602506B1 (ja) 2019-11-06
RU2710926C1 (ru) 2020-01-14
CN110301121B (zh) 2024-03-01
CA3052035C (en) 2021-03-30
KR20190089850A (ko) 2019-07-31
EP3544249A4 (en) 2020-04-08
JP6757448B2 (ja) 2020-09-16
US10666486B2 (en) 2020-05-26
EP3544249A1 (en) 2019-09-25
CA3070694A1 (en) 2019-07-25
CN110301121A (zh) 2019-10-01
CA3070694C (en) 2022-05-24
US11108615B2 (en) 2021-08-31
SG10202107059PA (en) 2021-07-29
AU2018338597B2 (en) 2019-09-12
JP2019176505A (ja) 2019-10-10
US20200204426A1 (en) 2020-06-25
WO2019142372A1 (ja) 2019-07-25
SG11201902729VA (en) 2019-08-27
US20190229972A1 (en) 2019-07-25
AU2018338597A1 (en) 2019-08-08
JPWO2019142372A1 (ja) 2020-01-23
CA3052035A1 (en) 2019-07-25
IL268326A (en) 2019-09-26
CA3065055A1 (en) 2019-12-13

Similar Documents

Publication Publication Date Title
KR102120021B1 (ko) 수신 방법, 수신 장치, 송신 방법, 송신 장치, 송수신 시스템
US8108438B2 (en) Finite harmonic oscillator
Yin et al. Integrated waveform for continuous active sonar detection and communication
CN104601512B (zh) 一种检测相位调制信号载波频偏的方法及***
CN107864020A (zh) 水下小目标单分量声散射回波的变换域提取方法
TWI730284B (zh) 接收方法、接收裝置、傳送方法、傳送裝置、傳送接收系統
JP5751553B2 (ja) 拡散装置、通信装置、送信装置、通信方法及びプログラム
JP5725458B2 (ja) 受信器、受信方法及びプログラム
Li et al. Wideband spectrum sensing based on modulated wideband converter with nested array
Tsihrintzis et al. Fractional lower-order statistics-based ambiguity functions for differential delay Doppler estimation
Alimosaymer et al. Systematic approach in designing wavelet packet modulation‐orthogonal frequency division multiplexing radar signal by applying the criterion of least‐squares
Khan et al. Iterative missing data recovery algorithm for non-stationary signals
Didkowsky et al. Optimal code sequences in problems of clock synchronization of communication systems with noise signal phase-shift keying
JP6982288B2 (ja) 信号処理システム、受信方法及びプログラム
Bose Waveform Synthesis for Active Sensing with Emerging Applications
West Wideband spectrum sensing with machine learning
Ipson Viability and Performance of RF Source Localization Using Autocorrelation-Based Fingerprinting
Popoff Fundamentals of Signal Processing in Generalized Metric Spaces: Algorithms and Applications
Zaheer et al. A Fourier-Mellin based Equalization Technique for Wideband Linear Time Varying Channels
Hong Advanced techniques for automatic classification of digitally modulated communication signals
Oberst et al. Binary Phase-shift Keyed Communication Systems
Weber Spread spectrum signal characteristic estimation using exponential averaging and an ad-hoc chip rate estimator
Bykhovskiy ìDigital Signal Processingî No. 2-2018

Legal Events

Date Code Title Description
A201 Request for examination
A302 Request for accelerated examination
E902 Notification of reason for refusal
AMND Amendment
E902 Notification of reason for refusal
AMND Amendment
E601 Decision to refuse application
AMND Amendment
X701 Decision to grant (after re-examination)
GRNT Written decision to grant