Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP6910641B2 - Small speaker design support device and speaker design support method - Google Patents
[go: Go Back, main page]

JP6910641B2 - Small speaker design support device and speaker design support method - Google Patents

Small speaker design support device and speaker design support method Download PDF

Info

Publication number
JP6910641B2
JP6910641B2 JP2017102449A JP2017102449A JP6910641B2 JP 6910641 B2 JP6910641 B2 JP 6910641B2 JP 2017102449 A JP2017102449 A JP 2017102449A JP 2017102449 A JP2017102449 A JP 2017102449A JP 6910641 B2 JP6910641 B2 JP 6910641B2
Authority
JP
Japan
Prior art keywords
acoustic
design
data
small speaker
characteristic data
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
JP2017102449A
Other languages
Japanese (ja)
Other versions
JP2018197965A (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.)
Kansai University
Original Assignee
Kansai University
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 Kansai University filed Critical Kansai University
Priority to JP2017102449A priority Critical patent/JP6910641B2/en
Publication of JP2018197965A publication Critical patent/JP2018197965A/en
Application granted granted Critical
Publication of JP6910641B2 publication Critical patent/JP6910641B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Description

本開示は、小型スピーカの設計支援装置及びスピーカの設計支援方法に関し、特に、スマートフォン等の小型機器に用いられるスピーカの設計を支援するための技術に関する。 The present disclosure relates to a small speaker design support device and a speaker design support method, and more particularly to a technique for supporting the design of a speaker used in a small device such as a smartphone.

近年、スマートフォンに代表される小型機器の需要が拡大し、それに伴ない、小型機器に内蔵されているスピーカの需要も拡大するとともに当該スピーカの小型化・複雑化が進んでいる。 In recent years, the demand for small devices such as smartphones has expanded, and along with this, the demand for speakers built into small devices has also increased, and the speakers have become smaller and more complex.

従来より、スピーカの設計手法としては、いわゆる音響等価回路解析法が知られている。音響等価回路解析法とは、振動板やダンパ等の機械振動系、及びキャビネットやホーン等の音響振動系を等価的な電気回路に置換えて、電気回路の知識を用いて対象の音響機器(スピーカ等)の音響特性を解析する手法である。このような音響等価回路解析法については、たとえば非特許文献1,2に記載されている。 Conventionally, a so-called acoustic equivalent circuit analysis method has been known as a speaker design method. The acoustic equivalent circuit analysis method replaces the mechanical vibration system such as a diaphragm and a damper and the acoustic vibration system such as a cabinet and a horn with an equivalent electric circuit, and uses the knowledge of the electric circuit to use the target acoustic equipment (speaker). Etc.) is a method for analyzing the acoustic characteristics. Such an acoustic equivalent circuit analysis method is described in, for example, Non-Patent Documents 1 and 2.

山本武夫編著、「スピーカ・システム(上)」、株式会社ラジオ技術社、1977年7月15日、p.35−44Edited by Takeo Yamamoto, "Speaker System (1)", Radio Technology Co., Ltd., July 15, 1977, p. 35-44 三井田惇郎著、「音響工学(21世紀を指向した電子・通信・情報カリキュラムシリーズ)」、昭晃堂、1987年3月、p.32−39"Acoustic Engineering (21st Century Oriented Electronics, Communication, Information Curriculum Series)" by Atsuro Mitsuida, Shokodo, March 1987, p. 32-39 豊田政弘編著、日本音響学会編、「FDTD法で視る音の世界(音響サイエンスシリーズ14)、コロナ社、2015年12月16日、p.1−46Edited by Masahiro Toyoda, edited by Acoustical Society of Japan, "The World of Sounds Seen by the FDTD Method (Acoustic Science Series 14), Corona Publishing Co., Ltd., December 16, 2015, p.1-46"

スマートフォン等の小型機器に内蔵されて用いられるような小型スピーカにおいては、音源となる振動板が筐体の内部に設けられ、筐体に設けられる音孔の寸法・位置や、振動板と音孔との間に設けられる気室の寸法等が重要な設計パラメータとなる。しかしながら、上述の音響等価回路解析法では、音孔の寸法や位置等を考慮することはできないので、上記解析法は、音孔の寸法や位置等について特段考慮する必要のない大型又は中型のスピーカの設計には有用であるけれども、音孔の寸法や位置等が重要な設計パラメータとなる小型スピーカの設計には適用することができない。このため、小型スピーカの設計については、十分な経験を積んだ技術者によるカット&トライ設計が行なわれており、設計に要する時間及びコストの削減が課題となっている。 In a small speaker that is built into a small device such as a smartphone, a diaphragm that serves as a sound source is provided inside the housing, and the dimensions and positions of the sound holes provided in the housing, as well as the diaphragm and sound hole. The dimensions of the air chamber provided between and the air chamber are important design parameters. However, since the above-mentioned acoustic equivalent circuit analysis method cannot consider the size and position of the sound hole, the above-mentioned analysis method does not need to consider the size and position of the sound hole in particular. However, it cannot be applied to the design of a small speaker in which the size and position of the sound hole are important design parameters. For this reason, in the design of a small speaker, a cut-and-try design is performed by an engineer who has sufficient experience, and reduction of time and cost required for the design is an issue.

本開示は、かかる課題を解決するためになされたものであり、その目的は、小型スピーカの設計において、設計に要する時間及びコストを削減可能な設計支援装置を提供することである。 The present disclosure has been made to solve such a problem, and an object of the present invention is to provide a design support device capable of reducing the time and cost required for designing a small speaker.

本開示の設計支援装置は、小型スピーカの設計支援装置であって、入力装置と、演算装置とを備える。入力装置は、小型スピーカの音響特性を示す音響特性データを利用者が入力するための装置である。演算装置は、音響特性データに基づいて小型スピーカの設計パラメータを出力するようにコンピュータを機能させるための学習済みモデルを用いて、入力装置により入力された音響特性データから設計パラメータを算出するように構成される。ここで、学習済みモデルは、音響FDTD(Finite-Difference Time-Domain)法を用いたコンピュータシミュレーションによって生成される、各々が設計パラメータと当該設計パラメータを有する小型スピーカの音響特性データとを含む複数の学習用データを用いて、音響特性データから設計パラメータを算出するように学習されたディープニューラルネットワーク(以下「DNN(Deep Neural Network)」と称する。)によって構成される。 The design support device of the present disclosure is a design support device for a small speaker, and includes an input device and an arithmetic unit. The input device is a device for the user to input acoustic characteristic data indicating the acoustic characteristics of the small speaker. The arithmetic unit now calculates the design parameters from the acoustic characteristic data input by the input device, using a trained model to make the computer function to output the design parameters of the small speaker based on the acoustic characteristic data. It is composed. Here, the trained model is generated by computer simulation using the acoustic FDTD (Finite-Difference Time-Domain) method, and includes a plurality of design parameters and acoustic characteristic data of a small speaker having the design parameters. It is composed of a deep neural network (hereinafter referred to as "DNN (Deep Neural Network)") trained to calculate design parameters from acoustic characteristic data using training data.

好ましくは、小型スピーカは、音源となる振動板と、小型スピーカの筐体の内部に形成される筐体内空間とを含んで構成される。そして、音響FDTD法を用いたコンピュータシミュレーションにおける分割グリッドの大きさは、振動板、筐体内空間、筐体の外部における音響伝搬空間の順に小さい。 Preferably, the small speaker is configured to include a diaphragm serving as a sound source and a space inside the housing formed inside the housing of the small speaker. The size of the divided grid in the computer simulation using the acoustic FDTD method is smaller in the order of the diaphragm, the space inside the housing, and the acoustic propagation space outside the housing.

好ましくは、各学習用データに含まれる設計パラメータは、予め定められた範囲内においてランダムに生成される。 Preferably, the design parameters included in each learning data are randomly generated within a predetermined range.

好ましくは、音響特性データは、所定の解析点における音圧の周波数特性を示すデータ含む。 Preferably, the acoustic characteristic data includes data indicating the frequency characteristic of sound pressure at a predetermined analysis point.

好ましくは、複数の学習用データとして用いられる音響特性データには、利用者の要求に従って、音圧の絶対値の周波数特性を示すデータ及び音圧レベルの周波数特性を示すデータのいずれかが選択される。 Preferably, as the acoustic characteristic data used as a plurality of learning data, either data showing the frequency characteristic of the absolute value of the sound pressure or data showing the frequency characteristic of the sound pressure level is selected according to the request of the user. NS.

好ましくは、音響特性データは、上記解析点における位相の周波数特性を示すデータをさらに含む。 Preferably, the acoustic characteristic data further includes data indicating the frequency characteristic of the phase at the analysis point.

好ましくは、設計パラメータは、小型スピーカに形成される音孔の寸法及び位置並びに気室の寸法を含む。 Preferably, the design parameters include the dimensions and location of the sound holes formed in the small speaker as well as the dimensions of the air chamber.

また、本開示の設計支援方法は、スピーカの設計支援方法であって、音響特性データに基づいてスピーカの設計パラメータを出力するようにコンピュータを機能させるための学習済みモデルを生成するステップと、音響特性データを利用者が入力するステップと、学習済みモデルを用いて、入力された音響特性データから設計パラメータを算出するステップとを含む。学習済みモデルを生成するステップは、音響FDTD法を用いたコンピュータシミュレーションによって、各々が設計パラメータと当該設計パラメータを有するスピーカの音響特性データとを含む複数の学習用データを生成するステップと、複数の学習用データを用いて、音響特性データから設計パラメータを算出するようにDNNを学習させることによって、学習済みモデルを生成するステップとを含む。 Further, the design support method of the present disclosure is a speaker design support method, which includes a step of generating a trained model for operating a computer to output a speaker design parameter based on acoustic characteristic data, and acoustics. It includes a step of inputting characteristic data by the user and a step of calculating design parameters from the input acoustic characteristic data using the trained model. The step of generating the trained model includes a step of generating a plurality of training data including a design parameter and acoustic characteristic data of a speaker having the design parameter, respectively, by computer simulation using the acoustic FDTD method, and a plurality of steps of generating the trained model. It includes a step of generating a trained model by training the DNN so as to calculate the design parameters from the acoustic characteristic data using the training data.

本開示の小型スピーカの設計支援装置においては、小型スピーカの音響特性データに基づいて小型スピーカの設計パラメータを出力するようにコンピュータを機能させるための学習済みモデルを用いて、入力装置により入力された音響特性データから設計パラメータが算出される。学習済みモデルは、音響特性データから設計パラメータを算出するように学習されたDNNによって構成される。これにより、技術者によるカット&トライ設計に頼ることなく小型スピーカの設計を行なうことができる。 In the small speaker design support device of the present disclosure, input is performed by an input device using a trained model for operating the computer to output the design parameters of the small speaker based on the acoustic characteristic data of the small speaker. Design parameters are calculated from the acoustic characteristic data. The trained model is composed of DNNs trained to calculate design parameters from acoustic characteristic data. This makes it possible to design a small speaker without relying on a cut-and-try design by an engineer.

ここで、各々が設計パラメータと当該設計パラメータを有する小型スピーカの音響特性データとを含む、DNNの学習に用いられる複数の学習用データは、音響FDTD法を用いたコンピュータシミュレーションによって生成される。これにより、複数の学習用データを用意するために実際に様々な音響構造を作成して実測を行なうことなく、DNNの学習に用いられる多数の学習用データを用意することができる。 Here, a plurality of learning data used for DNN learning, including design parameters and acoustic characteristic data of a small speaker having the design parameters, are generated by computer simulation using the acoustic FDTD method. As a result, it is possible to prepare a large number of learning data used for learning DNN without actually creating various acoustic structures to prepare a plurality of learning data and performing actual measurement.

このように、本開示の小型スピーカの設計支援装置によれば、技術者によるカット&トライ設計に頼ることなく、また、DNNの学習用データを用意するために実際に様々な音響構造を作成して実測を行なう必要もないので、小型スピーカの設計に要する時間及びコストを大幅に削減することができる。 As described above, according to the design support device for the small speaker of the present disclosure, various acoustic structures are actually created in order to prepare the learning data of DNN without relying on the cut and try design by the engineer. Since it is not necessary to actually measure the speaker, the time and cost required for designing the small speaker can be significantly reduced.

本開示の実施の形態に従う小型スピーカの設計支援装置の構成を示すブロック図である。It is a block diagram which shows the structure of the design support apparatus of a small speaker according to the embodiment of this disclosure. 図1に示す設計支援装置を用いて設計する小型スピーカの断面図である。FIG. 5 is a cross-sectional view of a small speaker designed by using the design support device shown in FIG. 図2に示す小型スピーカをx方向から視た平面図である。FIG. 5 is a plan view of the small speaker shown in FIG. 2 as viewed from the x direction. 図2,図3に示す小型スピーカの設計パラメータ、及び各パラメータのレンジの一例を示した図である。It is a figure which showed the design parameter of the small speaker shown in FIGS. 2 and 3 and an example of the range of each parameter. 図2に示す解析点P0における音圧(Pa)の周波数特性を示した図である。It is a figure which showed the frequency characteristic of the sound pressure (Pa) at the analysis point P0 shown in FIG. 図2に示す解析点P0における音波の位相の周波数特性を示した図である。It is a figure which showed the frequency characteristic of the phase of the sound wave at the analysis point P0 shown in FIG. 図1に示す演算装置において用いられる学習済みモデルを生成する学習システムの全体ブロック図である。It is an whole block diagram of the learning system which generates the trained model used in the arithmetic unit shown in FIG. 図7に示す学習システムに用いられる学習用データを生成する生成システムの構成図である。It is a block diagram of the generation system which generates the learning data used for the learning system shown in FIG. 7. 小型スピーカをx方向から視たときの外形図である。It is an outline drawing when a small speaker is viewed from the x direction. 小型スピーカをx方向から視たときの外形図である。It is an outline drawing when a small speaker is viewed from the x direction. 条件設定部によるグリッド設定の考え方の一例を説明する図である。It is a figure explaining an example of the concept of grid setting by a condition setting part. DNN学習部により生成される学習済みモデルに用いられるDNNの構成図である。It is a block diagram of the DNN used for the trained model generated by the DNN learning unit. 本実施の形態におけるDNNの入出力を示した図である。It is a figure which showed the input / output of DNN in this embodiment. 図1に示す演算装置において用いられる学習済みモデルの生成手順を説明するフローチャートである。It is a flowchart explaining the generation procedure of the trained model used in the arithmetic unit shown in FIG. 図1に示す演算装置において実行される処理の手順を説明するフローチャートである。It is a flowchart explaining the procedure of the process executed in the arithmetic unit shown in FIG. 音圧の絶対値(Pa)の周波数特性を示した図である。It is a figure which showed the frequency characteristic of the absolute value (Pa) of a sound pressure. 音圧レベル(dB)の周波数特性を示した図である。It is a figure which showed the frequency characteristic of the sound pressure level (dB).

以下、本開示の実施の形態について、図面を参照しながら詳細に説明する。なお、図中同一又は相当部分には同一符号を付してその説明は繰り返さない。 Hereinafter, embodiments of the present disclosure will be described in detail with reference to the drawings. The same or corresponding parts in the drawings are designated by the same reference numerals, and the description thereof will not be repeated.

<設計支援装置の全体構成>
図1は、本開示の実施の形態に従う小型スピーカの設計支援装置の構成を示すブロック図である。図1を参照して、この設計支援装置1は、スマートフォンに代表される小型機器に内蔵される小型スピーカの設計をコンピュータを用いて行なうものである。設計支援装置1は、入力装置10と、演算装置20と、表示装置30とを備える。
<Overall configuration of design support device>
FIG. 1 is a block diagram showing a configuration of a design support device for a small speaker according to an embodiment of the present disclosure. With reference to FIG. 1, the design support device 1 uses a computer to design a small speaker built in a small device represented by a smartphone. The design support device 1 includes an input device 10, an arithmetic unit 20, and a display device 30.

入力装置10は、小型スピーカの設計者である本設計支援装置1のユーザが、設計対象の小型スピーカの音響特性を示すデータ(音響特性データ)を入力するための装置である。小型スピーカの音響特性は、所定の解析点における音圧や位相の周波数特性等によって表すことができる。この実施の形態では、ユーザは、設計対象の小型スピーカについて、所望の音圧周波数特性及び位相周波数特性を入力装置10から入力することができる。音響特性データについては、後ほど改めて詳しく説明する。 The input device 10 is a device for a user of the design support device 1 who is a designer of the small speaker to input data (acoustic characteristic data) indicating the acoustic characteristics of the small speaker to be designed. The acoustic characteristics of the small speaker can be represented by the sound pressure at a predetermined analysis point, the frequency characteristics of the phase, and the like. In this embodiment, the user can input desired sound pressure frequency characteristics and phase frequency characteristics from the input device 10 for the small speaker to be designed. The acoustic characteristic data will be described in detail later.

演算装置20は、CPU(Central Processing Unit)、処理プログラム及び後述の学習済みモデルを記憶するROM(Read Only Memory)、データを一時的に記憶するRAM(Random Access Memory)等を含んで構成される。学習済みモデルについては、その生成方法を含めて後ほど詳細に説明するが、この学習済みモデルは、音響特性データに基づいて小型スピーカの設計パラメータ(音孔の寸法や位置等)を出力するようにコンピュータ(演算装置20)を機能させるためのモデルであって、音響特性データと設計パラメータとの因果関係が学習されたDNNによって構成される。演算装置20は、この学習済みモデルを用いて、入力装置10から受ける音響特性データから小型スピーカの設計パラメータを算出する。 The arithmetic unit 20 includes a CPU (Central Processing Unit), a ROM (Read Only Memory) for storing a processing program and a learned model described later, a RAM (Random Access Memory) for temporarily storing data, and the like. .. The trained model will be explained in detail later, including how to generate it, but this trained model should output the design parameters (sound hole dimensions, position, etc.) of the small speaker based on the acoustic characteristic data. It is a model for operating a computer (computing device 20), and is composed of a DNN in which a causal relationship between acoustic characteristic data and design parameters has been learned. Using this trained model, the arithmetic unit 20 calculates the design parameters of the small speaker from the acoustic characteristic data received from the input device 10.

表示装置30は、演算装置20によって算出された小型スピーカの設計パラメータをユーザに表示するための装置である。入力装置10から所望の音響特性データを入力したユーザは、演算装置20によって算出された設計パラメータ、或いはそのような設計パラメータを有する小型スピーカの構造を表示装置30において確認することができる。表示装置30は、たとえばタッチパネルディスプレイ等によって入力装置10と一体的に構成されてもよい。 The display device 30 is a device for displaying the design parameters of the small speaker calculated by the arithmetic unit 20 to the user. A user who has input desired acoustic characteristic data from the input device 10 can confirm the design parameters calculated by the arithmetic unit 20 or the structure of a small speaker having such design parameters on the display device 30. The display device 30 may be integrally configured with the input device 10 by, for example, a touch panel display or the like.

<設計パラメータ及び音響特性データの説明>
まず、図1に示した設計支援装置1を用いて設計される小型スピーカの構成及び設計パラメータ、並びにそのような小型スピーカの音響特性を示す音響特性データについて説明する。
<Explanation of design parameters and acoustic characteristic data>
First, the configuration and design parameters of the small speaker designed by using the design support device 1 shown in FIG. 1 and the acoustic characteristic data showing the acoustic characteristics of such a small speaker will be described.

図2及び図3は、図1に示した設計支援装置1を用いて設計する小型スピーカの構成図である。図2は、小型スピーカの断面図を示し、図3は、図2に示した小型スピーカをx方向から視た平面図である。この図2及び図3に示される小型スピーカの構成は、小型スピーカの設計において重要なパラメータ(音孔の寸法や位置等)に絞って小型スピーカを模式化したものである。 2 and 3 are block diagrams of a small speaker designed by using the design support device 1 shown in FIG. FIG. 2 shows a cross-sectional view of the small speaker, and FIG. 3 is a plan view of the small speaker shown in FIG. 2 as viewed from the x direction. The configuration of the small speaker shown in FIGS. 2 and 3 is a schematic representation of the small speaker focusing on important parameters (dimensions and positions of sound holes, etc.) in the design of the small speaker.

図2及び図3を参照して、設計支援装置1を用いて設計される小型スピーカ200は、筐体210と、振動板220と、音孔230と、気室240とを備える。音響伝搬空間250は、音孔230を通じて筐体210外に出力される音波が伝搬する空間である。なお、境界260は、後述の音響FDTD法によるコンピュータシミュレーションにおいて解析対象領域を有限にするために導入される吸収境界であって、たとえば完全吸収層(Perfectly matched Layer:PML)である。 The small speaker 200 designed by using the design support device 1 with reference to FIGS. 2 and 3 includes a housing 210, a diaphragm 220, a sound hole 230, and an air chamber 240. The acoustic propagation space 250 is a space in which sound waves output to the outside of the housing 210 propagate through the sound holes 230. The boundary 260 is an absorbing boundary introduced to make the analysis target region finite in a computer simulation by the acoustic FDTD method described later, and is, for example, a perfectly matched layer (PML).

振動板220は、小型スピーカ200の音源である。すなわち、振動板220が振動することによって音波が発生する。なお、特に図示しないが、振動板220の駆動源には、ボイスコイルを用いてもよいし(動電変換方式)、印加電圧に応じて伸縮する圧電素子を用いてもよい(圧電変換方式)。 The diaphragm 220 is a sound source of the small speaker 200. That is, sound waves are generated by the vibration of the diaphragm 220. Although not particularly shown, a voice coil may be used as the drive source of the diaphragm 220 (piezoelectric conversion method), or a piezoelectric element that expands and contracts according to the applied voltage may be used (piezoelectric conversion method). ..

音孔230は、筐体210に形成され、振動板220で発生した音波は、音孔230を通じて筐体210外の音響伝搬空間250へ出力される。気室240は、振動板220が設けられる空間と音孔230との間に設けられる。この実施の形態では、気室240は、一つだけ設けられるものとしているが、x方向に複数の気室を設けてもよい。 The sound hole 230 is formed in the housing 210, and the sound wave generated by the diaphragm 220 is output to the sound propagation space 250 outside the housing 210 through the sound hole 230. The air chamber 240 is provided between the space where the diaphragm 220 is provided and the sound hole 230. In this embodiment, only one air chamber 240 is provided, but a plurality of air chambers may be provided in the x direction.

LHx(図2)は、音孔230のx方向の長さを示す。また、LHy(図2,図3)は、音孔230のy方向の長さを示し、LHz(図3)は、音孔230のz方向の長さを示す。LCx(図2)は、気室240のx方向の長さを示す。Hymin,Hzminは、筐体210における音孔230の位置を示す。たとえば、Hymin=0,Hzmin=0であれば、音孔230の位置は、図3において筐体210の左下隅である。 LHx (FIG. 2) shows the length of the sound hole 230 in the x direction. Further, LHy (FIGS. 2 and 3) indicates the length of the sound hole 230 in the y direction, and LHz (FIG. 3) indicates the length of the sound hole 230 in the z direction. LCx (FIG. 2) shows the length of the air chamber 240 in the x direction. Hymin and Hzmin indicate the positions of the sound holes 230 in the housing 210. For example, if Hymin = 0 and Hzmin = 0, the position of the sound hole 230 is the lower left corner of the housing 210 in FIG.

音響伝搬空間250には、解析点P0が設定される(図2)。解析点P0は、小型スピーカ200の音響特性を測定・解析するポイントであり、小型スピーカ200の音響特性は、解析点P0における音響特性で代表されるものとし、解析点P0において所望の音響特性が得られるように小型スピーカ200の設計が行なわれる。なお、この実施の形態では、解析点P0は、筐体210の端面における音孔230の中心からx方向にL0だけ離れた位置に設定される。 An analysis point P0 is set in the acoustic propagation space 250 (FIG. 2). The analysis point P0 is a point for measuring and analyzing the acoustic characteristics of the small speaker 200. The acoustic characteristics of the small speaker 200 are represented by the acoustic characteristics at the analysis point P0, and the desired acoustic characteristics are obtained at the analysis point P0. The small speaker 200 is designed so as to be obtained. In this embodiment, the analysis point P0 is set at a position separated by L0 in the x direction from the center of the sound hole 230 on the end face of the housing 210.

図4は、図2,図3に示した小型スピーカ200の設計パラメータ、及び各パラメータのレンジの一例を示した図である。図4とともに図2,図3も参照して、この実施の形態では、設計パラメータは、LCx,LHx,LHy,LHz,Hymin,Hzmin(いずれも10-3m)である。すなわち、LCxは気室の寸法を示し、LHx,LHy,LHzは音孔230の寸法を示し、Hymin,Hzminは音孔230の位置を示す。なお、以下では、これらの設計パラメータは、設計変数ベクトルX=[LCx,LHx,LHy,LHz,Hymin,Hzmin]で表されることがある。 FIG. 4 is a diagram showing an example of the design parameters of the small speaker 200 shown in FIGS. 2 and 3 and the range of each parameter. With reference to FIGS. 2 and 3 as well as FIG. 4, in this embodiment, the design parameters are LCx, LHx, LHy, LHz, Hymin, Hzmin (all 10-3 m). That is, LCx indicates the size of the air chamber, LHx, LHy, and LHz indicate the size of the sound hole 230, and Hymin and Hzmin indicate the position of the sound hole 230. In the following, these design parameters may be represented by the design variable vector X = [LCx, LHx, LHy, LHz, Hymin, Hzmin].

各パラメータのレンジは、設計可能な範囲を示す。この実施の形態では、LCx,LHxの各レンジは、0.2〜2.0(10-3m)であり、LHy,LHzの各レンジは、0.5〜3.0(10-3m)であり、Hymin,Hzminの各レンジは、0.5〜11.5(10-3m)である。 The range of each parameter indicates the designable range. In this embodiment, the LCx and LHx ranges are 0.2 to 2.0 (10 -3 m), and the LHy and LHz ranges are 0.5 to 3.0 (10 -3 m). ), And each range of Hymin and Hzmin is 0.5 to 11.5 (10 -3 m).

なお、小型スピーカにおいては、気室に対する音孔の寸法及び位置が小型スピーカの音響特性に特に大きく影響するため、この実施の形態では、気室240の幅(y方向及びz方向の長さ)は大きく一定であるとし、気室240のx方向の幅、並びに音孔230の寸法及び位置を設計パラメータとしている。 In the small speaker, the size and position of the sound hole with respect to the air chamber have a particularly large effect on the acoustic characteristics of the small speaker. Therefore, in this embodiment, the width of the air chamber 240 (the length in the y direction and the z direction). Is large and constant, and the width of the air chamber 240 in the x direction and the dimensions and position of the sound hole 230 are set as design parameters.

図5及び図6は、小型スピーカ200の音響特性データの一例を示した図である。図5は、図2に示した解析点P0における音圧(Pa)の周波数特性を示し、図6は、解析点P0における音波の位相(rad)の周波数特性を示す。なお、振動板220(図2)には、ガウシアンパルスを入力として与えている(後述)。 5 and 6 are diagrams showing an example of acoustic characteristic data of the small speaker 200. FIG. 5 shows the frequency characteristic of the sound pressure (Pa) at the analysis point P0 shown in FIG. 2, and FIG. 6 shows the frequency characteristic of the phase (rad) of the sound wave at the analysis point P0. A Gaussian pulse is given as an input to the diaphragm 220 (FIG. 2) (described later).

図5及び図6を参照して、このような音圧の周波数特性のプロフィール(周波数波形)や、位相の周波数特性のプロフィールは、小型スピーカ200の音響特性を示す音響特性データの一例である。また、音圧がピークをとる周波数(位相が急変する周波数)は、小型スピーカ200の共振周波数を示し、このような共振周波数も、小型スピーカ200の音響特性データの一例である。 With reference to FIGS. 5 and 6, such a profile (frequency waveform) of the frequency characteristic of sound pressure and a profile of the frequency characteristic of the phase are examples of acoustic characteristic data showing the acoustic characteristic of the small speaker 200. Further, the frequency at which the sound pressure peaks (the frequency at which the phase changes suddenly) indicates the resonance frequency of the small speaker 200, and such a resonance frequency is also an example of the acoustic characteristic data of the small speaker 200.

再び図1を参照して、この設計支援装置1は、図4に示した設計パラメータで示される小型スピーカ200(図2,図3)の設計に用いられる。この設計支援装置1においては、所望の音響特性データ(図5,図6)がユーザによって入力装置10から入力され、そのような音響特性データを示す小型スピーカ200の設計パラメータ(図4)が演算装置20によって算出される。 With reference to FIG. 1 again, this design support device 1 is used for designing the small speaker 200 (FIGS. 2 and 3) shown by the design parameters shown in FIG. In the design support device 1, desired acoustic characteristic data (FIGS. 5 and 6) is input from the input device 10 by the user, and a design parameter (FIG. 4) of the small speaker 200 showing such acoustic characteristic data is calculated. Calculated by device 20.

上述のように、小型スピーカ200においては、音孔230の寸法・位置(LHx,LHy,LHz,Hymin,Hzmin)や、気室240の寸法(LCx)等が重要な設計パラメータである。しかしながら、上記の特許文献1,2に示される音響等価回路解析法では、音孔の寸法や位置等を考慮することはできないので、この解析法は、音孔の寸法や位置等について特段考慮する必要のない大型又は中型のスピーカの設計には有用であるけれども、音孔230の寸法や位置等が重要な設計パラメータとなる小型スピーカ200の設計には適用できない。 As described above, in the small speaker 200, the dimensions and positions of the sound holes 230 (LHx, LHy, LHz, Hymin, Hzmin), the dimensions of the air chamber 240 (LCx), and the like are important design parameters. However, since the acoustic equivalent circuit analysis method shown in Patent Documents 1 and 2 above cannot consider the size and position of the sound hole, this analysis method particularly considers the size and position of the sound hole. Although it is useful for designing a large or medium-sized speaker that does not need it, it cannot be applied to the design of a small speaker 200 in which the size and position of the sound hole 230 are important design parameters.

そこで、この実施の形態に従う設計支援装置1では、音響特性データに基づいて小型スピーカ200の設計パラメータを出力するようにコンピュータ(演算装置20)を機能させるための学習済みモデルを用いて、入力装置10により入力される音響特性データから設計パラメータを算出することとしたものである。学習済みモデルは、音響特性データから設計パラメータを算出するように学習されたDNNによって構成される。これにより、技術者によるカット&トライ設計に頼ることなく小型スピーカ200の設計を行なうことができる。 Therefore, in the design support device 1 according to this embodiment, an input device is used by using a trained model for operating the computer (arithmetic unit 20) so as to output the design parameters of the small speaker 200 based on the acoustic characteristic data. The design parameters are calculated from the acoustic characteristic data input by 10. The trained model is composed of DNNs trained to calculate design parameters from acoustic characteristic data. As a result, the small speaker 200 can be designed without relying on a cut-and-try design by an engineer.

ここで、DNNによる音響特性データと設計パラメータとの因果関係を精度よく学習するには、多数の学習用データを用意する必要がある。各学習用データは、設計パラメータのとり得るレンジ(図4)において適宜設定される設計パラメータと、当該設計パラメータを有する小型スピーカ200の音響特性データとから成るものである。そして、この実施の形態に従う設計支援装置1では、DNNの学習に用いられる複数の学習用データが、音響FDTD法を用いたコンピュータシミュレーションによって生成される。これにより、多数の学習用データを用意するために実際に様々な音響構造を作成して実測を行なうことなく、DNNの学習に用いられる多数の学習用データを用意することができる。 Here, in order to accurately learn the causal relationship between the acoustic characteristic data by DNN and the design parameters, it is necessary to prepare a large number of learning data. Each learning data is composed of design parameters appropriately set in a range (FIG. 4) in which design parameters can be taken, and acoustic characteristic data of a small speaker 200 having the design parameters. Then, in the design support device 1 according to this embodiment, a plurality of learning data used for learning the DNN are generated by a computer simulation using the acoustic FDTD method. As a result, it is possible to prepare a large number of learning data used for DNN learning without actually creating various acoustic structures to prepare a large number of learning data and performing actual measurement.

以下、この実施の形態に従う設計支援装置1の演算装置20において用いられる学習済みモデルについて説明する。 Hereinafter, the trained model used in the arithmetic unit 20 of the design support device 1 according to this embodiment will be described.

<学習済みモデルの説明>
図7は、図1に示した演算装置20において用いられる学習済みモデルを生成する学習システムの全体ブロック図である。図7を参照して、この学習システム100は、データベース120と、DNN学習部130と、学習条件設定部140とを含む。
<Explanation of trained model>
FIG. 7 is an overall block diagram of a learning system that generates a trained model used in the arithmetic unit 20 shown in FIG. With reference to FIG. 7, the learning system 100 includes a database 120, a DNN learning unit 130, and a learning condition setting unit 140.

データベース120は、DNN学習部130において学習済みモデルを生成するための学習用データ110−1,110−2,・・・,110−Nを格納する。学習用データ110−1〜110−Nの各々は、上述の設計パラメータ(設計変数ベクトルX)と、その設計パラメータを有する小型スピーカ200の音響特性データ(音圧周波数特性及び位相周波数特性)との対のデータである。学習用データ110−1〜110−Nは、図4に示した各パラメータのレンジ内で互いに異なる設計パラメータを有する。 The database 120 stores learning data 110-1, 110-2, ..., 110-N for generating a trained model in the DNN learning unit 130. Each of the learning data 110-1 to 110-N has the above-mentioned design parameter (design variable vector X) and the acoustic characteristic data (sound pressure frequency characteristic and phase frequency characteristic) of the small speaker 200 having the design parameter. It is a pair of data. The training data 110-1 to 110-N have different design parameters within the range of each parameter shown in FIG.

DNN学習部130において精度の高い学習済みモデルを生成するには、設計パラメータのとり得るレンジ内で様々な設計パラメータに対する多数の学習用データ110−1〜110−Nを用意することが重要である。多数の学習用データ110−1〜110−Nを用意するために実際に様々な音響構造を作成して実測を行なうことは、多大な時間とコストを有する。そこで、この実施の形態では、音響FDTD法を用いたコンピュータシミュレーションによって多数の学習用データ110−1〜110−Nが生成される。 In order to generate a highly accurate trained model in the DNN learning unit 130, it is important to prepare a large number of training data 110-1 to 110-N for various design parameters within the range that the design parameters can take. .. It takes a lot of time and cost to actually create various acoustic structures and actually measure them in order to prepare a large number of learning data 110-1 to 110-N. Therefore, in this embodiment, a large number of learning data 110-1 to 110-N are generated by computer simulation using the acoustic FDTD method.

図8は、図7に示した学習システム100に用いられる学習用データを生成する生成システムの構成図である。図8を参照して、この生成システム300は、設計パラメータ設定部310と、音響FDTD演算部320と、条件設定部330とを含む。 FIG. 8 is a configuration diagram of a generation system that generates learning data used in the learning system 100 shown in FIG. 7. With reference to FIG. 8, the generation system 300 includes a design parameter setting unit 310, an acoustic FDTD calculation unit 320, and a condition setting unit 330.

設計パラメータ設定部310は、多数の学習用データ110−1〜110−Nを生成するために、図4に示したレンジ内で互いに異なる多数の設計パラメータを生成する。たとえば、設計パラメータ設定部310は、数千個レベルの互いに異なる設計パラメータを生成する。 The design parameter setting unit 310 generates a large number of design parameters different from each other within the range shown in FIG. 4 in order to generate a large number of learning data 110-1 to 110-N. For example, the design parameter setting unit 310 generates thousands of different design parameters from each other.

多数の設計パラメータから多数の学習用データを生成するには、多大な負荷と時間がかかり得るところ、この実施の形態では、学習用データを生成する音響FDTD演算部320の演算負荷の低減を目的に、設計パラメータの対称性(線対称、点対称)が考慮され、重複する演算を回避している。さらに、設計パラメータ設定部310は、多数の学習用データを用いて生成される学習済みモデルの精度を高めるために、多数の設計パラメータをランダムに生成する。 It may take a large load and time to generate a large number of training data from a large number of design parameters. However, in this embodiment, the purpose is to reduce the calculation load of the acoustic FDTD calculation unit 320 that generates the training data. In addition, the symmetry of the design parameters (line symmetry, point symmetry) is taken into consideration, and duplicate operations are avoided. Further, the design parameter setting unit 310 randomly generates a large number of design parameters in order to improve the accuracy of the trained model generated by using the large number of training data.

図9,図10は、図2,図3に示した小型スピーカ200をx方向から視たときの外形図である。なお、図9では、音孔230については図示していない。図9を参照して、小型スピーカ200をx方向から視たときの小型スピーカ200の外形は、正方形であるものとする。この場合、形状の対称性を考慮すると、音孔230の位置については、音孔230の中心が斜線で示される領域270に含まれる設計パラメータを用意すれば十分であり、音孔230の中心がその他の領域に含まれる設計パラメータを生成する必要はない。 9 and 10 are external views of the small speaker 200 shown in FIGS. 2 and 3 when viewed from the x direction. Note that FIG. 9 does not show the sound hole 230. With reference to FIG. 9, it is assumed that the outer shape of the small speaker 200 when the small speaker 200 is viewed from the x direction is a square. In this case, considering the symmetry of the shape, it is sufficient to prepare the design parameters included in the region 270 in which the center of the sound hole 230 is shaded, and the center of the sound hole 230 is located. There is no need to generate design parameters contained in other areas.

そこで、設計パラメータ設定部310は、小型スピーカ200の形状の対称性を考慮して、音孔230の中心が領域270に含まれる設計パラメータについて、多数の設計パラメータを生成する。このように、設計パラメータの対称性を考慮することによって、生成された設計パラメータ毎に音響特性データを算出する音響FDTD演算部320の演算負荷及び演算時間を低減することができる。 Therefore, the design parameter setting unit 310 generates a large number of design parameters for the design parameters in which the center of the sound hole 230 is included in the region 270 in consideration of the symmetry of the shape of the small speaker 200. In this way, by considering the symmetry of the design parameters, it is possible to reduce the calculation load and the calculation time of the acoustic FDTD calculation unit 320 that calculates the acoustic characteristic data for each generated design parameter.

さらに、図10を参照して、設計パラメータ設定部310は、設計パラメータのとり得るレンジ内で(対称性を考慮して音孔230の中心が図9の領域270に含まれる範囲で)、ランダムに多数の設計パラメータを生成する。「ランダム」とは、決められた刻み幅で設計パラメータを変化させるのではなく、任意に設計パラメータを変化させるという意味である。ランダムに変化させた設計パラメータに基づく多数の学習用データを用いることによって、DNNにより学習済みモデルを生成する際に学習が局所解に陥ることを回避して学習済みモデルの精度を高めることができる。 Further, referring to FIG. 10, the design parameter setting unit 310 randomly sets the design parameter within the range that the design parameter can take (within the range in which the center of the sound hole 230 is included in the region 270 of FIG. 9 in consideration of symmetry). Generate a large number of design parameters. "Random" means that the design parameters are not changed at a fixed step size, but are arbitrarily changed. By using a large number of training data based on randomly changed design parameters, it is possible to improve the accuracy of the trained model by avoiding the learning from falling into a local solution when generating the trained model by DNN. ..

再び図8を参照して、音響FDTD演算部320は、CPU、GPU(Graphics Processing Unit)、処理プログラムを記憶するROM、データを一時的に記憶するRAM等を含んで構成される。そして、音響FDTD演算部320は、設計パラメータ設定部310において生成される多数の設計パラメータを受け、音響FDTD法を用いたコンピュータシミュレーションによって、各設計パラメータに対して小型スピーカ200が当該設計パラメータを有する場合の音響特性データを算出する。 With reference to FIG. 8 again, the acoustic FDTD calculation unit 320 includes a CPU, a GPU (Graphics Processing Unit), a ROM for storing a processing program, a RAM for temporarily storing data, and the like. Then, the acoustic FDTD calculation unit 320 receives a large number of design parameters generated by the design parameter setting unit 310, and the small speaker 200 has the design parameters for each design parameter by computer simulation using the acoustic FDTD method. Calculate the acoustic characteristic data of the case.

FDTD法とは、もともと電磁波の伝搬解析に用いられてきた手法であり、電磁波の支配式であるマクスウェル方程式を時間及び空間で差分化し、リープフロッグアルゴリズムを用いて解析空間の電磁界を時間的に更新することで、出力点の時間応答を得る手法である。音響FDTD法は、このFDTD法を音波の伝搬解析に適用したものであり、音圧と粒子速度とを用いて波動伝搬の計算が行なわれる。この音響FDTD法は、音孔及び気室の形状や位置関係を考慮した音響特性解析を行なうことが可能である。以下、音響FDTD法の概略を説明する(詳細は上記の特許文献3を参照)。 The FDTD method is a method originally used for propagation analysis of electromagnetic waves. Maxwell's equations, which are the dominant equations of electromagnetic waves, are differentiated in time and space, and the electromagnetic field in the analysis space is temporally differentiated by using the leapfrog algorithm. It is a method to obtain the time response of the output point by updating. The acoustic FDTD method is an application of this FDTD method to sound wave propagation analysis, and wave propagation is calculated using sound pressure and particle velocity. This acoustic FDTD method can perform acoustic characteristic analysis in consideration of the shape and positional relationship of the sound hole and the air chamber. The outline of the acoustic FDTD method will be described below (see Patent Document 3 above for details).

音響FDTD法では、スタガードグリッドと呼ばれる、音圧及び粒子速度の参照点が空間的・時間的に互い違いに配置された構造格子を用いて、解析空間における音圧及び粒子速度が空間的・時間的に離散化される。そして、リープフロッグアルゴリズムを用いて、空間的・時間的に離散化された音圧及び粒子速度が交互に計算される。 In the acoustic FDTD method, the sound pressure and particle velocity in the analysis space are spatially and temporally determined by using a structured grid in which the reference points of sound pressure and particle velocity are arranged in a spatially and temporally staggered manner, which is called a staggered grid. Is discreteized into. Then, using the leap frog algorithm, the spatially and temporally discretized sound pressure and particle velocity are alternately calculated.

音響FDTD法では、解析空間における空気粒子(格子状の微小体積)の運動及び変形は、以下に示される方程式によって支配される。 In the acoustic FDTD method, the motion and deformation of air particles (small grid-like volumes) in the analysis space are governed by the equations shown below.

Figure 0006910641
Figure 0006910641

式(1)〜(3)は、それぞれ空気粒子のx,y,z方向の運動方程式であり、式(4)は、音圧に関する連続方程式である。ρは空気密度を示し、u,v,wはそれぞれx,y,z方向の粒子速度を示す。また、pは音圧を示し、κは体積弾性率を示す。上記の式(1)〜(4)は、纏めて「支配方程式」と称される。 Equations (1) to (3) are equations of motion of air particles in the x, y, and z directions, respectively, and equation (4) is a continuity equation relating to sound pressure. ρ indicates the air density, and u, v, and w indicate the particle velocities in the x, y, and z directions, respectively. Further, p indicates sound pressure and κ indicates bulk modulus. The above equations (1) to (4) are collectively referred to as "governing equations".

上記の支配方程式を解くことによって解析空間内の音圧及び粒子速度を求めることができるが、これをコンピュータ上で実現するために、音響FDTD法では、上記の支配方程式を空間及び時間について離散化する。具体的には、音圧に関して、x,y,z方向の空間離散化幅をそれぞれΔx,Δy,Δzとし、時間離散化幅をΔtとした場合、粒子速度に関して、空間についてはΔx/2,Δy/2,Δz/2、時間についてはΔt/2だけ音圧の参照点からずらして離散化する。 The sound pressure and particle velocity in the analysis space can be obtained by solving the above governing equation, but in order to realize this on a computer, the acoustic FDTD method discretizes the above governing equation with respect to space and time. do. Specifically, when the spatial discretization widths in the x, y, and z directions are Δx, Δy, and Δz for sound pressure and the time discretization width is Δt, respectively, the particle velocity is Δx / 2, for space. Discretize Δy / 2, Δz / 2 and time by shifting from the sound pressure reference point by Δt / 2.

このような離散化に伴ない上記の支配方程式を差分近似すると、式(1)〜(4)の支配方程式はそれぞれ次式のように変形される。 When the above governing equations are differentially approximated with such discretization, the governing equations of equations (1) to (4) are transformed as shown in the following equations.

Figure 0006910641
Figure 0006910641

i,j,kは、解析空間においてそれぞれx,y,z方向に関して何番目の音圧参照点を示す空間ステップである。nは、時間に関して何番目の参照点であるかを示す時間ステップである。たとえば、un+1(i+1/2,j,k)は、空間ステップがi+1/2,j,kで時間ステップがn+1であるときのx方向の粒子速度の値を示し、pn+1/2(i,j,k)は、空間ステップがi,j,kで時間ステップがn+1/2であるときの音圧の値を示す。 i, j, and k are spatial steps indicating the number of sound pressure reference points in the x, y, and z directions in the analysis space, respectively. n is a time step indicating the number of reference points in terms of time. For example, un + 1 (i + 1/2, j, k) indicates the value of the particle velocity in the x direction when the spatial step is i + 1/2, j, k and the time step is n + 1, and p n + 1. / 2 (i, j, k) indicates the value of sound pressure when the spatial step is i, j, k and the time step is n + 1/2.

式(5)〜(7)で示される粒子速度と、式(8)で示される音圧とを解析空間(全空間ステップ)について交互に計算することによって、次々と新しい時間ステップにおける状態を算出することができる(リープフロッグアルゴリズム)。 By alternately calculating the particle velocity represented by the equations (5) to (7) and the sound pressure represented by the equation (8) for the analysis space (entire space step), the state in the new time step is calculated one after another. Can be (Leap Frog algorithm).

図8を参照して、条件設定部330は、音響FDTD演算部320における音響FDTD法を用いたコンピュータシミュレーションの演算条件を設定する。具体的には、条件設定部330は、支配方程式で用いられる音圧pや体積弾性率κ、空間離散化幅Δx,Δy,Δz、時間離散化幅Δt、音源である振動板220(図2)への入力、時間ステップ数、その他音速等の各種定数を音響FDTD演算部320に設定する。なお、筐体210(図2)については、剛体によって構成され、音波は全反射するものとして設定される。 With reference to FIG. 8, the condition setting unit 330 sets the calculation conditions of the computer simulation using the acoustic FDTD method in the acoustic FDTD calculation unit 320. Specifically, the condition setting unit 330 includes the sound pressure p and bulk modulus κ used in the governing equation, the spatial discretization width Δx, Δy, Δz, the time discretization width Δt, and the diaphragm 220 which is a sound source (FIG. 2). ), The number of time steps, and other constants such as the speed of sound are set in the acoustic FDTD calculation unit 320. The housing 210 (FIG. 2) is formed of a rigid body, and the sound wave is set to be totally reflected.

なお、振動板220への入力について、この実施の形態では、下式にて表されるガウシアンパルスを振動板220に与えるものとし、ガウシアンパルスによって振動板220を励振させたときの音響特性データが音響FDTD演算部320により算出される。 Regarding the input to the diaphragm 220, in this embodiment, the Gaussian pulse represented by the following equation is given to the diaphragm 220, and the acoustic characteristic data when the diaphragm 220 is excited by the Gaussian pulse is obtained. It is calculated by the acoustic FDTD calculation unit 320.

g(t)=A・exp[−{(nΔt−T)/0.29T}2] …(9)
T=0.646/fc …(10)
上記において、Aは最大振幅値であり、fcはガウシアンパルスの振幅の周波数スペクトルが3dB低下する周波数である。
g (t) = A · exp [− {(nΔt−T) /0.29T} 2 ]… (9)
T = 0.646 / fc ... (10)
In the above, A is the maximum amplitude value, and fc is the frequency at which the frequency spectrum of the amplitude of the Gaussian pulse is reduced by 3 dB.

音響FDTD演算部320は、条件設定部330により設定される演算条件に従って、設計パラメータ設定部310から受ける、ランダムに生成された多数の設計パラメータの各々に対して、音響FDTD法を用いたコンピュータシミュレーションにより解析空間(図2,図3に示した小型スピーカ200の境界260内の空間)における音圧及び粒子速度を逐次算出する。そして、音響FDTD演算部320は、各設計パラメータに対して、図2に示した解析点P0における音圧の時間応答データを取得し、その時間応答データに対して高速フーリエ変換(FFT)処理等を行なうことによって音響特性データ(解析点P0における音圧周波数特性及び位相周波数特性のデータ)を算出する。 The acoustic FDTD calculation unit 320 uses a computer simulation using the acoustic FDTD method for each of a large number of randomly generated design parameters received from the design parameter setting unit 310 according to the calculation conditions set by the condition setting unit 330. The sound pressure and particle velocity in the analysis space (the space within the boundary 260 of the small speaker 200 shown in FIGS. 2 and 3) are sequentially calculated. Then, the acoustic FDTD calculation unit 320 acquires the time response data of the sound pressure at the analysis point P0 shown in FIG. 2 for each design parameter, and performs a fast Fourier transform (FFT) process or the like on the time response data. The acoustic characteristic data (data of the sound pressure frequency characteristic and the phase frequency characteristic at the analysis point P0) is calculated by performing the above.

このようにして算出された音響特性データは、その音響特性データに対応する設計パラメータと対をなして学習用データを構成し、設計パラメータ設定部310により生成される多数の設計パラメータに対応する多数の学習用データ110−1〜110−Nがデータベース120(図7)に格納される。 The acoustic characteristic data calculated in this way constitutes learning data by pairing with the design parameters corresponding to the acoustic characteristic data, and corresponds to a large number of design parameters generated by the design parameter setting unit 310. The training data 110-1 to 110-N of the above is stored in the database 120 (FIG. 7).

なお、条件設定部330によって設定される空間離散化幅Δx,Δy,Δz(以下、この空間離散化幅を有する格子を「グリッド」と称する。)については、解析空間(図2,図3に示した小型スピーカ200の境界260内の空間)において、領域に応じてグリッドの大きさ(粗さ)を適宜変えるのが好ましい。 The spatial discretization widths Δx, Δy, Δz (hereinafter, the grid having this spatial discretization width is referred to as a “grid”) set by the condition setting unit 330 is described in the analysis space (FIGS. 2 and 3). In the space (space within the boundary 260 of the small speaker 200 shown), it is preferable to appropriately change the size (roughness) of the grid according to the region.

図11は、条件設定部330によるグリッド設定の考え方の一例を説明する図である。この図11は、小型スピーカ200の断面図であり、図2の断面図に模式的なグリッドを重ね合わせたものである。 FIG. 11 is a diagram illustrating an example of the concept of grid setting by the condition setting unit 330. FIG. 11 is a cross-sectional view of the small speaker 200, in which a schematic grid is superimposed on the cross-sectional view of FIG.

図11を参照して、条件設定部330は、音源である振動板220の近傍の領域410、筐体210において領域410を除く領域420、及び筐体210外の領域430(音響伝搬空間250)で、グリッドの大きさ(粗さ)を変えて設定する。具体的には、振動板220を含む領域410におけるグリッドの大きさが相対的に小さく、筐体210外の領域430(音響伝搬空間250)におけるグリッドの大きさが相対的に大きくなるように、グリッドが設定される。 With reference to FIG. 11, the condition setting unit 330 includes a region 410 in the vicinity of the diaphragm 220 which is a sound source, a region 420 in the housing 210 excluding the region 410, and a region 430 outside the housing 210 (acoustic propagation space 250). Then, change the size (roughness) of the grid and set it. Specifically, the size of the grid in the region 410 including the diaphragm 220 is relatively small, and the size of the grid in the region 430 (acoustic propagation space 250) outside the housing 210 is relatively large. The grid is set.

このように、解析空間においてグリッドの大きさを一律にするのではなく、領域に応じてグリッドの大きさを適宜変えることにより、音響FDTD法によるコンピュータシミュレーションの精度を高めることができるとともに、計算時間を短縮することができる。 In this way, the accuracy of computer simulation by the acoustic FDTD method can be improved and the calculation time can be improved by appropriately changing the size of the grid according to the region, instead of making the size of the grid uniform in the analysis space. Can be shortened.

なお、図11に示したような3段階のグリッド分割は一例であって、グリッドの分割方法は、図11に示されるものに限定されるものではない。また、図11に示した例は、3段階のグリッド間の相対的な大きさの違いを説明するために模式的に示したものであって、各グリッドの絶対的な大きさを示すものではない。 The three-stage grid division as shown in FIG. 11 is an example, and the grid division method is not limited to that shown in FIG. Further, the example shown in FIG. 11 is schematically shown for explaining the difference in relative size between the three stages of grids, and does not show the absolute size of each grid. do not have.

再び図7を参照して、DNN学習部130は、CPU、処理プログラムを記憶するROM、データを一時的に記憶するRAM等を含んで構成される。そして、DNN学習部130は、データベース120に格納された多数の学習用データ110−1〜110−Nを用いてDNNの学習を行ない、音響特性データから小型スピーカ200の設計パラメータを算出するようにコンピュータ(図1の演算装置20)を機能させるための学習済みモデルを生成する。すなわち、この学習済みモデルは、実際には、データベース120に格納された学習用データ110−1〜110−Nを用いて音響特性データと設計パラメータとの因果関係が学習されたDNNそのものである。 With reference to FIG. 7 again, the DNN learning unit 130 includes a CPU, a ROM for storing a processing program, a RAM for temporarily storing data, and the like. Then, the DNN learning unit 130 learns the DNN using a large number of learning data 110-1 to 110-N stored in the database 120, and calculates the design parameters of the small speaker 200 from the acoustic characteristic data. A trained model for operating a computer (computing device 20 in FIG. 1) is generated. That is, this trained model is actually the DNN itself in which the causal relationship between the acoustic characteristic data and the design parameters is learned using the training data 110-1 to 110-N stored in the database 120.

図12は、DNN学習部130により生成される学習済みモデルに用いられるDNNの構成図である。図12を参照して、DNNは、入力層と、出力層と、複数の中間層とを含んで構成される。DNNは、複数の中間層を備えることによって、従来のニューラルネットワーク(NN)に比べて高い汎化性能を有し、小型スピーカ200のような複雑な形状(設計パラメータ)を有するモデルには特に有用である。 FIG. 12 is a configuration diagram of a DNN used in the trained model generated by the DNN learning unit 130. With reference to FIG. 12, the DNN includes an input layer, an output layer, and a plurality of intermediate layers. By providing a plurality of intermediate layers, the DNN has higher generalization performance than the conventional neural network (NN), and is particularly useful for a model having a complicated shape (design parameter) such as a small speaker 200. Is.

入力層に入力されたデータは、第1の中間層との結合部分において重み付けされて第1の中間層に入力される。第1の中間層の各ユニット(ニューロン)には、活性化関数が設けられており、各ユニットの活性化関数の出力は、次の第2の中間層(図示せず)との結合部分において重み付けされて第2の中間層に入力される。このようなデータの伝播が複数の中間層において順次行なわれ、最終的に出力層にデータが出力される。 The data input to the input layer is weighted at the joint portion with the first intermediate layer and input to the first intermediate layer. Each unit (neuron) in the first intermediate layer is provided with an activation function, and the output of the activation function of each unit is in the connection portion with the next second intermediate layer (not shown). It is weighted and input to the second intermediate layer. Such data propagation is sequentially performed in the plurality of intermediate layers, and finally the data is output to the output layer.

出力層の出力データは、教師データと比較される。そして、出力層の出力データと教師データとの誤差に基づいて入力層までの勾配が算出され、算出された勾配を用いて各ユニット間の重み付け(パラメータ)が更新される(誤差逆伝播法)。このような処理を多数の学習用データで行なうことにより、入出力間の因果関係がモデル化される(学習済みモデルの生成)。 The output data of the output layer is compared with the teacher data. Then, the gradient to the input layer is calculated based on the error between the output data of the output layer and the teacher data, and the weighting (parameter) between each unit is updated using the calculated gradient (error backpropagation method). .. By performing such processing with a large amount of training data, the causal relationship between input and output is modeled (generation of trained model).

図13は、本実施の形態におけるDNNの入出力を示した図である。図13を参照して、DNNの入力層には、データベース120(図7)から受ける学習用データの音響特性データが入力される。この実施の形態では、音響特性データとして、音圧の絶対値(Pa)の周波数特性のプロフィールと、位相の周波数特性のプロフィールとが用いられる。一例として、音圧及び位相の周波数特性の各々のプロフィールを所定の周波数幅Δfで刻んで離散化し、DNNの入力層に与えられる。 FIG. 13 is a diagram showing input / output of DNN in the present embodiment. With reference to FIG. 13, the acoustic characteristic data of the learning data received from the database 120 (FIG. 7) is input to the input layer of the DNN. In this embodiment, the profile of the frequency characteristic of the absolute value (Pa) of the sound pressure and the profile of the frequency characteristic of the phase are used as the acoustic characteristic data. As an example, each profile of the frequency characteristics of sound pressure and phase is discretized in a predetermined frequency width Δf and given to the input layer of the DNN.

出力層には、設計パラメータに対応するデータが出力される。すなわち、出力層のユニット数は、設計変数ベクトルX=[LCx,LHx,LHy,LHz,Hymin,Hzmin]の変数の数と同じである(この実施の形態では6個)。そして、出力層の出力データは、入力層に入力された学習用データの音響特性データと対をなす設計パラメータ(教師データ)と比較され、その比較誤差に基づいて各ユニット間の重み付け(パラメータ)が更新される。DNN学習部130は、データベース120(図7)に格納される多数の学習用データに対して上記の処理を行なうことにより、学習済みモデルを生成する。 Data corresponding to the design parameters is output to the output layer. That is, the number of units in the output layer is the same as the number of variables of the design variable vector X = [LCx, LHx, LHy, LHz, Hymin, Hzmin] (6 in this embodiment). Then, the output data of the output layer is compared with the design parameter (teacher data) paired with the acoustic characteristic data of the learning data input to the input layer, and the weighting (parameter) between each unit is based on the comparison error. Is updated. The DNN learning unit 130 generates a trained model by performing the above processing on a large number of learning data stored in the database 120 (FIG. 7).

再び図7を参照して、学習条件設定部140は、DNN学習部130において学習済みモデルを生成する際の学習条件を設定する。たとえば、学習条件設定部140は、DNNの中間層の数(この実施の形態では4層で構成)、並びに入力層、各中間層、及び出力層の各々のユニット数(ニューロン数)を設定する。また、たとえば、学習条件設定部140は、各学習用データによる学習を終了する条件を設定する。具体的には、学習用データによるDNNの過学習を防ぐために、学習用データとは別のテストデータ(設計パラメータ及び音響特性データ)が準備される。そして、各学習用データによる学習において、学習が行なわれる毎にDNNの入力層にテストデータ(音響特性データ)を与え、出力層の出力データとテストデータ(設計パラメータ)との誤差(テスト誤差)の最小値が所定回数連続して更新されなければ当該学習用データの学習が終了するように(早期終了)、各学習用データによる学習を終了する条件として上記の所定回数が設定される。 With reference to FIG. 7 again, the learning condition setting unit 140 sets the learning conditions for generating the trained model in the DNN learning unit 130. For example, the learning condition setting unit 140 sets the number of intermediate layers of the DNN (composed of four layers in this embodiment) and the number of units (number of neurons) of each of the input layer, each intermediate layer, and the output layer. .. Further, for example, the learning condition setting unit 140 sets a condition for ending learning by each learning data. Specifically, in order to prevent overfitting of DNN by the training data, test data (design parameters and acoustic characteristic data) different from the training data are prepared. Then, in learning with each learning data, test data (acoustic characteristic data) is given to the input layer of the DNN each time learning is performed, and an error (test error) between the output data of the output layer and the test data (design parameter). The above-mentioned predetermined number of times is set as a condition for ending the learning by each learning data so that the learning of the learning data is completed if the minimum value of is not updated continuously a predetermined number of times (early end).

図14は、図1に示した演算装置20において用いられる学習済みモデルの生成手順を説明するフローチャートである。図14とともに図7,図8を参照して、最初に、ステップS10からS40までの処理において、多数の学習用データ110−1〜110−Nが生成される。すなわち、生成システム300において、設計パラメータ設定部310により、学習用データを生成するための設計パラメータが設定される(ステップS10)。なお、設計パラメータ設定部310は、図9で説明したように設計パラメータの対称性を考慮し、かつ、図10で説明したようにランダムに、設計パラメータを設定する。 FIG. 14 is a flowchart illustrating a procedure for generating a trained model used in the arithmetic unit 20 shown in FIG. With reference to FIGS. 7 and 8 together with FIG. 14, first, a large number of learning data 110-1 to 110-N are generated in the processes of steps S10 to S40. That is, in the generation system 300, the design parameter setting unit 310 sets the design parameters for generating the learning data (step S10). The design parameter setting unit 310 considers the symmetry of the design parameters as described in FIG. 9, and randomly sets the design parameters as described in FIG.

次いで、音響FDTD演算部320において、上述の音響FDTD法を用いたコンピュータシミュレーションにより、ステップS10において設定された設計パラメータに対応する音響特性データが算出される(ステップS20)。そして、ステップS10において設定された設計パラメータと、当該設計パラメータに基づきステップS20において算出された音響特性データとを対とする学習用データが、学習システム100のデータベース120に格納される(ステップS30)。 Next, the acoustic FDTD calculation unit 320 calculates acoustic characteristic data corresponding to the design parameters set in step S10 by computer simulation using the above-mentioned acoustic FDTD method (step S20). Then, the learning data for pairing the design parameter set in step S10 and the acoustic characteristic data calculated in step S20 based on the design parameter is stored in the database 120 of the learning system 100 (step S30). ..

次いで、音響FDTD演算部320において、データベース120に格納した学習用データの数が所定数N以上であるか否かが判定される(ステップS40)。なお、この所定数Nは、たとえば数千レベルの値である。ステップS40において学習用データ数が所定数Nよりも少ないと判定されると(ステップS40においてNO)、ステップS10へ処理が戻され、設計パラメータの対称性を考慮し、かつ、ランダムに、次の設計パラメータが設定される。 Next, the acoustic FDTD calculation unit 320 determines whether or not the number of learning data stored in the database 120 is a predetermined number N or more (step S40). The predetermined number N is, for example, a value of several thousand levels. When it is determined in step S40 that the number of training data is less than the predetermined number N (NO in step S40), the process is returned to step S10, considering the symmetry of the design parameters, and randomly, the next Design parameters are set.

ステップS40において学習用データ数が所定数N以上であると判定されると(ステップS40においてYES)、続いて、ステップS50からS70までの処理において、学習済みモデルが生成される。すなわち、学習システム100において、学習条件設定部140により、DNN学習部130におけるDNNの学習条件が設定される(ステップS50)。上述のように、学習条件は、たとえば、DNNの構成(層数やユニット数(ニューロンの数)等)や、各学習用データによる学習を終了する条件等が設定される。 When it is determined in step S40 that the number of training data is a predetermined number N or more (YES in step S40), a trained model is subsequently generated in the processes from steps S50 to S70. That is, in the learning system 100, the learning condition setting unit 140 sets the DNN learning conditions in the DNN learning unit 130 (step S50). As described above, as the learning conditions, for example, the configuration of the DNN (the number of layers, the number of units (the number of neurons), etc.), the conditions for terminating the learning with each learning data, and the like are set.

次いで、DNN学習部130において、データベース120に格納された多数の学習用データ110−1〜110−Nを用いて、音響特性データと設計パラメータとの間の因果関係がDNNにより学習される(ステップS60)。具体的には、DNNの入力層に学習用データの音響特性データが与えられ、当該音響特性データと対をなす学習用データの設計パラメータを教師データとして、多数の学習用データ110−1〜110−Nを用いてDNNの学習が行なわれる。 Next, in the DNN learning unit 130, the causal relationship between the acoustic characteristic data and the design parameters is learned by the DNN using a large number of learning data 110-1 to 110-N stored in the database 120 (step). S60). Specifically, the acoustic characteristic data of the learning data is given to the input layer of the DNN, and a large number of learning data 110-1 to 110 are used as the teacher data using the design parameters of the learning data paired with the acoustic characteristic data. DNN learning is performed using −N.

そして、学習用データ110−1〜110−NによるDNNの学習が終了すると、DNN学習部130から学習済みモデル(実際には、DNNの構造と、各ユニット間の重み付けを示すパラメータ)が出力される(ステップS70)。 Then, when the learning of the DNN by the learning data 110-1 to 110-N is completed, the learned model (actually, the structure of the DNN and the parameter indicating the weighting between the units) is output from the DNN learning unit 130. (Step S70).

図15は、図1に示した演算装置20において実行される処理の手順を説明するフローチャートである。図15を参照して、演算装置20は、ユーザにより入力装置10から入力された音響特性データ(図5及び図6に示されるような音圧(Pa)及び位相の各周波数特性のプロフィール)を入力装置10から取得する(ステップS110)。なお、この音響特性データは、学習システム100のデータベース120に格納される学習用データの音響特性データではなく、小型スピーカの設計をこれから行なうユーザが入力装置10から入力した所望の音響特性データである。 FIG. 15 is a flowchart illustrating a procedure of processing executed by the arithmetic unit 20 shown in FIG. With reference to FIG. 15, the arithmetic apparatus 20 obtains acoustic characteristic data (profiles of sound pressure (Pa) and phase frequency characteristics as shown in FIGS. 5 and 6) input from the input device 10 by the user. Obtained from the input device 10 (step S110). The acoustic characteristic data is not the acoustic characteristic data of the learning data stored in the database 120 of the learning system 100, but the desired acoustic characteristic data input from the input device 10 by the user who is about to design the small speaker. ..

次いで、演算装置20は、ステップS10において取得された音響特性データを、図14に示した手順に従って予め生成された学習済みモデル(学習後DNN)に入力し、学習済みモデルを用いて設計パラメータ(図2,図3に示したLCx,LHx,LHy,LHz,Hymin,Hzmin)を算出する(ステップS120)。 Next, the arithmetic unit 20 inputs the acoustic characteristic data acquired in step S10 into the trained model (post-learning DNN) generated in advance according to the procedure shown in FIG. 14, and uses the trained model to design parameters ( LCx, LHx, LHy, LHz, Hymin, Hzmin) shown in FIGS. 2 and 3) are calculated (step S120).

そして、演算装置20は、算出された設計パラメータであるLCx,LHx,LHy,LHz,Hymin,Hzminを表示装置30へ出力する(ステップS130)。なお、表示装置30においては、演算装置20から受ける設計パラメータの数値をそのまま表示してもよいし、そのような設計パラメータを有する小型スピーカ200の構造を図示してもよい。 Then, the arithmetic unit 20 outputs the calculated design parameters LCx, LHx, LHy, LHz, Hymin, and Hzmin to the display device 30 (step S130). In the display device 30, the numerical value of the design parameter received from the arithmetic unit 20 may be displayed as it is, or the structure of the small speaker 200 having such a design parameter may be illustrated.

以上のように、この実施の形態においては、小型スピーカの音響特性データに基づいて小型スピーカの設計パラメータを出力するようにコンピュータ(演算装置20)を機能させるための学習済みモデルを用いて、入力装置10により入力された所望の音響特性データから設計パラメータが算出される。学習済みモデルは、音響特性データから設計パラメータを算出するように学習されたDNNによって構成される。これにより、技術者によるカット&トライ設計に頼ることなく小型スピーカの設計を行なうことができる。 As described above, in this embodiment, input is performed using a trained model for operating the computer (computing device 20) so as to output the design parameters of the small speaker based on the acoustic characteristic data of the small speaker. Design parameters are calculated from the desired acoustic characteristic data input by the device 10. The trained model is composed of DNNs trained to calculate design parameters from acoustic characteristic data. This makes it possible to design a small speaker without relying on a cut-and-try design by an engineer.

そして、各々が設計パラメータと当該設計パラメータを有する小型スピーカの音響特性データとを含む、DNNの学習に用いられる複数の学習用データ110−1〜110−Nは、音響FDTD法を用いたコンピュータシミュレーションによって生成される。これにより、複数の学習用データを用意するために実際に様々な音響構造を作成して実測を行なうことなく、DNNの学習に用いられる多数の学習用データ110−1〜110−Nを用意することができる。 A plurality of learning data 110-1 to 110-N used for DNN learning, each including a design parameter and acoustic characteristic data of a small speaker having the design parameter, are computer simulations using the acoustic FDTD method. Generated by. As a result, a large number of learning data 110-1 to 110-N used for learning DNN are prepared without actually creating various acoustic structures and actually measuring them in order to prepare a plurality of learning data. be able to.

このように、本開示の小型スピーカの設計支援装置1によれば、技術者によるカット&トライ設計に頼ることなく、また、DNNの学習用データを用意するために実際に様々な音響構造を作成して実測を行なう必要もないので、小型スピーカの設計に要する時間及びコストを大幅に削減することができる。 As described above, according to the design support device 1 for the small speaker of the present disclosure, various acoustic structures are actually created in order to prepare the learning data of DNN without relying on the cut and try design by the engineer. Since it is not necessary to carry out actual measurement, the time and cost required for designing a small speaker can be significantly reduced.

また、この実施の形態においては、音響FDTD法を用いたコンピュータシミュレーションにおけるグリッドの大きさ(粗さ)を領域によって適宜変えている。具体的には、音源である振動板220の近傍の領域におけるグリッドを小さく、筐体外の音響伝搬空間250におけるグリッドを大きくしている(図11)。これにより、グリッドの大きさが一律である場合に比べて、音響FDTD法の精度を高めることができるとともに、計算時間を短縮することができる。 Further, in this embodiment, the size (roughness) of the grid in the computer simulation using the acoustic FDTD method is appropriately changed depending on the region. Specifically, the grid in the region near the diaphragm 220, which is a sound source, is made small, and the grid in the acoustic propagation space 250 outside the housing is made large (FIG. 11). As a result, the accuracy of the acoustic FDTD method can be improved and the calculation time can be shortened as compared with the case where the grid size is uniform.

また、この実施の形態によれば、設計対象(小型スピーカ200)の形状の対称性を考慮して解析空間を限定することによって、設計パラメータ毎に音響特性データを算出する音響FDTD演算部320の演算負荷及び演算時間を低減することができる。 Further, according to this embodiment, the acoustic FDTD calculation unit 320 calculates acoustic characteristic data for each design parameter by limiting the analysis space in consideration of the symmetry of the shape of the design target (small speaker 200). The calculation load and calculation time can be reduced.

さらに、この実施の形態によれば、ランダムに変化させた設計パラメータに基づく多数の学習用データを用いてDNNの学習が行なわれるので、DNNにより学習済みモデルを生成する際に学習が局所解に陥ることを回避して学習済みモデルの精度を高めることができる。 Further, according to this embodiment, since the DNN is trained using a large number of training data based on the randomly changed design parameters, the training becomes a local solution when the trained model is generated by the DNN. It is possible to avoid falling and improve the accuracy of the trained model.

なお、学習済みモデルの生成に用いられる学習用データ110−1〜110−Nについて、上記の実施の形態においては、DNN学習部130による学習時にDNNの入力層に与えられる音響特性データのうちの音圧周波数特性は、図16に示されるような音圧の絶対値(Pa)の周波数特性のプロフィールとしたが、音圧の絶対値(Pa)の周波数特性のプロフィールに代えて、図17に示されるような音圧レベル(dB)の周波数特性のプロフィールをDNNの入力層に与えてもよい。音圧レベル(dB)の周波数特性のプロフィールは、音圧の絶対値(Pa)の周波数特性のプロフィールよりも、共振周波数fr以外の周波数における特徴がよく表れており、周波数特性のプロフィールを重視する設計を行なう場合には、学習済みモデルの精度を高めることができる。 Regarding the training data 110-1 to 110-N used for generating the trained model, in the above embodiment, of the acoustic characteristic data given to the input layer of the DNN at the time of training by the DNN learning unit 130. The sound pressure frequency characteristic is a profile of the frequency characteristic of the absolute value (Pa) of the sound pressure as shown in FIG. 16, but instead of the profile of the frequency characteristic of the absolute value (Pa) of the sound pressure, FIG. 17 shows. A profile of the frequency characteristics of the sound pressure level (dB) as shown may be given to the input layer of the DNN. The profile of the frequency characteristic of the sound pressure level (dB) shows the characteristics at frequencies other than the resonance frequency fr better than the profile of the frequency characteristic of the absolute value (Pa) of the sound pressure, and emphasizes the profile of the frequency characteristic. When designing, the accuracy of the trained model can be improved.

また、学習済みモデルの生成の際に、音圧の絶対値(Pa)の周波数特性のプロフィールをDNNの入力層に与えて学習するか、それとも音圧レベル(dB)の周波数特性のプロフィールをDNNの入力層に与えて学習するかを、ユーザが選択可能としてもよい。たとえば、再び図7を参照して、ユーザの選択結果を示す情報を学習条件設定部140からDNN学習部130に設定するようにし、音圧レベル(dB)の周波数特性のプロフィールを用いて学習することがユーザにより選択された場合には、データベース120から受ける音圧の絶対値(Pa)の周波数特性データをデシベル(dB)に換算してDNNの入力層に与えればよい。 Further, when the trained model is generated, the profile of the frequency characteristic of the absolute value (Pa) of the sound pressure is given to the input layer of the DNN for learning, or the profile of the frequency characteristic of the sound pressure level (dB) is given to the DNN. The user may be able to select whether to give the sound to the input layer for learning. For example, referring to FIG. 7 again, information indicating the user's selection result is set in the learning condition setting unit 140 to the DNN learning unit 130, and learning is performed using the profile of the frequency characteristic of the sound pressure level (dB). When this is selected by the user, the frequency characteristic data of the absolute value (Pa) of the sound pressure received from the database 120 may be converted into decibels (dB) and given to the input layer of the DNN.

なお、音圧レベル(dB)の周波数特性のプロフィールをDNNの入力層に与えて学習済みモデルが生成された場合には、設計支援装置1の演算装置20において学習済みモデルに与える所望の音響特性データのうちの音圧周波数特性についても、音圧レベル(dB)の周波数特性のプロフィールとする必要がある。 When a trained model is generated by giving a profile of the frequency characteristic of the sound pressure level (dB) to the input layer of the DNN, the desired acoustic characteristic given to the trained model by the arithmetic unit 20 of the design support device 1. The sound pressure frequency characteristic in the data also needs to be a profile of the frequency characteristic of the sound pressure level (dB).

また、上記の実施の形態においては、音響特性データは、音圧の周波数特性を示すデータと、位相の周波数特性を示すデータとによって構成されるものとしたが、音圧の周波数特性を示すデータのみとしてもよい。なお、上述の実施の形態のように、音圧の周波数特性だけでなく位相の周波数特性を併用することによって、DNNの学習精度、及び学習済みモデルを用いた小型スピーカの設計精度の向上が期待できる。 Further, in the above embodiment, the acoustic characteristic data is composed of data indicating the frequency characteristic of the sound pressure and data indicating the frequency characteristic of the phase, but data indicating the frequency characteristic of the sound pressure. May be only. By using not only the frequency characteristics of sound pressure but also the frequency characteristics of phase as in the above-described embodiment, it is expected that the learning accuracy of DNN and the design accuracy of a small speaker using the trained model will be improved. can.

今回開示された実施の形態は、すべての点で例示であって制限的なものではないと考えられるべきである。本開示の範囲は、上記した実施の形態の説明ではなくて特許請求の範囲によって示され、特許請求の範囲と均等の意味及び範囲内でのすべての変更が含まれることが意図される。 The embodiments disclosed this time should be considered to be exemplary in all respects and not restrictive. The scope of the present disclosure is shown by the claims rather than the description of the embodiments described above, and is intended to include all modifications within the meaning and scope equivalent to the claims.

1 設計支援装置、10 入力装置、20 演算装置、30 表示装置、100 学習システム、110−1,110−2,・・・,110−N 学習用データ、120 データベース、130 DNN学習部、140 学習条件設定部、200 小型スピーカ、210 筐体、220 振動板、230 音孔、240 気室、250 音響伝搬空間、260 境界、270,410〜430 領域、300 生成システム、310 設計パラメータ設定部、320 音響FDTD演算部、330 条件設定部。 1 Design support device, 10 input device, 20 arithmetic device, 30 display device, 100 learning system, 110-1, 110-2, ..., 110-N learning data, 120 database, 130 DNN learning unit, 140 learning Condition setting unit, 200 small speaker, 210 housing, 220 diaphragm, 230 sound hole, 240 air chamber, 250 acoustic propagation space, 260 boundary, 270, 410-430 area, 300 generation system, 310 design parameter setting unit, 320 Acoustic FDTD calculation unit, 330 condition setting unit.

Claims (8)

小型スピーカの設計支援装置であって、
前記小型スピーカの音響特性を示す音響特性データを利用者が入力するための入力装置と、
前記音響特性データに基づいて前記小型スピーカの設計パラメータを出力するようにコンピュータを機能させるための学習済みモデルを用いて、前記入力装置により入力された音響特性データから前記設計パラメータを算出するように構成された演算装置とを備え、
前記学習済みモデルは、音響FDTD法を用いたコンピュータシミュレーションによって生成される、各々が前記設計パラメータと当該設計パラメータを有する小型スピーカの前記音響特性データとを含む複数の学習用データを用いて、前記音響特性データから前記設計パラメータを算出するように学習されたディープニューラルネットワークによって構成される、小型スピーカの設計支援装置。
It is a design support device for small speakers.
An input device for the user to input acoustic characteristic data indicating the acoustic characteristics of the small speaker, and
The design parameters are calculated from the acoustic characteristic data input by the input device using a trained model for operating the computer to output the design parameters of the small speaker based on the acoustic characteristic data. Equipped with a configured arithmetic unit
The trained model uses a plurality of training data, each including the design parameter and the acoustic characteristic data of a small speaker having the design parameter, generated by computer simulation using the acoustic FDTD method. A design support device for a small speaker configured by a deep neural network learned to calculate the design parameters from acoustic characteristic data.
前記小型スピーカは、音源となる振動板と、前記小型スピーカの筐体の内部に形成される筐体内空間とを含んで構成され、
音響FDTD法を用いた前記コンピュータシミュレーションにおける分割グリッドの大きさは、前記振動板、前記筐体内空間、前記筐体の外部における音響伝搬空間の順に小さい、請求項1に記載の小型スピーカの設計支援装置。
The small speaker is configured to include a diaphragm serving as a sound source and a space inside the housing formed inside the housing of the small speaker.
The design support for a small speaker according to claim 1, wherein the size of the divided grid in the computer simulation using the acoustic FDTD method is smaller in the order of the diaphragm, the space inside the housing, and the acoustic propagation space outside the housing. Device.
前記複数の学習用データの各々に含まれる前記設計パラメータは、予め定められた範囲内においてランダムに生成される、請求項1又は請求項2に記載の小型スピーカの設計支援装置。 The design support device for a small speaker according to claim 1 or 2, wherein the design parameters included in each of the plurality of learning data are randomly generated within a predetermined range. 前記音響特性データは、所定の解析点における音圧の周波数特性を示すデータを含む、請求項1から請求項3のいずれか1項に記載の小型スピーカの設計支援装置。 The small speaker design support device according to any one of claims 1 to 3, wherein the acoustic characteristic data includes data indicating a frequency characteristic of sound pressure at a predetermined analysis point. 前記複数の学習用データとして用いられる前記音響特性データには、利用者の要求に従って、音圧の絶対値の周波数特性を示すデータ及び音圧レベルの周波数特性を示すデータのいずれかが選択される、請求項4に記載の小型スピーカの設計支援装置。 For the acoustic characteristic data used as the plurality of learning data, either data indicating the frequency characteristic of the absolute value of the sound pressure or data indicating the frequency characteristic of the sound pressure level is selected according to the request of the user. , The design support device for a small speaker according to claim 4. 前記音響特性データは、前記解析点における位相の周波数特性を示すデータをさらに含む、請求項4又は請求項5に記載の小型スピーカの設計支援装置。 The design support device for a small speaker according to claim 4 or 5, wherein the acoustic characteristic data further includes data indicating a frequency characteristic of the phase at the analysis point. 前記設計パラメータは、前記小型スピーカに形成される音孔の寸法及び位置並びに気室の寸法を含む、請求項1から請求項6のいずれか1項に記載の小型スピーカの設計支援装置。 The design support device for a small speaker according to any one of claims 1 to 6, wherein the design parameter includes the size and position of a sound hole formed in the small speaker and the size of an air chamber. スピーカの設計支援方法であって、
前記スピーカの音響特性を示す音響特性データに基づいて前記スピーカの設計パラメータを出力するようにコンピュータを機能させるための学習済みモデルを生成するステップと、
前記音響特性データを利用者が入力するステップと、
前記学習済みモデルを用いて、入力された前記音響特性データから前記設計パラメータを算出するステップとを含み、
前記学習済みモデルを生成するステップは、
音響FDTD法を用いたコンピュータシミュレーションによって、各々が前記設計パラメータと当該設計パラメータを有するスピーカの前記音響特性データとを含む複数の学習用データを生成するステップと、
前記複数の学習用データを用いて、前記音響特性データから前記設計パラメータを算出するようにディープニューラルネットワークを学習させることによって、前記学習済みモデルを生成するステップとを含む、スピーカの設計支援方法。
It is a speaker design support method.
A step of generating a trained model for the computer to function to output the design parameters of the speaker based on the acoustic characteristic data indicating the acoustic characteristics of the speaker.
The step of inputting the acoustic characteristic data by the user and
Including the step of calculating the design parameter from the input acoustic characteristic data using the trained model.
The step of generating the trained model is
A step of generating a plurality of learning data including the design parameter and the acoustic characteristic data of the speaker having the design parameter, respectively, by computer simulation using the acoustic FDTD method.
A speaker design support method including a step of generating the trained model by training a deep neural network so as to calculate the design parameter from the acoustic characteristic data using the plurality of learning data.
JP2017102449A 2017-05-24 2017-05-24 Small speaker design support device and speaker design support method Active JP6910641B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2017102449A JP6910641B2 (en) 2017-05-24 2017-05-24 Small speaker design support device and speaker design support method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2017102449A JP6910641B2 (en) 2017-05-24 2017-05-24 Small speaker design support device and speaker design support method

Publications (2)

Publication Number Publication Date
JP2018197965A JP2018197965A (en) 2018-12-13
JP6910641B2 true JP6910641B2 (en) 2021-07-28

Family

ID=64663184

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017102449A Active JP6910641B2 (en) 2017-05-24 2017-05-24 Small speaker design support device and speaker design support method

Country Status (1)

Country Link
JP (1) JP6910641B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI694252B (en) * 2019-01-03 2020-05-21 逢甲大學 A predicting method of porous material and prediction system thereof
CN110321619B (en) * 2019-06-26 2020-09-15 深圳技术大学 Parameterized custom model generation method based on sound data

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4200989B2 (en) * 2005-07-19 2008-12-24 ヤマハ株式会社 Acoustic design support device and acoustic design support program
CN107996028A (en) * 2015-03-10 2018-05-04 Ossic公司 Calibrate listening devices

Also Published As

Publication number Publication date
JP2018197965A (en) 2018-12-13

Similar Documents

Publication Publication Date Title
US9041662B2 (en) Touch sensitive device
Wang et al. Design of a Bézier-profile horn for high displacement amplification
Senturia et al. Simulating the behavior of MEMS devices: computational methods and needs
JP2019530849A (en) Contribution of flow-induced noise sources
CN111159928A (en) A method and system for calculating transformer noise based on multi-line sound source model
Aucejo et al. Experimental simulation of turbulent boundary layer induced vibrations by using a synthetic array
JP6910641B2 (en) Small speaker design support device and speaker design support method
Gatica et al. Classifying acoustic cavitation with machine learning trained on multiple physical models
JP4973296B2 (en) Vibration countermeasure support system, vibration countermeasure support method, vibration countermeasure support program
CN104239676B (en) In wide frequency range to method and apparatus of the system sound emission with propagating modeling
Wang et al. Investigation and analysis of the influence of excitation signal on radiation characteristics of capacitive micromachined ultrasonic transducer
Langlois et al. Partition of unity finite element method for 2d vibro-acoustic modeling
Yu et al. A novel version of hierarchical genetic algorithm and its application for hyperparameters optimization in CNN models for structural delamination identification
JP7113168B2 (en) Sound pressure calculation method and sound pressure calculation program
Anderson Hybrid frequency-time analysis and numerical methods for time-dependent wave propagation
An et al. A new method for predicting the acoustic transfer function inside vehicles using physics-informed neural networks
CN113221409B (en) A two-dimensional acoustic wave numerical simulation method and device coupled with finite element and boundary element
Darus et al. Finite difference simulation of a flexible plate structure
JP2021018221A (en) Sound field analysis device, sound field analysis method, and program
JP7285513B2 (en) SOUND FIELD ANALYZER, SOUND FIELD ANALYSIS METHOD AND PROGRAM
MAJERCSIK On the semiactive control of a Stewart platform with flexible joints
JP2022112625A (en) Sound field analysis apparatus, sound analysis method, and program
JP2012014246A (en) Sound field analyzing method
JP7470248B2 (en) Tire external noise simulation method
Wagner et al. An acoustic finite-difference time-domain algorithm with isotropic dispersion

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200519

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20210517

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: 20210601

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210630

R150 Certificate of patent or registration of utility model

Ref document number: 6910641

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250