JP6543318B2 - Particle size distribution measuring device - Google Patents
Particle size distribution measuring device Download PDFInfo
- Publication number
- JP6543318B2 JP6543318B2 JP2017214108A JP2017214108A JP6543318B2 JP 6543318 B2 JP6543318 B2 JP 6543318B2 JP 2017214108 A JP2017214108 A JP 2017214108A JP 2017214108 A JP2017214108 A JP 2017214108A JP 6543318 B2 JP6543318 B2 JP 6543318B2
- Authority
- JP
- Japan
- Prior art keywords
- vector
- light
- coefficient matrix
- particle size
- size distribution
- 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
Links
- 239000002245 particle Substances 0.000 title claims description 84
- 238000009826 distribution Methods 0.000 title claims description 51
- 239000011159 matrix material Substances 0.000 claims description 55
- 230000006870 function Effects 0.000 claims description 5
- 230000001678 irradiating effect Effects 0.000 claims description 5
- 238000003860 storage Methods 0.000 description 15
- 238000000034 method Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 230000010287 polarization Effects 0.000 description 2
- 238000003892 spreading Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Description
本発明は、測定対象である粒子に光を照射し、その照射により生じる回折光や散乱光（以下、回折／散乱光とも言う）の光強度分布に基づいて粒子の粒子径分布を算出する粒子径分布測定装置に関するものである。 The present invention irradiates light to a particle to be measured, and calculates the particle size distribution of the particle based on the light intensity distribution of diffracted light or scattered light (hereinafter also referred to as diffracted / scattered light) generated by the irradiation. The present invention relates to a diameter distribution measuring device.
この種の粒子径分布測定装置は、特許文献１に示されるように、回折／散乱光を複数の光検出器で検出し、各光検出器から出力された光強度信号の値により得られる光強度分布ベクトルを用いて、下記の式（１）に基づき粒子径分布を算出するようにしている。 As shown in Patent Document 1, this type of particle size distribution measuring apparatus detects diffracted / scattered light with a plurality of photodetectors, and obtains light obtained by the value of the light intensity signal output from each photodetector. The particle size distribution is calculated based on the following equation (1) using the intensity distribution vector.
ところで、係数行列Ｋは、粒子の屈折率等の物性や粒子径や光検出器の配置位置によって定まるものであるところ、上記の式（１）に基づき粒子径分布を算出するためには、まず係数行列Ｋを算出しておく必要がある。 By the way, the coefficient matrix K is determined by the physical properties such as the refractive index of the particles, the particle diameter, and the arrangement position of the light detector, but in order to calculate the particle diameter distribution based on the above equation (1) It is necessary to calculate the coefficient matrix K.
しかしながら、従来は、係数行列Ｋの各要素を順に、ＭＩＥ散乱理論等に基づいて１つずつ求めているので、係数行列Ｋの算出にはかなりの時間を要する。 However, conventionally, since each element of the coefficient matrix K is sequentially obtained one by one based on the MIE scattering theory or the like, it takes a considerable amount of time to calculate the coefficient matrix K.
そこで本発明は、かかる課題を鑑みてなされたものであって、係数行列Ｋの算出時間を短くすることによって、粒子径分布の算出時間を短縮することをその主たる目的とするものである。 Therefore, the present invention has been made in view of such problems, and its main object is to shorten the calculation time of the particle size distribution by shortening the calculation time of the coefficient matrix K.
すなわち、本発明に係る粒子径分布測定装置は、測定対象である粒子に光を照射する光源と、前記光の照射により生じる回折／散乱光の光強度を検出する複数の光検出器と、前記各光検出器から出力された光強度信号を受信し、ベクトルｓがベクトルｑと係数行列Ｋとの積を含む所定式で表されることに基づいて前記粒子の粒子径分布を算出する演算部とを具備し、前記演算部が、前記係数行列Ｋの要素のうち少なくとも２つの要素を並列計算して算出することを特徴とする粒子径分布測定装置。
ここで、ベクトルｓは、各光検出器から出力された光強度信号の値から得られる回折／散乱光の拡がり角度毎の光強度分布を示すベクトル、ベクトルｑは、測定対象である粒子の粒子径分布を示すベクトル、係数行列Ｋは、ベクトルｑをベクトルｓに変換するための行列である。
なお、所定式としては、例えば上述した式（１）や、この式に例えばノイズ等を表す項を加えた式が挙げられる。
That is, the particle size distribution measuring apparatus according to the present invention comprises: a light source for irradiating light to particles to be measured; a plurality of photodetectors for detecting light intensity of diffracted / scattered light generated by the irradiation of the light; A computing unit that receives the light intensity signal output from each photodetector and calculates the particle size distribution of the particles based on that the vector s is represented by a predetermined equation including the product of the vector q and the coefficient matrix K And the operation unit calculates in parallel at least two of the elements of the coefficient matrix K to calculate.
Here, the vector s is a vector indicating the light intensity distribution for each spread angle of the diffracted / scattered light obtained from the value of the light intensity signal output from each light detector, and the vector q is a particle of the particle to be measured A vector representing a radius distribution, the coefficient matrix K is a matrix for converting the vector q into a vector s.
Examples of the predetermined formula include the above-mentioned formula (1) and a formula obtained by adding a term representing, for example, noise to this formula.
このようなものであれば、係数行列Ｋの要素を並列計算して算出しているので、各要素を１つずつ算出する場合に比べて、係数行列Ｋの算出時間を短くすることができ、ひいては粒子径分布を算出するまでの時間を短縮することができる。 In such a case, since the elements of the coefficient matrix K are calculated in parallel, calculation time of the coefficient matrix K can be shortened as compared with the case where each element is calculated one by one, As a result, the time to calculate the particle size distribution can be shortened.
具体的な実施態様としては、前記演算部が、前記係数行列Ｋの要素のうち１つの行に含まれる２以上の要素又は１つの列に含まれる２以上の要素を並列計算して算出するものが挙げられる。 As a specific embodiment, the operation unit calculates two or more elements included in one row of the elements of the coefficient matrix K or two or more elements included in one column in parallel calculation. Can be mentioned.
また、係数行列Ｋの算出時間をより短くするためには、前記演算部が、前記係数行列Ｋの１つの要素を算出するために用いられる複数のパラメータの値を算出するとともに、これらの値のうち少なくとも１つを記憶して、前記係数行列Ｋの別の要素を算出するために用いることが好ましい。 In order to further shorten the calculation time of the coefficient matrix K, the operation unit calculates values of a plurality of parameters used to calculate one element of the coefficient matrix K, and Preferably, at least one of them is stored and used to calculate another element of the coefficient matrix K.
ここで、係数行列Ｋは、下記の式（２）に表されるように、
この係数行列Ｋのうち、同一行に属する要素を算出する際には、回折／散乱光の拡がり角に依存するパラメータを共通して用いることができ、同一列に属する要素を算出する際には、粒子の屈折率及び粒子径に依存するパラメータを共通して用いることができる。
Here, the coefficient matrix K is expressed by the following equation (2),
When calculating an element belonging to the same row in the coefficient matrix K, a parameter depending on the spread angle of the diffracted / scattered light can be used in common, and when calculating an element belonging to the same column Parameters that depend on the refractive index of the particles and the particle size can be used in common.
そこで、具体的な実施態様としては、前記演算部が、前記係数行列Ｋの要素のうちある行に含まれる要素を算出する際に記憶した複数の前記パラメータの値と、ある列に含まれる要素を算出する際に記憶した複数の前記パラメータの値とを用いて、これらの行及び列が交わる位置の要素を算出するものが挙げられる。 Therefore, as a specific embodiment, values of a plurality of the parameters stored when the operation unit calculates an element included in a certain row of the elements of the coefficient matrix K, and an element included in a certain column And calculating the element at the position where these rows and columns intersect, using the values of a plurality of the parameters stored when calculating.
また、本発明に係る粒子径分布測定装置は、測定対象である粒子に光を照射する光源と、前記光の照射により生じる回折／散乱光の光強度を検出する複数の光検出器と、前記各光検出器から出力された光強度信号を受信し、ベクトルｓがベクトルｑと係数行列Ｋとの積を含む所定式で表されることに基づいて前記粒子の粒子径分布を算出する演算部とを具備し、前記演算部が、前記係数行列Ｋの要素のうち１つの要素を算出するために用いられる前記粒子の粒子径に依存する複数の第１パラメータ及び前記回折／散乱光の拡がり角に依存する複数の第２パラメータのうち少なくとも２つのパラメータを並列計算して算出することを特徴とするものである。 A particle size distribution measuring device according to the present invention comprises a light source for irradiating light to particles to be measured, a plurality of light detectors for detecting light intensity of diffracted / scattered light generated by the light irradiation, and A computing unit that receives the light intensity signal output from each photodetector and calculates the particle size distribution of the particles based on that the vector s is represented by a predetermined equation including the product of the vector q and the coefficient matrix K A plurality of first parameters depending on the particle diameter of the particles used to calculate one of the elements of the coefficient matrix K, and the spread angle of the diffracted / scattered light. And calculating at least two of the plurality of second parameters depending on.
このようなものであれば、係数行列Ｋの要素のうち１つの要素を算出するために用いられる複数の第１パラメータ及び複数の第２パラメータのうち少なくとも２つのパラメータを並列計算しているので、従来に比べて係数行列Ｋを算出する時間を短くすることができ、粒子径分布を算出するまでの時間を短縮することができる。 In such a case, since at least two of the plurality of first parameters and the plurality of second parameters used to calculate one of the elements of the coefficient matrix K are calculated in parallel, The time to calculate the coefficient matrix K can be shortened as compared with the conventional case, and the time to calculate the particle size distribution can be shortened.
係数行列Ｋの算出時間をより短くするためには、前記演算部が、２種類の前記第１パラメータ及び２種類の前記第２パラメータを並列計算して算出し、下記の式（３）に基づき係数行列Ｋの各要素を算出するものが好ましい。
このように構成した本発明によれば、係数行列Ｋの算出時間を従来に比べて大幅に短くすることが可能となり、ひいては測定全体にかかる時間を短縮することができる。 According to the present invention configured as described above, the calculation time of the coefficient matrix K can be significantly shortened as compared with the prior art, and hence the time taken for the entire measurement can be shortened.
以下に本発明に係る粒子径分布測定装置１の一実施形態について図面を参照して説明する。 Hereinafter, an embodiment of a particle size distribution measuring device 1 according to the present invention will be described with reference to the drawings.
本実施形態に係る粒子径分布測定装置１は、粒子に光を照射した際に生じる回折／散乱光の拡がり角に応じる光強度分布が、ＭＩＥ散乱理論から粒子径によって定まることを利用し、前記回折／散乱光を検出することによって粒子径分布を測定するようにしたものである。 The particle size distribution measuring apparatus 1 according to the present embodiment utilizes the fact that the light intensity distribution according to the spread angle of the diffracted / scattered light generated when the particles are irradiated with light is determined by the particle size from the MIE scattering theory. The particle size distribution is measured by detecting the diffracted / scattered light.
粒子径分布測定装置１は、図１に模式的に示すように、装置本体２と演算部３とを備えている。 The particle size distribution measuring apparatus 1 includes an apparatus main body 2 and a calculating unit 3 as schematically shown in FIG.
装置本体２は、粒子を分散させたサンプルを収容するセル２１と、そのセル２１内の粒子にレンズ２２を介してレーザ光を照射する光源２３たるレーザ装置と、レーザ光の照射により生じる回折／散乱光の光強度を拡がり角に応じて検出する複数の光検出器２４（Ａ）、２４（Ｂ）とを備えたものである。
なお、セル２１は、本実施形態では、バッチ式セルを用いているが、循環式セルを用いても構わない。
The device body 2 includes a cell 21 for storing a sample in which particles are dispersed, a laser device as a light source 23 for irradiating the particles in the cell 21 with a laser beam through a lens 22, and diffraction / A plurality of light detectors 24 (A) and 24 (B) for detecting the light intensity of the scattered light according to the spread angle is provided.
In addition, although the cell 21 uses a batch type cell in the present embodiment, a circulation type cell may be used.
演算部３は、物理的に言えば、ＣＰＵ、メモリ、入出力インターフェース等を備えた汎用乃至専用のコンピュータであり、各光検出器２４（Ａ）、２４（Ｂ）から出力された光強度信号を受信して式（１）に基づき粒子径分布を算出するものである。 The arithmetic unit 3 is physically or generally a general-purpose or dedicated computer provided with a CPU, a memory, an input / output interface and the like, and the light intensity signal output from each of the light detectors 24 (A) and 24 (B) To calculate the particle size distribution based on the equation (1).
この演算部３は、前記メモリの所定領域に記憶させた所定プログラムにしたがって、ＣＰＵや周辺機器を協働させることにより、図２に示すように、係数行列Ｋを算出する係数行列算出部３１、各光強度信号の値から求められる光強度分布ベクトルｓと前記係数行列Ｋとを用いて式（１）に基づき粒子径分布を算出する粒子径分布算出部３２としての機能を少なくとも発揮するようにしたものである。 The calculation unit 3 calculates a coefficient matrix K as shown in FIG. 2 by causing the CPU and peripheral devices to cooperate in accordance with a predetermined program stored in a predetermined area of the memory, In order to exhibit at least the function as the particle size distribution calculating unit 32 that calculates the particle size distribution based on the equation (1) using the light intensity distribution vector s obtained from the values of the respective light intensity signals and the coefficient matrix K It is
本実施形態では、係数行列算出部３１が特徴的であるので、以下に詳述する。 In the present embodiment, since the coefficient matrix calculation unit 31 is characteristic, it will be described in detail below.
係数行列算出部３１は、測定対象である粒子の物性、粒子径、及び光検出器２４の配置位置に応じた係数行列Ｋを算出するものであり、本実施形態では、オペレータが入力した粒子の例えば屈折率の値を受け付けて、粒子径及び光検出器２４の配置位置に応じた係数行列Ｋの各要素を算出するように構成されている。 The coefficient matrix calculation unit 31 calculates the coefficient matrix K according to the physical properties of the particles to be measured, the particle diameter, and the arrangement position of the light detector 24. In the present embodiment, the particles input by the operator For example, the value of the refractive index is received, and each element of the coefficient matrix K according to the particle diameter and the arrangement position of the light detector 24 is calculated.
より詳細には、この係数行列算出部３１は、図３に示すように、式（３）に基づいて要素を算出する複数の要素算出部４と、前記メモリの所定領域に設定され、前記各要素算出部４の計算過程で算出される２種類の第１パラメータａ、ｂ及び２種類の第２パラメータπ、τの値を記憶する記憶部５としての機能を備えている。
なお、記憶部５は、本実施形態では、キャッシュメモリの所定領域に設定されたものである。
More specifically, as shown in FIG. 3, the coefficient matrix calculation unit 31 is set to a plurality of element calculation units 4 that calculate elements based on the equation (3), and in each of the predetermined areas of the memory. It has a function as a storage unit 5 that stores values of two types of first parameters a and b and two types of second parameters π and τ calculated in the calculation process of the element calculation unit 4.
In the present embodiment, the storage unit 5 is set in a predetermined area of the cache memory.
各要素算出部４は、図３に示すように、式（３）に基づく演算を実行する演算実行部４１と、演算実行部４１が演算に用いるパラメータを算出するパラメータ算出部４２としての機能を具備している。 As shown in FIG. 3, each element calculation unit 4 has functions as an operation execution unit 41 that executes an operation based on Equation (3) and a parameter calculation unit 42 that calculates a parameter that the operation execution unit 41 uses for operation. It is equipped.
しかして、本実施形態では、係数行列算出部３１が、上述した要素算出部４を係数行列Ｋの行の数（ｘ）と同数有しており、図４に示すように、これらの要素算出部４が、１つの列に含まれる要素を並列して計算し、この並列計算を第１列から第ｙ列まで順次進めていくことで全ての要素を算出するように構成されている。
なお、ここで言う並列計算とは、計算を開始するタイミングや終了するタイミング等が必ずしも一致している必要はなく、各要素算出部４が、多少のタイミングのずれを有して要素を算出している計算状態も含む。
また、これらの要素算出部４は、複数のパソコンにまたがらせて構成されていても良い。
Thus, in the present embodiment, the coefficient matrix calculation unit 31 includes the above-described element calculation unit 4 in the same number as the number (x) of the rows of the coefficient matrix K, and as shown in FIG. The unit 4 is configured to calculate the elements included in one column in parallel and calculate all the elements by sequentially advancing the parallel calculation from the first column to the y-th column.
In addition, the parallel calculation said here does not necessarily need to coincide the timing which starts calculation, the timing which complete | finishes, etc., and each element calculation part 4 has a deviation of some timing, and calculates an element. Calculation state is included.
Further, these element calculation units 4 may be configured to extend over a plurality of personal computers.
続いて、要素算出部４が、例えばｉ行ｊ列に位置する要素ｋｉｊを算出する手順について、図３及び図５を参照しながら、要素算出部４における各部の動作説明を兼ねて詳述する。ただし、ｉは、１≦ｉ≦ｘの整数とし、ｊは、１≦ｊ≦ｙの整数とする。
なお、下記の説明では、オペレータにより係数行列算出部３１に粒子の屈折率ｍが入力されているものとする。
Subsequently, the procedure for calculating, for example, the element k ij located at the i-th row and j-th column by the element calculation unit 4 will be described in detail with the explanation of the operation of each part in the element calculation unit 4 referring to FIGS. Do. Here, i is an integer of 1 ≦ i ≦ x, and j is an integer of 1 ≦ j ≦ y.
In the following description, it is assumed that the refractive index m of the particle is input to the coefficient matrix calculation unit 31 by the operator.
ここで、要素ｋｉｊは、測定対象である粒子径範囲を分割したうちｊ番目の範囲に属する単位粒子量の粒子によって生じる回折／散乱光が、複数の光検出器２４のうちｉ番目の光検出器２４に入射したときに検出される光強度を示している。したがって、この要素ｋｉｊを求めるうえで、式（３）における粒子径に関連する値α及び回折／散乱光の拡がり角θは定数となり、要素算出部４は、ｎのみを変数とした下記の式（３）’に基づき要素ｋｉｊを算出する。
ここで、説明の便宜上、
上述した式（３）’及び（４）において、本実施形態では、１種目の第１パラメータは、複数あり、それらは、それぞれａｊ（１）、ａｊ（２）、・・・、ａｊ（Ｎ）である。２種目の第１パラメータも同様に複数あり、それらは、それぞれｂｊ（１）、ｂｊ（２）、・・・、ｂｊ（Ｎ）である。
上述した式（３）’及び（４）において、本実施形態では、１種目の第２パラメータは、複数あり、それらは、それぞれπｉ（１）、πｉ（２）、・・・、πｉ（Ｎ）である。２種目の第２パラメータも同様に複数あり、それらは、それぞれτｉ（１）、τｉ（２）、・・・、τｉ（Ｎ）である。
In the equations (3) ′ and (4) described above, in the present embodiment, there are a plurality of first parameters of the first type, and they are respectively a j (1), a j (2),. j (N). Similarly, there are a plurality of second parameters of the first parameter, which are b j (1), b j (2),..., B j (N), respectively.
In the above-described formulas (3) ′ and (4), in the present embodiment, there are a plurality of first parameters of the second type, and they are respectively π i (1), π i (2),. i (N). Similarly, there are a plurality of second parameters of the second type, which are respectively τ i (1), τ i (2),..., Τ i (N).
始めに、演算実行部４１は、記憶部５にアクセスして（Ｓ１）、ｊ列又はｉ行に含まれる別の要素を算出する要素算出部４がパラメータａｊ（ｍ）、πｉ（ｍ）、ｂｊ（ｍ）、τｉ（ｍ）を既に算出しており、その値を記憶部５が記憶しているか確認する（Ｓ２１〜２４）。 First, the operation execution unit 41 accesses the storage unit 5 (S1), and the element calculation unit 4 that calculates another element included in the j-th column or i-th row receives the parameters a j (m), π i (m ), B j (m), τ i (m) are already calculated, and it is confirmed whether the storage unit 5 stores the values (S21 to 24).
Ｓ２１〜２４において、前述した別の要素を算出する要素算出部４が既に算出しているパラメータに関しては、演算実行部４１がそのパラメータの値を記憶部５から取得する。 In S21 to 24, with respect to the parameter already calculated by the element calculation unit 4 that calculates the above-described other element, the operation execution unit 41 acquires the value of the parameter from the storage unit 5.
Ｓ２１〜２４において、前述した別の要素を算出する要素算出部４がいずれも算出していないパラメータに関しては、演算実行部４１が、そのパラメータを算出するための算出信号をパラメータ算出部４２へ送信する。 In S21 to 24, with respect to parameters that are not calculated by the element calculation unit 4 that calculates the other elements described above, the calculation execution unit 41 transmits a calculation signal for calculating the parameters to the parameter calculation unit 42. Do.
算出信号を受信したパラメータ算出部４２は、その算出信号に対応するパラメータの値を算出して演算実行部４１へ送信するとともに、その値を記憶部５にも送信する（Ｓ３１〜３４）。 The parameter calculating unit 42 that has received the calculation signal calculates the value of the parameter corresponding to the calculation signal and transmits it to the calculation execution unit 41, and also transmits the value to the storage unit 5 (S31 to 34).
演算実行部４１は、記憶部５及びパラメータ算出部４２から取得した各パラメータの値を用いてｆ（ｍ）を算出し、その値をメモリの所定領域に記憶させる（Ｓ４）。 The operation execution unit 41 calculates f (m) using the values of the parameters acquired from the storage unit 5 and the parameter calculation unit 42, and stores the values in a predetermined area of the memory (S4).
続いて、演算実行部４１は、ｆ（１）からｆ（ｍ）までの値を取得するとともに、それらの値を足し合わせ（Ｓ５）、その総和が一定の値に収束しているかを判断する（Ｓ６）。
なお、収束しているかの判断は、例えば、ｆ（ｍ−１）までの総和とｆ（ｍ）までの総和との比率が一定の値以下となるかによって判断する。
Subsequently, the operation execution unit 41 obtains values from f (1) to f (m), adds those values (S5), and determines whether the sum converges to a constant value. (S6).
It is to be noted that the determination as to convergence is made, for example, based on whether the ratio of the sum up to f (m-1) and the sum up to f (m) is equal to or less than a certain value.
Ｓ６において、前記総和が収束していないと判断された場合は、ｎ＝ｍ＋１として（Ｓ７）、再びＳ１へ戻り、Ｓ６において収束されると判断されるまでステッ１からＳ５までを繰り返す。 If it is determined in S6 that the sum does not converge, n = m + 1 is established (S7), and the process returns to S1 again, and steps 1 to S5 are repeated until it is determined in S6 that convergence is performed.
Ｓ６において、前記総和が収束していると判断された場合は、演算実行部４１は、その収束値を要素ｋｉｊの値として粒子径算出部に出力する（Ｓ８）。 If it is determined in S6 that the sum has converged, the operation execution unit 41 outputs the convergence value to the particle diameter calculation unit as the value of the element k ij (S8).
以上のように構成された本実施形態に係る粒子径分布測定装置１によれば、以下のような効果を得ることができる。 According to the particle size distribution measuring device 1 of the present embodiment configured as described above, the following effects can be obtained.
複数の要素算出部４が、ある列に含まれている全ての要素を並列計算して算出するとともに、ある要素算出部４が、計算過程で算出したパラメータの値を記憶部５に送信し、別の要素算出部４が、このパラメータの値を記憶部５から取得して要素を算出するために用いるので、従来のように１つずつ要素を算出する場合に比べて、係数行列Ｋの算出時間を大幅に短くすることができる。 The plurality of element calculation units 4 calculate in parallel all the elements included in a certain column and calculate them, and the certain element calculation unit 4 transmits the values of the parameters calculated in the calculation process to the storage unit 5, Since another element calculation unit 4 obtains the value of this parameter from the storage unit 5 and uses it to calculate an element, calculation of the coefficient matrix K is performed as compared to the case where elements are calculated one by one as in the conventional case. The time can be significantly shortened.
以下にこの効果を詳述する。 This effect is described in detail below.
複数の要素算出部４が、例えば第１列の要素を並列計算する場合、図６に示すように、２種類の第１パラメータａ１（ｎ）及びｂ１（ｎ）は、各要素を算出するために共通して用いられるパラメータである。したがって、複数の要素算出部４のうちいずれか１つが第１パラメータａ１（ｎ）又はｂ１（ｎ）を算出すれば、別の要素算出部４は、記憶部５からこの値を取得して要素の算出に用いることができる。 When the plurality of element calculation units 4 calculate, for example, the elements of the first column in parallel, as shown in FIG. 6, the two types of first parameters a 1 (n) and b 1 (n) calculate each element Parameters commonly used to Therefore, if any one of the plurality of element calculation units 4 calculates the first parameter a 1 (n) or b 1 (n), another element calculation unit 4 acquires this value from the storage unit 5. Can be used to calculate the
次に、複数の要素算出部４が、例えば第２列の要素を並列計算する場合、図６に示すように、各要素を算出するために用いられる２種類の第２パラメータπｉ（ｎ）及びτｉ（ｎ）は、第１列に含まれる要素を算出する過程で既に算出されている。したがって、各要素算出部４は、既に算出されている２種類の第２パラメータπｉ（ｎ）及びτｉ（ｎ）の値を記憶部５から取得し、各要素の算出に用いることができる。
このことは、第３列から第ｙ列の要素を算出する場合においても同様である。
Next, when the plurality of element calculation units 4 calculate the elements of the second column in parallel, for example, as shown in FIG. 6, two types of second parameters π i (n) used to calculate each element And τ i (n) are already calculated in the process of calculating the elements included in the first column. Therefore, each element calculation unit 4 can obtain the values of the two types of second parameters π i (n) and τ i (n) already calculated from the storage unit 5 and can use them for calculation of each element. .
The same applies to the case of calculating the elements of the third to y-th columns.
したがって、前述したように、本実施形態に係る粒子径分布測定装置１によれば、従来に比べて、各要素を算出する時間が短くなり、係数行列Ｋの算出時間を大幅に短くすることができ、粒子径分布の算出時間を短縮することができる。 Therefore, as described above, according to the particle size distribution measuring device 1 according to the present embodiment, the time for calculating each element is shorter than in the conventional case, and the time for calculating the coefficient matrix K is significantly shortened. It is possible to shorten the calculation time of the particle size distribution.
なお、本発明は前記実施形態に限定されるものではない。 The present invention is not limited to the above embodiment.
例えば、前記実施形態では、１つの列に含まれる複数の要素を並列計算するように構成されていたが、図７上段に示すように、１つの行に含まれる複数の要素を並列計算するように構成されていても良い。
また、図７下段に示すように、斜めに位置する複数の要素を並列計算するように構成されていても良い。この場合も、複数の要素を並列計算して算出されたパラメータの値を、別の要素の算出に用いることが可能である。
For example, in the above embodiment, a plurality of elements included in one column are calculated in parallel, but as shown in the upper part of FIG. 7, a plurality of elements included in one row is calculated in parallel. It may be configured.
Further, as shown in the lower part of FIG. 7, a plurality of obliquely positioned elements may be calculated in parallel. Also in this case, it is possible to use the value of the parameter calculated by parallel calculation of a plurality of elements for calculation of another element.
さらに、各列の計算順序は自由に選択することができ、第ｙ列から第１列まで順次計算するように構成しても良く、このように構成することで、後方又は側方への回折／散乱光に対応する要素であって、Ｓ６における収束に比較的長い時間を要する要素を優先して計算することができる。 Furthermore, the calculation order of each column can be freely selected, and may be configured to sequentially calculate from the y-th column to the first column, and by this configuration, backward or sideward diffraction can be performed. The element corresponding to / scattered light, which requires a relatively long time for convergence in S6, can be preferentially calculated.
また、前記実施形態では、パラメータ算出部４２が、ａ、π、ｂ、τを順次算出するように構成されているが、これらのパラメータのうち少なくとも２種類のパラメータを並列計算して算出するように構成されていても良い。さらに、例えば１種目の第１パラメータａｊ（１）及び１種目の第１パラメータａｊ（２）のように、同じ種類のパラメータを並列計算して算出するように構成されていても良い。
このように構成することで、係数行列Ｋの算出時間をより短くすることができる。
Further, in the above embodiment, the parameter calculating unit 42 is configured to sequentially calculate a, π, b, and τ, but at least two of these parameters are calculated in parallel and calculated. It may be configured. Furthermore, for example, the first type of first parameter a j (1) and the first type of first parameter a j (2) may be configured to calculate the parameters of the same type in parallel.
With this configuration, the calculation time of the coefficient matrix K can be further shortened.
さらに、前記実施形態では、パラメータ算出部４２が、算出したパラメータの値を全て記憶部５に送信していたが、予め定めた項数（例えばｎ＝１００）までに算出した値を記憶部５に送信し、それ以降（例えばｎ＝１０１以降）に算出した値は記憶部５に送信しないように構成しても良い。 Furthermore, in the above-described embodiment, although the parameter calculation unit 42 transmits all the calculated parameter values to the storage unit 5, the storage unit 5 may calculate values calculated up to a predetermined number of terms (for example, n = 100). The value calculated after that (for example, after n = 101) may be configured not to be transmitted to the storage unit 5.
粒子径分布算出部３２は、前記実施形態では、式（１）に基づき粒子径分布を算出するものであったが、この式（１）に例えばノイズ等を表す項を加えた式に基づいて粒子径分布を算出するものであっても良い。 The particle size distribution calculating unit 32 calculates the particle size distribution based on the equation (1) in the above embodiment, but based on an equation obtained by adding a term representing noise or the like to the equation (1), for example. The particle size distribution may be calculated.
また、式（１）のベクトルｓは、各光検出器２４から出力された光強度信号の値から得られる回折／散乱光の光強度分布を示すベクトルであったが、ベクトルｓを前記光強度信号の値そのものを要素とするベクトルとしても良い。
この場合、粒子径分布を示すベクトルｑをベクトルｓに変換すべく、係数行列Ｋの各要素が、式（２）から得られる値を、例えば各光検出器２４からセル２１までの距離等を用いてさらに変換した値であれば良い。
Moreover, although vector s of Formula (1) was vector which shows light intensity distribution of the diffracted / scattered light obtained from the value of the light intensity signal output from each light detector 24, vector s is the said light intensity It may be a vector having the signal value itself as an element.
In this case, in order to convert the vector q indicating the particle size distribution into a vector s, each element of the coefficient matrix K has a value obtained from the equation (2), for example, the distance from each photodetector 24 to the cell 21 It may be a value further converted by using.
装置本体２に関して言えば、前記実施形態では、光源２３としてレーザ装置を用いていたが、光源２３として例えば無偏向な光を発するＬＥＤ素子を用いても良い。
この場合、要素算出部４は、光強度の垂直な偏向成分と水平な偏向成分とを表している下記の２つの式（５）及び（６）により得られるｋ１及びｋ２の平均値に基づいて各要素を算出するように構成されていれば良い。
In this case, the element calculation unit 4 calculates the average value of k 1 and k 2 obtained by the following two equations (5) and (6) representing the vertical deflection component and the horizontal deflection component of the light intensity. It may be configured to calculate each element based on it.
その他、本発明は前記各実施形態に限られないし、その各部分構成を組み合わせて
も良く、その趣旨を逸脱しない範囲で種々の変形が可能であるのは言うまでもない。
In addition, the present invention is not limited to the above embodiments, and the respective partial configurations thereof may be combined, and it goes without saying that various modifications can be made without departing from the scope of the invention.
１ ・・・粒子径分布測定装置
２ ・・・装置本体
３ ・・・演算部
３１ ・・・係数行列算出部
４ ・・・要素算出部
４１ ・・・演算実行部
４２ ・・・パラメータ算出部
５ ・・・記憶部
1 ... particle size distribution measuring device 2 ... device main body 3 ... computing unit 31 ... coefficient matrix computing unit 4 ... element computing unit 41 ... computation executing unit 42 ... parameter computing unit 5 ··· Storage unit
Claims (4)
前記演算部が、前記粒子の粒子径に依存する複数の第１パラメータ及び前記回折／散乱光の拡がり角に依存する複数の第２パラメータを用いて前記係数行列Ｋの要素を算出するものであり、
前記演算部が、前記係数行列Ｋの要素のうち少なくとも２つの要素を並列計算して算出することを特徴とする粒子径分布測定装置。
ここで、ベクトルｓは、各光検出器から出力された光強度信号の値から得られる回折／散乱光の拡がり角度毎の光強度分布を示すベクトル、ベクトルｑは、測定対象である粒子の粒子径分布を示すベクトル、係数行列Ｋは、ベクトルｑをベクトルｓに変換するための行列である。 A light source for irradiating light to particles to be measured, a plurality of light detectors for detecting light intensity of diffracted / scattered light generated by the light irradiation, and light intensity signals output from the light detectors are received And an operation unit for calculating the particle size distribution of the particles based on the fact that the vector s is represented by a predetermined equation including the product of the vector q and the coefficient matrix K,
The calculation unit calculates the elements of the coefficient matrix K using a plurality of first parameters depending on the particle diameter of the particles and a plurality of second parameters depending on the spread angle of the diffracted / scattered light. ,
The particle diameter distribution measuring apparatus according to claim 1, wherein the calculation unit calculates at least two of the elements of the coefficient matrix K in parallel.
Here, the vector s is a vector indicating the light intensity distribution for each spread angle of the diffracted / scattered light obtained from the value of the light intensity signal output from each light detector, and the vector q is a particle of the particle to be measured A vector representing a radius distribution, the coefficient matrix K is a matrix for converting the vector q into a vector s.
第１の演算実行部が、互いに異なる複数の行それぞれに含まれる要素又は互いに異なる複数の列それぞれに含まれる要素を算出し、
第２の演算実行部が、前記互いに異なる複数の行又は前記互いに異なる複数の列それぞれに含まれ、前記第１の演算実行部により算出された要素とは別の要素を算出することを特徴とする請求項１記載の粒子径分布測定装置。 The operation unit has a plurality of operation execution units that calculate the same number of elements in parallel and calculate them in parallel;
A first operation execution unit calculates an element included in each of a plurality of different rows or an element included in each of a plurality of different columns;
A second operation execution unit is included in each of the plurality of different rows or each of the plurality of different columns, and calculates an element different from the element calculated by the first operation execution unit. The particle size distribution measuring device according to claim 1.
前記第２の演算実行部が、前記第１の演算実行部により算出された行とは別の１つの行又は前記第１の演算実行部により算出された列とは別の１つの列に含まれるそれぞれの要素の算出することを特徴とする請求項２記載の粒子径分布測定装置。 The first operation execution unit calculates each element included in one row or one column,
The second operation execution unit is included in one row different from the row calculated by the first operation execution unit or in one column different from the column calculated by the first operation execution unit The particle size distribution measuring apparatus according to claim 2, wherein each of the elements to be calculated is calculated.
前記各光検出器から出力された光強度信号を受信し、ベクトルｓがベクトルｑと係数行列Ｋとの積を含む所定式で表されることに基づいて前記粒子の粒子径分布を算出する演算部としての機能をコンピュータに発揮させるものであり、
前記演算部が、前記粒子の粒子径に依存する複数の第１パラメータ及び前記回折／散乱光の拡がり角に依存する複数の第２パラメータを用いて前記係数行列Ｋの要素を算出するように構成されるとともに、
前記演算部が、前記係数行列Ｋの要素のうち少なくとも２つの要素を並列計算して算出するように構成されていることを特徴とするプログラム。
ここで、ベクトルｓは、各光検出器から出力された光強度信号の値から得られる回折／散乱光の拡がり角度毎の光強度分布を示すベクトル、ベクトルｑは、測定対象である粒子の粒子径分布を示すベクトル、係数行列Ｋは、ベクトルｑをベクトルｓに変換するための行列である。 A program mounted on a particle size distribution measuring apparatus, comprising: a light source for irradiating light to particles to be measured; and a plurality of light detectors for detecting light intensity of diffracted / scattered light generated by the light irradiation. ,
Calculation to calculate the particle size distribution of the particles based on the fact that the light intensity signal output from each of the light detectors is received and the vector s is expressed by a predetermined equation including the product of the vector q and the coefficient matrix K To make the computer perform its functions as
The calculation unit is configured to calculate the elements of the coefficient matrix K using a plurality of first parameters depending on the particle diameter of the particles and a plurality of second parameters depending on the spread angle of the diffracted / scattered light. As well as
A program characterized in that the calculation unit calculates in parallel calculation of at least two elements of the elements of the coefficient matrix K.
Here, the vector s is a vector indicating the light intensity distribution for each spread angle of the diffracted / scattered light obtained from the value of the light intensity signal output from each light detector, and the vector q is a particle of the particle to be measured A vector representing a radius distribution, the coefficient matrix K is a matrix for converting the vector q into a vector s.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2017214108A JP6543318B2 (en) | 2017-11-06 | 2017-11-06 | Particle size distribution measuring device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2017214108A JP6543318B2 (en) | 2017-11-06 | 2017-11-06 | Particle size distribution measuring device |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2013132140A Division JP6240416B2 (en) | 2013-06-24 | 2013-06-24 | Particle size distribution measuring device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2018036277A JP2018036277A (en) | 2018-03-08 |
| JP6543318B2 true JP6543318B2 (en) | 2019-07-10 |
Family
ID=61565732
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2017214108A Active JP6543318B2 (en) | 2017-11-06 | 2017-11-06 | Particle size distribution measuring device |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6543318B2 (en) |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0224533A (en) * | 1988-07-13 | 1990-01-26 | Shimadzu Corp | Apparatus for measuring distribution of grain size |
| JPH05324700A (en) * | 1992-05-19 | 1993-12-07 | N T T Data Tsushin Kk | Matrix multiplication device |
| JP2862077B2 (en) * | 1996-08-28 | 1999-02-24 | 株式会社島津製作所 | Particle size distribution analyzer |
| JP3704931B2 (en) * | 1997-12-25 | 2005-10-12 | 株式会社島津製作所 | Laser diffraction / scattering particle size distribution analyzer |
| US20010035954A1 (en) * | 2000-03-10 | 2001-11-01 | Rahn John Richard | Method and apparatus for measuring particle size distributions using light scattering |
| JP2008164539A (en) * | 2006-12-29 | 2008-07-17 | Horiba Ltd | Particle size distribution measuring apparatus |
-
2017
- 2017-11-06 JP JP2017214108A patent/JP6543318B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2018036277A (en) | 2018-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6240416B2 (en) | Particle size distribution measuring device | |
| US9538970B2 (en) | Generating attenuation image data and phase image data in an X-ray system | |
| JP2014526805A5 (en) | ||
| CN106815878A (en) | A kind of image rebuilding method based on the multispectral tomography of various visual angles | |
| JP6218453B2 (en) | Particle size distribution measuring device | |
| Kuramoto et al. | The nature of progression in Parkinson’s disease: an application of non-linear, multivariate, longitudinal random effects modelling | |
| Stephany et al. | A pre-regularization scheme for the reconstruction of a spatial dependent scattering albedo using a hybrid ant colony optimization implementation | |
| US9863764B2 (en) | Storage medium, shape calculation device, and shape measurement method | |
| Sarkar et al. | Perturbing the perturbed: Stability of quasinormal modes in presence of a positive cosmological constant | |
| JP6218449B2 (en) | Particle size distribution measuring device | |
| JPWO2019077955A1 (en) | Spectrum analyzer and spectrum analysis method | |
| JPWO2018092462A1 (en) | Particle size distribution measuring device and program for particle size distribution measuring device | |
| JP2015519649A5 (en) | ||
| JP6543318B2 (en) | Particle size distribution measuring device | |
| JP2018189516A (en) | Device, method, and program for calculating surface contamination density distribution | |
| Fritz-Popovski | Two-dimensional indirect Fourier transformation for evaluation of small-angle scattering data of oriented samples | |
| JP6555164B2 (en) | Particle size distribution measuring apparatus, data processing method, and data processing program | |
| Harder et al. | Emulating aerosol microphysics with machine learning | |
| US10712251B2 (en) | System for estimating planar drop sizes | |
| JP6112025B2 (en) | Data processing device for particle size distribution measurement, particle size distribution measuring device equipped with the same, data processing method for particle size distribution measurement, and data processing program for particle size distribution measurement | |
| US20220253508A1 (en) | Consecutive approximation calculation method, consecutive approximation calculation device, and program | |
| JP2012100961A (en) | Biological light measuring instrument and method for estimating absorption coefficient changes in living body | |
| JP7045910B2 (en) | Particle size distribution measuring device, particle size distribution measuring method and particle size distribution measuring program | |
| Konieczka et al. | Transformation of PET raw data into images for event classification using convolutional neural networks | |
| JP5443411B2 (en) | Method and apparatus for detecting the size of particles in a liquid |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20171205 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20181128 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20181213 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20190111 |
|
| 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: 20190604 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20190614 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6543318 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |