JP7156666B2 - Parameter search device, parameter search method, and program - Google Patents

Parameter search device, parameter search method, and program Download PDF

Info

Publication number
JP7156666B2
JP7156666B2 JP2018116286A JP2018116286A JP7156666B2 JP 7156666 B2 JP7156666 B2 JP 7156666B2 JP 2018116286 A JP2018116286 A JP 2018116286A JP 2018116286 A JP2018116286 A JP 2018116286A JP 7156666 B2 JP7156666 B2 JP 7156666B2
Authority
JP
Japan
Prior art keywords
search
square
objective function
squares
value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2018116286A
Other languages
Japanese (ja)
Other versions
JP2019219869A (en
Inventor
雅次 伊藤
安規男 宇田
孝行 天羽
哲也 常友
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NEC Solutions Innovators Ltd
Original Assignee
NEC Solutions Innovators Ltd
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 NEC Solutions Innovators Ltd filed Critical NEC Solutions Innovators Ltd
Priority to JP2018116286A priority Critical patent/JP7156666B2/en
Publication of JP2019219869A publication Critical patent/JP2019219869A/en
Application granted granted Critical
Publication of JP7156666B2 publication Critical patent/JP7156666B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Complex Calculations (AREA)

Description

本開示は、パラメータ探索装置、パラメータ探索方法、及びプログラムに関する。 The present disclosure relates to a parameter search device, parameter search method, and program.

図1を参照して、物流センター等における商品の在庫数の時系列的な変化について説明する。図1は、物流センター等における商品の在庫数の時系列データの例を示す図である。
図1において、横軸は時間、縦軸は商品の在庫数である。
また、図1の縦軸において、最大在庫数は、商品の必要最大限の在庫数である。平均在庫数は、特定期間中の商品の在庫数の平均値である。発注点数は、商品を発注するタイミングの基準となる在庫数であり、商品の在庫数が発注点数まで減少すると、商品を発注する。安全在庫数は、商品の品切れや欠品を防ぐために必要となる在庫数である。
With reference to FIG. 1, chronological changes in the number of products in stock at a distribution center or the like will be described. FIG. 1 is a diagram showing an example of time-series data of the number of products in stock at a distribution center or the like.
In FIG. 1, the horizontal axis is time, and the vertical axis is the number of products in stock.
In addition, on the vertical axis of FIG. 1, the maximum inventory quantity is the maximum required inventory quantity of the product. The average number of inventories is the average value of the number of inventories of products during a specific period. The number of items to be ordered is the number of items in stock that serves as a reference for the timing of ordering products. The number of safety stocks is the number of stocks required to prevent out-of-stock or shortage of products.

また、図1の横軸において、最大在庫日数(=保有日数)は、最大在庫数分の商品が減少して0になるまでの日数である。発注リードタイム日数は、商品の調達に要する日数である。安全在庫日数は、安全在庫数分の商品が減少して0になるまでの日数である。安全在庫日数は、安全在庫数を平均出荷数で除算することで求められる。 In addition, on the horizontal axis of FIG. 1, the maximum inventory days (=holding days) is the number of days until the number of products corresponding to the maximum inventory decreases to zero. The number of days of order lead time is the number of days required to procure the product. The number of safety stock days is the number of days until the number of safety stock items decreases to 0. The safety stock days is obtained by dividing the number of safety stocks by the average number of shipments.

図1の例では、商品の在庫数が発注点数まで減少したタイミングで、商品の発注を行う。発注点数は、最大在庫日数及び安全在庫日数を基に求められる。従って、適切なタイミングで商品の発注を行うには、発注点数を決定するパラメータとしての最大在庫日数及び安全在庫日数を適切に設定する必要がある。 In the example of FIG. 1, the product is ordered at the timing when the inventory quantity of the product has decreased to the order quantity. The number of items to be ordered is calculated based on the maximum inventory days and safety inventory days. Therefore, in order to order products at appropriate timing, it is necessary to appropriately set the maximum number of days in stock and the number of days in safety stock as parameters for determining the number of items to be ordered.

従来、最大在庫日数及び安全在庫日数は、ユーザが過去の経験に基づき手動で設定していた。しかし、最大在庫日数及び安全在庫日数は、平均在庫日数、発注回数といった複数の指標に依存して決定されるものであるため、経験豊富なユーザであっても、最大在庫日数及び安全在庫日数を適切に設定することは困難であった。 Conventionally, the maximum number of days in stock and the number of days in safety stock were manually set by the user based on past experience. However, since the maximum days in stock and the safety days in stock are determined depending on multiple indicators such as the average days in stock and the number of orders, even an experienced user cannot determine the maximum days in stock and the safety days in stock. It was difficult to set properly.

そのため、最近は、AI(Artificial Intelligence)技術を利用して、パラメータの最適解を探索する技術が提案されている。例えば、特許文献1,2には、勾配降下法を利用して、パラメータの最適解を探索する技術が開示されている。 Therefore, recently, a technique of searching for the optimum solution of parameters using AI (Artificial Intelligence) technique has been proposed. For example, Patent Literatures 1 and 2 disclose techniques for searching for optimal solutions for parameters using the gradient descent method.

特許第4349271号公報Japanese Patent No. 4349271 特表2005-527034号公報Japanese Patent Publication No. 2005-527034

しかし、特許文献1,2に開示されるような一般的な勾配降下法は、パラメータの最適解の探索に微分を利用するものであるため、探索に複雑な計算が必要であり、探索に時間が掛かるといった問題がある。 However, the general gradient descent method disclosed in Patent Documents 1 and 2 uses differentiation to search for the optimum solution of parameters, so the search requires complicated calculations, and the search takes a long time. There is a problem that it hangs.

そこで、本開示の目的は、上述した課題を解決し、複雑な計算を行うことなく、パラメータの最適解を探索することができるパラメータ探索装置、パラメータ探索方法、及びプログラムを提供することにある。 Therefore, an object of the present disclosure is to solve the above-described problems and to provide a parameter search device, parameter search method, and program capable of searching for the optimum solution of parameters without performing complicated calculations.

一態様によるパラメータ探索装置は、
N(Nは2以上の自然数)個のパラメータを軸とする所定のN次元空間をマス目状に分割する分割部と、
前記所定のN次元空間上のマスの中から、当該マスに対応するN個のパラメータのそれぞれの値を用いて得られる目的関数の値が最小値となるマスを探索する探索部と、
前記目的関数の値が最小値となるマスに対応するN個のパラメータのそれぞれの値を最適解に決定する決定部と、
を備える。
A parameter search device according to one aspect includes:
a dividing unit that divides a predetermined N-dimensional space with N (N is a natural number equal to or greater than 2) parameters as axes;
a searching unit for searching, from among the squares in the predetermined N-dimensional space, a square having a minimum value of an objective function obtained using the respective values of N parameters corresponding to the square;
a determination unit that determines the values of the N parameters corresponding to the mass at which the value of the objective function is the minimum value as an optimal solution;
Prepare.

一態様によるパラメータ探索方法は、
パラメータ探索装置によるパラメータ探索方法であって、
N(Nは2以上の自然数)個のパラメータを軸とする所定のN次元空間をマス目状に分割するステップと、
前記所定のN次元空間上のマスの中から、当該マスに対応するN個のパラメータのそれぞれの値を用いて得られる目的関数の値が最小値となるマスを探索するステップと、
前記目的関数の値が最小値となるマスに対応するN個のパラメータのそれぞれの値を最適解に決定するステップと、
を含む。
A parameter search method according to one aspect comprises:
A parameter search method by a parameter search device,
a step of dividing a predetermined N-dimensional space with N (N is a natural number equal to or greater than 2) parameters as axes;
a step of searching from among the squares in the predetermined N-dimensional space for a square having the minimum value of an objective function obtained using the respective values of N parameters corresponding to the square;
a step of determining the values of each of the N parameters corresponding to the mass at which the value of the objective function is the minimum value as an optimal solution;
including.

一態様によるプログラムは、
コンピュータを、
N(Nは2以上の自然数)個のパラメータを軸とする所定のN次元空間をマス目状に分割する分割手段と、
前記所定のN次元空間上のマスの中から、当該マスに対応するN個のパラメータのそれぞれの値を用いて得られる目的関数の値が最小値となるマスを探索する探索手段と、
前記目的関数の値が最小値となるマスに対応するN個のパラメータのそれぞれの値を最適解に決定する決定手段と、
として機能させる。
A program according to one aspect comprises:
the computer,
dividing means for dividing a predetermined N-dimensional space with N (N is a natural number equal to or greater than 2) parameters as axes;
search means for searching, from among the squares in the predetermined N-dimensional space, a square having a minimum value of an objective function obtained using respective values of N parameters corresponding to the square;
Determination means for determining the values of each of the N parameters corresponding to the mass at which the value of the objective function is the minimum value as an optimal solution;
function as

上記の態様によれば、複雑な計算を行うことなく、パラメータを探索することができるパラメータ探索装置、パラメータ探索方法、及びプログラムを提供できるという効果が得られる。 According to the above aspect, it is possible to provide a parameter search device, a parameter search method, and a program capable of searching for parameters without performing complicated calculations.

物流センター等における商品の在庫数の時系列データの例を示す図である。FIG. 4 is a diagram showing an example of time-series data of the number of products in stock in a distribution center or the like; 実施の形態に係るパラメータ探索装置の構成例を示すブロック図である。1 is a block diagram showing a configuration example of a parameter search device according to an embodiment; FIG. 実施の形態に係る探索部によるマスの探索方法のイメージ例を示す図である。FIG. 7 is a diagram showing an image example of a method of searching for a square by a searching unit according to the embodiment; 実施の形態に係る探索部によるマスの探索方法の例を示す図である。It is a figure which shows the example of the search method of the square by the search part which concerns on embodiment. 実施の形態に係る探索部によるマスの探索方法の例を示す図である。It is a figure which shows the example of the search method of the square by the search part which concerns on embodiment. 実施の形態に係る探索部が図5の方法でマスを探索する場合の探索経路の例を示す図である。6 is a diagram showing an example of a search path when the searching unit according to the embodiment searches for squares by the method of FIG. 5; FIG. 実施の形態に係るパラメータ探索装置によるパラメータ探索方法の例を示すフロー図である。FIG. 4 is a flow chart showing an example of a parameter searching method by the parameter searching device according to the embodiment; 図7のステップS102の処理の例を示すフロー図である。FIG. 8 is a flowchart showing an example of processing in step S102 of FIG. 7; 実施の形態に係る探索部が図8の方法でマスを探索する場合の探索経路の例を示す図である。FIG. 9 is a diagram showing an example of a search path when the searching unit according to the embodiment searches for squares by the method of FIG. 8; 実施の形態に係る探索部によるマスの探索方法の例を示す図である。It is a figure which shows the example of the search method of the square by the search part which concerns on embodiment. 実施の形態に係るパラメータ探索装置を概念的に示したパラメータ探索装置の構成例を示すブロック図である。1 is a block diagram showing a configuration example of a parameter search device conceptually showing a parameter search device according to an embodiment; FIG.

以下、図面を参照して本開示の実施の形態について説明する。なお、以下の記載及び図面は、説明の明確化のため、適宜、省略及び簡略化がなされている。また、以下の各図面において、同一の要素には同一の符号が付されており、必要に応じて重複説明は省略されている。 Embodiments of the present disclosure will be described below with reference to the drawings. Note that the following descriptions and drawings are appropriately omitted and simplified for clarity of explanation. Further, in each drawing below, the same elements are denoted by the same reference numerals, and redundant description is omitted as necessary.

<実施の形態>
本実施の形態においては、物流センター等における商品の最大在庫日数及び安全在庫日数という2個のパラメータの最適解を探索するものとする。また、商品の最大在庫日数及び安全在庫日数を決定する指標は、商品の平均在庫日数、平均在庫金額、発注回数、及び欠品率の4個であるものとする。
<Embodiment>
In this embodiment, it is assumed that an optimum solution is searched for two parameters, ie, the maximum number of days in stock and the number of days in safe stock of products in a distribution center or the like. It is also assumed that the four indices for determining the maximum number of days in stock and the number of safe days in stock of a product are the average number of days in stock, the average amount of inventory, the number of orders placed, and the out-of-stock rate.

ここで、平均在庫日数は、平均在庫数分の商品が減少して0になるまでの日数である。平均在庫金額は、特定期間中の商品の在庫金額累積を特定期間中の日数で除算することで求められる金額である。発注回数は、特定期間中の商品を発注した回数である。欠品率は、特定期間中の商品の注文に対して、商品を提供できなかった割合である。欠品率は、欠品数(受注数-出荷数)を受注数で除算することで求められる。 Here, the average number of days in stock is the number of days until the average number of stocked items decreases to 0. The average inventory value is an amount obtained by dividing the accumulated inventory value of the product during the specified period by the number of days during the specified period. The number of orders is the number of times a product is ordered during a specific period. The out-of-stock rate is the ratio of products that could not be supplied to orders for products during a specific period. The out-of-stock rate is obtained by dividing the number of out-of-stock items (number of orders received - number of shipments) by the number of orders received.

まず、図2を参照して、本実施の形態に係るパラメータ探索装置10の構成について説明する。図2は、本実施の形態に係るパラメータ探索装置10の構成例を示すブロック図である。
図2に示されるように、本実施の形態に係るパラメータ探索装置10は、分割部11、探索部12、決定部13、及び入力部14を備えている。
First, the configuration of the parameter search device 10 according to the present embodiment will be described with reference to FIG. FIG. 2 is a block diagram showing a configuration example of the parameter search device 10 according to this embodiment.
As shown in FIG. 2, the parameter search device 10 according to the present embodiment includes a division section 11, a search section 12, a determination section 13, and an input section .

分割部11は、最大在庫日数及び安全在庫日数を軸とする所定の2次元空間をマス目状に分割する。 The dividing unit 11 divides a predetermined two-dimensional space with the maximum inventory days and the safety inventory days as axes into squares.

探索部12は、所定の2次元空間上のマスについて、目的関数fの値を計算する機能を備えている。目的関数fは、例えば、次の数式(1)で表される。

Figure 0007156666000001
ここで、k(i=1~4)は、各指標の桁合わせをするための除数であり、指標毎に計算される。また、ωは、各指標の重要度に応じて各指標の重み付けをするための重み係数であり、指標毎に設定される。 The search unit 12 has a function of calculating the value of the objective function f for a predetermined two-dimensional space. The objective function f is represented, for example, by the following formula (1).
Figure 0007156666000001
Here, k i (i=1 to 4) is a divisor for adjusting the digits of each index and is calculated for each index. ω i is a weighting factor for weighting each index according to the importance of each index, and is set for each index.

除数kは、例えば、探索部12により、後述の探索開始点となるマスにおける指標の値を基に、次の数式(2)で計算される。

Figure 0007156666000002
ここで、xは、指標である。また、floor(t)は、t以下の整数の中で最大の整数を求める床関数である。
これにより、後述の探索開始点となるマスでは、数式(1)におけるx/kの値は、0.5以上1.0未満になる。 The divisor k i is calculated, for example, by the searching unit 12 by the following formula (2) based on the value of the index in the cell serving as the search starting point, which will be described later.
Figure 0007156666000002
where x i is the index. Also, floor(t) is a floor function that obtains the largest integer among integers less than or equal to t.
As a result, the value of x i /k i in the equation (1) is 0.5 or more and less than 1.0 in the squares that serve as search starting points, which will be described later.

重み係数ωは、例えば、全指標を等価とみなす場合は全て1.0に設定する。また、ωは、例えば、欠品率の重要度が他の指標の重要度よりも高い場合は、欠品率のωだけを2.0に設定し、他の指標のωを1.0に設定するなどする。
重み係数ωは、例えば、ユーザにより入力部14を介して入力される。
For example, all weighting factors ω i are set to 1.0 when all indices are considered equivalent. Also, with respect to ω i , for example, if the importance of the out-of-stock rate is higher than the importance of other indices, only the out-of-stock rate ω i is set to 2.0, and the other indices ω i are set to 1 .0 for example.
The weighting factor ω i is input by the user through the input unit 14, for example.

探索部12は、所定の2次元空間上のマスについて、目的関数fの値を計算する場合、以下のようにして、計算を行う。
まず、探索部12は、該当するマスに対応する最大在庫日数及び安全在庫日数のそれぞれの値となるように、平均在庫日数、平均在庫金額、発注回数、及び欠品率のそれぞれの値を(シミュレーションなどで)計算する。
次に、探索部12は、上記で計算した平均在庫日数、平均在庫金額、発注回数、及び欠品率の値を、目的関数fに代入する。
これにより、該当するマスの目的関数fの値が計算される。
When the search unit 12 calculates the value of the objective function f for a predetermined two-dimensional space, the calculation is performed as follows.
First, the search unit 12 sets the average number of days in stock, the average number of days in stock, the number of orders, and the out-of-stock rate so that the values of the maximum number of days in stock and the number of days in safety stock correspond to the corresponding square ( (by simulation, etc.).
Next, the search unit 12 substitutes the values of the average number of days in stock, the average amount of inventory, the number of orders, and the out-of-stock rate calculated above into the objective function f.
Thereby, the value of the objective function f of the corresponding mass is calculated.

探索部12は、所定の2次元空間上のマスの中から、目的関数fの値が最小値となるマスを探索する。
決定部13は、目的関数fの値が最小値となるマスに対応する最大在庫日数及び安全在庫日数のそれぞれの値を最適解に決定する。
入力部14は、ユーザにより各種入力が行われる部分であり、例えば、上述のように、指標x毎に設定される重み係数ωが入力される。
The search unit 12 searches for a square having the minimum value of the objective function f from among the squares in a predetermined two-dimensional space.
The determining unit 13 determines the values of the maximum number of days in stock and the number of safety days in stock corresponding to the mass having the minimum value of the objective function f as the optimum solution.
The input unit 14 is a part where various inputs are performed by the user, and for example, the weighting coefficient ω i set for each index x i is input as described above.

ここで、図3を参照して、本実施の形態に係る探索部12によるマスの探索方法のイメージについて説明する。図3は、本実施の形態に係る探索部12によるマスの探索方法のイメージ例を示す図である。なお、図3においては、所定の2次元空間の2軸を最大在庫日数及び安全在庫日数とし、この所定の2次元空間に対して垂直な軸を目的関数fとしている。 Here, with reference to FIG. 3, an image of a method of searching for a square by the searching unit 12 according to the present embodiment will be described. FIG. 3 is a diagram showing an image example of a square search method by the search unit 12 according to the present embodiment. In FIG. 3, the two axes of a predetermined two-dimensional space are the maximum inventory days and the safe inventory days, and the axis perpendicular to this predetermined two-dimensional space is the objective function f.

図3に示されるように、探索部12は、所定の2次元空間上のマスの1つを探索開始点とし、目的関数fの値が最も大きく減少する方向(すなわち、傾斜が最も急な方向)にマスを移動させながら、目的関数fの値が最小値となるマスを探索する。 As shown in FIG. 3, the search unit 12 uses one of the squares in a predetermined two-dimensional space as a search starting point, and uses the direction in which the value of the objective function f decreases the most (that is, the direction in which the slope is steepest). ), searching for the square with the minimum value of the objective function f.

ここで、図4及び図5を参照して、本実施の形態に係る探索部12によるマスの探索方法について詳細に説明する。図4及び図5は、本実施の形態に係る探索部12によるマスの探索方法の例を示す図である。 Here, with reference to FIGS. 4 and 5, a method of searching for a square by the search unit 12 according to the present embodiment will be described in detail. 4 and 5 are diagrams showing an example of a method of searching for a square by the searching unit 12 according to the present embodiment.

図4の例では、まず、探索部12は、所定の2次元空間上のマスの1つを探索開始点に設定する。探索開始点のマスは、例えば、所定の2次元空間上のマスのうち、目的関数fの値が最大値となると推測されるマスに設定される。 In the example of FIG. 4, the searching unit 12 first sets one of the squares in a predetermined two-dimensional space as a search starting point. The square of the search starting point is set, for example, to a square that is estimated to have the maximum value of the objective function f among squares in a predetermined two-dimensional space.

続いて、探索部12は、探索開始点のマス(現マスとする)及びその周囲8方向(すなわち、全周囲)のマスの目的関数fの値を計算し、周囲8方向のマスのうち目的関数fの値が最小であるマスに移動し、移動先のマスを新たな現マスとする。 Subsequently, the search unit 12 calculates the values of the objective function f for the square at the search start point (assumed to be the current square) and the squares in eight directions around it (that is, all around). Move to the square with the smallest value of the function f, and set the destination square as the new current square.

以降、探索部12は、現マス及びその周囲8方向のマスの目的関数fの値を計算し、現マスの周囲8方向のマスの中から、現マスよりも目的関数fの値が小さいマスを探索する探索処理を再帰的に実行する。ただし、探索範囲から外れるマス(日数が0未満になるマスなど)は探索対象外とする。 After that, the search unit 12 calculates the values of the objective function f of the current square and the squares in the eight directions around it, and selects the squares with smaller objective function f values than the current square from among the squares in the eight directions around the current square. recursively search for . However, squares outside the search range (such as squares where the number of days is less than 0) are excluded from the search.

探索部12は、探索処理で探索対象であった周囲のマスを見つけられなかった場合、現マスが所定の2次元空間上のマスの中で目的関数fの値が最小値となるマスであると判断し、マスの探索を終了する。 When the searching unit 12 fails to find the surrounding squares to be searched in the search processing, the current square is the square having the minimum value of the objective function f among the squares in the predetermined two-dimensional space. , and ends the search for the square.

ただし、図4の例では、探索部12は、マスを移動する度に、現マスの周囲8方向のマスの目的関数fの値を全て計算するため、マスの探索時間が増大してしまうことが懸念される。 However, in the example of FIG. 4, the searching unit 12 calculates all the values of the objective function f of the squares in the eight directions around the current square every time the square is moved, so the search time for the square increases. is concerned.

図5の例は、図4の例と比較して、マスの探索時間の短縮化を図るものである。図5の例では、探索開始点のマスから移動した後に行う探索処理は、斜め探索モードと、上下左右探索モードと、の2個の探索モードを含んでいる。 The example in FIG. 5 is intended to shorten the search time for cells compared to the example in FIG. In the example of FIG. 5, the search processing performed after moving from the search start point square includes two search modes, an oblique search mode and an up/down/left/right search mode.

斜め探索モードは、現マス及びその周囲の斜め4方向のマスの目的関数fの値を計算し、現マスの周囲の斜め4方向のマスの中から、現マスよりも目的関数fの値が小さいマスを探索する処理である。 In the diagonal search mode, the values of the objective function f are calculated for the current square and its surrounding squares in four diagonal directions. This is the process of searching for small squares.

上下左右探索モードは、現マス及びその周囲の上下左右4方向のマスの目的関数fの値を計算し、現マスの周囲の上下左右4方向のマスの中から、現マスよりも目的関数fの値が小さいマスを探索する処理である。 In the up, down, left, and right search mode, the values of the objective function f of the current square and the surrounding squares in four directions, up, down, left, and right, are calculated. This is a process of searching for a square with a small value of .

図5の例は、探索開始点のマスから移動した移動方向が、斜め4方向のいずれかである場合の例である。この場合、探索開始点のマスの周囲のマスのうち目的関数fの値が最小であるのは、斜め方向のマスであることになる。そのため、探索開始点の近傍領域で最も傾斜が急な方向は、探索開始点のマスから斜め方向である確率が高いと考えられる。そのため、探索部12は、斜め探索モードを実行する。 The example of FIG. 5 is an example in which the direction of movement from the square of the search start point is one of four oblique directions. In this case, among the squares around the square of the search start point, the squares in the oblique direction have the smallest value of the objective function f. Therefore, it is highly probable that the steepest direction in the vicinity of the search start point is the diagonal direction from the search start point square. Therefore, the search unit 12 executes the diagonal search mode.

探索部12は、斜め探索モードで探索対象であった周囲のマスを見つけた場合は、斜め探索モードで見つけた周囲のマスのうち目的関数fの値が最小である周囲のマスに移動し、移動先のマスを新たな現マスとして、斜め探索モードを再帰的に実行する。このように、斜め探索モードで探索対象であった周囲のマスを見つけた場合は、上下左右4方向のマスの計算を省略できるため、マスの探索時間の短縮化を図ることが可能となる。 When the search unit 12 finds the surrounding squares that were the search target in the diagonal search mode, the search unit 12 moves to the surrounding squares having the smallest value of the objective function f among the surrounding squares found in the diagonal search mode, The diagonal search mode is recursively executed with the destination square as the new current square. In this way, when the surrounding squares to be searched in the oblique search mode are found, calculation of squares in four directions, up, down, left, and right, can be omitted, so that the square search time can be shortened.

一方、探索部12は、斜め探索モードで探索対象であった周囲のマスを見つけられなかった場合は、斜め探索モードを上下左右探索モードに切り替え、上下左右探索モードを実行する。探索部12は、上下左右探索モードで探索対象であった周囲のマスを見つけた場合は、上下左右探索モードで見つけた周囲のマスのうち目的関数fの値が最小である周囲のマスに移動し、移動先のマスを新たな現マスとして、上下左右探索モードを再帰的に実行する。逆に、探索部12は、上下左右探索モードでも探索対象であった周囲のマスを見つけられなかった場合は、現マスが所定の2次元空間上のマスの中で目的関数fの値が最小値となるマスであると判断し、マスの探索を終了する。 On the other hand, if the search unit 12 cannot find the surrounding squares that were the search target in the oblique search mode, it switches the oblique search mode to the up, down, left, and right search mode, and executes the up, down, left, and right search mode. When the search unit 12 finds the surrounding squares that were the search target in the up/down/left/right search mode, the search unit 12 moves to the surrounding squares having the smallest value of the objective function f among the surrounding squares found in the up/down/left/right search mode. Then, the up/down/left/right search mode is recursively executed with the destination square as the new current square. Conversely, if the search unit 12 cannot find the surrounding squares that were the search target even in the up, down, left, and right search mode, the search unit 12 determines that the current square is the smallest value of the objective function f among the squares in the predetermined two-dimensional space. It is determined that the square is a value, and the search for the square ends.

このように、図5の例では、探索部12は、斜め探索モードと上下左右探索モードとを切り替えながら、探索処理を実行する。そして、探索部12は、斜め探索モード及び上下左右探索モードの両方で探索対象であった周囲のマスを見つけられなかった場合、現マスが所定の2次元空間上のマスの中で目的関数fの値が最小値となるマスであると判断する。 In this manner, in the example of FIG. 5, the search unit 12 performs the search process while switching between the diagonal search mode and the up/down/left/right search mode. Then, if the search unit 12 fails to find the surrounding squares that were the search target in both the diagonal search mode and the up/down/left/right search mode, the search unit 12 determines that the current square is a square in a predetermined two-dimensional space with the objective function f value is the minimum value.

図6は、本実施の形態に係る探索部12が図5の方法でマスを探索する場合の探索経路の例を示す図である。
図6の例では、左上のマスを探索開始点とし、まずは、斜め探索モードを実行し、続いて、上下左右探索モードへの切り替えを行い、その後、再度、斜め探索モードへの切り替えを行って、右下のマスで探索を終了している。すなわち、図6の例では、所定の2次元空間上のマスの中で、右下のマスが、目的関数fの値が最小値となるマスであると判断されている。
FIG. 6 is a diagram showing an example of a search path when the searching unit 12 according to this embodiment searches for a square by the method of FIG.
In the example of FIG. 6, the upper left square is set as the search starting point, and the diagonal search mode is first executed, followed by switching to the up/down/left/right search mode, and then switching to the diagonal search mode again. , the search ends in the lower right square. That is, in the example of FIG. 6, among the squares in the predetermined two-dimensional space, the lower right square is determined to be the square having the minimum value of the objective function f.

以下では、図7を参照して、本実施の形態に係るパラメータ探索装置10によるパラメータ探索方法について説明する。図7は、本実施の形態に係るパラメータ探索装置10によるパラメータ探索方法の例を示すフロー図である。 A parameter searching method by the parameter searching device 10 according to the present embodiment will be described below with reference to FIG. FIG. 7 is a flowchart showing an example of a parameter searching method by the parameter searching device 10 according to this embodiment.

図7に示されるように、分割部11は、最大在庫日数及び安全在庫日数を軸とする所定の2次元空間をマス目状に分割する(ステップS101)。例えば、1日単位で分割を行う。
探索部12は、所定の2次元空間上のマスの中から、目的関数fの値が最小値となるマスを探索する(ステップS102)。
決定部13は、目的関数fの値が最小値となるマスに対応する最大在庫日数及び安全在庫日数のそれぞれの値を最適解に決定する(ステップS103)。
As shown in FIG. 7, the dividing unit 11 divides a predetermined two-dimensional space around the maximum number of days in stock and the number of days in safety stock into squares (step S101). For example, division is performed in units of one day.
The searching unit 12 searches for a square having the minimum value of the objective function f from among the squares in a predetermined two-dimensional space (step S102).
The determination unit 13 determines the values of the maximum number of days in stock and the number of days in safety corresponding to the mass having the minimum value of the objective function f as the optimum solution (step S103).

続いて、図8を参照して、図7のステップS102の処理について詳細に説明する。図8は、図7のステップS102の処理の例を示すフロー図である。なお、図8は、図5の方法によりマスの探索を行う場合の処理を示している。 Next, referring to FIG. 8, the process of step S102 in FIG. 7 will be described in detail. FIG. 8 is a flow chart showing an example of the processing in step S102 of FIG. Note that FIG. 8 shows the processing when searching for a square by the method of FIG.

図8に示されるように、まず、探索部12は、所定の2次元空間上のマスの1つを探索開始点に決定する(ステップS201)。続いて、探索部12は、探索開始点のマス及びその周囲8方向のマスの目的関数fの値を計算し(ステップS202)、周囲8方向のマスのうち目的関数fの値が最小であるマスに移動し、移動先のマスを新たな現マスとする(ステップS203)。 As shown in FIG. 8, the search unit 12 first determines one of the squares in a predetermined two-dimensional space as a search start point (step S201). Subsequently, the search unit 12 calculates the value of the objective function f for the square at the search start point and the squares in the eight directions around it (step S202). Move to a square and set the destination square as a new current square (step S203).

続いて、探索部12は、探索開始点のマスから移動した移動方向(上下左右4方向のいずれか、又は、斜め4方向のいずれか)を判断する(ステップS204)。 Subsequently, the search unit 12 determines the direction of movement (one of the four directions up, down, left, and right, or one of the four diagonal directions) from the square of the search start point (step S204).

探索開始点のマスから移動した移動方向が上下左右4方向のいずれかであった場合、探索部12は、上下左右探索モードを実行する。すなわち、探索部12は、現マス及びその周囲の上下左右4方向のマスの目的関数fの値を計算し(ステップS205)、現マスの周囲の上下左右4方向のマスの中から、現マスよりも目的関数fの値が小さいマスを探索する(ステップS206)。 If the direction of movement from the square of the search start point is one of the four directions of up, down, left, and right, the search unit 12 executes the up, down, left, and right search mode. That is, the search unit 12 calculates the values of the objective function f of the current square and the squares in four directions, up, down, left, and right (step S205). A square with a smaller value of the objective function f than is searched for (step S206).

探索部12は、上下左右探索モードで探索対象であった周囲のマスを見つけた場合は(ステップS206のYes)、上下左右探索モードで見つけた周囲のマスのうち目的関数fの値が最小である周囲のマスに移動し(ステップS207)、移動先のマスを新たな現マスとする。そして、探索部12は、ステップS205に移行し、上下左右探索モードを再帰的に実行する。 When the search unit 12 finds the surrounding squares that were the search target in the up/down/left/right search mode (Yes in step S206), the search unit 12 finds the smallest value of the objective function f among the surrounding squares found in the up/down/left/right search mode. Move to a surrounding square (step S207), and set the new square as the new current square. Then, the search unit 12 proceeds to step S205 and recursively executes the up/down/left/right search mode.

探索部12は、上下左右探索モードで探索対象であった周囲のマスを見つけられなかった場合は(ステップS206のNo)、続いて、斜め探索モードを実行する。すなわち、探索部12は、現マス及びその周囲の斜め4方向のマスの目的関数fの値を計算し(ステップS208)、現マスの周囲の斜め4方向のマスの中から、現マスよりも目的関数fの値が小さいマスを探索する(ステップS209)。 If the search unit 12 cannot find the surrounding squares that were the search target in the up/down/left/right search mode (No in step S206), the search unit 12 then executes the diagonal search mode. That is, the search unit 12 calculates the values of the objective function f of the current square and its surrounding squares in four diagonal directions (step S208), A square with a small value of the objective function f is searched (step S209).

探索部12は、斜め探索モードで探索対象であった周囲のマスを見つけた場合は(ステップS209のYes)、斜め探索モードで見つけた周囲のマスのうち目的関数fの値が最小である周囲のマスに移動し(ステップS210)、移動先のマスを新たな現マスとする。そして、探索部12は、ステップS211に移行し、斜め探索モードを再帰的に実行する。 When the search unit 12 finds the surrounding squares that were the search target in the diagonal search mode (Yes in step S209), the searching unit 12 searches the surrounding squares having the smallest value of the objective function f among the surrounding squares found in the diagonal search mode. (step S210), and the destination square is set as a new current square. Then, the search unit 12 proceeds to step S211 and recursively executes the diagonal search mode.

探索部12は、斜め探索モードでも探索対象であった周囲のマスを見つけられなかった場合は(ステップS209のNo)、現マスが所定の2次元空間上のマスの中で目的関数fの値が最小値となるマスであると判断し、マスの探索を終了する。 If the search unit 12 cannot find the surrounding squares that were the search target even in the diagonal search mode (No in step S209), the search unit 12 determines that the current square is the value of the objective function f among the squares in the predetermined two-dimensional space. is determined to be the square with the minimum value, and the search for the square ends.

一方、探索開始点のマスから移動した移動方向が斜め4方向のいずれかであった場合、探索部12は、上記のステップS209,S210と同様の斜め探索モードを実行する(ステップS211,S212)。 On the other hand, if the movement direction from the search start point square is one of the four diagonal directions, the search unit 12 executes the diagonal search mode similar to steps S209 and S210 (steps S211 and S212). .

探索部12は、斜め探索モードで探索対象であった周囲のマスを見つけた場合は(ステップS212のYes)、斜め探索モードで見つけた周囲のマスのうち目的関数fの値が最小である周囲のマスに移動し(ステップS213)、移動先のマスを新たな現マスとする。そして、探索部12は、ステップS211に移行し、斜め探索モードを再帰的に実行する。 If the search unit 12 finds the surrounding squares that were the search target in the diagonal search mode (Yes in step S212), the searching unit 12 searches the surrounding squares that have the smallest value of the objective function f among the surrounding squares found in the diagonal search mode. (step S213), and the destination square is set as a new current square. Then, the search unit 12 proceeds to step S211 and recursively executes the diagonal search mode.

探索部12は、斜め探索モードで探索対象であった周囲のマスを見つけられなかった場合は(ステップS212のNo)、続いて、上記のステップS206,S207と同様の上下左右探索モードを実行する(ステップS214,S215)。 If the search unit 12 cannot find the surrounding squares to be searched in the oblique search mode (No in step S212), then the search unit 12 executes the up/down/left/right search mode similar to steps S206 and S207 described above. (Steps S214, S215).

探索部12は、上下左右探索モードで探索対象であった周囲のマスを見つけた場合は(ステップS215のYes)、上下左右探索モードで見つけた周囲のマスのうち目的関数fの値が最小である周囲のマスに移動し(ステップS216)、移動先のマスを新たな現マスとする。そして、探索部12は、ステップS205に移行し、上下左右探索モードを再帰的に実行する。 If the search unit 12 finds the surrounding squares that were the search target in the up/down/left/right search mode (Yes in step S215), the search unit 12 finds the smallest value of the objective function f among the surrounding squares found in the up/down/left/right search mode. Move to a surrounding square (step S216), and set the destination square as a new current square. Then, the search unit 12 proceeds to step S205 and recursively executes the up/down/left/right search mode.

探索部12は、上下左右探索モードでも探索対象であった周囲のマスを見つけられなかった場合は(ステップS215のNo)、現マスが所定の2次元空間上のマスの中で目的関数fの値が最小値となるマスであると判断し、マスの探索を終了する。 If the search unit 12 cannot find the surrounding squares that were the search target even in the up/down/left/right search mode (No in step S215), the search unit 12 determines that the current square is the target function f among the squares in the predetermined two-dimensional space. It is determined that the square has the minimum value, and the search for the square ends.

図9は、本実施の形態に係る探索部12が図8の方法でマスを探索する場合の探索経路の例を示す図である。
図9の例では、探索部12は、最大在庫日数が20日で安全在庫日数が10日の右下のマスを探索開始点に決定している。
また、探索開始点のマスでは、その周囲8方向のマスのうち目的関数fの値が最小であるマスが、左方向のマスである。
FIG. 9 is a diagram showing an example of a search path when the searching unit 12 according to the present embodiment searches for squares by the method of FIG.
In the example of FIG. 9, the search unit 12 has determined the lower right square with the maximum inventory days of 20 days and the safety inventory days of 10 days as the search start point.
In addition, in the search starting point, the cell with the smallest value of the objective function f among the eight surrounding cells is the cell in the left direction.

そのため、探索部12は、探索開始点のマスから左方向のマスに移動し、以降、上下左右探索モードを再帰的に実行し、最大在庫日数が16日で安全在庫日数が10日のマスまで移動している。 Therefore, the search unit 12 moves from the search start point to the leftward square, and thereafter recursively executes the up, down, left, and right search mode until the maximum inventory days is 16 days and the safety inventory days is 10 days. moving.

ただし、最大在庫日数が16日で安全在庫日数が10日の現マスでは、その周囲の上下左右方向のマスの中に、現マスよりも目的関数fの値が小さいマスはない。一方で、現マスでは、その周囲の斜め方向のマスの中に、現マスよりも目的関数fの値が小さいマスがあり、また、これらのマスのうち目的関数fの値が最小である周囲のマスが、左上斜め方向のマスである。 However, in the current cell with the maximum inventory days of 16 days and the safety inventory days of 10 days, none of the cells in the vertical and horizontal directions surrounding it has a smaller value of the objective function f than the current cell. On the other hand, in the current cell, among the cells in the oblique direction around it, there is a cell with a smaller objective function f value than the current cell. is the square in the upper left diagonal direction.

そのため、探索部12は、現マスから左上斜めのマスに移動し、以降、斜め探索モードを再帰的に実行し、最大在庫日数が9日で安全在庫日数が3日のマスまで移動している。 Therefore, the search unit 12 moves from the current square to the upper left diagonal square, recursively executes the diagonal search mode, and moves to the square with the maximum inventory days of 9 days and the safety inventory days of 3 days. .

ただし、最大在庫日数が9日で安全在庫日数が3日の現マスでは、その周囲の斜め方向のマスの中に、現マスよりも目的関数fの値が小さいマスはない。また、現マスでは、その周囲の上下左右方向のマスの中にも、現マスよりも目的関数fの値が小さいマスはない。
そのため、探索部12は、現マスが、目的関数fの値が最小値となるマスであると判断し、決定部13は、現マスに対応する最大在庫日数の値(9日)及び安全在庫日数の値(3日)を、それぞれ、最大在庫日数及び安全在庫日数の最適解に決定する。
However, for the current cell with the maximum inventory days of 9 days and the safety inventory days of 3 days, there is no cell with a smaller value of the objective function f than the current cell among the cells in the oblique direction around it. Further, in the current cell, there is no cell having a smaller value of the objective function f than the current cell in the cells in the vertical and horizontal directions around it.
Therefore, the search unit 12 determines that the current mass is the mass with the minimum value of the objective function f, and the determination unit 13 determines the value of the maximum inventory days (9 days) corresponding to the current mass and the safety stock The value of days (3 days) is determined to be the optimal solution for maximum days of supply and safety days of supply, respectively.

上述したように本実施の形態によれば、パラメータ探索装置10は、最大在庫日数及び安全在庫日数を軸とする所定の2次元空間をマス目状に分割し、所定の2次元空間上のマスの中から、目的関数fの値が最小値となるマスを探索し、見つけたマスに対応する最大在庫日数及び安全在庫日数のそれぞれの値を最適解に決定する。そのため、複雑な計算を行うことなく、最大在庫日数及び安全在庫日数という2個のパラメータのそれぞれの最適解を探索することができる。 As described above, according to the present embodiment, the parameter search device 10 divides a predetermined two-dimensional space with the maximum inventory days and the safety inventory days as axes into squares, and divides the predetermined two-dimensional space into squares. Among them, the mass with the minimum value of the objective function f is searched, and the respective values of the maximum number of days in stock and the number of days in safety stock corresponding to the found mass are determined as optimum solutions. Therefore, it is possible to search for the optimum solution for each of the two parameters, the maximum number of days in stock and the number of safety days in stock, without performing complicated calculations.

したがって、経験が不足しているユーザでも、最大在庫日数及び安全在庫日数を適切に設定することができるため、適切に発注点数を決定し、決定した発注点数まで商品の在庫数が減少した適切なタイミングで、商品の発注を行うことができる。これにより、商品の在庫数が不足して、商品の品切れや欠品が発生することを防ぐことに寄与し得る。 Therefore, even an inexperienced user can set the maximum number of days in stock and the number of days in safety stock appropriately. You can place an order for the product in time. This can contribute to preventing the occurrence of out-of-stock or shortage of products due to shortage of the number of products in stock.

<実施の形態の変形例>
上記の実施の形態では、目的関数fの指標x毎の桁合わせ用の除数kは、探索開始点のマスにおける指標xの値を基に計算し、以降、マスが移動しても、除数kは固定されていた。しかし、この除数kは、探索開始点のマスにおける指標xの値を基に計算したものであるため、マスの移動が進むにしたがって、適切な値にならない可能性がある。
そのため、本開示では、除数kを、マスを移動する度に計算しても良い。この場合、移動先のマスで周囲を探索するのに使用する目的関数fの指標x毎の除数kを、移動先のマスにおける指標xの値とすれば良い。
<Modified example of the embodiment>
In the above embodiment, the divisor k i for digit alignment for each index x i of the objective function f is calculated based on the value of the index x i in the square at the search start point, and thereafter, even if the square moves, , the divisor k i was fixed. However, since this divisor k i is calculated based on the value of the index x i in the search starting point square, it may not become an appropriate value as the square moves.
Therefore, in this disclosure, the divisor k i may be calculated each time a square is moved. In this case, the divisor k i for each index x i of the objective function f used for searching the surroundings in the destination cell should be the value of the index x i in the destination cell.

また、上記の実施の形態では、最大在庫日数及び安全在庫日数を決定する指標が、平均在庫日数、平均在庫金額、発注回数、及び欠品率の4個である例について説明したが、本開示はこれには限定されない。本開示では、最大在庫日数及び安全在庫日数を決定する指標は、これら4個の指標のうち少なくとも1個を含む、複数個の指標であれば良い。 In addition, in the above embodiment, an example in which the indicators for determining the maximum number of days in stock and the number of days in safety stock are the average number of days in stock, the average amount of inventory, the number of orders, and the out-of-stock rate has been described, but the present disclosure is not limited to this. In the present disclosure, the indices for determining the maximum number of days in stock and the number of safety days in stock may be a plurality of indices including at least one of these four indices.

また、上記の実施の形態では、最適解を探索するパラメータが、商品の最大在庫日数及び安全在庫日数の2個である例について説明したが、本開示はこれには限定されない。本開示では、最適解を探索するパラメータが、物流センター等における在庫管理の分野以外の他分野のパラメータであっても良い。この場合、他分野のパラメータを決定する指標は、該当する分野に応じて適宜決定すれば良い。また、本開示は、N(Nは2以上の自然数)個のパラメータの最適解を探索する場合に適用可能であり、Nは2に限定されない。この場合、N個のパラメータを軸とする所定のN次元空間をマス目状に分割し、所定のN次元空間上のマスの中から、目的関数の値が最小値となるマスを探索し、目的関数の値が最小値となるマスに対応するN個のパラメータのそれぞれの値を最適解に決定する。 Further, in the above-described embodiment, an example has been described in which two parameters for searching for the optimum solution are the maximum inventory days and the safety inventory days of the product, but the present disclosure is not limited to this. In the present disclosure, the parameters for searching for the optimum solution may be parameters in fields other than the field of inventory management in distribution centers and the like. In this case, the indices for determining the parameters of other fields may be appropriately determined according to the relevant fields. In addition, the present disclosure is applicable to searching for optimal solutions for N (N is a natural number equal to or greater than 2) parameters, and N is not limited to 2. In this case, a predetermined N-dimensional space with N parameters as an axis is divided into squares, and a square having the minimum value of the objective function is searched from among the squares in the predetermined N-dimensional space, The value of each of the N parameters corresponding to the mass with the minimum value of the objective function is determined as the optimum solution.

<Nが2以外になる事例>
以下では、Nが2以外になる具体的な事例について説明する。
本事例は、本開示を朝食の栄養管理の分野に適用したものである。本事例では、例えば、寮の朝食が幾つかのメニューから選択する方式である場合に、塩分、糖分などを抑え、かつ、食費も安くなるように栄養管理を行う。
<Case where N is other than 2>
A specific case where N is other than 2 will be described below.
This example is an application of the present disclosure to the field of breakfast nutrition management. In this case, for example, when the breakfast in the dormitory is a system in which a choice is made from several menus, nutritional management is performed so as to reduce salt and sugar content and reduce food costs.

本事例では、パラメータは、各メニューを選択する日数とする。具体的には、パラメータは、以下の4個とする(すなわち、Nが4)。
納豆ごはん …… N日/月
鮭の塩焼き …… N日/月
ハムエッグ …… N日/月
シリアル …… N日/月
ここで、N~Nの拘束条件は以下とする。
+N+N+N=30(又は、31)
また、上記の各メニューの日数を決定する指標は、食費、塩分摂取量、糖分摂取量、及びコレステロール摂取量の4個とする。
In this example, the parameter is the number of days to select each menu. Specifically, the following four parameters are used (that is, N is 4).
Rice with natto …… N 1 day/month Salmon grilled with salt …… N 2 days/month Ham and eggs …… N 3 days/month Cereal …… N 4 days/month Here, the constraint conditions for N 1 to N 4 are as follows. .
N1 + N2 +N3+ N4 = 30 (or 31)
In addition, four indices for determining the number of days for each menu are food expenses, salt intake, sugar intake, and cholesterol intake.

本事例の場合、所定の4次元空間上のマスの中から、目的関数fの値が最小値となるマスを探索していく。
本事例の目的関数fは、上記の数式(1)に準じたもので、食費、塩分摂取量、糖分摂取量、及びコレステロール摂取量の4個の指標を用いた関数とする。
In the case of this example, a square having the minimum value of the objective function f is searched from squares in a predetermined four-dimensional space.
The objective function f in this example conforms to the above formula (1), and is a function using the four indices of food expenses, salt intake, sugar intake, and cholesterol intake.

本事例において、所定の4次元空間上のマスについて、目的関数fの値を計算する場合、まず、該当するマスに対応する4個のメニューの日数のそれぞれの値となるように、食費、塩分摂取量、糖分摂取量、及びコレステロール摂取量のそれぞれの値を計算する。次に、上記で計算した食費、塩分摂取量、糖分摂取量、及びコレステロール摂取量の値を、目的関数fに代入する。これにより、該当するマスの目的関数fの値が計算される。 In this example, when calculating the value of the objective function f for a predetermined square in a four-dimensional space, first, food expenses, salt content Values for intake, sugar intake, and cholesterol intake are calculated. Next, the values of food expenses, salt intake, sugar intake, and cholesterol intake calculated above are substituted into the objective function f. As a result, the value of the objective function f of the corresponding mass is calculated.

<N次元空間上のマスの探索処理>
以下では、所定のN次元空間上のマスの探索処理について、詳細に説明する。
(A)2次元空間の場合(すなわち、Nが2の場合)
上記の実施の形態は、2次元空間の場合の例であった。
上記の実施の形態のように2次元空間の場合、探索開始点のマスから移動した後に行う探索処理においては、図5を用いて説明したように、上下左右探索モード及び斜め探索モードという、探索方向が互いに異なる2個の探索モードを切り替えて実行する。この2個の探索モードは、言い換えれば、以下のようになる。
1番目の探索モード(上下左右探索モード):
2軸から1軸を取出し、取出した1軸の軸方向に現マスから±1した位置にあるマスを探索する(他の軸は0)。探索方向は、×2=4方向となる。
2番目の探索モード(斜め探索モード):
2軸から2軸を取出し、取出した2軸のそれぞれの軸方向に現マスから±1した位置にあるマスを探索する。探索方向は、×2=4方向となる。
なお、2次元空間の場合、現マスの全周囲は、3-1=8方向となる。
<Processing to search for squares in N-dimensional space>
In the following, a detailed description will be given of the process of searching for a square in a predetermined N-dimensional space.
(A) For two-dimensional space (i.e., when N is 2)
The above embodiment is an example of a two-dimensional space.
In the case of a two-dimensional space as in the above embodiment, in the search processing performed after moving from the search start point square, as described with reference to FIG. Two search modes with different directions are switched and executed. These two search modes are, in other words, as follows.
First search mode (up/down/left/right search mode):
One axis is extracted from the two axes, and a square located at a position ±1 from the current square in the axial direction of the extracted one axis is searched (other axes are 0). The search directions are 2 C 1 ×2 1 =4 directions.
Second search mode (diagonal search mode):
2 axes are extracted from the 2 axes, and a square located at a position ±1 from the current square in each axial direction of the extracted 2 axes is searched. The search directions are 2 C 2 ×2 2 =4 directions.
In the case of a two-dimensional space, the total circumference of the current square is 3 2 -1=8 directions.

(B)3次元空間の場合(すなわち、Nが3の場合)
3次元空間の場合、探索開始点のマスから移動した後に行う探索処理においては、図10に示されるように、探索方向が互いに異なる3個の探索モードを切り替えて実行する。この3個の探索モードは、以下のようになる。
1番目の探索モード:
3軸から1軸を取出し、取出した1軸の軸方向に現マスから±1した位置にあるマスを探索する(他の軸は0)。探索方向は、×2=6方向となる。
2番目の探索モード:
3軸から2軸を取出し、取出した2軸のそれぞれの軸方向に現マスから±1した位置にあるマスを探索する(他の軸は0)。探索方向は、×2=12方向となる。
3番目の探索モード:
3軸から3軸を取出し、取出した3軸のそれぞれの軸方向に現マスから±1した位置にあるマスを探索する。探索方向は、×2=8方向となる。
なお、3次元空間の場合、現マスの全周囲は、3-1=26方向となる。
(B) For 3-dimensional space (i.e., when N is 3)
In the case of a three-dimensional space, in the search processing performed after moving from the search start point square, as shown in FIG. 10, three search modes with mutually different search directions are switched and executed. The three search modes are as follows.
First search mode:
One axis is taken out of the three axes, and a mass located at a position ±1 from the current mass in the axial direction of the taken out one axis is searched (other axes are 0). The search directions are 3 C 1 ×2 1 =6 directions.
Second search mode:
2 axes are extracted from the 3 axes, and a square located at a position ±1 from the current square in each axial direction of the extracted 2 axes is searched (other axes are 0). The search directions are 3 C 2 ×2 2 =12 directions.
Third search mode:
3 axes are taken out from the 3 axes, and a square located at a position ±1 from the current square in each axial direction of the taken out 3 axes is searched. The search directions are 3 C 3 ×2 3 =8 directions.
In the case of a three-dimensional space, there are 3 3 −1=26 directions around the current square.

(C)N次元空間の場合
2次元空間及び3次元空間の場合の探索処理は、上記の通りである。
これを踏まえると、N次元空間の場合、探索開始点のマスから移動した後に行う探索処理においては、探索方向が互いに異なるN個の探索モードを切り替えて実行する。このN個の探索モードは、以下のようになる。
1番目の探索モード:
N軸から1軸を取出し、取出した1軸の軸方向に現マスから±1した位置にあるマスを探索する(他の軸は0)。探索方向は、×2方向となる。
2番目の探索モード:
N軸から2軸を取出し、取出した2軸のそれぞれの軸方向に現マスから±1した位置にあるマスを探索する(他の軸は0)。探索方向は、×2方向となる。
N番目の探索モード:
N軸からN軸を取出し、取出したN軸のそれぞれの軸方向に現マスから±1した位置にあるマスを探索する。探索方向は、×2方向となる。
なお、N次元空間の場合、現マスの全周囲の方向は、次の数式(3)で表される。

Figure 0007156666000003
(C) Case of N-Dimensional Space Search processing for two-dimensional space and three-dimensional space is as described above.
Based on this, in the case of the N-dimensional space, in the search processing performed after moving from the search start point square, N search modes with different search directions are switched and executed. The N search modes are as follows.
First search mode:
One axis is taken out from the N axis, and a mass located at a position ±1 from the current mass in the axial direction of the taken out one axis is searched (other axes are 0). The search direction is the NC 1 ×2 1 direction.
Second search mode:
Two axes are extracted from the N axis, and a mass located at a position ±1 from the current mass is searched for in each axial direction of the extracted two axes (other axes are 0). The search direction is NC 2 ×2 2 directions.
Nth search mode:
The N-axis is extracted from the N-axis, and a mass located at a position ±1 from the current mass is searched for in each axial direction of the extracted N-axis. The search direction is the NCN×2N direction .
In the case of an N-dimensional space, the directions around the current square are represented by the following formula (3).
Figure 0007156666000003

<実施の形態の概念>
続いて、上記の実施の形態に係るパラメータ探索装置10を概念的に示した構成について説明する。図11は、上記の実施の形態に係るパラメータ探索装置10を概念的に示したパラメータ探索装置100の構成例を示すブロック図である。
<Concept of Embodiment>
Next, a configuration conceptually showing the parameter search device 10 according to the above embodiment will be described. FIG. 11 is a block diagram showing a configuration example of a parameter searching device 100 conceptually showing the parameter searching device 10 according to the above embodiment.

図11に示されるように、パラメータ探索装置100は、分割部101、探索部102、及び決定部103を備えている。分割部101は、分割部11に相当し、探索部102は、探索部12に相当し、決定部103は、決定部13に相当する。 As shown in FIG. 11, the parameter searching device 100 comprises a dividing section 101, a searching section 102, and a determining section 103. FIG. The dividing section 101 corresponds to the dividing section 11 , the searching section 102 corresponds to the searching section 12 , and the determining section 103 corresponds to the determining section 13 .

分割部101は、所定のN次元空間をマス目状に分割する。
探索部102は、所定のN次元空間上のマスの中から、目的関数の値が最小値となるマスを探索する。
決定部103は、目的関数の値が最小値となるマスに対応するN個のパラメータのそれぞれの値を最適解に決定する。
A dividing unit 101 divides a predetermined N-dimensional space into grids.
The searching unit 102 searches for a square having the minimum value of the objective function from squares in a predetermined N-dimensional space.
The determination unit 103 determines the values of the N parameters corresponding to the mass with the minimum value of the objective function as the optimum solution.

以下では、探索部102について詳細に説明する。
まず、探索部102は、所定のN次元空間上のマスの1つを探索開始点とする。探索開始点のマスは、例えば、所定のN次元空間上のマスのうち、目的関数の値が最大値となると推測されるマスとする。
The search unit 102 will be described in detail below.
First, the search unit 102 sets one of the squares on a predetermined N-dimensional space as a search start point. The search starting point is, for example, the square in which the value of the objective function is assumed to be the maximum among the squares in the predetermined N-dimensional space.

続いて、探索部102は、探索開始点のマス及びその周囲のマスの目的関数の値を計算し、その周囲のマスのうち目的関数の値が最小である周囲のマスに移動し、移動先のマスを新たな現マスとする。そして、探索部102は、現マス及びその周囲のマスの目的関数の値を計算し、現マスの周囲のマスの中から、現マスよりも目的関数の値が小さいマスを探索する探索処理を再帰的に実行する。 Next, the search unit 102 calculates the objective function values of the search start point square and the squares around it, moves to the squares around it that have the smallest objective function value, and moves to the surrounding squares. is the new current cell. Then, the search unit 102 calculates the objective function values of the current square and its surrounding squares, and performs a search process of searching for a square having a smaller objective function value than the current square from among the squares surrounding the current square. Run recursively.

探索部102は、探索処理で探索対象であった周囲のマスを見つけられなかった場合、現マスが所定のN次元空間上のマスの中で目的関数の値が最小値となるマスであると判断し、探索を終了する。 If the searching unit 102 fails to find the surrounding squares to be searched in the search processing, the search unit 102 determines that the current square is the square having the minimum value of the objective function among the squares in the predetermined N-dimensional space. Make a decision and end the search.

ここで、探索処理は、マスの探索方向が互いに異なるN個の探索モードを含み、探索部102は、N個の探索モードを切り替えながら、探索処理を実行する。Nが2の場合の探索処理は、上述のように、上下左右探索モードと斜め探索モードとの2個の探索モードを含むことになる。 Here, the search processing includes N search modes in which the square search directions are different from each other, and the search unit 102 performs the search processing while switching between the N search modes. The search process when N is 2 includes two search modes, the up/down/left/right search mode and the oblique search mode, as described above.

また、目的関数は、N個のパラメータのそれぞれの値を決定する複数の指標のそれぞれの値を、指標毎に計算された除数で除算し、かつ、指標毎に設定された重み係数で乗算した上で、加算した関数である。探索部102は、該当マスに対応するN個のパラメータのそれぞれの値となるように複数の指標のそれぞれの値を計算し、計算した複数の指標のそれぞれの値を目的関数に代入することで、該当マスの目的関数の値を計算する。 Also, the objective function is obtained by dividing each value of a plurality of indices that determine the values of each of the N parameters by a divisor calculated for each index and multiplying by a weighting factor set for each index. Above is the added function. The search unit 102 calculates the values of each of the plurality of indices so as to be the values of the N parameters corresponding to the square, and substitutes the calculated values of the plurality of indices into the objective function. , to calculate the value of the objective function of the corresponding cell.

なお、探索部102は、目的関数の指標毎の除数を、探索開始点のマスにおけるその指標の値を基に計算しても良い。又は、探索部102は、目的関数の指標毎の除数を、マスを移動する度に計算しても良い。マスを移動する度に計算する場合、探索部102は、移動先のマスで周囲を探索するのに使用する目的関数の指標毎の除数を、移動先のマスにおけるその指標の値とすれば良い。 Note that the search unit 102 may calculate the divisor for each index of the objective function based on the value of the index in the search start point cell. Alternatively, the search unit 102 may calculate the divisor for each index of the objective function each time the square is moved. When calculating each time a square is moved, the search unit 102 may set the divisor for each index of the objective function used to search the surroundings of the destination square as the value of the index in the destination square. .

以上、実施の形態を参照して本開示を説明したが、本開示は上記の実施の形態に限定されるものではない。本開示の構成や詳細には、本開示のスコープ内で当業者が理解し得る様々な変更をすることができる。 Although the present disclosure has been described above with reference to the embodiments, the present disclosure is not limited to the above embodiments. Various changes that can be understood by those skilled in the art can be made to the configuration and details of the present disclosure within the scope of the present disclosure.

例えば、上記の実施の形態では、本開示のパラメータ探索装置をハードウェアの構成として説明したが、本開示は、これに限定されるものではない。本開示は、パラメータ探索装置が、プロセッサ及びメモリを内蔵し、プロセッサが、メモリに格納されたコンピュータプログラムを読み出し実行することにより、パラメータ探索装置の任意の処理を、実現することも可能である。 For example, in the above embodiment, the parameter search device of the present disclosure has been described as a hardware configuration, but the present disclosure is not limited to this. According to the present disclosure, the parameter search device incorporates a processor and memory, and the processor reads and executes a computer program stored in the memory, thereby realizing arbitrary processing of the parameter search device.

上記の例において、プログラムは、様々なタイプの非一時的なコンピュータ可読媒体(non-transitory computer readable medium)を用いて格納され、コンピュータに供給することができる。非一時的なコンピュータ可読媒体は、様々なタイプの実体のある記録媒体(tangible storage medium)を含む。非一時的なコンピュータ可読媒体の例は、磁気記録媒体(例えばフレキシブルディスク、磁気テープ、ハードディスクドライブ)、光磁気記録媒体(例えば光磁気ディスク)、CD-ROM(Compact Disc-Read Only Memory)、CD-R(CD-Recordable)、CD-R/W(CD-ReWritable)、半導体メモリ(例えば、マスクROM、PROM(Programmable ROM)、EPROM(Erasable PROM)、フラッシュROM、RAM(Random Access Memory))を含む。また、プログラムは、様々なタイプの一時的なコンピュータ可読媒体(transitory computer readable medium)によってコンピュータに供給されても良い。一時的なコンピュータ可読媒体の例は、電気信号、光信号、及び電磁波を含む。一時的なコンピュータ可読媒体は、電線及び光ファイバ等の有線通信路、又は無線通信路を介して、プログラムをコンピュータに供給できる。 In the above examples, the programs can be stored and delivered to the computer using various types of non-transitory computer readable media. Non-transitory computer-readable media include various types of tangible storage media. Examples of non-transitory computer-readable media include magnetic recording media (eg, flexible discs, magnetic tapes, hard disk drives), magneto-optical recording media (eg, magneto-optical discs), CD-ROMs (Compact Disc-Read Only Memory), CDs - R (CD-Recordable), CD-R/W (CD-ReWritable), semiconductor memory (e.g. mask ROM, PROM (Programmable ROM), EPROM (Erasable PROM), flash ROM, RAM (Random Access Memory)) include. The program may also be delivered to the computer by various types of transitory computer readable medium. Examples of transitory computer-readable media include electrical signals, optical signals, and electromagnetic waves. Transitory computer-readable media can deliver the program to the computer via wired channels, such as wires and optical fibers, or wireless channels.

上記の実施の形態の一部又は全部は、以下の付記のようにも記載されうるが、以下には限られない。
(付記1)
N(Nは2以上の自然数)個のパラメータを軸とする所定のN次元空間をマス目状に分割する分割部と、
前記所定のN次元空間上のマスの中から、当該マスに対応するN個のパラメータのそれぞれの値を用いて得られる目的関数の値が最小値となるマスを探索する探索部と、
前記目的関数の値が最小値となるマスに対応するN個のパラメータのそれぞれの値を最適解に決定する決定部と、
を備えるパラメータ探索装置。
(付記2)
前記探索部は、
前記所定のN次元空間上のマスの1つを開始点とし、
開始点とする現マス及びその周囲のマスの前記目的関数の値を計算し、開始点のマスの周囲のマスのうち前記目的関数の値が最小である周囲のマスに移動し、移動先のマスを新たな現マスとする探索処理を再帰的に実行し、
前記探索処理で探索対象であった周囲のマスを見つけられなかった場合、現マスが前記所定のN次元空間上のマスの中で前記目的関数の値が最小値となるマスであると判断する、
付記1に記載のパラメータ探索装置。
(付記3)
前記開始点のマスは、前記所定のN次元空間上のマスのうち、前記目的関数の値が最大値となると推測されるマスである、
付記2に記載のパラメータ探索装置。
(付記4)
前記探索処理は、
マスの探索方向が互いに異なるN個の探索モードを含み、
前記探索部は、
N個の探索モードを切り替えながら、前記探索処理を実行する、
付記2又は3に記載のパラメータ探索装置。
(付記5)
Nは2であり、
前記探索処理は、
現マス及びその周囲の第1探索方向のマスの前記目的関数の値を計算し、現マスの周囲の前記第1探索方向のマスの中から、現マスよりも前記目的関数の値が小さいマスを探索する第1探索モードと、
現マス及びその周囲の第2探索方向のマスの前記目的関数の値を計算し、現マスの周囲の前記第2探索方向のマスの中から、現マスよりも前記目的関数の値が小さいマスを探索する第2探索モードと、を含み、
前記探索部は、
前記開始点とする現マスから移動した移動先のマスの方向が前記第1探索方向である場合、前記第1探索モードを実行し、前記第1探索モードで探索対象であった周囲のマスを見つけた場合、前記第1探索モードで見つけた周囲のマスのうち前記目的関数の値が最小である周囲のマスに移動し、移動先のマスを新たな現マスとして、前記第1探索モードを再帰的に実行し、
前記開始点とする現マスから移動した移動先のマスの方向が前記第2探索方向である場合、前記第2探索モードを実行し、前記第2探索モードで探索対象であった周囲のマスを見つけた場合、前記第2探索モードで見つけた周囲のマスのうち前記目的関数の値が最小である周囲のマスに移動し、移動先のマスを新たな現マスとして、前記第2探索モードを再帰的に実行する、
付記4に記載のパラメータ探索装置。
(付記6)
前記探索部は、
前記第1探索モードで探索対象であった周囲のマスを見つけられなかった場合、続いて前記第2探索モードを実行し、前記第2探索モードでも探索対象であった周囲のマスを見つけられなかった場合、現マスが前記所定の2次元空間上のマスの中で前記目的関数の値が最小値となるマスであると判断し、前記第2探索モードで探索対象であった周囲のマスを見つけた場合、前記第2探索モードで見つけた周囲のマスのうち前記目的関数の値が最小である周囲のマスに移動し、移動先のマスを新たな現マスとして、前記第2探索モードを再帰的に実行し、
前記第2探索モードで探索対象であった周囲のマスを見つけられなかった場合、続いて前記第1探索モードを実行し、前記第1探索モードでも探索対象であった周囲のマスを見つけられなかった場合、現マスが前記所定の2次元空間上のマスの中で前記目的関数の値が最小値となるマスであると判断し、前記第1探索モードで探索対象であった周囲のマスを見つけた場合、前記第1探索モードで見つけた周囲のマスのうち前記目的関数の値が最小である周囲のマスに移動し、移動先のマスを新たな現マスとして、前記第1探索モードを再帰的に実行する、
付記5に記載のパラメータ探索装置。
(付記7)
前記目的関数は、N個のパラメータのそれぞれの値を決定する複数の指標のそれぞれの値を、前記指標毎に計算された除数で除算し、かつ、前記指標毎に設定された重み係数で乗算した上で、加算した関数であり、
前記探索部は、
前記所定のN次元空間上のマスの前記目的関数の値を計算する場合、当該マスに対応するN個のパラメータのそれぞれの値となるように前記複数の指標のそれぞれの値を計算し、計算した前記複数の指標のそれぞれの値を前記目的関数に代入する、
付記2から4のいずれか1項に記載のパラメータ探索装置。
(付記8)
前記探索部は、
前記指標毎の前記除数を、前記開始点のマスにおける当該指標の値を基に計算する、
付記7に記載のパラメータ探索装置。
(付記9)
前記探索部は、
前記指標毎の前記除数を、マスを移動する度に計算し、
移動先のマスで周囲を探索するのに使用する前記目的関数の前記指標毎の前記除数を、移動先のマスにおける当該指標の値とする、
付記7に記載のパラメータ探索装置。
(付記10)
Nは2であり、
前記パラメータは、商品の最大在庫日数及び安全在庫日数であり、
前記指標は、商品の平均在庫日数、平均在庫金額、発注回数、及び欠品率のうち少なくとも1個を含む、
付記7から9のいずれか1項に記載のパラメータ探索装置。
(付記11)
パラメータ探索装置によるパラメータ探索方法であって、
N(Nは2以上の自然数)個のパラメータを軸とする所定のN次元空間をマス目状に分割するステップと、
前記所定のN次元空間上のマスの中から、当該マスに対応するN個のパラメータのそれぞれの値を用いて得られる目的関数の値が最小値となるマスを探索するステップと、
前記目的関数の値が最小値となるマスに対応するN個のパラメータのそれぞれの値を最適解に決定するステップと、
を含むパラメータ探索方法。
(付記12)
コンピュータを、
N(Nは2以上の自然数)個のパラメータを軸とする所定のN次元空間をマス目状に分割する分割手段と、
前記所定のN次元空間上のマスの中から、当該マスに対応するN個のパラメータのそれぞれの値を用いて得られる目的関数の値が最小値となるマスを探索する探索手段と、
前記目的関数の値が最小値となるマスに対応するN個のパラメータのそれぞれの値を最適解に決定する決定手段と、
として機能させるためのプログラム。
Some or all of the above embodiments can also be described in the following additional remarks, but are not limited to the following.
(Appendix 1)
a dividing unit that divides a predetermined N-dimensional space with N (N is a natural number equal to or greater than 2) parameters as axes;
a searching unit for searching, from among the squares in the predetermined N-dimensional space, a square having a minimum value of an objective function obtained using the respective values of N parameters corresponding to the square;
a determination unit that determines the values of the N parameters corresponding to the mass at which the value of the objective function is the minimum value as an optimal solution;
A parameter search device comprising:
(Appendix 2)
The search unit is
With one of the squares on the predetermined N-dimensional space as a starting point,
Calculate the value of the objective function for the current square serving as the starting point and surrounding squares, move to the surrounding square with the smallest objective function value among the squares around the starting square, and recursively execute search processing with the square as the new current square,
If the search process fails to find the surrounding squares to be searched, it is determined that the current square is the square having the minimum value of the objective function among the squares in the predetermined N-dimensional space. ,
A parameter search device according to appendix 1.
(Appendix 3)
The square of the starting point is, among the squares in the predetermined N-dimensional space, the square that is estimated to have the maximum value of the objective function.
The parameter search device according to appendix 2.
(Appendix 4)
The search process includes
including N search modes in which the search directions of the squares are different from each other;
The search unit is
performing the search process while switching between N search modes;
A parameter search device according to appendix 2 or 3.
(Appendix 5)
N is 2;
The search process includes
calculating the value of the objective function of the current square and its surrounding squares in the first search direction; a first search mode for searching for
calculating the value of the objective function of the current square and its surrounding squares in the second search direction; a second search mode that searches for
The search unit is
If the direction of the destination square to which the current square has been moved from the starting point is the first search direction, the first search mode is executed, and the surrounding squares that were searched in the first search mode are searched. If it is found, it moves to the surrounding square having the smallest value of the objective function among the surrounding squares found in the first search mode, sets the destination square as a new current square, and resumes the first search mode. run recursively,
If the direction of the destination square to which the current square as the starting point has been moved is the second search direction, the second search mode is executed, and the surrounding squares that were searched in the second search mode are searched. If it is found, it moves to the surrounding square having the smallest value of the objective function among the surrounding squares found in the second search mode, sets the destination square as a new current square, and resumes the second search mode. run recursively,
The parameter search device according to appendix 4.
(Appendix 6)
The search unit is
When the search target surrounding squares cannot be found in the first search mode, the second search mode is subsequently executed, and the search target surrounding squares cannot be found even in the second search mode. In this case, the current square is determined to be the square having the minimum value of the objective function among the squares in the predetermined two-dimensional space, and the surrounding squares that were searched in the second search mode are searched. If it is found, it moves to the surrounding square having the smallest value of the objective function among the surrounding squares found in the second search mode, sets the destination square as a new current square, and resumes the second search mode. run recursively,
If the search target surrounding squares cannot be found in the second search mode, then the first search mode is executed, and the search target surrounding squares cannot be found even in the first search mode. In this case, the current square is determined to be the square having the minimum value of the objective function among the squares in the predetermined two-dimensional space, and the surrounding squares that were searched in the first search mode are searched. If it is found, it moves to the surrounding square having the smallest value of the objective function among the surrounding squares found in the first search mode, sets the destination square as a new current square, and resumes the first search mode. run recursively,
A parameter search device according to appendix 5.
(Appendix 7)
The objective function divides each value of a plurality of indices that determine the values of each of the N parameters by a divisor calculated for each of the indices, and multiplies by a weighting factor set for each of the indices. and then add the function,
The search unit is
When calculating the value of the objective function of the mass on the predetermined N-dimensional space, calculating the values of each of the plurality of indices so as to be the respective values of N parameters corresponding to the mass, and calculating substituting the value of each of the plurality of indices obtained into the objective function;
5. The parameter search device according to any one of appendices 2 to 4.
(Appendix 8)
The search unit is
calculating the divisor for each index based on the value of the index in the starting point square;
A parameter search device according to appendix 7.
(Appendix 9)
The search unit is
calculating the divisor for each index each time a square is moved;
Let the divisor for each index of the objective function used to search the surroundings in the destination square be the value of the index in the destination square;
A parameter search device according to appendix 7.
(Appendix 10)
N is 2;
The parameters are the maximum inventory days and safety inventory days of the product,
The indicator includes at least one of the average inventory days, average inventory value, number of orders, and out-of-stock rate of the product,
The parameter search device according to any one of appendices 7 to 9.
(Appendix 11)
A parameter search method by a parameter search device,
a step of dividing a predetermined N-dimensional space with N (N is a natural number equal to or greater than 2) parameters as axes;
a step of searching from among the squares in the predetermined N-dimensional space for a square having the minimum value of an objective function obtained using the respective values of N parameters corresponding to the square;
a step of determining the values of each of the N parameters corresponding to the mass at which the value of the objective function is the minimum value as an optimal solution;
Parameter search method, including
(Appendix 12)
the computer,
dividing means for dividing a predetermined N-dimensional space with N (N is a natural number equal to or greater than 2) parameters as axes;
search means for searching, from among the squares in the predetermined N-dimensional space, a square having a minimum value of an objective function obtained using respective values of N parameters corresponding to the square;
Determination means for determining the values of each of the N parameters corresponding to the mass at which the value of the objective function is the minimum value as an optimal solution;
A program to function as

10 パラメータ探索装置
11 分割部
12 探索部
13 決定部
14 入力部
100 パラメータ探索装置
101 分割部
102 探索部
103 決定部
REFERENCE SIGNS LIST 10 parameter search device 11 division unit 12 search unit 13 determination unit 14 input unit 100 parameter search device 101 division unit 102 search unit 103 determination unit

Claims (6)

N(Nは2以上の自然数)個のパラメータを軸とする所定のN次元空間をマス目状に分割する分割部と、
前記所定のN次元空間上のマスの中から、当該マスに対応するN個のパラメータのそれぞれの値を用いて得られる目的関数の値が最小値となるマスを探索する探索部と、
前記目的関数の値が最小値となるマスに対応するN個のパラメータのそれぞれの値を最適解に決定する決定部と、
を備え
前記探索部は、
前記所定のN次元空間上のマスの1つを開始点とし、
開始点とする現マス及びその周囲のマスの前記目的関数の値を計算し、開始点のマスの周囲のマスのうち前記目的関数の値が最小である周囲のマスに移動し、移動先のマスを新たな現マスとする探索処理を再帰的に実行し、
前記探索処理で探索対象であった周囲のマスを見つけられなかった場合、現マスが前記所定のN次元空間上のマスの中で前記目的関数の値が最小値となるマスであると判断し、
前記目的関数は、N個のパラメータのそれぞれの値を決定する複数の指標のそれぞれの値を、前記指標毎に計算された除数で除算し、かつ、前記指標毎に設定された重み係数で乗算した上で、加算した関数であり、
前記探索部は、
前記所定のN次元空間上のマスの前記目的関数の値を計算する場合、当該マスに対応するN個のパラメータのそれぞれの値となるように前記複数の指標のそれぞれの値を計算し、計算した前記複数の指標のそれぞれの値を前記目的関数に代入する、
パラメータ探索装置。
a dividing unit that divides a predetermined N-dimensional space with N (N is a natural number equal to or greater than 2) parameters as axes;
a searching unit for searching, from among the squares in the predetermined N-dimensional space, a square having a minimum value of an objective function obtained using the respective values of N parameters corresponding to the square;
a determination unit that determines the values of the N parameters corresponding to the mass at which the value of the objective function is the minimum value as an optimal solution;
with
The search unit is
With one of the squares on the predetermined N-dimensional space as a starting point,
Calculate the value of the objective function for the current square serving as the starting point and surrounding squares, move to the surrounding square with the smallest objective function value among the squares around the starting square, and recursively execute search processing with the square as the new current square,
If the search processing fails to find the surrounding squares to be searched, it is determined that the current square is the square having the minimum value of the objective function among the squares in the predetermined N-dimensional space. ,
The objective function divides each value of a plurality of indices that determine the values of each of the N parameters by a divisor calculated for each of the indices, and multiplies by a weighting factor set for each of the indices. and then add the function,
The search unit is
When calculating the value of the objective function of the mass on the predetermined N-dimensional space, calculating the values of each of the plurality of indices so as to be the respective values of N parameters corresponding to the mass, and calculating substituting the value of each of the plurality of indices obtained into the objective function;
Parameter searcher.
前記探索処理は、
マスの探索方向が互いに異なるN個の探索モードを含み、
前記探索部は、
N個の探索モードを切り替えながら、前記探索処理を実行する、
請求項に記載のパラメータ探索装置。
The search process includes
including N search modes in which the search directions of the squares are different from each other;
The search unit is
performing the search process while switching between N search modes;
The parameter searching device according to claim 1 .
前記探索部は、
前記指標毎の前記除数を、前記開始点のマスにおける当該指標の値を基に計算する、
請求項に記載のパラメータ探索装置。
The search unit is
calculating the divisor for each index based on the value of the index in the starting point square;
The parameter searching device according to claim 1 .
前記探索部は、
前記指標毎の前記除数を、マスを移動する度に計算し、
移動先のマスで周囲を探索するのに使用する前記目的関数の前記指標毎の前記除数を、移動先のマスにおける当該指標の値とする、
請求項に記載のパラメータ探索装置。
The search unit is
calculating the divisor for each index each time a square is moved;
Let the divisor for each index of the objective function used to search the surroundings in the destination square be the value of the index in the destination square;
The parameter searching device according to claim 1 .
パラメータ探索装置によるパラメータ探索方法であって、
N(Nは2以上の自然数)個のパラメータを軸とする所定のN次元空間をマス目状に分割するステップと、
前記所定のN次元空間上のマスの中から、当該マスに対応するN個のパラメータのそれぞれの値を用いて得られる目的関数の値が最小値となるマスを探索する探索ステップと、
前記目的関数の値が最小値となるマスに対応するN個のパラメータのそれぞれの値を最適解に決定するステップと、
を含み、
前記探索ステップでは、
前記所定のN次元空間上のマスの1つを開始点とし、
開始点とする現マス及びその周囲のマスの前記目的関数の値を計算し、開始点のマスの周囲のマスのうち前記目的関数の値が最小である周囲のマスに移動し、移動先のマスを新たな現マスとする探索処理を再帰的に実行し、
前記探索処理で探索対象であった周囲のマスを見つけられなかった場合、現マスが前記所定のN次元空間上のマスの中で前記目的関数の値が最小値となるマスであると判断し、
前記目的関数は、N個のパラメータのそれぞれの値を決定する複数の指標のそれぞれの値を、前記指標毎に計算された除数で除算し、かつ、前記指標毎に設定された重み係数で乗算した上で、加算した関数であり、
前記探索ステップでは、
前記所定のN次元空間上のマスの前記目的関数の値を計算する場合、当該マスに対応するN個のパラメータのそれぞれの値となるように前記複数の指標のそれぞれの値を計算し、計算した前記複数の指標のそれぞれの値を前記目的関数に代入する、
パラメータ探索方法。
A parameter search method by a parameter search device,
a step of dividing a predetermined N-dimensional space with N (N is a natural number equal to or greater than 2) parameters as axes;
a searching step of searching, from among the squares in the predetermined N-dimensional space, the square having the minimum value of the objective function obtained using the respective values of the N parameters corresponding to the square;
a step of determining the values of each of the N parameters corresponding to the mass at which the value of the objective function is the minimum value as an optimal solution;
including
In the search step,
With one of the squares on the predetermined N-dimensional space as a starting point,
Calculate the value of the objective function for the current square serving as the starting point and surrounding squares, move to the surrounding square with the smallest objective function value among the squares around the starting square, and recursively execute search processing with the square as the new current square,
If the search processing fails to find the surrounding squares to be searched, it is determined that the current square is the square having the minimum value of the objective function among the squares in the predetermined N-dimensional space. ,
The objective function divides each value of a plurality of indices that determine the values of each of the N parameters by a divisor calculated for each of the indices, and multiplies by a weighting factor set for each of the indices. and then add the function,
In the search step,
When calculating the value of the objective function of the mass on the predetermined N-dimensional space, calculating the values of each of the plurality of indices so as to be the respective values of N parameters corresponding to the mass, and calculating substituting the value of each of the plurality of indices obtained into the objective function;
Parameter search method.
コンピュータを、
N(Nは2以上の自然数)個のパラメータを軸とする所定のN次元空間をマス目状に分割する分割手段と、
前記所定のN次元空間上のマスの中から、当該マスに対応するN個のパラメータのそれぞれの値を用いて得られる目的関数の値が最小値となるマスを探索する探索手段と、
前記目的関数の値が最小値となるマスに対応するN個のパラメータのそれぞれの値を最適解に決定する決定手段と、
として機能させるためのプログラムであって、
前記探索手段は、
前記所定のN次元空間上のマスの1つを開始点とし、
開始点とする現マス及びその周囲のマスの前記目的関数の値を計算し、開始点のマスの周囲のマスのうち前記目的関数の値が最小である周囲のマスに移動し、移動先のマスを新たな現マスとする探索処理を再帰的に実行し、
前記探索処理で探索対象であった周囲のマスを見つけられなかった場合、現マスが前記所定のN次元空間上のマスの中で前記目的関数の値が最小値となるマスであると判断し、
前記目的関数は、N個のパラメータのそれぞれの値を決定する複数の指標のそれぞれの値を、前記指標毎に計算された除数で除算し、かつ、前記指標毎に設定された重み係数で乗算した上で、加算した関数であり、
前記探索手段は、
前記所定のN次元空間上のマスの前記目的関数の値を計算する場合、当該マスに対応するN個のパラメータのそれぞれの値となるように前記複数の指標のそれぞれの値を計算し、計算した前記複数の指標のそれぞれの値を前記目的関数に代入する、
プログラム
the computer,
dividing means for dividing a predetermined N-dimensional space with N (N is a natural number equal to or greater than 2) parameters as axes;
search means for searching, from among the squares in the predetermined N-dimensional space, a square having a minimum value of an objective function obtained using respective values of N parameters corresponding to the square;
Determination means for determining the values of each of the N parameters corresponding to the mass at which the value of the objective function is the minimum value as an optimal solution;
A program for functioning as
The search means is
With one of the squares on the predetermined N-dimensional space as a starting point,
Calculate the value of the objective function for the current square serving as the starting point and surrounding squares, move to the surrounding square with the smallest objective function value among the squares around the starting square, and recursively execute search processing with the square as the new current square,
If the search processing fails to find the surrounding squares to be searched, it is determined that the current square is the square having the minimum value of the objective function among the squares in the predetermined N-dimensional space. ,
The objective function divides each value of a plurality of indices that determine the values of each of the N parameters by a divisor calculated for each of the indices, and multiplies by a weighting factor set for each of the indices. and then add the function,
The search means is
When calculating the value of the objective function of the mass on the predetermined N-dimensional space, calculating the values of each of the plurality of indices so as to be the respective values of N parameters corresponding to the mass, and calculating substituting the value of each of the plurality of indices obtained into the objective function;
program .
JP2018116286A 2018-06-19 2018-06-19 Parameter search device, parameter search method, and program Active JP7156666B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2018116286A JP7156666B2 (en) 2018-06-19 2018-06-19 Parameter search device, parameter search method, and program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2018116286A JP7156666B2 (en) 2018-06-19 2018-06-19 Parameter search device, parameter search method, and program

Publications (2)

Publication Number Publication Date
JP2019219869A JP2019219869A (en) 2019-12-26
JP7156666B2 true JP7156666B2 (en) 2022-10-19

Family

ID=69096694

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018116286A Active JP7156666B2 (en) 2018-06-19 2018-06-19 Parameter search device, parameter search method, and program

Country Status (1)

Country Link
JP (1) JP7156666B2 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109751994B (en) * 2019-01-26 2022-11-29 西安邮电大学 Autonomous homing underwater vehicle geomagnetic bionic navigation method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005527034A (en) 2002-05-22 2005-09-08 ビッグウッド テクノロジー,インク. A dynamic method for solving large-scale discrete and continuous optimization problems.
JP2005304118A (en) 2004-04-07 2005-10-27 Nippon Telegr & Teleph Corp <Ntt> Controller and control method of distributed energy system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005527034A (en) 2002-05-22 2005-09-08 ビッグウッド テクノロジー,インク. A dynamic method for solving large-scale discrete and continuous optimization problems.
JP2005304118A (en) 2004-04-07 2005-10-27 Nippon Telegr & Teleph Corp <Ntt> Controller and control method of distributed energy system

Also Published As

Publication number Publication date
JP2019219869A (en) 2019-12-26

Similar Documents

Publication Publication Date Title
CN108009016B (en) Resource load balancing control method and cluster scheduler
US10839338B2 (en) Order sourcing with asynchronous communication and using optimization for large sourcing networks
US10936983B2 (en) Optimized menu planning
US9104497B2 (en) Method and system for work load balancing
Wang et al. Task scheduling for multi-tier shuttle warehousing systems
TW202036393A (en) Partitioning of deep learning inference with dynamic offloading
US9384053B2 (en) Task allocation optimization system, task allocation optimization method, and non-transitory computer readable medium storing task allocation optimization program
US10839420B2 (en) Constrained large-data markdown optimizations based upon markdown budget
US10397315B2 (en) Information processing apparatus and load distribution control method
CN106202092A (en) The method and system that data process
JP7156666B2 (en) Parameter search device, parameter search method, and program
Javadi et al. A hybrid electromagnetism-like algorithm for dynamic inter/intra-cell layout problem
JP5987720B2 (en) Binary decision graph processing system and method
US11887017B2 (en) Automatically predicting transaction likelihood information and transaction-related temporal information using machine learning techniques
CN110456921A (en) Predict the method and device of user&#39;s keyboard operation behavior
CN113962540A (en) Maintenance task scheduling method and related equipment
Nirmal et al. Issues of K means clustering while migrating to map reduce paradigm with big data: A survey
CN108536759B (en) Sample playback data access method and device
JP6963465B2 (en) Computer system and data processing control method
KR101919698B1 (en) Group search optimization data clustering method and system using silhouette
CN115061811A (en) Resource scheduling method, device, equipment and storage medium
Gkiotsalitis Multi-objective optimization
Malakooti Independent, convergent, and divergent decision behaviour for interactive multiple objectives linear programming
CN111967682A (en) Equipment shelving scheme optimizing method and device and terminal equipment
CN114581223B (en) Distribution task processing method, equipment, distributed computing system and storage medium

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20210512

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20220331

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20220419

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20220614

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20220906

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20220929

R151 Written notification of patent or utility model registration

Ref document number: 7156666

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151