KR20220064154A - Apparatus and method for estimating channel of wireless power transmission system using an intelligent reflective surface - Google Patents

Apparatus and method for estimating channel of wireless power transmission system using an intelligent reflective surface Download PDF

Info

Publication number
KR20220064154A
KR20220064154A KR1020200150352A KR20200150352A KR20220064154A KR 20220064154 A KR20220064154 A KR 20220064154A KR 1020200150352 A KR1020200150352 A KR 1020200150352A KR 20200150352 A KR20200150352 A KR 20200150352A KR 20220064154 A KR20220064154 A KR 20220064154A
Authority
KR
South Korea
Prior art keywords
channel
reflective surface
updated
phase variable
measurement
Prior art date
Application number
KR1020200150352A
Other languages
Korean (ko)
Other versions
KR102457507B1 (en
KR102457507B9 (en
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 연세대학교 산학협력단
Priority to KR1020200150352A priority Critical patent/KR102457507B1/en
Publication of KR20220064154A publication Critical patent/KR20220064154A/en
Application granted granted Critical
Publication of KR102457507B1 publication Critical patent/KR102457507B1/en
Publication of KR102457507B9 publication Critical patent/KR102457507B9/en

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J50/00—Circuit arrangements or systems for wireless supply or distribution of electric power
    • H02J50/20—Circuit arrangements or systems for wireless supply or distribution of electric power using microwaves or radio frequency waves
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/04013—Intelligent reflective surfaces
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J50/00—Circuit arrangements or systems for wireless supply or distribution of electric power
    • H02J50/80—Circuit arrangements or systems for wireless supply or distribution of electric power involving the exchange of data, concerning supply or distribution of electric power, between transmitting devices and receiving devices
    • H04B5/0037—
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B5/00—Near-field transmission systems, e.g. inductive or capacitive transmission systems
    • H04B5/70—Near-field transmission systems, e.g. inductive or capacitive transmission systems specially adapted for specific purposes
    • H04B5/79—Near-field transmission systems, e.g. inductive or capacitive transmission systems specially adapted for specific purposes for data transfer in combination with power transfer
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/06—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
    • H04B7/0613—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission
    • H04B7/0615—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal
    • H04B7/0617—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal for beam forming
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

The present invention can provide a channel estimating device of a wireless power transmitting system using an intelligent reflective surface and a method thereof which check the strength of a received signal received by a terminal for a measurement signal transmitted from a transmitter regardless of the number of reflective elements of an IRS or the number of terminals and can estimate a channel between the transmitter and the IRS and a channel between the IRS and the terminal based on the checked strength of the received signal so as to efficiently obtain channel information even in a wireless power transmitting system including a low power terminal.

Description

지능형 반사 표면을 이용한 무선 전력 전송 시스템의 채널 추정 장치 및 방법{APPARATUS AND METHOD FOR ESTIMATING CHANNEL OF WIRELESS POWER TRANSMISSION SYSTEM USING AN INTELLIGENT REFLECTIVE SURFACE}Apparatus and method for estimating channels in a wireless power transmission system using intelligent reflective surfaces

본 발명은 무선 전력 전송 시스템의 채널 추정 장치 및 방법에 관한 것으로, 지능형 반사 표면을 이용한 무선 전력 전송 시스템의 채널 추정 장치 및 방법에 관한 것이다.The present invention relates to an apparatus and method for estimating a channel in a wireless power transmission system, and to an apparatus and method for estimating a channel in a wireless power transmission system using an intelligent reflective surface.

최근 무선 네트워크 기술의 발전과 함께 사물 인터넷(Internet of Things: 이하 IoT) 장치의 수가 급격하게 증가하는 추세이다. 대규모 무선 네트워크를 효과적으로 관리하기 위해서는 IoT 장치를 위한 안정적인 전력 공급 시스템이 필요하다. 이를 위해 최근에는 무선주파수(Radio Frequency: RF) 신호 기반의 무선전력 전송(Wireless Energy Transfer: WET) 시스템이 주목받고 있다.Recently, with the development of wireless network technology, the number of Internet of Things (IoT) devices is rapidly increasing. A reliable power supply system for IoT devices is required to effectively manage large-scale wireless networks. To this end, recently, a wireless energy transfer (WET) system based on a radio frequency (RF) signal is receiving attention.

특히 무선전력 전송의 효율을 높이기 위해 지능형 반사 표면(Intelligent Reflecting Surfaces: 이하 IRS) 사용에 대한 연구가 활발히 이루어지고 있다In particular, studies on the use of Intelligent Reflecting Surfaces (IRS) are being actively conducted to increase the efficiency of wireless power transmission.

도 1은 지능형 반사 표면을 활용한 무선 전력 전송 시스템을 개념을 설명하기 위한 도면이다.1 is a diagram for explaining the concept of a wireless power transmission system utilizing an intelligent reflective surface.

도 1을 참조하면, IRS를 활용한 무선 전력 전송 시스템에서는 송신기(100)와 컨트롤러(200), IRS(300) 및 K개의 단말(401 ~ 40K)이 포함될 수 있다. 송신기(100)는 K개의 단말(401 ~ 40K)로 전력 신호를 송신하고, K개의 단말(401 ~ 40K) 각각은 송신기(100)에서 전송된 전력 신호로부터 동작에 필요한 무선 전력을 수확(wireless energy harvesting)한다.Referring to FIG. 1 , a wireless power transmission system using IRS may include a transmitter 100 , a controller 200 , an IRS 300 , and K terminals 401 to 40K. The transmitter 100 transmits a power signal to the K terminals 401 to 40K, and each of the K terminals 401 to 40K harvests wireless power required for operation from the power signal transmitted from the transmitter 100 (wireless energy). harvesting).

여기서 송신기(100)는 M개의 안테나를 구비하여 전력 신호를 빔포밍하여 전송함으로써 전력 전송 효율성을 향상시킬 수 있으며, K개의 단말(401 ~ 40K)은 각각 하나의 안테나를 구비하는 것으로 가정한다.Here, it is assumed that the transmitter 100 is provided with M antennas to improve power transmission efficiency by beamforming and transmitting a power signal, and each of the K terminals 401 to 40K is provided with one antenna.

IRS(300)는 N개의 반사 소자를 포함하고, 컨트롤러(200)의 제어에 따라 N개의 반사 소자가 송신기(100)에서 전송된 전력 신호를 반사하면서 위상을 변화시켜, 송신기(100)에서 방사된 전력 신호가 K개의 단말(401 ~ 40K)로 더 효율적으로 전송되도록 한다.The IRS 300 includes N reflective elements, and the N reflective elements change the phase while reflecting the power signal transmitted from the transmitter 100 under the control of the controller 200, It allows the power signal to be transmitted more efficiently to the K terminals (401 to 40K).

IoT 시스템에서 K개의 단말(401 ~ 40K)은 다양한 임의의 위치에 분산 배치될 수 있다. 따라서 다수의 안테나를 구비하는 송신기(100)가 빔포밍 기법으로 전력 신호를 송신하더라도, K개의 단말(401 ~ 40K)의 배치 위치에 따라 요구되는 전력을 충분하게 공급하기 어려운 경우가 많다. 특히 송신기(100)와 K개의 단말(401 ~ 40K) 사이에 장애물 등이 위치하는 경우, 빔포밍 기법을 이용하더라도 각 단말에 안정적으로 전력을 공급하기 매우 어렵다.In the IoT system, K terminals 401 to 40K may be distributed in various arbitrary locations. Therefore, even if the transmitter 100 having a plurality of antennas transmits a power signal using the beamforming technique, it is often difficult to sufficiently supply the required power according to the arrangement positions of the K terminals 401 to 40K. In particular, when an obstacle is located between the transmitter 100 and the K terminals 401 to 40K, it is very difficult to stably supply power to each terminal even using the beamforming technique.

이러한 한계를 극복하기 위해 IRS(300)는 송신기(100)에서 전송된 전력 신호를 반사하여 K개의 단말(401 ~ 40K) 방향으로 안정적으로 전송되도록 한다. 이때, IRS(300)는 송신기(100)에서 전송된 전력 신호를 위상을 변화시켜 반사하여 전력 전송 효율성을 높일 수 있다.In order to overcome this limitation, the IRS 300 reflects the power signal transmitted from the transmitter 100 so as to be stably transmitted in the direction of the K terminals 401 to 40K. In this case, the IRS 300 may change the phase of the power signal transmitted from the transmitter 100 and reflect it to increase power transmission efficiency.

컨트롤러(200)는 송신기(100)에서 전송되는 전력 신호의 빔포밍 벡터에 따라 IRS(300)의 반사 위상을 조절하는 한편, IRS(300)의 반사 위상 정보를 송신기(100)로 전달하여, 송신기(100)가 빔포밍 벡터를 변경하도록 한다.The controller 200 adjusts the reflection phase of the IRS 300 according to the beamforming vector of the power signal transmitted from the transmitter 100, while transmitting the reflection phase information of the IRS 300 to the transmitter 100, the transmitter Let (100) change the beamforming vector.

다만 전력 신호가 효율적으로 K개의 단말(401 ~ 40K)로 전송되도록 송신기(100)와 컨트롤러(200)가 빔포밍 벡터와 반사 위상을 조절하기 위해서는 송신기(100)와 K개의 단말(401 ~ 40K) 사이의 채널 상태 정보(channel state information: 이하 CSI)를 미리 알고 있어야만 한다.However, in order for the transmitter 100 and the controller 200 to adjust the beamforming vector and the reflection phase so that the power signal is efficiently transmitted to the K terminals 401 to 40K, the transmitter 100 and the K terminals 401 to 40K) Channel state information (hereinafter referred to as CSI) must be known in advance.

도 1과 같은 IRS(300)를 활용한 무선 전력 전송 시스템에서는 채널은 송신기(100)와 IRS(300) 사이의 채널(G(t))과 IRS(300)와 K개의 단말(401 ~ 40K) 사이의 채널(h(t))로 구분된다. 또한 IRS(300)와 K개의 단말(401 ~ 40K) 사이의 채널(h(t))은 IRS(300)와 K개의 단말(401 ~ 40K) 각각에 대한 K개의 채널(h(t) = {h1(t), h2(t), …, hk(t)})로 구분될 수 있다. 따라서 IRS를 활용한 무선 전력 전송 시스템에서 송신기(100)와 K개의 단말(401 ~ 40K) 사이의 CSI를 획득하는 것은 용이하지 않다.In the wireless power transmission system using the IRS 300 as shown in FIG. 1 , the channel is a channel (G(t)) between the transmitter 100 and the IRS 300 and the IRS 300 and K terminals (401 to 40K). It is separated by a channel h(t) between In addition, the channel (h(t)) between the IRS 300 and the K terminals 401 to 40K is K channels for each of the IRS 300 and the K terminals 401 to 40K (h(t) = { h 1 (t), h 2 (t), …, h k (t)}). Therefore, it is not easy to obtain CSI between the transmitter 100 and the K terminals 401 to 40K in a wireless power transmission system using IRS.

가장 단순한 방법으로는 송신기(100)가 IRS(300)의 N개의 반사 소자와 K개의 단말(401 ~ 40K) 모두에 대한 파일럿 신호를 개별 전송하여 채널을 추정하는 방법이 있다. 이 경우, 파일럿 신호의 길이가 반사 소자의 개수(N)와 단말의 개수(K)의 곱(NK) 이상이어야 하며, 이에 K개의 단말(401 ~ 40K) 각각은 NK 이상의 길이를 갖는 파일럿 신호에 대한 응신을 송신기(100)로 전송해야 한다. As the simplest method, there is a method in which the transmitter 100 individually transmits pilot signals for all of the N reflective elements of the IRS 300 and the K terminals 401 to 40K to estimate the channel. In this case, the length of the pilot signal must be equal to or greater than the product (NK) of the number of reflective elements (N) and the number of terminals (K). An acknowledgment should be transmitted to the transmitter 100 .

일반적으로 IRS(300)의 반사 소자 개수(N)는 16개 이상이므로 단말의 개수에 따라 파일럿 신호의 길이는 기하급수적으로 증가하게 된다. 따라서 K개의 단말(401 ~ 40K)가 파일럿 신호를 수신하고 이에 대한 응신을 송신기(100)로 전송하는 경우 전력 소모가 크게 발생한다.In general, since the number N of reflective elements of the IRS 300 is 16 or more, the length of the pilot signal increases exponentially according to the number of terminals. Accordingly, when the K terminals 401 to 40K receive the pilot signal and transmit a response thereto to the transmitter 100, power consumption is greatly increased.

전력 신호를 통해 무선으로 전력을 공급받는 K개의 단말(401 ~ 40K)은 전력 소비에 취약할 뿐만 아니라, 무선 전력 공급 효용성을 높이기 위해 길이가 긴 파일럿 신호를 송수신하는 것은 매우 비효율적이다.The K terminals 401 to 40K, which are wirelessly supplied with power through a power signal, are not only vulnerable to power consumption, but it is very inefficient to transmit and receive a long pilot signal to increase the effectiveness of wireless power supply.

따라서 IRS(300)를 활용한 무선 전력 전송 시스템에서 효율적으로 CSI를 추정할 수 있는 기술이 요구되고 있다.Therefore, there is a need for a technology capable of efficiently estimating CSI in a wireless power transmission system using the IRS 300 .

한국 등록 특허 제10-1678640호 (2016.11.16 등록)Korean Patent Registration No. 10-1678640 (Registered on November 16, 2016)

본 발명의 목적은 지능형 반사 표면을 이용한 무선 전력 전송 시스템에서 IRS의 반사 소자의 개수나 단말의 개수에 무관하게 수신 신호의 세기를 기반으로 용이하게 채널을 추정할 수 있는 무선 전력 전송 시스템의 채널 추정 장치 및 방법을 제공하는데 있다.It is an object of the present invention to estimate a channel of a wireless power transmission system that can easily estimate a channel based on the strength of a received signal regardless of the number of reflective elements or the number of terminals of the IRS in a wireless power transmission system using an intelligent reflective surface To provide an apparatus and method.

본 발명의 다른 목적은 지능형 반사 표면을 이용한 무선 전력 전송 시스템에서 저전력 IoT 단말의 작은 전력량으로도 채널을 추정할 수 있는 무선 전력 전송 시스템의 채널 추정 장치 및 방법을 제공하는데 있다.Another object of the present invention is to provide an apparatus and method for estimating a channel of a wireless power transmission system capable of estimating a channel even with a small amount of power of a low-power IoT terminal in a wireless power transmission system using an intelligent reflective surface.

상기 목적을 달성하기 위한 본 발명의 일 실시예에 따른 지능형 반사 표면을 이용한 무선 전력 전송 시스템의 채널 추정 장치는 무선 전력 전송 시스템의 송신기가 측정 신호를 전송하기 위해 생성한 측정 빔포밍 벡터와 상기 송신기에서 전송된 신호를 반사하여 다수의 단말로 전달하는 지능형 반사표면(이하 IRS)이 측정 신호를 반사하면서 조절하는 위상을 지정하기 위해 생성된 측정 위상 변화 행렬 및 상기 측정 신호를 수신한 다수의 단말 각각에서 피드백된 수신 신호 세기(이하 RSS)를 인가받고, 최대 우도 추정 기법에 따라 추정되어야 하는 상기 다수의 단말 각각에 수신된 측정 신호의 위상을 나타내는 위상 변수와 상기 IRS와 상기 다수의 단말 사이의 채널인 반사표면 채널 및 상기 송신기와 상기 IRS 사이의 채널인 송신 채널 중 반사표면 채널과 송신 채널을 고정하고, 고정된 반사표면 채널과 송신 채널과 인가된 측정 빔포밍 벡터와 측정 위상 변화 행렬 및 RSS를 기반으로 위상 변수를 계산하여 업데이트하는 위상 변수 업데이트부; 상기 송신 채널과 업데이트된 위상 변수를 기반으로 반사표면 채널을 계산하여 업데이트하는 반사표면 채널 업데이트부; 업데이트된 위상 변수와 반사표면 채널을 기반으로 연속 컨벡스 근사 기법에 따라 송신 채널을 획득하여 업데이트하는 송신 채널 업데이트부; 및 업데이트된 위상 변수와 반사표면 채널 및 송신 채널의 수렴 여부를 판별하여, 수렴되면 최종 업데이트된 반사표면 채널 및 송신 채널을 추정된 채널 정보로 획득하는 수렴 판별부를 포함한다.In order to achieve the above object, an apparatus for estimating a channel of a wireless power transmission system using an intelligent reflective surface according to an embodiment of the present invention includes a measurement beamforming vector generated by a transmitter of a wireless power transmission system to transmit a measurement signal and the transmitter A measurement phase change matrix generated to designate a phase that an intelligent reflective surface (hereinafter IRS) that reflects the signal transmitted from the and transmits it to a plurality of terminals reflects the measurement signal and adjusts the phase, and each of the plurality of terminals that received the measurement signal A channel between the IRS and the plurality of terminals and a phase variable indicating the phase of the measured signal received by each of the plurality of terminals to be applied with received signal strength (hereinafter referred to as RSS) fed back from and estimated according to the maximum likelihood estimation technique of the reflective surface channel and the transmission channel that is the channel between the transmitter and the IRS, the reflective surface channel and the transmission channel are fixed, and the fixed reflective surface channel and the transmission channel, the applied measurement beamforming vector, the measurement phase change matrix, and the RSS a phase variable update unit that calculates and updates the phase variable based on the calculation; a reflective surface channel update unit for calculating and updating a reflective surface channel based on the transmission channel and the updated phase variable; a transmission channel updater for acquiring and updating a transmission channel according to a continuous convex approximation technique based on the updated phase variable and the reflective surface channel; and a convergence determining unit for determining whether the updated phase variable and the reflective surface channel and the transmission channel converge, and obtaining the finally updated reflective surface channel and the transmission channel as estimated channel information when convergence.

상기 위상 변수 업데이트부는 측정 빔포밍 벡터와 측정 위상 변화 행렬 및 RSS를 인가받고, 반사표면 채널과 송신 채널과 인가된 측정 빔포밍 벡터와 측정 위상 변화 행렬 및 RSS를 기반으로 수학식 The phase variable update unit receives the measurement beamforming vector, the measurement phase change matrix, and the RSS, and an equation based on the reflective surface channel and the transmission channel, the applied measurement beamforming vector, the measurement phase change matrix, and the RSS

Figure pat00001
Figure pat00001

(여기서 hk 는 이전 기지정된 L개의 타임 슬롯 구간 동안 IRS와 다수의 단말 중 제k 단말 사이의 제k 반사표면 채널이고, φ k는 이전 L개의 타임 슬롯 구간 동안 K개의 단말에서 피드백된 RSS에 위상(

Figure pat00002
) 정보를 추가하기 위해 적용되는 위상 변수이며, M G 는 현지 타임 슬롯(T)로부터 이전 L개의 타임 슬롯 구간 동안 송신기로부터 IRS가 반사할 때까지의 측정 빔포밍 벡터(w)와 송신 채널(G) 및 측정 위상 변화 행렬(Θ)에 의한 송신 반사 채널로서 수학식
Figure pat00003
으로 계산된다.)에 따라 업데이트되는 위상 변수를 계산할 수 있다.(here, h k is the k-th reflective surface channel between the IRS and the k-th terminal among the plurality of terminals during the previously predetermined L time slot period, and φ k is the RSS fed back from the K terminals during the previous L time slot period. Phase(
Figure pat00002
) is a phase variable applied to add information, and M G is the measurement beamforming vector (w) and the transmission channel (G) from the local time slot (T) until the IRS reflects from the transmitter during the previous L time slot sections ) and the measured phase shift matrix ( Θ ) as the transmit reflection channel by the equation
Figure pat00003
Calculated as .), the updated phase variable can be calculated.

상기 반사표면 채널 업데이트부는 송신 채널과 업데이트된 위상 변수를 기반으로 수학식 The reflective surface channel update unit is an equation based on the transmission channel and the updated phase variable.

Figure pat00004
Figure pat00004

(여기서 βk는 이전 L개의 타임 슬롯 구간 동안 제k 단말의 속도에 따라 결정되는 채널 상관 계수이고, r k 는 이전 L개의 타임 슬롯 구간 동안 제k 단말로부터 피드백된 RSS이며, ㅀ는 하다마드 곱셈 기호이다. 그리고

Figure pat00005
는 이전 L개의 타임 슬롯 구간 동안 송신기와 제k 단말 사이의 가우시안 노이즈 벡터에 대한 공분산 행렬(R k)을 근사하여 획득되는 근사 공분산 행렬을 나타낸다.)에 따라 업데이트되는 반사표면 채널을 계산할 수 있다.(where β k is a channel correlation coefficient determined according to the speed of the k-th terminal during the previous L time slot sections, r k is the RSS fed back from the k-th terminal during the previous L time slot sections, and ㅀ is the Hadamard multiplication It is a sign, and
Figure pat00005
denotes an approximate covariance matrix obtained by approximating the covariance matrix ( R k ) for the Gaussian noise vector between the transmitter and the k-th terminal during the previous L time slot periods.) It is possible to calculate the updated reflective surface channel.

상기 반사표면 채널 업데이트부는 상기 근사 공분산 행렬을 수학식 The reflective surface channel updater calculates the approximate covariance matrix by Equation

Figure pat00006
Figure pat00006

(여기서, p와 q는 각각 공분산 행렬(R k)의 행과 열을 나타낸다.)에 따라 계산하는 가우시안 노이즈 벡터에 대한 공분산 행렬(R k)을 수학식 (Here, p and q represent the row and column of the covariance matrix ( R k ) , respectively. )

Figure pat00007
Figure pat00007

에 따라 근사하여 획득할 수 있다.It can be obtained by approximation according to

상기 송신 채널 업데이트부는 업데이트된 위상 변수와 업데이트된 반사표면 채널을 기반으로 수학식 The transmission channel update unit is based on the updated phase variable and the updated reflective surface channel Equation

Figure pat00008
Figure pat00008

(여기서

Figure pat00009
는 업데이트되기 이전 송신 채널이고,
Figure pat00010
는 역추적 라인 탐색법으로 획득 스텝 사이즈이며,
Figure pat00011
는 송신 채널(G)에 따른 반사표면 채널(h)의 기울기를 나타낸다.)에 따라 업데이트되는 송신 채널을 계산할 수 있다.(here
Figure pat00009
is the transmission channel before being updated,
Figure pat00010
is the step size acquired by the trace line search method,
Figure pat00011
denotes the slope of the reflective surface channel h according to the transmission channel G).

상기 송신 채널 업데이트부는 반사표면 채널(h)의 기울기(

Figure pat00012
)를 수학식 The transmission channel update unit slope (h) of the reflective surface channel (h)
Figure pat00012
) to the formula

Figure pat00013
Figure pat00013

(여기서 송신기와 제k 단말 사이의 가우시안 노이즈의 표준 편차를 나타낸다.)에 따라 계산할 수 있다.(Here, the standard deviation of Gaussian noise between the transmitter and the k-th terminal is indicated.).

상기 수렴 판별부는 업데이트된 위상 변수와 반사표면 채널 및 송신 채널의 누적 변화를 분석하여 위상 변수와 반사표면 채널 및 송신 채널과 이전 업데이트된 위상 변수와 반사표면 채널 및 송신 채널 사이의 편차가 기지정된 범위 이내로 수렴되는 것으로 판별되면, 최종 업데이트된 반사표면 채널 및 송신 채널을 추정된 채널 정보로 획득하고, 수렴되지 않는 것으로 판별되면, 위상 변수와 반사표면 채널 및 송신 채널이 반복 업데이트되도록 업데이트된 위상 변수와 반사표면 채널 및 송신 채널을 상기 위상 변수 업데이트부로 전달할 수 있다.The convergence determining unit analyzes the updated phase variable, the cumulative change of the reflective surface channel and the transmission channel, and the deviation between the phase variable, the reflective surface channel, and the transmission channel and the previously updated phase variable, the reflective surface channel, and the transmission channel is a predetermined range If it is determined to converge within A reflective surface channel and a transmission channel may be transmitted to the phase variable updater.

상기 위상 변수 업데이트부는 현재 타임 슬롯에서 가장 최근 업데이트된 위상 변수와 반사표면 채널 및 송신 채널을 이용하여 업데이트되는 위상 변수를 계산하고, 현재 타임 슬롯에서 이전 업데이트된 위상 변수와 반사표면 채널 및 송신 채널이 존재하지 않으면, 이전 타임 슬롯에서 최종 업데이트된 위상 변수와 반사표면 채널 및 송신 채널을 이용하여 업데이트되는 위상 변수를 계산할 수 있다.The phase variable update unit calculates the updated phase variable using the most recently updated phase variable, the reflective surface channel and the transmission channel in the current time slot, and the phase variable, the reflective surface channel and the transmission channel updated previously in the current time slot are If it does not exist, it is possible to calculate the updated phase variable using the last updated phase variable in the previous time slot and the reflective surface channel and transmit channel.

상기 측정 빔포밍 벡터는 상기 송신기에 의해 평균이 0이고, 평균 측정 전력과 송신기 안테나 개수의 비에 따른 분산을 갖는 복소 가우시안 분포를 따르도록 임의로 생성될 수 있다.The measurement beamforming vector may be arbitrarily generated by the transmitter to follow a complex Gaussian distribution having an average of 0 and a variance according to a ratio of the average measured power and the number of antennas of the transmitter.

상기 측정 위상 변화 행렬은 상기 IRS를 제어하는 컨트롤러에 의해 IRS의 N개의 반사 소자에 따른 측정 신호의 위상 변화가 시간에 독립적으로 0 ~ 2π의 범위에서 균일하게 분포하도록 임의로 생성될 수 있다.The measurement phase change matrix may be arbitrarily generated so that the phase change of the measurement signal according to the N reflective elements of the IRS is uniformly distributed in the range of 0 to 2π independently of time by the controller controlling the IRS.

상기 목적을 달성하기 위한 본 발명의 다른 실시예에 따른 무선 전력 전송 시스템의 채널 추정 방법은 컨트롤러의 제어에 따라 송신기에서 전송된 무선 전력 신호를 반사하여 다수의 단말로 전달하는 지능형 반사표면(이하 IRS)을 포함하는 무선 전력 전송 시스템에 있어서, 상기 송신기가 측정 신호를 전송하기 위해 생성한 측정 빔포밍 벡터와 상기 IRS가 측정 신호를 반사하면서 조절하는 위상을 지정하기 위해 생성된 측정 위상 변화 행렬을 획득하는 단계; 상기 측정 신호를 수신한 다수의 단말 각각에서 피드백된 RSS를 인가받는 단계; 최대 우도 추정 기법에 따라 추정되어야 하는 상기 다수의 단말 각각에 수신된 측정 신호의 위상을 나타내는 위상 변수와 상기 IRS와 상기 다수의 단말 사이의 채널인 반사표면 채널 및 상기 송신기와 상기 IRS 사이의 채널인 송신 채널 중 반사표면 채널과 송신 채널을 고정하고, 고정된 반사표면 채널과 송신 채널과 인가된 측정 빔포밍 벡터와 측정 위상 변화 행렬 및 RSS를 기반으로 위상 변수를 계산하여 업데이트하는 단계; 상기 송신 채널과 업데이트된 위상 변수를 기반으로 반사표면 채널을 계산하여 업데이트하는 단계; 업데이트된 위상 변수와 반사표면 채널을 기반으로 연속 컨벡스 근사 기법에 따라 송신 채널을 획득하여 업데이트하는 단계; 및 업데이트된 위상 변수와 반사표면 채널 및 송신 채널의 수렴 여부를 판별하여, 수렴되면 최종 업데이트된 반사표면 채널 및 송신 채널을 추정된 채널 정보로 획득하는 단계를 포함한다.The channel estimation method of the wireless power transmission system according to another embodiment of the present invention for achieving the above object is an intelligent reflective surface (hereinafter IRS) that reflects the wireless power signal transmitted from the transmitter under the control of the controller and transmits it to a plurality of terminals ) In a wireless power transmission system comprising: a measurement beamforming vector generated by the transmitter to transmit a measurement signal and a measurement phase change matrix generated to specify a phase that the IRS adjusts while reflecting the measurement signal to do; receiving RSS fed back from each of a plurality of terminals receiving the measurement signal; A phase variable indicating the phase of the measurement signal received in each of the plurality of terminals to be estimated according to the maximum likelihood estimation technique, a reflection surface channel that is a channel between the IRS and the plurality of terminals, and a channel between the transmitter and the IRS fixing a reflective surface channel and a transmission channel among transmission channels, calculating and updating a phase variable based on the fixed reflective surface channel and transmission channel, an applied measurement beamforming vector, a measurement phase change matrix, and RSS; calculating and updating a reflective surface channel based on the transmission channel and the updated phase variable; acquiring and updating a transmission channel according to a continuous convex approximation technique based on the updated phase variable and the reflective surface channel; and determining whether the updated phase variable and the reflective surface channel and the transmission channel converge, and when the convergence is achieved, the finally updated reflective surface channel and the transmission channel are acquired as estimated channel information.

따라서, 본 발명의 실시예에 따른 지능형 반사 표면을 이용한 무선 전력 전송 시스템의 채널 추정 장치 및 방법은 IRS의 반사 소자의 개수나 단말의 개수에 무관하게 송신기에서 전송되는 측정 신호에 대해 단말이 수신한 수신 신호의 세기를 확인하고, 확인된 수신 신호의 세기를 기반으로 송신기와 IRS 사이의 채널과 IRS와 단말 사이의 채널을 추정할 수 있으므로 저전력 단말을 포함하는 무선 전력 전송 시스템에서도 효율적으로 채널 정보를 획득할 수 있다.Therefore, the channel estimation apparatus and method of a wireless power transmission system using an intelligent reflective surface according to an embodiment of the present invention is a measurement signal transmitted from a transmitter regardless of the number of reflective elements of the IRS or the number of terminals. Since the strength of the received signal can be checked and the channel between the transmitter and the IRS and the channel between the IRS and the terminal can be estimated based on the strength of the received signal, it is possible to efficiently transmit channel information even in a wireless power transmission system including a low-power terminal. can be obtained

도 1은 지능형 반사 표면을 활용한 무선 전력 전송 시스템을 개념을 설명하기 위한 도면이다.
도 2는 본 발명의 일 실시예에 따른 지능형 반사 표면을 이용한 무선 전력 전송 시스템의 채널 추정 장치의 개략적 구조를 나타낸다.
도 3은 도 2에 도시된 무선 전력 전송 시스템의 채널 추정 장치가 채널 추정을 수행하는 동작을 설명하기 위한 도면이다.
도 4는 본 발명의 일 실시예에 따른 지능형 반사 표면을 이용한 무선 전력 전송 시스템의 채널 추정 방법을 나타낸다.1 is a diagram for explaining the concept of a wireless power transmission system utilizing an intelligent reflective surface.
2 shows a schematic structure of an apparatus for estimating a channel of a wireless power transmission system using an intelligent reflective surface according to an embodiment of the present invention.
FIG. 3 is a diagram for explaining an operation in which the channel estimator of the wireless power transmission system shown in FIG. 2 performs channel estimation.
4 shows a channel estimation method of a wireless power transmission system using an intelligent reflective surface according to an embodiment of the present invention.

본 발명과 본 발명의 동작상의 이점 및 본 발명의 실시에 의하여 달성되는 목적을 충분히 이해하기 위해서는 본 발명의 바람직한 실시예를 예시하는 첨부 도면 및 첨부 도면에 기재된 내용을 참조하여야만 한다. In order to fully understand the present invention, the operational advantages of the present invention, and the objects achieved by the practice of the present invention, reference should be made to the accompanying drawings illustrating preferred embodiments of the present invention and the contents described in the accompanying drawings.

이하, 첨부한 도면을 참조하여 본 발명의 바람직한 실시예를 설명함으로써, 본 발명을 상세히 설명한다. 그러나, 본 발명은 여러 가지 상이한 형태로 구현될 수 있으며, 설명하는 실시예에 한정되는 것이 아니다. 그리고, 본 발명을 명확하게 설명하기 위하여 설명과 관계없는 부분은 생략되며, 도면의 동일한 참조부호는 동일한 부재임을 나타낸다. Hereinafter, the present invention will be described in detail by describing preferred embodiments of the present invention with reference to the accompanying drawings. However, the present invention may be embodied in various different forms, and is not limited to the described embodiments. In addition, in order to clearly explain the present invention, parts irrelevant to the description are omitted, and the same reference numerals in the drawings indicate the same members.

명세서 전체에서, 어떤 부분이 어떤 구성요소를 "포함"한다고 할 때, 이는 특별히 반대되는 기재가 없는 한 다른 구성요소를 제외하는 것이 아니라, 다른 구성요소를 더 포함할 수 있는 것을 의미한다. 또한, 명세서에 기재된 "...부", "...기", "모듈", "블록" 등의 용어는 적어도 하나의 기능이나 동작을 처리하는 단위를 의미하며, 이는 하드웨어나 소프트웨어 또는 하드웨어 및 소프트웨어의 결합으로 구현될 수 있다. Throughout the specification, when a part "includes" a certain component, it does not exclude other components, unless otherwise stated, meaning that other components may be further included. In addition, terms such as "...unit", "...group", "module", and "block" described in the specification mean a unit that processes at least one function or operation, which is hardware, software, or hardware. and a combination of software.

도 2는 본 발명의 일 실시예에 따른 지능형 반사 표면을 이용한 무선 전력 전송 시스템의 채널 추정 장치의 개략적 구조를 나타내고, 도 3은 도 2에 도시된 무선 전력 전송 시스템의 채널 추정 장치가 채널 추정을 수행하는 동작을 설명하기 위한 도면이다.2 shows a schematic structure of a channel estimating apparatus of a wireless power transmission system using an intelligent reflective surface according to an embodiment of the present invention, and FIG. 3 is a channel estimation apparatus of the wireless power transmission system shown in FIG. It is a diagram for explaining the operation to be performed.

본 실시예에 따른 채널 추정 장치는 도 1의 무선 전력 전송 시스템에서 송신기(100)에 포함되어 구성될 수 있으며, 경우에 따라서는 컨트롤러(200)에 포함되어 구성될 수도 있다. 여기서는 채널 추정 장치가 송신기(100)에 포함된 경우를 가정하여 설명한다.The channel estimation apparatus according to the present embodiment may be configured to be included in the transmitter 100 in the wireless power transmission system of FIG. 1 , and in some cases may be included and configured to be included in the controller 200 . Here, it is assumed that the channel estimation apparatus is included in the transmitter 100 .

도 1 내지 도 3을 참조하면, 본 실시예에 따른 채널 추정 장치는 측정 신호 생성부(110), 수신 세기 획득부(120) 및 채널 추정부(130)를 포함할 수 있다.1 to 3 , the channel estimator according to the present embodiment may include a measurement signal generator 110 , a reception strength acquirer 120 , and a channel estimator 130 .

측정 신호 생성부(110)는 채널 추정을 위해 송신기(100)로부터 K개의 단말(401 ~ 40K)로 전송할 측정 신호를 생성한다. 여기서 측정 신호는 K개의 단말(401 ~ 40K)에서 수신되는 수신 신호의 세기를 측정하기 위한 기지정된 신호로서 단위 신호 또는 파일럿 신호일 수 있다.The measurement signal generator 110 generates measurement signals to be transmitted from the transmitter 100 to the K terminals 401 to 40K for channel estimation. Here, the measurement signal is a predetermined signal for measuring the strength of the received signal received by the K terminals 401 to 40K, and may be a unit signal or a pilot signal.

이때 측정 신호 생성부(110)는 M개의 안테나를 구비하는 송신기(100)가 제t 타임 슬롯(t)에서 측정 신호를 전송하기 위한 측정 빔포밍 벡터(w(t))를 임의로 생성한다. 여기서 측정 빔포밍 벡터(w(t))는 평균이 0이고, 송신기(100)의 평균 측정 전력(p0)과 안테나 개수(M)의 비(p0/M)에 따른 분산을 갖는 복소 가우시안 분포를 따르도록 생성할 수 있다.In this case, the measurement signal generator 110 arbitrarily generates a measurement beamforming vector w(t) for the transmitter 100 having M antennas to transmit a measurement signal in the t-th time slot t. Here, the measurement beamforming vector w(t) has an average of 0, and is a complex Gaussian having a variance according to the ratio (p 0 /M) of the average measured power (p 0 ) of the transmitter 100 to the number of antennas (M). It can be created to follow a distribution.

또한 측정 신호 생성부(110)는 제t 타임 슬롯(t)에서 IRS(300)의 N개의 반사 소자가 송신기(100)에서 전송된 측정 신호를 반사하여 위상 변화시키는 수준을 나타내는 측정 위상 변화 행렬(Θ(t))을 임의로 생성한다. 여기서 N개의 반사 소자를 구비하는 IRS(300)의 위상 변화 행렬(Θ(t))은 N개의 반사 소자 각각의 위상 변화(

Figure pat00014
)를 원소로 포함하는 대각 행렬(
Figure pat00015
)로 표현될 수 있으며, 측정 신호 생성부(110)는 측정 위상 변화(
Figure pat00016
)가 시간에 독립적으로 0 ~ 2π의 범위에서 균일하게 분포하도록 임의로 생성할 수 있다.In addition, the measurement signal generation unit 110 is a measurement phase change matrix ( Θ(t)) is randomly generated. Here, the phase change matrix Θ(t) of the IRS 300 including N reflective elements is the phase change (
Figure pat00014
) as an element of a diagonal matrix (
Figure pat00015
), and the measurement signal generating unit 110 measures the phase change (
Figure pat00016
) can be arbitrarily generated so that it is uniformly distributed in the range of 0 to 2π independently of time.

송신기(100)는 생성된 빔포밍 벡터(w(t))와 위상 변화 행렬(Θ(t))을 기반으로 측정 신호를 전송한다.The transmitter 100 transmits a measurement signal based on the generated beamforming vector w(t) and the phase change matrix Θ(t).

송신기(100)에서 전송된 측정 신호는 송신기(100)와 IRS(300) 사이의 채널(G(t))을 통해 IRS(300)에 도달하고, IRS(300)의 N개의 반사 소자에 의해 반사되어 IRS(300)와 K개의 단말(401 ~ 40K) 사이의 채널(h(t))을 통해 K개의 단말(401 ~ 40K)로 전송된다.The measurement signal transmitted from the transmitter 100 arrives at the IRS 300 through the channel G(t)) between the transmitter 100 and the IRS 300 , and is reflected by N reflective elements of the IRS 300 . and transmitted to the K terminals 401 to 40K through the channel h(t) between the IRS 300 and the K terminals 401 to 40K.

이하에서는 설명의 편의를 위하여, 송신기(100)와 IRS(300) 사이의 채널(G(t))을 송신 채널이라 하고, IRS(300)와 K개의 단말(401 ~ 40K) 사이의 채널(h(t))을 반사표면 채널이라 한다.Hereinafter, for convenience of description, the channel G(t) between the transmitter 100 and the IRS 300 is referred to as a transmission channel, and the channel h between the IRS 300 and the K terminals 401 to 40K (t)) is called the reflection surface channel.

이에 K개의 단말(401 ~ 40K) 중 반사표면 채널(h(t))을 통해 전송된 측정 신호를 수신 신호로 수신한 단말은 수신 신호 세기(Received Signal Strength: 이하 RSS)를 측정하고, 측정된 RSS를 다시 송신기(100)로 피드백한다. 본 실시예에서 K개의 단말(401 ~ 40K)은 전송된 측정 신호로부터 측정되는 RSS만을 피드백할 뿐, 다른 정보는 전송하지 않는다. 즉 스칼라 값인 RSS만을 회신하므로, 추가 정보를 전송하기 위한 전력 소비가 없어 전력 소모를 저감할 수 있다. 이는 제한된 전력으로 구동되는 IoT 장치에서는 매우 중요한 이슈이다.Accordingly, the terminal receiving the measurement signal transmitted through the reflective surface channel (h(t)) as the received signal among the K terminals (401 ~ 40K) measures the received signal strength (hereinafter referred to as RSS), and the measured The RSS is fed back to the transmitter 100 . In this embodiment, the K terminals 401 to 40K feed back only the RSS measured from the transmitted measurement signal, and do not transmit other information. That is, since only RSS, which is a scalar value, is returned, there is no power consumption for transmitting additional information, so power consumption can be reduced. This is a very important issue for IoT devices powered by limited power.

수신 세기 획득부(120)는 K개의 단말(401 ~ 40K)에서 피드백 되는 RSS를 인가받아 채널 추정부(130)로 전달한다.The reception strength obtaining unit 120 receives RSS fed back from the K terminals 401 to 40K and transmits it to the channel estimator 130 .

채널 추정부(130)는 수신 세기 획득부(120)에서 전달된 RSS와 측정 신호 생성부(110)에서 생성된 측정 빔포밍 벡터(w(t))와 측정 위상 변화 행렬(Θ(t))을 인가받아 송신 채널(G(t))과 반사표면 채널(h(t))을 추정한다.The channel estimator 130 includes the RSS transmitted from the reception intensity acquisition unit 120 , the measurement beamforming vector w(t) generated by the measurement signal generator 110 , and the measurement phase change matrix Θ(t)) is applied and the transmission channel (G(t)) and the reflection surface channel (h(t)) are estimated.

채널 추정부(130)가 채널을 추정하기 위해서는 우선 RSS를 살펴보아야 한다.In order for the channel estimator 130 to estimate the channel, first, RSS should be looked at.

송신기(100)에서 측정 빔포밍 벡터(w(t))와 측정 위상 변화 행렬(Θ(t))에 따라 전송된 측정 신호를 K개의 단말(401 ~ 40K) 중 제k 단말(40k)이 수신한 경우를 고려하면, 제k 단말(40k)의 수신 신호는 송신 채널(G(t))을 통해 IRS(300)에 전달된 후, IRS(300)과 제k 단말(40k) 사이의 제k 반사표면 채널(hk(t))을 통해 전달된 측정 신호이다.The k-th terminal 40k among the K terminals 401 to 40K receives the measurement signal transmitted from the transmitter 100 according to the measurement beamforming vector w(t) and the measurement phase change matrix Θ(t)) Considering one case, the received signal of the k-th terminal 40k is transmitted to the IRS 300 through the transmission channel G(t), and then the k-th terminal between the IRS 300 and the k-th terminal 40k. It is the measurement signal transmitted through the reflective surface channel h k (t).

따라서 제t 타임 슬롯(t)에서 제k 단말(40k)은 RSS(rk(t))를 수신된 수신 신호 전력의 제곱근으로 수학식 1로 계산하여 피드백할 수 있다.Accordingly, in the t-th time slot t, the k-th terminal 40k may calculate RSS(r k (t)) as the square root of the received received signal power by Equation 1 and feed it back.

Figure pat00017
Figure pat00017

여기서 nk(t)는 평균이 0이고 분산이 σk 2 인 가우시안 노이즈이며, H는 에르미트 행렬(Hermitian matrix)을 의미한다.Here, n k (t) is Gaussian noise with an average of 0 and a variance of σ k 2 , and H denotes a Hermitian matrix.

수학식 1에 나타난 바와 같이, 제k 단말(40k)의 무선 전력 수신 성능은 송신 채널(G(t))과 제k 반사표면 채널(hk(t)), 빔포밍 벡터(w(t)) 및 위상 변화 행렬(Θ(t))에 의해 결정됨을 알 수 있다.As shown in Equation 1, the wireless power reception performance of the k-th terminal 40k is a transmission channel (G(t)), a k-th reflective surface channel (h k (t)), and a beamforming vector (w(t)). ) and the phase change matrix Θ(t).

이중 빔포밍 벡터(w(t))와 위상 변화 행렬(Θ(t))은 각각 송신기(100)와 컨트롤러(200)에 의해 조절될 수 있으므로, 결과적으로 무선 전력 수신 성능은 송신 채널(G(t))과 제k 반사표면 채널(hk(t))에 의해 결정된다. 따라서 IRS(300)을 이용하여 요구되는 무선 전력 수신 성능을 획득하기 위해서는 송신 채널(G(t))과 제k 반사표면 채널(hk(t))에 대한 정보를 획득해야만 한다.Since the double beamforming vector w(t) and the phase change matrix Θ(t) can be adjusted by the transmitter 100 and the controller 200, respectively, as a result, the wireless power reception performance is reduced to the transmission channel G( t)) and the kth reflective surface channel h k (t). Therefore, in order to obtain the required wireless power reception performance using the IRS 300, information on the transmission channel G(t) and the k-th reflective surface channel h k (t) must be acquired.

한편, 본 실시예에서는 제k 반사표면 채널(hk(t))이 시간에 따라 상관 관계를 가지고 변화하는 1차 자기 회귀(first order auto-regressive) 채널 모델을 따르는 것으로 가정한다.Meanwhile, in the present embodiment, it is assumed that the k-th reflective surface channel h k (t) follows a first order auto-regressive channel model that has a correlation with time and changes.

제k 반사표면 채널(hk(t))이 1차 자기 회귀 채널 모델을 따르는 경우, 제t 타임 슬롯(t)에서의 제k 반사표면 채널(hk(t))은 이전 제t-1 타임 슬롯(t-1)에서의 제k 반사표면 채널(hk(t))과의 상관 관계 형식으로 수학식 2와 같이 표현될 수 있다.If the k th reflective surface channel h k (t) follows the first-order autoregressive channel model, then the k th reflective surface channel h k (t) in the t th time slot t is the previous t-1 th It can be expressed as Equation 2 in the form of correlation with the k-th reflective surface channel h k (t) in the time slot t-1.

Figure pat00018
Figure pat00018

여기서 βk는 βk ∈ (0, 1] 인 제k 반사표면 채널(hk(t))에 대한 시간적 페이딩 계수로서 제k 단말(40k)의 속도에 따라 결정되는 채널 상관 계수이며, ek(t)는 공간적 상관 관계없이 평균이 0이고 분산이 1-βk 2 인 가우시안 랜덤 벡터로 모델링되는 시변 잡음 항목이다.Here, β k is a temporal fading coefficient for the k-th reflective surface channel h k (t) with β k ∈ (0, 1], and is a channel correlation coefficient determined according to the speed of the k-th terminal 40k, e k (t) is a time-varying noise term modeled as a Gaussian random vector with mean 0 and variance 1-β k 2 , regardless of spatial correlation.

한편, 본 실시예에서 송신기(100)와 IRS(300)는 이동하지 않고 고정된 위치를 갖는 것으로 가정한다. 이에 송신기(100)와 IRS(300)의 상대 속도가 0이므로 송신 채널(G(t))에 대한 시간적 페이딩 계수(βk)는 1이고, 따라서 제t 타임 슬롯(t)에서의 송신 채널(G(t))과 제t-1 타임 슬롯(t-1)에서의 송신 채널(G(t-1))은 수학식 3과 같이 동일하게 나타난다.Meanwhile, in this embodiment, it is assumed that the transmitter 100 and the IRS 300 do not move and have fixed positions. Accordingly, since the relative speed of the transmitter 100 and the IRS 300 is 0, the temporal fading coefficient β k for the transmission channel G(t) is 1, and thus the transmission channel ( G(t)) and the transmission channel G(t-1) in the t-1 th time slot t-1 are the same as in Equation (3).

Figure pat00019
Figure pat00019

본 실시예에서 채널 추정부(130)는 최근 획득된 L개의 RSS(rk(t))를 이용하여 송신 채널(G(t))과 제k 반사표면 채널(hk(t))을 추정한다. 즉 현재 타임 슬롯이 T라고 할 때, 채널 추정부(130)는 슬라이딩 윈도우의 길이(L)에 따라 제T-L+1 타임 슬롯(T-L+1)부터 제T 타임 슬롯(T)까지 L개의 타임 슬롯 구간 동안 획득된 RSS(rk(T-L+1), …, rk(t))를 이용하여 송신 채널(G(t))과 제k 반사표면 채널(hk(t))을 추정한다.In this embodiment, the channel estimator 130 estimates the transmission channel G(t) and the k-th reflective surface channel h k (t) using the recently acquired L RSS(r k (t)). do. That is, when the current time slot is T, the channel estimator 130 extends from the T-L+1 time slot (T-L+1) to the T-th time slot (T) according to the length L of the sliding window. The transmission channel (G(t)) and the kth reflective surface channel (h k (t)) using RSS(r k (T-L+1), ..., r k (t)) obtained during L time slot sections )) is estimated.

제T 타임 슬롯(T)(여기서 t < T)의 제k 반사표면 채널(hk(t))은 수학식 2로부터 수학식 4와 같이 표현될 수 있다.The k-th reflective surface channel h k (t) of the T-th time slot T (where t < T) may be expressed as in Equation 4 from Equation (2).

Figure pat00020
Figure pat00020

그리고 수학식 4를 수학식 1에 대입하면, RSS(rk(t))는 수학식 5로 다시 표현될 수 있다.And, by substituting Equation 4 into Equation 1, RSS(r k (t)) can be expressed again as Equation 5.

Figure pat00021
Figure pat00021

여기서

Figure pat00022
이고, 수학식 3에 따라 G(t) = G(t-1)이므로, 수학식 5에서 타임 슬롯에 대한 색인은 생략하였다.here
Figure pat00022
, and since G(t) = G(t-1) according to Equation 3, an index for a time slot is omitted in Equation 5.

수학식 5에서는 절대값 함수(|·|)를 이용하여 위상에 대한 정보가 없는 RSS(rk(t))를 표현하였으므로, 수학식 5의 RSS(rk(t))에 위상 변수(φk(t))를 추가로 반영하면, 수학식 6과 같이 표현될 수 있다.In Equation 5, RSS(r k (t)) without phase information is expressed using the absolute value function (|·|), so the phase variable (φ) is added to RSS(r k (t)) of Equation 5 If k (t)) is additionally reflected, it can be expressed as Equation (6).

Figure pat00023
Figure pat00023

수학식 6으로부터 상기한 바와 같이, L개의 RSS(rk(t))로 이루어지는 RSS 벡터를 굵은 글씨(r k(t))로 RSS(rk(t))와 구분하여 표현하면, RSS 벡터(r k(t) = [rk(t-L+1), rk(t-L+2), …, rk(t)]T)는 수학식 7로 간략화될 수 있다.As described above from Equation 6, if an RSS vector composed of L RSS(r k (t)) is expressed in bold text ( r k (t)) separately from RSS(r k (t)), the RSS vector ( r k (t) = [r k (t-L+1), r k (t-L+2), ..., r k (t)] T ) may be simplified to Equation 7.

Figure pat00024
Figure pat00024

여기서 ㅀ는 하다마드 곱셈(Hadamard product) 기호이고, 굵은 글씨는 벡터를 나타내며,

Figure pat00025
는 평균이 0이고, 공분산 행렬이 R k(T)인 가우시안 노이즈 벡터로 볼 수 있다.where ㅀ is the Hadamard product symbol, bold text indicates a vector,
Figure pat00025
can be viewed as a Gaussian noise vector with a mean of 0 and a covariance matrix R k (T).

수학식 7에서 M G(T)는 송신기(100)의 빔포밍 벡터(w)와 송신 채널(G) 및 측정 위상 변화 행렬(Θ(t))이 결합되어, 송신기(100)로부터 전송된 측정 신호를 IRS(300)가 반사할 때까지의 채널 특성을 일괄하여 나타내는 송신 반사 채널 정보로 볼 수 있다.In Equation 7, M G (T) is the beamforming vector ( w ) of the transmitter 100 , the transmission channel ( G ), and the measurement phase change matrix (Θ(t)) are combined, and the measurement transmitted from the transmitter 100 It can be seen as transmission reflection channel information that collectively indicates the channel characteristics until the signal is reflected by the IRS 300 .

그리고 노이즈 벡터(

Figure pat00026
)의 공분산 행렬(R k(T))은 수학식 8과 같이 나타날 수 있다.and the noise vector (
Figure pat00026
) of the covariance matrix R k (T) can be expressed as in Equation 8.

Figure pat00027
Figure pat00027

여기서 p와 q는 각각 공분산 행렬(R k(T))의 행과 열을 나타낸다.where p and q denote the row and column of the covariance matrix R k (T), respectively.

한편 본 실시예에서는 송신 채널(G)과 K개의 단말(401 ~ 40K)에 대한 반사표면 채널(h =

Figure pat00028
)을 추정하기 위해 최대 우도 추정(maximum likelihood estimation: MLE) 기법을 사용한다. 다만 최대 우도 추정 기법을 적용하기 위해서는 추정되어야 하는 송신 채널(G)과 반사표면 채널(
Figure pat00029
) 사이의 확률 모델을 정의하는 과정이 필요하다.Meanwhile, in this embodiment, the transmission channel (G) and the reflection surface channel (h =
Figure pat00028
), a maximum likelihood estimation (MLE) technique is used. However, in order to apply the maximum likelihood estimation method, the transmission channel (G) and the reflective surface channel (
Figure pat00029
), a process of defining a probabilistic model between

이에 RSS 벡터(r k(t))와 송신 채널(G) 및 반사표면 채널(

Figure pat00030
) 사이의 결합 우도 함수(Joint likelihood function)를 수학식 7로부터 수학식 9와 같이 도출할 수 있다.Accordingly, the RSS vector ( r k (t)) and the transmit channel (G) and the reflective surface channel (
Figure pat00030
) can be derived from Equation 7 to Equation 9 as in Equation 9.

Figure pat00031
Figure pat00031

여기서 Ψ는 추정되어야 하는 추정 변수 집합으로 여기서는 송신 채널(G)과 반사표면 채널(

Figure pat00032
) 및 위상 변수(
Figure pat00033
)의 집합(
Figure pat00034
)을 나타낸다.where Ψ is the set of estimation variables to be estimated, here the transmission channel (G) and the reflective surface channel (
Figure pat00032
) and the phase variable (
Figure pat00033
) of the set (
Figure pat00034
) is indicated.

수학식 9를 이용하여 최대 우도 추정 기법을 적용하면 채널 추정 문제는 수학식 10과 같이 정리된다.When the maximum likelihood estimation technique is applied using Equation 9, the channel estimation problem is summarized as Equation 10.

Figure pat00035
Figure pat00035

수학식 10으로 정리된 채널 추정 문제를 여기서는 최대 우도 위상 탐색(Maximum Likelihood Phase Retrieval: 이하 MLPR) 문제라 한다. 수학식 10의 채널 추정 문제는 해를 도출하기 어려운 형태이므로, 우선 노이즈 벡터(

Figure pat00036
)의 공분산 행렬(R k(T))을 수학식 11과 같이 근사한다.The channel estimation problem summarized in Equation (10) is referred to herein as a Maximum Likelihood Phase Retrieval (MLPR) problem. Since the channel estimation problem of Equation 10 is difficult to derive a solution, first of all, the noise vector (
Figure pat00036
) of the covariance matrix R k (T) is approximated as in Equation 11.

Figure pat00037
Figure pat00037

수학식 11로 근사된 공분산 행렬(

Figure pat00038
)을 이용하면, 최대 우도 추정 기법에 따른 수학식 10의 MLPR 문제는 수학식 12와 같이 근사된 MLPR(Approximated MLPR: AMLPR) 문제 형태로 재정리될 수 있다.The covariance matrix approximated by Equation 11 (
Figure pat00038
), the MLPR problem of Equation 10 according to the maximum likelihood estimation technique can be rearranged in the form of an Approximated MLPR (AMLPR) problem as shown in Equation 12.

Figure pat00039
Figure pat00039

수학식 12와 같이, 정리된 AMLPR 문제는 추정 변수 집합(

Figure pat00040
)의 변수들인 송신 채널(G)과 반사표면 채널(
Figure pat00041
) 및 위상 변수(
Figure pat00042
) 중 두 변수를 고정하고 나머지 하나의 변수를 교대로 업데이트하는 과정을 반복하는 반복 알고리즘(block coordinate descent algorithm)을 이용하여 해결할 수 있다.As shown in Equation 12, the sorted AMLPR problem is a set of estimated variables (
Figure pat00040
), the transmission channel (G) and the reflective surface channel (
Figure pat00041
) and the phase variable (
Figure pat00042
), it can be solved by using a block coordinate descent algorithm that repeats the process of fixing two variables and updating the other variable alternately.

채널 추정부(130)는 우선 송신 채널(G)과 반사표면 채널(

Figure pat00043
) 및 위상 변수(
Figure pat00044
)를 초기값을 설정한다. 여기서 초기값은 이전 타임 슬롯(T-1)에서 추정된 추정 변수 집합(
Figure pat00045
)으로 설정될 수 있다. 이전 타임 슬롯(T-1)에서 추정된 추정 변수 집합(
Figure pat00046
)이 존재하지 않는 경우에는 추정 변수 집합(
Figure pat00047
)의 송신 채널(G)과 반사표면 채널(
Figure pat00048
) 및 위상 변수(
Figure pat00049
)의 초기값은 임의의 값으로 설정될 수 있다.The channel estimator 130 first comprises a transmission channel (G) and a reflective surface channel (
Figure pat00043
) and the phase variable (
Figure pat00044
) to set the initial value. Here, the initial value is the set of estimated variables estimated in the previous time slot (T-1) (
Figure pat00045
) can be set. The set of estimated variables estimated in the previous time slot (T-1) (
Figure pat00046
) does not exist, the set of estimated variables (
Figure pat00047
) of the transmit channel (G) and the reflective surface channel (
Figure pat00048
) and the phase variable (
Figure pat00049
) may be set to an arbitrary value.

추정 변수 집합(

Figure pat00050
)의 초기값이 설정되면, 송신 채널(G)과 반사표면 채널(
Figure pat00051
)의 값을 고정하고, 고정된 송신 채널(G)과 반사표면 채널(
Figure pat00052
)에 대응하는 최적의 위상 변수(
Figure pat00053
)를 우선 수학식 13에 따라 계산하여 업데이트한다.set of estimated variables (
Figure pat00050
) is set, the transmission channel (G) and the reflective surface channel (
Figure pat00051
) is fixed, and the fixed transmission channel (G) and the reflective surface channel (
Figure pat00052
) corresponding to the optimal phase variable (
Figure pat00053
) is first calculated and updated according to Equation 13.

Figure pat00054
Figure pat00054

여기서 arg(ㅇ) 함수는 벡터의 각 성분의 위상을 나타낸다.Here, the arg(o) function represents the phase of each component of the vector.

수학식 13에 따라 고정된 송신 채널(G)과 반사표면 채널(

Figure pat00055
)에 대응하는 위상 변수(
Figure pat00056
)가 계산되면, 송신 채널(G)과 계산된 위상 변수(
Figure pat00057
)가 고정된 상태에서 최적의 반사표면 채널(
Figure pat00058
)을 수학식 14에 따라 계산하여 업데이트한다.According to Equation 13, the fixed transmission channel (G) and the reflective surface channel (
Figure pat00055
) corresponding to the phase variable (
Figure pat00056
) is calculated, the transmission channel (G) and the calculated phase variable (
Figure pat00057
) is fixed and the optimal reflective surface channel (
Figure pat00058
) is calculated and updated according to Equation 14.

Figure pat00059
Figure pat00059

위상 변수(

Figure pat00060
)와 반사표면 채널(
Figure pat00061
)이 각각 수학식 13 및 수학식 14에 의해 계산되면, 계산된 위상 변수(
Figure pat00062
)와 반사표면 채널(
Figure pat00063
)을 기반으로 최적의 송신 채널(G)을 계산한다. 다만 최적의 송신 채널(G)을 획득하는 문제는 논컨벡스(non-convex) 문제로 연속 컨벡스 근사(Successive Convex Approximation) 기법을 적용하여 수학식 15에 따라 획득되어 업데이트 될 수 있다.phase variable (
Figure pat00060
) and the reflective surface channel (
Figure pat00061
) is calculated by Equations 13 and 14, respectively, then the calculated phase variable (
Figure pat00062
) and the reflective surface channel (
Figure pat00063
) to calculate the optimal transmission channel (G). However, the problem of obtaining the optimal transmission channel G is a non-convex problem and may be obtained and updated according to Equation 15 by applying a Successive Convex Approximation technique.

Figure pat00064
Figure pat00064

여기서

Figure pat00065
는 업데이트되기 이전 송신 채널을 나타내고,
Figure pat00066
는 스텝 사이즈로서 역추적 라인 탐색법(backtracking line search method)으로 획득될 수 있다. 그리고
Figure pat00067
는 송신 채널(G)에 따른 반사표면 채널(h)의 기울기(gradient)로서 수학식 16에 따라 계산될 수 있다.here
Figure pat00065
represents the transmission channel before being updated,
Figure pat00066
may be obtained by a backtracking line search method as a step size. And
Figure pat00067
is a gradient of the reflective surface channel h according to the transmission channel G, and may be calculated according to Equation 16.

Figure pat00068
Figure pat00068

채널 추정부(130)는 추정 변수 집합(

Figure pat00069
)에서 송신 채널(G)과 반사표면 채널(
Figure pat00070
) 및 위상 변수(
Figure pat00071
)를 하나씩 교대로 선택하여 업데이트를 수행한다.The channel estimator 130 is a set of estimated variables (
Figure pat00069
) in the transmit channel (G) and the reflective surface channel (
Figure pat00070
) and the phase variable (
Figure pat00071
) is selected one by one to perform the update.

그리고 반복하여 업데이트되는 추정 변수 집합(

Figure pat00072
)의 변화를 분석하여, 반복 업데이트되는 추정 변수 집합(
Figure pat00073
)이 기지정된 범위 이내의 편차로 수렴되는지 판별하고, 수렴되지 않으면 다시 반복하여 송신 채널(G)과 반사표면 채널(
Figure pat00074
) 및 위상 변수(
Figure pat00075
)를 하나씩 교대로 선택하여 수학식 13 내지 15에 따라 업데이트한다. 그러나 수렴하는 것으로 판별되면, 업데이트를 종료하고, 마지막으로 업데이트된 추정 변수 집합(
Figure pat00076
)의 송신 채널(G)과 반사표면 채널(
Figure pat00077
)를 추정 채널 정보로 설정한다.and an iteratively updated set of estimated variables (
Figure pat00072
) by analyzing the change in the iteratively updated estimated variable set (
Figure pat00073
) converges with a deviation within a predetermined range, and if it does not converge, it is repeated again for the transmission channel (G) and the reflective surface channel (
Figure pat00074
) and the phase variable (
Figure pat00075
) are alternately selected one by one and updated according to Equations 13 to 15. However, if it is determined to converge, the update ends, and the last updated set of estimated variables (
Figure pat00076
) of the transmit channel (G) and the reflective surface channel (
Figure pat00077
) is set as estimated channel information.

이에 송신기(100)는 추정 채널 정보에 따라 빔포밍 벡터(w(t))를 생성하여 무선 전력 신호를 전송한다. 그리고 IRS(300)를 제어하는 컨트롤러(200)는 추정 채널 정보에 따라 위상 변화 행렬(Θ(t))을 생성하여, IRS(300)의 N개의 반사 소자의 반사 위상을 조절함으로써, 송신기(100)에서 전송된 무선 전력 신호가 효율적으로 K개의 단말(401 ~ 40K)로 전달되도록 한다.Accordingly, the transmitter 100 transmits a wireless power signal by generating a beamforming vector w(t) according to the estimated channel information. In addition, the controller 200 controlling the IRS 300 generates a phase change matrix Θ(t) according to the estimated channel information, and adjusts the reflection phases of the N reflective elements of the IRS 300, so that the transmitter 100 ) so that the transmitted wireless power signal is efficiently transmitted to the K terminals (401 ~ 40K).

다시 도 2를 참조하면, 채널 추정부(130)는 일 예로 위상 변수 업데이트부(131), 반사표면 채널 업데이트부(132), 송신 채널 업데이트부(133) 및 수렴 판별부(134)를 포함할 수 있다.Referring back to FIG. 2 , the channel estimator 130 may include, for example, a phase variable update unit 131 , a reflective surface channel update unit 132 , a transmission channel update unit 133 , and a convergence determining unit 134 . can

위상 변수 업데이트부(131)는 측정 신호 생성부(110)로부터 빔포밍 벡터(w(t))와 위상 변화 행렬(Θ(t))을 인가받고, 수신 신호 세기 획득부(120)로부터 RSS(rk(t))를 인가받아, 이전 획득된 추정 변수 집합(

Figure pat00078
)의 송신 채널(G)과 반사표면 채널(
Figure pat00079
)를 이용하여 위상 변수(
Figure pat00080
)를 수학식 13에 따라 계산하여 업데이트한다.The phase variable update unit 131 receives the beamforming vector w(t) and the phase change matrix Θ(t) from the measurement signal generation unit 110 , and receives the RSS ( r k (t)) is authorized, and the previously obtained set of estimated variables (
Figure pat00078
) of the transmit channel (G) and the reflective surface channel (
Figure pat00079
) using the phase variable (
Figure pat00080
) is calculated and updated according to Equation 13.

그리고 반사표면 채널 업데이트부(132)는 추정 변수 집합(

Figure pat00081
)의 송신 채널(G)과 위상 변수 업데이트부(131)에서 업데이트된 위상 변수(
Figure pat00082
)를 이용하여 수학식 14에 따라 반사표면 채널(
Figure pat00083
)을 계산하여 업데이트 한다.And the reflective surface channel updater 132 is an estimated variable set (
Figure pat00081
) of the transmission channel (G) and the phase variable (
Figure pat00082
) using the reflective surface channel (
Figure pat00083
) is calculated and updated.

한편 송신 채널 업데이트부(133)는 위상 변수 업데이트부(131)와 반사표면 채널 업데이트부(132) 각각에서 업데이트된 위상 변수(

Figure pat00084
)와 반사표면 채널(
Figure pat00085
)을 이용하여 수학식 15에 따라 송신 채널(G)을 획득하여 업데이트 한다.On the other hand, the transmission channel updater 133 updates the phase variable (
Figure pat00084
) and the reflective surface channel (
Figure pat00085
) to obtain and update the transmission channel G according to Equation 15.

추정 변수 집합(

Figure pat00086
)의 모든 변수가 업데이트되면, 수렴 판별부(134)는 업데이트된 추정 변수 집합(
Figure pat00087
)과 이전 추정 변수 집합을 비교하여, 추정 변수 집합(
Figure pat00088
)의 각 변수값의 변화가 이전 변화에 비해 점차 감소하는지 판별한다.set of estimated variables (
Figure pat00086
When all the variables of ) are updated, the convergence determining unit 134 sets the updated estimated variable (
Figure pat00087
) and the previous set of estimated variables,
Figure pat00088
), it is determined whether the change in the value of each variable is gradually decreasing compared to the previous change.

수렴 판별부(134)는 변수값의 변화로부터 각 변수의 값이 기지정된 범위 이내로 수렴되는 것으로 판별되면, 최종 업데이트된 추정 변수 집합(

Figure pat00089
)의 송신 채널(G)과 반사표면 채널(
Figure pat00090
)를 추정 채널 정보로 설정한다. 그러나 변수값의 변화가 감소되지 않거나, 감소되더라도 기지정된 범위 이내로 수렴되지 않는 것으로 판별되면, 업데이트된 추정 변수 집합(
Figure pat00091
)을 위상 변수 업데이트부(131)로 전달하여, 추정 변수 집합(
Figure pat00092
)이 다시 반복하여 업데이트 되도록 한다. 즉 수렴 판별부(134)는 반복 업데이트되는 추정 변수 집합(
Figure pat00093
)의 누적 변화를 분석하여 특정값으로 수렴되는 것으로 판별되는 경우, 업데이트를 중지하고 최종 업데이트된 추정 변수 집합(
Figure pat00094
)을 추정된 채널 정보로 획득한다.When it is determined that the value of each variable converges within a predetermined range from the change in the variable value, the convergence determining unit 134 determines the final updated estimated variable set (
Figure pat00089
) of the transmit channel (G) and the reflective surface channel (
Figure pat00090
) is set as estimated channel information. However, if it is determined that the change in the variable value does not decrease or does not converge within a predetermined range even if it is reduced, the updated estimated variable set (
Figure pat00091
) is transmitted to the phase variable update unit 131, and the estimated variable set (
Figure pat00092
) is repeated again to be updated. That is, the convergence determining unit 134 is an iteratively updated estimated variable set (
Figure pat00093
), if it is determined to converge to a specific value by analyzing the cumulative change in
Figure pat00094
) as estimated channel information.

결과적으로 본 실시예의 지능형 반사 표면을 이용한 무선 전력 전송 시스템의 채널 추정 장치는 K개의 단말(401 ~ 40K)로부터 스칼라값이 RSS(rk(t))만을 피드백 받아도 지능형 반사 표면으로 인해 구분되는 송신 채널(G)과 반사표면 채널(

Figure pat00095
)을 모두 추정할 수 있으므로, 제한된 전력으로 동작하는 단말(401 ~ 40K)의 전력 소비를 저감시킬 수 있다.As a result, the channel estimation apparatus of the wireless power transmission system using the intelligent reflective surface of this embodiment receives only the RSS (r k (t)) of the scalar value from the K terminals 401 to 40K, and the transmission is differentiated due to the intelligent reflective surface. Channel (G) and the reflective surface channel (
Figure pat00095
) can be estimated, so power consumption of the terminals 401 to 40K operating with limited power can be reduced.

상기에서는 설명의 편의를 위해 측정 신호 생성부(110)와 수신 세기 획득부(120)를 포함하는 것으로 설명하였으나, 이는 채널 추정 장치가 송신기(100)에 포함되는 경우를 가정한 것으로서, 본 실시예에 따른 채널 추정 장치는 채널 추정부(130)로 구성될 수 있다.In the above description, it has been described that the measurement signal generator 110 and the reception strength acquirer 120 are included for convenience of explanation, but it is assumed that the channel estimator is included in the transmitter 100 , and this embodiment The channel estimating apparatus according to , may include the channel estimator 130 .

도 4는 본 발명의 일 실시예에 따른 지능형 반사 표면을 이용한 무선 전력 전송 시스템의 채널 추정 방법을 나타낸다.4 shows a channel estimation method of a wireless power transmission system using an intelligent reflective surface according to an embodiment of the present invention.

도 2 및 도 3을 참조하여, 도 4의 채널 추정 방법을 설명하면, 우선 송신기(100)가 측정 신호를 송신하기 위한 측정 빔포밍 벡터(w(t))와 IRS(300)가 전송된 측정 신호의 반사 시에 변경하는 위상를 설정하기 위한 측정 위상 변화 행렬(Θ(t))을 생성한다(S10). 여기서 빔포밍 벡터(w(t))와 위상 변화 행렬(Θ(t))은 임의로 생성될 수 있다.Referring to FIGS. 2 and 3 , the channel estimation method of FIG. 4 will be described. First, the transmitter 100 transmits the measurement beamforming vector w(t) for transmitting the measurement signal and the IRS 300 is transmitted. A measurement phase change matrix Θ(t) for setting a phase to change upon reflection of a signal is generated (S10). Here, the beamforming vector w(t) and the phase change matrix Θ(t) may be arbitrarily generated.

측정 빔포밍 벡터(w(t))와 측정 위상 변화 행렬(Θ(t))이 생성되면, 송신기(100)는 생성된 측정 빔포밍 벡터(w(t))에 따라 측정 신호를 빔포밍하여 전송하고, IRS(300)는 측정 위상 변화 행렬(Θ(t))에 의해 지정된 반사 위상으로 송신기(100)에서 전송된 측정 신호의 위상을 조절하여 반사한다(S20).When the measurement beamforming vector w(t) and the measurement phase change matrix Θ(t) are generated, the transmitter 100 beamforms the measurement signal according to the generated measurement beamforming vector w(t). and the IRS 300 adjusts the phase of the measurement signal transmitted from the transmitter 100 to the reflection phase specified by the measurement phase change matrix Θ(t) and reflects it (S20).

그리고 무선 전력 전송 시스템의 다수의 단말 중 IRS(300)에 반사되어 전송된 측정 신호를 수신한 단말로부터 피드백되는 RSS(rk(t))를 확인한다(S30). 여기서 RSS(rk(t))는 측정 신호를 수신한 각 단말이 수학식 1에 따라 수신 신호 전력의 제곱근으로 계산하여 획득하여 송신기(100)로 피드백 할 수 있다.And check the RSS (r k (t)) fed back from the terminal receiving the transmitted measurement signal reflected by the IRS 300 among a plurality of terminals of the wireless power transmission system (S30). Here, RSS(r k (t)) may be obtained by each terminal receiving the measurement signal calculated as the square root of the received signal power according to Equation 1, and may be fed back to the transmitter 100 .

이후, 추정되어야 하는 추정 변수인 송신 채널(G)과 반사표면 채널(

Figure pat00096
) 및 위상 변수(
Figure pat00097
)의 초기값을 설정한다(S40). 여기서 초기값은 이전 마지막으로 추정된 추정 변수들의 값으로 설정될 수 있으며, 이전 설정된 추정 변수들이 존재하지 않는다면 임의의 값으로 설정될 수 있다.Then, the transmission channel (G) and the reflective surface channel (
Figure pat00096
) and the phase variable (
Figure pat00097
) to set an initial value (S40). Here, the initial value may be set as the values of the previously estimated estimated variables, or may be set to an arbitrary value if the previously set estimated variables do not exist.

측정 빔포밍 벡터(w(t))와 측정 위상 변화 행렬(Θ(t)) 및 RSS(rk(t))가 확인되고, 추정 변수인 송신 채널(G)과 반사표면 채널(

Figure pat00098
) 및 위상 변수(
Figure pat00099
)의 초기값이 설정되면, 3개의 추정 변수 중 우선 송신 채널(G)과 반사표면 채널(
Figure pat00100
)을 현재 설정된 값으로 고정한 상태에서의 최적의 위상 변수(
Figure pat00101
)를 수학식 13에 따라 계산하여 업데이트한다(S50).The measurement beamforming vector (w(t)), the measurement phase change matrix (Θ(t)), and the RSS(r k (t)) are confirmed, and the transmission channel (G) and the reflective surface channel (
Figure pat00098
) and the phase variable (
Figure pat00099
) is set, the preferred transmission channel (G) and the reflective surface channel (
Figure pat00100
) at the current set value, the optimal phase variable (
Figure pat00101
) is calculated and updated according to Equation 13 (S50).

위상 변수(

Figure pat00102
)가 업데이트되면, 송신 채널(G)과 업데이트된 위상 변수(
Figure pat00103
)가 고정된 상태에서 최적의 반사표면 채널(
Figure pat00104
)을 수학식 14에 따라 계산하여 업데이트한다(S60).phase variable (
Figure pat00102
) is updated, the transmit channel (G) and the updated phase variable (
Figure pat00103
) is fixed and the optimal reflective surface channel (
Figure pat00104
) is calculated and updated according to Equation 14 (S60).

이후, 업데이트된 위상 변수(

Figure pat00105
)와 반사표면 채널(
Figure pat00106
)을 고정된 상태에서 최적의 송신 채널(G)을 연속 컨벡스 근사(Successive Convex Approximation) 기법을 적용하여 수학식 15 및 16에 따라 획득하여 업데이트한다(S70).After that, the updated phase variable (
Figure pat00105
) and the reflective surface channel (
Figure pat00106
) is obtained and updated according to Equations 15 and 16 by applying a Successive Convex Approximation technique to an optimal transmission channel G in a fixed state (S70).

추정되어야 하는 3개의 추정 변수가 모두 업데이트되면, 이전 업데이트 추정 변수들과의 누적된 변화를 분석하여, 추정 변수들이 수렴되는지 여부를 판판별한다(S80).When all three estimated variables to be estimated are updated, it is determined whether the estimated variables converge or not by analyzing the accumulated changes with the previously updated estimated variables ( S80 ).

만일 수렴되는 것으로 판별되면, 추정 변수들에 대한 추가적인 업데이트를 중지하고 최종 업데이트된 추정 변수들 중 송신 채널(G)과 반사표면 채널(

Figure pat00107
)을 추정된 채널 정보로 획득하고, 무선 전력 전송 시스템의 송신기(100)와 컨트롤러(200)가 획득된 추정 채널 정보에 따라 빔포밍 벡터(w(t))와 위상 변화 행렬(Θ(t))을 생성하여 무선 전력 신호를 송신 및 반사하도록 한다.If it is determined that convergence is determined, additional update of the estimated variables is stopped, and among the last updated estimated variables, the transmission channel (G) and the reflective surface channel (
Figure pat00107
) as the estimated channel information, and the transmitter 100 and the controller 200 of the wireless power transmission system obtain the beamforming vector w(t) and the phase change matrix Θ(t) according to the obtained estimated channel information. ) to transmit and reflect wireless power signals.

그러나 수렴되지 않는 것으로 판별되면, 현재 업데이트된 추정 변수들을 기반으로 다시 추정 변수들을 계산하여 업데이트한다(S50 ~ S70).However, if it is determined that there is no convergence, the estimated variables are calculated and updated again based on the currently updated estimated variables (S50 to S70).

본 발명에 따른 방법은 컴퓨터에서 실행시키기 위한 매체에 저장된 컴퓨터 프로그램으로 구현될 수 있다. 여기서 컴퓨터 판독가능 매체는 컴퓨터에 의해 액세스 될 수 있는 임의의 가용 매체일 수 있고, 또한 컴퓨터 저장 매체를 모두 포함할 수 있다. 컴퓨터 저장 매체는 컴퓨터 판독가능 명령어, 데이터 구조, 프로그램 모듈 또는 기타 데이터와 같은 정보의 저장을 위한 임의의 방법 또는 기술로 구현된 휘발성 및 비휘발성, 분리형 및 비분리형 매체를 모두 포함하며, ROM(판독 전용 메모리), RAM(랜덤 액세스 메모리), CD(컴팩트 디스크)-ROM, DVD(디지털 비디오 디스크)-ROM, 자기 테이프, 플로피 디스크, 광데이터 저장장치 등을 포함할 수 있다.The method according to the present invention may be implemented as a computer program stored in a medium for execution by a computer. Here, the computer-readable medium may be any available medium that can be accessed by a computer, and may include all computer storage media. Computer storage media includes both volatile and nonvolatile, removable and non-removable media implemented in any method or technology for storage of information such as computer readable instructions, data structures, program modules or other data, and read dedicated memory), RAM (Random Access Memory), CD (Compact Disk)-ROM, DVD (Digital Video Disk)-ROM, magnetic tape, floppy disk, optical data storage, and the like.

본 발명은 도면에 도시된 실시예를 참고로 설명되었으나 이는 예시적인 것에 불과하며, 본 기술 분야의 통상의 지식을 가진 자라면 이로부터 다양한 변형 및 균등한 타 실시예가 가능하다는 점을 이해할 것이다.Although the present invention has been described with reference to the embodiment shown in the drawings, which is only exemplary, those skilled in the art will understand that various modifications and equivalent other embodiments are possible therefrom.

따라서, 본 발명의 진정한 기술적 보호 범위는 첨부된 청구범위의 기술적 사상에 의해 정해져야 할 것이다.Accordingly, the true technical protection scope of the present invention should be defined by the technical spirit of the appended claims.

100: 송신기 200: 컨트롤러
300: 지능형 반사 표면(IRS) 401 ~ 40K: 단말
110: 측정 신호 생성부 120: 수신 세기 획득부
130: 채널 추정부 131: 위상 변수 업데이트부
132: 반사표면 채널 업데이트부 133: 송신 채널 업데이트부
134: 수렴 판별부100: transmitter 200: controller
300: Intelligent Reflective Surface (IRS) 401 to 40K: Terminal
110: measurement signal generating unit 120: receiving strength obtaining unit
130: channel estimator 131: phase variable update unit
132: reflective surface channel update unit 133: transmit channel update unit
134: convergence discrimination unit

Claims (20)

무선 전력 전송 시스템의 송신기가 측정 신호를 전송하기 위해 생성한 측정 빔포밍 벡터와 상기 송신기에서 전송된 신호를 반사하여 다수의 단말로 전달하는 지능형 반사표면(이하 IRS)이 측정 신호를 반사하면서 조절하는 위상을 지정하기 위해 생성된 측정 위상 변화 행렬 및 상기 측정 신호를 수신한 다수의 단말 각각에서 피드백된 수신 신호 세기(이하 RSS)를 인가받고, 최대 우도 추정 기법에 따라 추정되어야 하는 상기 다수의 단말 각각에 수신된 측정 신호의 위상을 나타내는 위상 변수와 상기 IRS와 상기 다수의 단말 사이의 채널인 반사표면 채널 및 상기 송신기와 상기 IRS 사이의 채널인 송신 채널 중 반사표면 채널과 송신 채널을 고정하고, 고정된 반사표면 채널과 송신 채널과 인가된 측정 빔포밍 벡터와 측정 위상 변화 행렬 및 RSS를 기반으로 위상 변수를 계산하여 업데이트하는 위상 변수 업데이트부;
상기 송신 채널과 업데이트된 위상 변수를 기반으로 반사표면 채널을 계산하여 업데이트하는 반사표면 채널 업데이트부;
업데이트된 위상 변수와 반사표면 채널을 기반으로 연속 컨벡스 근사 기법에 따라 송신 채널을 획득하여 업데이트하는 송신 채널 업데이트부; 및
업데이트된 위상 변수와 반사표면 채널 및 송신 채널의 수렴 여부를 판별하여, 수렴되면 최종 업데이트된 반사표면 채널 및 송신 채널을 추정된 채널 정보로 획득하는 수렴 판별부를 포함하는 채널 추정 장치.The measurement beamforming vector generated by the transmitter of the wireless power transmission system to transmit the measurement signal and the intelligent reflective surface (hereinafter IRS) that reflects the signal transmitted from the transmitter and transmits it to a plurality of terminals reflects and adjusts the measurement signal The received signal strength (hereinafter, RSS) fed back from each of the plurality of terminals receiving the measurement phase change matrix generated to designate the phase and the measurement signal is applied, and each of the plurality of terminals to be estimated according to the maximum likelihood estimation technique Fixing a reflective surface channel and a transmission channel among a phase variable indicating the phase of the measured signal received in the reflective surface channel, which is a channel between the IRS and the plurality of terminals, and a transmission channel, which is a channel between the transmitter and the IRS, and fixed a phase variable update unit for calculating and updating a phase variable based on the applied reflective surface channel and transmission channel, the applied measurement beamforming vector, the measurement phase change matrix, and the RSS;
a reflective surface channel update unit for calculating and updating a reflective surface channel based on the transmission channel and the updated phase variable;
a transmission channel updater for acquiring and updating a transmission channel according to a continuous convex approximation technique based on the updated phase variable and the reflective surface channel; and
A channel estimator comprising a convergence determining unit for determining whether the updated phase variable, the reflective surface channel, and the transmission channel converge, and obtaining the finally updated reflective surface channel and the transmission channel as estimated channel information when convergence. 제1항에 있어서, 상기 위상 변수 업데이트부는
측정 빔포밍 벡터와 측정 위상 변화 행렬 및 RSS를 인가받고, 반사표면 채널과 송신 채널과 인가된 측정 빔포밍 벡터와 측정 위상 변화 행렬 및 RSS를 기반으로 수학식
Figure pat00108

(여기서 hk 는 이전 기지정된 L개의 타임 슬롯 구간 동안 IRS와 다수의 단말 중 제k 단말 사이의 제k 반사표면 채널이고, φ k는 이전 L개의 타임 슬롯 구간 동안 K개의 단말에서 피드백된 RSS에 위상(
Figure pat00109
) 정보를 추가하기 위해 적용되는 위상 변수이며, M G 는 현지 타임 슬롯(T)로부터 이전 L개의 타임 슬롯 구간 동안 송신기로부터 IRS가 반사할 때까지의 측정 빔포밍 벡터(w)와 송신 채널(G) 및 측정 위상 변화 행렬(Θ)에 의한 송신 반사 채널로서 수학식
Figure pat00110
으로 계산된다.)
에 따라 업데이트되는 위상 변수를 계산하는 채널 추정 장치.The method of claim 1, wherein the phase variable update unit
The measurement beamforming vector, the measurement phase change matrix, and the RSS are applied, and the reflective surface channel, the transmission channel, and the applied measurement beamforming vector, the measurement phase change matrix, and the RSS are based on an equation
Figure pat00108

(here, h k is the k-th reflective surface channel between the IRS and the k-th terminal among the plurality of terminals during the previously predetermined L time slot period, and φ k is the RSS fed back from the K terminals during the previous L time slot period. Phase(
Figure pat00109
) is a phase variable applied to add information, and M G is the measurement beamforming vector (w) and the transmission channel (G) from the local time slot (T) until the IRS reflects from the transmitter during the previous L time slot sections ) and the measured phase shift matrix ( Θ ) as the transmit reflection channel by the equation
Figure pat00110
is calculated.)
A channel estimator that calculates a phase variable that is updated according to . 제2항에 있어서, 상기 반사표면 채널 업데이트부는
송신 채널과 업데이트된 위상 변수를 기반으로 수학식
Figure pat00111

(여기서 βk는 이전 L개의 타임 슬롯 구간 동안 제k 단말의 속도에 따라 결정되는 채널 상관 계수이고, r k 는 이전 L개의 타임 슬롯 구간 동안 제k 단말로부터 피드백된 RSS이며, ㅀ는 하다마드 곱셈 기호이다. 그리고
Figure pat00112
는 이전 L개의 타임 슬롯 구간 동안 송신기와 제k 단말 사이의 가우시안 노이즈 벡터에 대한 공분산 행렬(R k)을 근사하여 획득되는 근사 공분산 행렬을 나타낸다.)
에 따라 업데이트되는 반사표면 채널을 계산하는 채널 추정 장치.The method of claim 2, wherein the reflective surface channel update unit
Equation based on transmit channel and updated phase variable
Figure pat00111

(where β k is a channel correlation coefficient determined according to the speed of the k-th terminal during the previous L time slot sections, r k is the RSS fed back from the k-th terminal during the previous L time slot sections, and ㅀ is the Hadamard multiplication It is a sign, and
Figure pat00112
denotes an approximate covariance matrix obtained by approximating a covariance matrix ( R k ) for a Gaussian noise vector between a transmitter and a k-th terminal during the previous L time slot periods.)
A channel estimator that computes a reflective surface channel that is updated according to 제3항에 있어서, 상기 반사표면 채널 업데이트부는
상기 근사 공분산 행렬을
수학식
Figure pat00113

(여기서, p와 q는 각각 공분산 행렬(R k)의 행과 열을 나타낸다.)
에 따라 계산하는 가우시안 노이즈 벡터에 대한 공분산 행렬(R k)을
수학식
Figure pat00114

에 따라 근사하여 획득하는 채널 추정 장치.The method of claim 3, wherein the reflective surface channel update unit
The approximate covariance matrix
formula
Figure pat00113

(Where p and q represent the row and column of the covariance matrix ( R k ), respectively.)
The covariance matrix ( R k ) for the Gaussian noise vector computed according to
formula
Figure pat00114

A channel estimation device obtained by approximation according to 제4항에 있어서, 상기 송신 채널 업데이트부는
업데이트된 위상 변수와 업데이트된 반사표면 채널을 기반으로 수학식
Figure pat00115

(여기서
Figure pat00116
는 업데이트되기 이전 송신 채널이고,
Figure pat00117
는 역추적 라인 탐색법으로 획득 스텝 사이즈이며,
Figure pat00118
는 송신 채널(G)에 따른 반사표면 채널(h)의 기울기를 나타낸다.)
에 따라 업데이트되는 송신 채널을 계산하는 채널 추정 장치.5. The method of claim 4, wherein the transmission channel update unit
Equation based on the updated phase variable and the updated reflective surface channel
Figure pat00115

(here
Figure pat00116
is the transmission channel before being updated,
Figure pat00117
is the step size acquired by the trace line search method,
Figure pat00118
denotes the slope of the reflection surface channel (h) according to the transmission channel (G).)
A channel estimation device for calculating a transmission channel updated according to . 제5항에 있어서, 상기 송신 채널 업데이트부는
반사표면 채널(h)의 기울기(
Figure pat00119
)를 수학식
Figure pat00120

(여기서 송신기와 제k 단말 사이의 가우시안 노이즈의 표준 편차를 나타낸다.)
에 따라 계산하는 채널 추정 장치.The method of claim 5, wherein the transmission channel update unit
The slope of the reflective surface channel (h) (
Figure pat00119
) to the formula
Figure pat00120

(Here, it represents the standard deviation of Gaussian noise between the transmitter and the kth terminal.)
A channel estimation device that calculates according to 제6항에 있어서, 상기 수렴 판별부는
업데이트된 위상 변수와 반사표면 채널 및 송신 채널의 누적 변화를 분석하여 위상 변수와 반사표면 채널 및 송신 채널과 이전 업데이트된 위상 변수와 반사표면 채널 및 송신 채널 사이의 편차가 기지정된 범위 이내로 수렴되는 것으로 판별되면, 최종 업데이트된 반사표면 채널 및 송신 채널을 추정된 채널 정보로 획득하고, 수렴되지 않는 것으로 판별되면, 위상 변수와 반사표면 채널 및 송신 채널이 반복 업데이트되도록 업데이트된 위상 변수와 반사표면 채널 및 송신 채널을 상기 위상 변수 업데이트부로 전달하는 채널 추정 장치.The method of claim 6, wherein the convergence determining unit
By analyzing the cumulative change of the updated phase variable and the reflective surface channel and the transmit channel, it is determined that the deviation between the phase variable, the reflective surface channel and the transmit channel and the previously updated phase variable and the reflective surface channel and the transmit channel converge within a predetermined range. If it is determined, the last updated reflective surface channel and transmission channel are acquired as estimated channel information, and if it is determined that there is no convergence, the updated phase variable and the reflective surface channel and the transmitted channel are repeatedly updated. A channel estimator for transferring a transmission channel to the phase variable updater. 제7항에 있어서, 상기 위상 변수 업데이트부는
현재 타임 슬롯에서 가장 최근 업데이트된 위상 변수와 반사표면 채널 및 송신 채널을 이용하여 업데이트되는 위상 변수를 계산하고, 현재 타임 슬롯에서 이전 업데이트된 위상 변수와 반사표면 채널 및 송신 채널이 존재하지 않으면, 이전 타임 슬롯에서 최종 업데이트된 위상 변수와 반사표면 채널 및 송신 채널을 이용하여 업데이트되는 위상 변수를 계산하는 채널 추정 장치.The method of claim 7, wherein the phase variable update unit
Calculates the updated phase variable using the most recently updated phase variable, reflective surface channel and transmit channel in the current time slot, and if there is no previous updated phase variable, reflective surface channel and transmit channel in the current time slot, A channel estimator that calculates a phase variable that is updated using the last updated phase variable in a time slot, a reflective surface channel, and a transmission channel. 제1항에 있어서, 상기 측정 빔포밍 벡터는
상기 송신기에 의해 평균이 0이고, 평균 측정 전력과 송신기 안테나 개수의 비에 따른 분산을 갖는 복소 가우시안 분포를 따르도록 임의로 생성되는 채널 추정 장치.The method of claim 1, wherein the measurement beamforming vector is
A channel estimation apparatus arbitrarily generated by the transmitter to follow a complex Gaussian distribution with an average of 0 and a variance according to a ratio of the average measured power and the number of antennas of the transmitter. 제1항에 있어서, 상기 측정 위상 변화 행렬은
상기 IRS를 제어하는 컨트롤러에 의해 IRS의 N개의 반사 소자에 따른 측정 신호의 위상 변화가 시간에 독립적으로 0 ~ 2π의 범위에서 균일하게 분포하도록 임의로 생성되는 채널 추정 장치.The method of claim 1, wherein the measured phase change matrix is
A channel estimation device arbitrarily generated by a controller for controlling the IRS so that the phase change of the measurement signal according to the N reflective elements of the IRS is uniformly distributed in the range of 0 to 2π independently of time. 컨트롤러의 제어에 따라 송신기에서 전송된 무선 전력 신호를 반사하여 다수의 단말로 전달하는 지능형 반사표면(이하 IRS)을 포함하는 무선 전력 전송 시스템의 채널 추정 방법에 있어서,
상기 송신기가 측정 신호를 전송하기 위해 생성한 측정 빔포밍 벡터와 상기 IRS가 측정 신호를 반사하면서 조절하는 위상을 지정하기 위해 생성된 측정 위상 변화 행렬을 획득하는 단계;
상기 측정 신호를 수신한 다수의 단말 각각에서 피드백된 수신 신호 세기(이하 RSS)를 인가받는 단계;
최대 우도 추정 기법에 따라 추정되어야 하는 상기 다수의 단말 각각에 수신된 측정 신호의 위상을 나타내는 위상 변수와 상기 IRS와 상기 다수의 단말 사이의 채널인 반사표면 채널 및 상기 송신기와 상기 IRS 사이의 채널인 송신 채널 중 반사표면 채널과 송신 채널을 고정하고, 고정된 반사표면 채널과 송신 채널과 인가된 측정 빔포밍 벡터와 측정 위상 변화 행렬 및 RSS를 기반으로 위상 변수를 계산하여 업데이트하는 단계;
상기 송신 채널과 업데이트된 위상 변수를 기반으로 반사표면 채널을 계산하여 업데이트하는 단계;
업데이트된 위상 변수와 반사표면 채널을 기반으로 연속 컨벡스 근사 기법에 따라 송신 채널을 획득하여 업데이트하는 단계; 및
업데이트된 위상 변수와 반사표면 채널 및 송신 채널의 수렴 여부를 판별하여, 수렴되면 최종 업데이트된 반사표면 채널 및 송신 채널을 추정된 채널 정보로 획득하는 단계를 포함하는 채널 추정 방법.In the channel estimation method of a wireless power transmission system comprising an intelligent reflective surface (hereinafter IRS) that reflects the wireless power signal transmitted from the transmitter under the control of a controller and transmits it to a plurality of terminals,
obtaining a measurement beamforming vector generated by the transmitter to transmit a measurement signal and a measurement phase change matrix generated to designate a phase that the IRS adjusts while reflecting the measurement signal;
receiving received signal strength (hereinafter, RSS) fed back from each of a plurality of terminals receiving the measurement signal;
A phase variable indicating the phase of the measurement signal received in each of the plurality of terminals to be estimated according to the maximum likelihood estimation technique, a reflection surface channel that is a channel between the IRS and the plurality of terminals, and a channel between the transmitter and the IRS fixing a reflective surface channel and a transmission channel among transmission channels, calculating and updating a phase variable based on the fixed reflective surface channel and transmission channel, an applied measurement beamforming vector, a measurement phase change matrix, and RSS;
calculating and updating a reflective surface channel based on the transmission channel and the updated phase variable;
acquiring and updating a transmission channel according to a continuous convex approximation technique based on the updated phase variable and the reflective surface channel; and
A method for estimating a channel, comprising: determining whether an updated phase variable, a reflective surface channel, and a transmission channel converge; and acquiring a final updated reflective surface channel and a transmission channel as estimated channel information when convergence. 제11항에 있어서, 상기 위상 변수를 계산하여 업데이트하는 단계는
측정 빔포밍 벡터와 측정 위상 변화 행렬 및 RSS를 인가받는 단계; 및
반사표면 채널과 송신 채널과 인가된 측정 빔포밍 벡터와 측정 위상 변화 행렬 및 RSS를 기반으로 수학식
Figure pat00121

(여기서 hk 는 이전 기지정된 L개의 타임 슬롯 구간 동안 IRS와 다수의 단말 중 제k 단말 사이의 제k 반사표면 채널이고, φ k는 이전 L개의 타임 슬롯 구간 동안 K개의 단말에서 피드백된 RSS에 위상(
Figure pat00122
) 정보를 추가하기 위해 적용되는 위상 변수이며, M G 는 현지 타임 슬롯(T)로부터 이전 L개의 타임 슬롯 구간 동안 송신기로부터 IRS가 반사할 때까지의 측정 빔포밍 벡터(w)와 송신 채널(G) 및 측정 위상 변화 행렬(Θ)에 의한 송신 반사 채널로서 수학식
Figure pat00123
으로 계산된다.)
에 따라 업데이트되는 위상 변수를 계산하는 단계를 포함하는 채널 추정 방법.12. The method of claim 11, wherein the calculating and updating the phase variable comprises:
receiving a measurement beamforming vector, a measurement phase change matrix, and RSS; and
Equation based on reflective surface channel and transmission channel, applied measurement beamforming vector, measurement phase change matrix, and RSS
Figure pat00121

(here, h k is the k-th reflective surface channel between the IRS and the k-th terminal among the plurality of terminals during the previously predetermined L time slot period, and φ k is the RSS fed back from the K terminals during the previous L time slot period. Phase(
Figure pat00122
) is a phase variable applied to add information, and M G is the measurement beamforming vector (w) and the transmission channel (G) from the local time slot (T) until the IRS reflects from the transmitter during the previous L time slot sections ) and the measured phase shift matrix ( Θ ) as the transmit reflection channel by the equation
Figure pat00123
is calculated.)
A method of estimating a channel comprising calculating a phase variable that is updated according to 제12항에 있어서, 상기 반사표면 채널을 계산하여 업데이트하는 단계는
송신 채널과 업데이트된 위상 변수를 기반으로 수학식
Figure pat00124

(여기서 βk는 이전 L개의 타임 슬롯 구간 동안 제k 단말의 속도에 따라 결정되는 채널 상관 계수이고, r k 는 이전 L개의 타임 슬롯 구간 동안 제k 단말로부터 피드백된 RSS이며, ㅀ는 하다마드 곱셈 기호이다. 그리고
Figure pat00125
는 이전 L개의 타임 슬롯 구간 동안 송신기와 제k 단말 사이의 가우시안 노이즈 벡터에 대한 공분산 행렬(R k)을 근사하여 획득되는 근사 공분산 행렬을 나타낸다.)
에 따라 업데이트되는 반사표면 채널을 계산하는 채널 추정 방법.13. The method of claim 12, wherein calculating and updating the reflective surface channel comprises:
Equation based on transmit channel and updated phase variable
Figure pat00124

(where β k is a channel correlation coefficient determined according to the speed of the k-th terminal during the previous L time slot sections, r k is the RSS fed back from the k-th terminal during the previous L time slot sections, and ㅀ is the Hadamard multiplication It is a sign, and
Figure pat00125
denotes an approximate covariance matrix obtained by approximating a covariance matrix ( R k ) for a Gaussian noise vector between a transmitter and a k-th terminal during the previous L time slot periods.)
A channel estimation method that computes a reflective surface channel that is updated according to 제13항에 있어서, 상기 근사 공분산 행렬은
수학식
Figure pat00126

(여기서, p와 q는 각각 공분산 행렬(R k)의 행과 열을 나타낸다.)
에 따라 계산하되 가우시안 노이즈 벡터에 대한 공분산 행렬(R k)을
수학식
Figure pat00127

에 따라 근사하여 획득되는 채널 추정 방법.14. The method of claim 13, wherein the approximate covariance matrix is
formula
Figure pat00126

(Where p and q represent the row and column of the covariance matrix ( R k ), respectively.)
Calculate the covariance matrix ( R k ) for the Gaussian noise vector according to
formula
Figure pat00127

A channel estimation method obtained by approximation according to 제14항에 있어서, 상기 송신 채널을 획득하여 업데이트하는 단계는
업데이트된 위상 변수와 업데이트된 반사표면 채널을 기반으로 수학식
Figure pat00128

(여기서
Figure pat00129
는 업데이트되기 이전 송신 채널이고,
Figure pat00130
는 역추적 라인 탐색법으로 획득 스텝 사이즈이며,
Figure pat00131
는 송신 채널(G)에 따른 반사표면 채널(h)의 기울기를 나타낸다.)
에 따라 업데이트되는 송신 채널을 계산하는 채널 추정 방법.15. The method of claim 14, wherein the step of obtaining and updating the transmission channel comprises:
Equation based on the updated phase variable and the updated reflective surface channel
Figure pat00128

(here
Figure pat00129
is the transmission channel before being updated,
Figure pat00130
is the step size acquired by the trace line search method,
Figure pat00131
denotes the slope of the reflection surface channel (h) according to the transmission channel (G).)
A channel estimation method for calculating a transmission channel that is updated according to . 제15항에 있어서, 상기 반사표면 채널(h)의 기울기(
Figure pat00132
)는
수학식
Figure pat00133

(여기서 송신기와 제k 단말 사이의 가우시안 노이즈의 표준 편차를 나타낸다.)
에 따라 계산되는 채널 추정 방법.16. The method of claim 15, wherein the slope of the reflective surface channel (h) (h)
Figure pat00132
)Is
formula
Figure pat00133

(Here, it represents the standard deviation of Gaussian noise between the transmitter and the kth terminal.)
Channel estimation method calculated according to . 제16항에 있어서, 상기 추정된 채널 정보로 획득하는 단계는
업데이트된 위상 변수와 반사표면 채널 및 송신 채널의 누적 변화를 분석하여 위상 변수와 반사표면 채널 및 송신 채널과 이전 업데이트된 위상 변수와 반사표면 채널 및 송신 채널 사이의 편차가 기지정된 범위 이내로 수렴되는지 판별하는 단계;
수렴되는 것으로 판별되면, 최종 업데이트된 반사표면 채널 및 송신 채널을 추정된 채널 정보로 획득하는 단계; 및
수렴되지 않는 것으로 판별되면, 위상 변수와 반사표면 채널 및 송신 채널을 반복 업데이트하는 단계를 포함하는 채널 추정 방법.The method of claim 16, wherein the acquiring with the estimated channel information comprises:
Analyzes the cumulative change of the updated phase variable and the reflective surface channel and the transmit channel to determine whether the deviation between the phase variable and the reflective surface channel and transmit channel and the previously updated phase variable and the reflective surface channel and transmit channel converges within a predetermined range to do;
if it is determined that the convergence is determined, acquiring the last updated reflective surface channel and transmission channel as estimated channel information; and
and if it is determined that there is no convergence, iteratively updating a phase variable, a reflective surface channel, and a transmission channel. 제17항에 있어서, 상기 위상 변수를 계산하는 단계는
현재 타임 슬롯에서 가장 최근 업데이트된 위상 변수와 반사표면 채널 및 송신 채널을 이용하여 업데이트되는 위상 변수를 계산하고, 현재 타임 슬롯에서 이전 업데이트된 위상 변수와 반사표면 채널 및 송신 채널이 존재하지 않으면, 이전 타임 슬롯에서 최종 업데이트된 위상 변수와 반사표면 채널 및 송신 채널을 이용하여 업데이트되는 위상 변수를 계산하는 채널 추정 방법.18. The method of claim 17, wherein calculating the phase variable comprises:
Calculates the updated phase variable using the most recently updated phase variable, reflective surface channel and transmit channel in the current time slot, and if there is no previous updated phase variable, reflective surface channel and transmit channel in the current time slot, A channel estimation method for calculating a phase variable that is updated using the last updated phase variable in a time slot and a reflective surface channel and a transmit channel. 제11항에 있어서, 상기 측정 빔포밍 벡터는
상기 송신기에 의해 평균이 0이고, 평균 측정 전력과 송신기 안테나 개수의 비에 따른 분산을 갖는 복소 가우시안 분포를 따르도록 임의로 생성되는 채널 추정 방법.The method of claim 11, wherein the measurement beamforming vector is
A channel estimation method arbitrarily generated by the transmitter to follow a complex Gaussian distribution with an average of 0 and a variance according to the ratio of the average measured power and the number of antennas of the transmitter. 제11항에 있어서, 상기 측정 위상 변화 행렬은
상기 컨트롤러에 의해 상기 IRS의 N개의 반사 소자에 따른 측정 신호의 위상 변화가 시간에 독립적으로 0 ~ 2π의 범위에서 균일하게 분포하도록 임의로 생성되는 채널 추정 방법.12. The method of claim 11, wherein the measured phase change matrix is
A channel estimation method arbitrarily generated by the controller so that the phase change of the measurement signal according to the N reflective elements of the IRS is uniformly distributed in the range of 0 to 2π independently of time.
KR1020200150352A 2020-11-11 2020-11-11 Apparatus and method for estimating channel of wireless power transmission system using an intelligent reflective surface KR102457507B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020200150352A KR102457507B1 (en) 2020-11-11 2020-11-11 Apparatus and method for estimating channel of wireless power transmission system using an intelligent reflective surface

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020200150352A KR102457507B1 (en) 2020-11-11 2020-11-11 Apparatus and method for estimating channel of wireless power transmission system using an intelligent reflective surface

Publications (3)

Publication Number Publication Date
KR20220064154A true KR20220064154A (en) 2022-05-18
KR102457507B1 KR102457507B1 (en) 2022-10-20
KR102457507B9 KR102457507B9 (en) 2023-03-13

Family

ID=81800812

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020200150352A KR102457507B1 (en) 2020-11-11 2020-11-11 Apparatus and method for estimating channel of wireless power transmission system using an intelligent reflective surface

Country Status (1)

Country Link
KR (1) KR102457507B1 (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115037406A (en) * 2022-05-27 2022-09-09 南京理工大学 Interference covariance matrix estimation method based on gradient descent method
CN115483957A (en) * 2022-09-20 2022-12-16 华工未来科技(江苏)有限公司 Phase regulation and control method and device of RIS auxiliary communication system, electronic equipment and medium
CN115694578A (en) * 2022-09-19 2023-02-03 华工未来科技(江苏)有限公司 Multi-intelligent reflecting surface synchronous control method, system, device and storage medium
KR102546184B1 (en) * 2022-09-19 2023-06-22 전남대학교산학협력단 Multi-task learning-based secure transmission system and method
CN117221978A (en) * 2023-09-11 2023-12-12 国网江苏省电力有限公司电力科学研究院 Passive array-assisted data transmission method and device between nodes of electric power Internet of things
KR20240002313A (en) 2022-06-29 2024-01-05 경희대학교 산학협력단 Intelligent reflecting surface assisted uplink nonorthogonal multiple access system, and optimization method to improve minimum rate using the same
WO2024010112A1 (en) * 2022-07-06 2024-01-11 엘지전자 주식회사 Device and method for estimating channel in wireless communication system
GB2624191A (en) * 2022-11-09 2024-05-15 Nokia Technologies Oy RIS assisted energy steering for energy harvesting devices
KR20240150923A (en) 2023-04-10 2024-10-17 국립금오공과대학교 산학협력단 RSMA-based IRS-supported OAM wireless communication system for non-line-of-sight users
KR20240150926A (en) 2023-04-10 2024-10-17 국립금오공과대학교 산학협력단 RIS-based OAM wireless communication system using unmanned aerial vehicles
KR20240150924A (en) 2023-04-10 2024-10-17 국립금오공과대학교 산학협력단 Split RIS-based OAM wireless communication system for multi-user

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101678640B1 (en) 2014-10-01 2016-11-23 한국과학기술원 Wireless Power Transfer Apparatus Using Reflector and Separator
KR20170100649A (en) * 2014-12-29 2017-09-04 에너저스 코포레이션 System and method for wireless power transmission
KR20180074425A (en) * 2016-12-23 2018-07-03 삼성전자주식회사 Wireless power transmitter, electronic device and method for controlling thereof
KR20180090893A (en) * 2015-12-24 2018-08-13 에너저스 코포레이션 Object detection system and method in wireless power charging system
KR20200076263A (en) * 2018-12-19 2020-06-29 한국과학기술원 Method and method for multi-antenna transmission to minimize charging time of users in wireless powered communications networks

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101678640B1 (en) 2014-10-01 2016-11-23 한국과학기술원 Wireless Power Transfer Apparatus Using Reflector and Separator
KR20170100649A (en) * 2014-12-29 2017-09-04 에너저스 코포레이션 System and method for wireless power transmission
KR20180090893A (en) * 2015-12-24 2018-08-13 에너저스 코포레이션 Object detection system and method in wireless power charging system
KR20180074425A (en) * 2016-12-23 2018-07-03 삼성전자주식회사 Wireless power transmitter, electronic device and method for controlling thereof
KR20200076263A (en) * 2018-12-19 2020-06-29 한국과학기술원 Method and method for multi-antenna transmission to minimize charging time of users in wireless powered communications networks

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115037406A (en) * 2022-05-27 2022-09-09 南京理工大学 Interference covariance matrix estimation method based on gradient descent method
KR20240002313A (en) 2022-06-29 2024-01-05 경희대학교 산학협력단 Intelligent reflecting surface assisted uplink nonorthogonal multiple access system, and optimization method to improve minimum rate using the same
WO2024010112A1 (en) * 2022-07-06 2024-01-11 엘지전자 주식회사 Device and method for estimating channel in wireless communication system
KR102546184B1 (en) * 2022-09-19 2023-06-22 전남대학교산학협력단 Multi-task learning-based secure transmission system and method
CN115694578B (en) * 2022-09-19 2023-11-24 华工未来科技(江苏)有限公司 Multi-intelligent reflecting surface synchronous control method, system, device and storage medium
CN115694578A (en) * 2022-09-19 2023-02-03 华工未来科技(江苏)有限公司 Multi-intelligent reflecting surface synchronous control method, system, device and storage medium
CN115483957B (en) * 2022-09-20 2024-01-09 华工未来科技(江苏)有限公司 Phase regulation method and device for RIS auxiliary communication system, electronic equipment and medium
CN115483957A (en) * 2022-09-20 2022-12-16 华工未来科技(江苏)有限公司 Phase regulation and control method and device of RIS auxiliary communication system, electronic equipment and medium
GB2624191A (en) * 2022-11-09 2024-05-15 Nokia Technologies Oy RIS assisted energy steering for energy harvesting devices
KR20240150923A (en) 2023-04-10 2024-10-17 국립금오공과대학교 산학협력단 RSMA-based IRS-supported OAM wireless communication system for non-line-of-sight users
KR20240150926A (en) 2023-04-10 2024-10-17 국립금오공과대학교 산학협력단 RIS-based OAM wireless communication system using unmanned aerial vehicles
KR20240150924A (en) 2023-04-10 2024-10-17 국립금오공과대학교 산학협력단 Split RIS-based OAM wireless communication system for multi-user
CN117221978A (en) * 2023-09-11 2023-12-12 国网江苏省电力有限公司电力科学研究院 Passive array-assisted data transmission method and device between nodes of electric power Internet of things
CN117221978B (en) * 2023-09-11 2024-04-30 国网江苏省电力有限公司电力科学研究院 Passive array-assisted data transmission method and device between nodes of electric power Internet of things

Also Published As

Publication number Publication date
KR102457507B1 (en) 2022-10-20
KR102457507B9 (en) 2023-03-13

Similar Documents

Publication Publication Date Title
KR102457507B1 (en) Apparatus and method for estimating channel of wireless power transmission system using an intelligent reflective surface
CN110365375B (en) Beam alignment and tracking method in millimeter wave communication system and computer equipment
US11967829B2 (en) Apparatus and method for estimating channel of wireless power transmission system using intelligent reflecting surface
Chiriac et al. Ziv–Zakai bound for joint parameter estimation in MIMO radar systems
WO2015192696A1 (en) Gossip algorithm based multiple target doa estimating system and estimating method in distributed network
KR102709853B1 (en) Method for determinig optimal beam and an electronic device thereof
US12019181B2 (en) Iterative focused millimeter wave integrated communication and sensing method
CN114024586B (en) Intelligent beam prediction method, device, equipment and medium for nonlinear track
KR20190043080A (en) Estimation method and compensation method for rf chain imbalance in uca oam radio system
US20240256738A1 (en) Environmental sensing method based on model evolution
Feriani et al. On the robustness of deep reinforcement learning in IRS-aided wireless communications systems
KR20220144679A (en) Method and Apparatus for Simultaneous Optimization of Transmitted Rate and The Harvested Energy in Intelligent Reflecting Surface-aided MIMO System
KR102027674B1 (en) Method and system for estimating i/q imbalance parameter of transceiver based on reinforcement learning
KR20170065252A (en) Method and apparatus for increasing channel capacity in line of sight circumstance
Sun et al. Adaptive source localization by a mobile robot using signal power gradient in sensor networks
KR20220017161A (en) Transmitting and receiving apparatus and method of wireless power transmission system to improve wireless power transmission efficiency
Seo et al. Pilot beam sequence design for channel estimation in millimeter-wave MIMO systems: A POMDP framework
CN107426748B (en) Method for estimating performance of multiple sensors in wireless network control system
CN115361258A (en) Large-scale MIMO sparse channel estimation method and related equipment
KR20200029646A (en) an internet of Things system
CN104023396B (en) Single goal DOA estimating systems based on gossip algorithms and method of estimation in distributed network
CN116233732A (en) Co-location method, apparatus, computer device and computer readable storage medium
WO2021199335A1 (en) Information processing device, information processing method, and communication device
KR101606159B1 (en) Pilot signal design method for massive antenna systems
Jiang et al. Cancellation strategy in dynamic framed slotted ALOHA for RFID system

Legal Events

Date Code Title Description
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
G170 Re-publication after modification of scope of protection [patent]