(57)【要約】【課題】アンテナの設計において搭載される実機を考慮して試作を行う場合には、過去の設計指針が適用できないため試作を繰り返すことになり、設計期間が長期化し、コストが増加するという問題があった。また試作の代わりに電磁界シミュレーションを用いた場合にも、実機を正確に考慮した電磁界シミュレーションを行うと長時間かかってしまうため現実的でなかった。【解決手段】あらかじめ用意したデータベースと設計パラメータと特性パラメータから得られる関係式を用いたアンテナの設計方法を用いることで、所望のアンテナ特性を短時間で得ることができる。【選択図】図１
2005158019.tif 000002
【特許請求の範囲】
【請求項１】
設計パラメータとしてアンテナを構成する材料と各部寸法を決定し、次に前記アンテナの設計パラメータの値をデータベースと比較し、前記比較工程で設計パラメータが一致するデータベースがある場合、そのデータベースの特性を取得し、前記比較工程で設計パラメータが一致するデータベースがない場合、データベース内のデータを内挿、または外挿して得られる関係式より予測された特性を取得するアンテナの設計方法。
【請求項２】
データベースは、属性としてアンテナの設計パラメータである材料とアンテナの各部寸法、及びこれらの設計パラメータで決定されるアンテナ特性が対になって蓄積されている請求項１に記載のアンテナの設計方法。
【請求項３】
データベース内のアンテナ設計パラメータに、アンテナ形状と、携帯電話の筐体内側に施されたメッキ形状と、アンテナとその周辺部分との距離と、アンテナを取り付ける回路基板の形状を選択する請求項１に記載のアンテナの設計方法。
【請求項４】
データベース内のアンテナ特性に、共振周波数、帯域幅、効率を選択する請求項１に記載のアンテナの設計方法。
【請求項５】
データベース内のデータは、電磁界シミュレーションにより作成される請求項１に記載のアンテナの設計方法。
【請求項６】
電磁界シミュレーションに用いるシミュレーションモデルが、実際のアンテナと携帯電話の筐体と、回路基板と、アンテナの周辺部品を、３次元レーザースキャナで取り込んだ電子データを基に作成される請求項５に記載のアンテナの設計方法。
【請求項７】
関係式は、アンテナ特性がアンテナ設計パラメータの線形和で表される請求項１に記載のアンテナの設計方法。
【請求項８】
関係式は、アンテナ特性とアンテナ設計パラメータからニューラルネットワークを用いて求められる請求項１に記載のアンテナの設計方法。
【請求項９】
関係式は、アンテナ特性とアンテナ設計パラメータから単一の補間多項式を用いて求められる請求項１に記載のアンテナの設計方法。
【請求項１０】
関係式は、アンテナ特性とアンテナ設計パラメータからBe' zier（ベジェ）曲線を用いて求められる請求項１に記載のアンテナの設計方法。
【請求項１１】
関係式は、アンテナ特性とアンテナ設計パラメータから区分多項式を用いて求められる請求項１に記載のアンテナの設計方法。
【請求項１２】
関係式は、アンテナ特性とアンテナ設計パラメータから有理Be' zier（ベジェ）曲線を用いて求められる請求項１に記載のアンテナの設計方法。
【請求項１３】
関係式は、アンテナ特性とアンテナ設計パラメータから三角関数による補間を用いて求められる請求項１に記載のアンテナの設計方法。
【請求項１４】
請求項１〜１３に記載のアンテナの設計方法を用いて、設計されることを特徴とするアンテナ。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、携帯電話等に用いられるアンテナの設計方法、及びこの設計方法を用いて設計されるアンテナに関するものである。
【背景技術】
【０００２】
近年の情報関連機器は小型化の傾向にあり、それに伴い各種電子部品にも小型化、低背化の波が押し寄せている。携帯電話等に搭載されるアンテナもその例外ではなく、小型化が要求されている。しかし、アンテナは一般的にサイズが小さくなると電磁波の放射効率が低下し、周辺部品に対する感度も大きくなる。よって実筐体やアンテナの周辺部品の影響も考慮に入れたアンテナの設計が必要となっている。従来のアンテナ設計は、実筐体やアンテナの周辺部品を用いて試作を繰り返すことで、アンテナ特性の最適化を行っている。一般に、携帯電話の筐体やアンテナが取り付けられる回路基板の形状、及びアンテナの周辺部品の形状や配置は多種多様（カスタマ、品種によって異なる）であり、過去の設計指針は適用できないため、実筐体やアンテナの周辺部品を用いた試作を何度も繰り返す必要があり、設計期間が長期化し、さらにコストが増加するといった課題があった。また、特許文献１に示すように、アンテナ設計方法として、電磁界シミュレーションを利用することもなされている。
【特許文献１】特開平１１−１６１６９０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
しかしながら、従来の電磁界シミュレーションでは、前記筐体や回路基板、アンテナの周辺部品を正確に考慮した場合には、シミュレーションモデルが複雑かつ大規模になるため、シミュレーションに数日を要していた。
【課題を解決するための手段】
【０００４】
この課題を解決するために本発明は、設計パラメータとしてアンテナを構成する材料と各部寸法を決定し、次に前記アンテナの設計パラメータの値をデータベースと比較し、前記比較工程で設計パラメータが一致するデータベースがある場合、そのデータベースの特性を取得し、前記比較工程で設計パラメータが一致するデータベースがない場合、データベース内のデータを内挿、または外挿して得られる関係式より予測された特性を取得することを特徴とするアンテナの設計方法およびこれを用いたアンテナであり、特性パラメータのデータベース、またはアンテナ設計パラメータと特性パラメータの関係式を利用することにより、電磁界シミュレーションを行わずに所望の特性パラメータの結果を得ることができるので、アンテナを構成する材料と各部寸法を決定してからアンテナの特性パラメータを取得するまでの時間が数秒程度に短縮されるという作用を有する。
【発明の効果】
【０００５】
本発明は、携帯電話等の実機を考慮してアンテナの設計を行うのに際し、電磁界シミュレーションを用いてあらかじめ用意したデータベースと、設計パラメータと特性パラメータから得られる関係式を用いて設計を行うことで、従来１回のシミュレーションに数日を要していたものが数秒程度で済むようになり、結果としてアンテナ設計時間を短縮し、アンテナ開発のコストダウンに大きく貢献することができる。
【発明を実施するための最良の形態】
【０００６】
（実施の形態１）
以下、本発明の実施の形態１について、本発明の特に請求項１〜７，１４の発明について、図面を参照しながら説明する。
【０００７】
まず本発明のアンテナを設計するにあたり、データベースを作成する工程、関係式を作成する工程の２つの工程を要する。前記の２つの工程を含めた本発明の実施の形態による携帯電話用アンテナ特性解析のフローチャートを図１に示す。
【０００８】
本発明の実施に当たって、まずアンテナを構成する設計パラメータである材料と各部寸法をセットする（ステップＳ１）。入力した設計パラメータの値をデータベースと比較する（ステップＳ２）。ここで、この比較工程で設計パラメータが一致するデータがある場合はそのデータベースの特性を取得し（ステップＳ３）、ステップ２の比較工程で設計パラメータが一致するデータパラメータがない場合、データベース内のデータを内挿、または外挿して得られる関係式より予測された特性を取得する（ステップＳ４）。ここで関係式で表される範囲外の入力があった場合には電磁界シミュレータを使用して数値計算を行い（ステップＳ５）、特性を取得する。この時に得られた計算結果は新たにデータベースに追加する。
【０００９】
本発明の実施に先立つシステムの構築例について説明する。図２は回路基板にモノポールアンテナを取り付けた時の正面図であり、図３は回路基板にモノポールアンテナを取り付けた時の側面図である。ここで図２において、材質が純銅であるアンテナ素子１の長手方向の長さをＸ１、アンテナ端部から回路基板３上に取り付けられたシールドケース２の端までの距離をＸ２、回路基板３の長手方向の長さをＸ３とする。
【００１０】
次に、データベースを作成するために、例としてＸ１の長さを７０ｍｍ、７５ｍｍ、８０ｍｍの３種類、Ｘ２の長さを５ｍｍ、７ｍｍ、９ｍｍの時の３種類、Ｘ３の長さを９５ｍｍ、１００ｍｍ、１０５ｍｍの３種類の計２７種類のデータについて電磁界シミュレータを用いて数値計算を行い、それぞれの時の共振周波数、ＶＳＷＲ（Voltage Standing Wave Ratio）＜３となる帯域幅、効率を求める。ここで実物のモデルから設計パラメータを得る場合には、３次元スキャナを利用することで、簡単に設計パラメータの値を得ることができる。ＶＳＷＲとは電圧定在波比のことであり、ＶＳＷＲが大きくなるほど、インピーダンス不整合によるリターンロスが大きくなることを意味する。ここでは、実用上のアンテナにおいて、電源側のインピーダンスとアンテナのインピーダンスの不整合に起因するロスとして許容できる範囲で、ＶＳＷＲ＜３とする。
【００１１】
また、インピーダンス特性を評価する範囲は、０．５［ＧＨｚ］〜１．５［ＧＨｚ］とし、放射効率を評価する周波数を１［ＧＨｚ］とする。ここで、共振が１つの場合、および２つの場合について説明する。ここで図４は０．５［ＧＨｚ］〜１．５［ＧＨｚ］の間に１つの共振を有する時のＶＳＷＲの周波数特性図、図５は０．５［ＧＨｚ］〜１．５［ＧＨｚ］の間に２つの共振を有する時のＶＳＷＲの周波数特性図である。本実施の形態に示すように、アンテナの構造によっては、指定した範囲内に図５のように複数の共振を持つ場合もある。図４に示すように、範囲内に１つの共振しか存在しない場合の共振における共振周波数をＹ１、ＶＳＷＲ＜３となる帯域幅をＹ２、１［ＧＨｚ］における効率をＹ３とする。また図５に示すように、範囲内に複数の共振が存在する場合の、１つ目の共振における共振周波数をＹ１、１つ目の共振におけるＶＳＷＲ＜３となる帯域幅をＹ２、１［ＧＨｚ］における放射効率をＹ３（図６）、２つ目の共振における共振周波数をＹ４、２つ目の共振におけるＶＳＷＲ＜３となる帯域幅をＹ５とする。
【００１２】
ここで、本実施の形態の例において、設計パラメータによって特性パラメータがどのように変化するかについて簡単に説明する。基板に取り付けられたモノポールアンテナは、基板もアンテナの一部として働かせることで、アンテナ長がλ／４（λ：波長）の時に共振することが知られている。ここで、周波数と波長は反比例の関係にあるため、アンテナ長Ｘ１が長くなるにつれて共振周波数Ｙ１は低下するという結果が得られる。また、アンテナに金属導体が近づいた場合、すなわちＸ２の長さが短くなると、アンテナの特性インピーダンスが変化するため、共振周波数Ｙ１、帯域幅Ｙ２、放射効率Ｙ３が変化する。また前述したように、基板に取り付けられたモノポールアンテナでは、基板もアンテナの一部として働くため、理想的な共振長から基板の長さがずれた場合には、共振周波数Ｙ１、帯域幅Ｙ２、放射効率Ｙ３の変化がみられる。
【００１３】
以上のような設計パラメータと数値計算により得られた特性パラメータを用いて、図７に示すようなデータベースを作成する。ここで図７より設計パラメータが３変数、特性パラメータが３変数であるので、設計パラメータ３変数により１つの特性パラメータが決定されるテーブルが３つ作成されることになる。また、インピーダンス特性を評価する範囲内において２つの共振を有する場合には、設計パラメータが３変数、特性パラメータが５変数となり、設計パラメータ３変数により１つの特性パラメータが決定されるテーブルが５つ作成されることになる。
【００１４】
すなわち設計パラメータがＭ個の変数、特性パラメータがＮ個の変数とした場合には、Ｍ個の変数により１つの特性が決定されるテーブルがＮ個作成されることになる。
【００１５】
次にこのデータベースのデータから関係式を導く。
【００１６】
ここで設計パラメータＸと、特性パラメータＹとの間に線形和の関係が成り立つとすると、Ｙ１，Ｙ２，Ｙ３はＸ１，Ｘ２，Ｘ３を用いて（数１）〜（数３）のように表すことができる。
【００１７】
【数１】
2005158019.tif 000003
【００１８】
【数２】
2005158019.tif 000004
【００１９】
【数３】
2005158019.tif 000005
【００２０】
ここでｋ11〜ｋ33は未知変数である。
【００２１】
データベースより（Ｘ２，Ｙ２）＝［（3、950），（5、990），（10、1000），（15、1005），（30、1007）］が得られるとする。ここで求める関数ｆ（Ｘ２）を（数４）のように設定する。
【００２２】
【数４】
2005158019.tif 000006
【００２３】
全ての点において最適な関係式を導出するために、（数５）に示すような２乗誤差の総和が最小となるようにｋ22とＡ2を決定する。上記５点を（数５）に代入することによりｋ22＝１．４８、Ａ2＝９７１．７２が得られる。
【００２４】
【数５】
2005158019.tif 000007
【００２５】
本関係式を用いることにより、新たに未知のＸ２が与えられた場合にもＹ２を得ることができる。ここではＹ２とＸ２のみの関係式を作成したが、同様の手順を行うことによりデータベースに存在しない設計パラメータＸの組み合わせが与えられた場合にもそれに対応した各特性パラメータＹを算出することが可能となる。
【００２６】
次に、作成された関係式を用いた本発明の実施例を具体的に示す。
【００２７】
まず、データベース内に、入力した設計パラメータと同じものが存在した場合について説明する。Ｘ１＝７５ｍｍ、Ｘ２＝７ｍｍ、Ｘ３＝１００ｍｍを設計パラメータとして入力して、この時の共振周波数Ｙ１、帯域幅Ｙ２、放射効率Ｙ３を求める。入力した設計パラメータとデータベース内の設計パラメータを比較すると、共振周波数Ｙ１、帯域幅Ｙ２、放射効率Ｙ３はデータベース内に存在するので、データベースより、この共振周波数Ｙ１、帯域幅Ｙ２、放射効率Ｙ３を特性パラメータとして得る。
【００２８】
以上の方法により、実筐体に搭載された際において、最適な特性を有するアンテナを設計することが可能となる。
【００２９】
次に、データベース内に、入力した設計パラメータと同じものが存在しないが、関係式より特性パラメータを導くことができる場合について説明する。Ｘ１＝７５ｍｍ、Ｘ２＝７ｍｍ、Ｘ３＝９８ｍｍを設計パラメータとして入力して、この時の共振周波数Ｙ１、帯域幅Ｙ２、放射効率Ｙ３を求める。ここで、入力した設計パラメータＸ１，Ｘ２，Ｘ３は図７のデータベース内に存在しないので、設計パラメータと特性パラメータより構成される関係式を用いて共振周波数Ｙ１、帯域幅Ｙ２、放射効率Ｙ３を得る。
【００３０】
次に、上記で求めた関係式で表される範囲外のパラメータの入力があった場合について説明する。
【００３１】
図２は回路基板にモノポールアンテナを取り付けた時の正面図であり、図３は回路基板にモノポールアンテナを取り付けた時の側面図である。図２において材質が純銅であるアンテナ素子１の長手方向の長さをＸ１、アンテナ端部から回路基板３上に取り付けられたシールドケース２の端までの距離をＸ２、回路基板３の長手方向の長さをＸ３、アンテナ端からスピーカー４までの距離をＸ４とする。次にＸ１＝７５ｍｍ、Ｘ２＝７ｍｍ、Ｘ３＝９５ｍｍ、Ｘ４＝３ｍｍを設計パラメータとして入力して、この時の共振周波数Ｙ１、帯域幅Ｙ２、放射効率Ｙ３を求める。ここで前記Ｘ１，Ｘ２，Ｘ３，Ｘ４はデータベースには存在しないので、共振周波数Ｙ１、帯域幅Ｙ２、放射効率Ｙ３はデータベース内に存在しない。また設計パラメータと特性パラメータより構成される前記（数１）〜（数３）の関係式においてもＸ４という変数は存在しないため、関係式からもＹ１，Ｙ２，Ｙ３を得ることができない。この場合には電磁界シミュレータを使用して数値計算を行い、共振周波数Ｙ１、帯域幅Ｙ２、放射効率Ｙ３を得る。ここで電磁界シミュレータを用いて得られたデータは新たにデータベースに蓄積し、次に本システムを利用する時のデータとすることができる。
【００３２】
以上のように、本実施の形態によれば、特性パラメータのデータベースと、アンテナ設計パラメータと特性パラメータの関係式を利用することにより、短期間でアンテナを設計することができる。
【００３３】
（実施の形態２）
以下、本発明の実施の形態２について、本発明の特に請求項８の発明について、図面を参照しながら説明する。なお、実施の形態１と同様の構成を有するものについては、その説明を省略する。
【００３４】
データベースの作成については実施の形態１と同様の手順を用いる。データベース作成後、データベースのデータから関係式を導く。本手法を用いるにあたり、用いるデータをニューラルネットワークに適したデータ形式に変化する必要がある。ニューラルネットワークに入力として与えるデータは学習精度・学習速度の観点から０〜１の間の数であることが望ましい。
【００３５】
ニューラルネットワークの入力パラメータである設計パラメータには寸法などの数値データとアンテナの種類などの非数値データの２種類のデータ形式が存在すると考えられる。ここで寸法などの数値データは、学習精度・速度を考慮すると、０〜１の間に正規化する、もしくは２進数により０と１のみで表される形に変換するのが望ましい。一方アンテナの種類などの非数値データについては、まず各種類に番号を割り振り、０〜１の間に正規化する、もしくは２進数により０と１のみで表される形に変換するのが望ましい。
【００３６】
次に実際のシステムについて説明する。図８に示すように、システムの構成は入力層、中間層、出力層の３層からなる。入力層、中間層、出力層の数は任意に決定することができるが、図８では例として入力層が４、中間層が３、出力層が１の場合を例に挙げている。入力層と中間層、中間層と出力層はそれぞれ重ね付けされた結線により接続されている。ここで中間層の出力を算出するために、（数６）に示すように入力層の値に結線の重みを掛け合わせたものを、（数７）に示すようなシグモイド関数に入力値として与え、これによって出力を得る。
【００３７】
【数６】
2005158019.tif 000008
【００３８】
【数７】
2005158019.tif 000009
【００３９】
次に中間層の出力を入力として同様の手順を繰り返す。これにより出力層の出力の値が得られる。ここでこの出力層での出力を出力信号、設計パラメータを入力として電磁界解析を行うことにより得られた特性パラメータを教師信号と呼ぶこととする。この出力信号と教師信号の値を（数８）に示すような誤差評価尺度を用いて比較を行う。（数８）の値が大きい場合にはシステム内部の結線の重みを修正し（数８）の右辺が最小となるように設定する。
【００４０】
【数８】
2005158019.tif 000010
【００４１】
データベースより（Ｘ１，Ｘ２，Ｙ２）＝［（75、10、1000），（70、3、1020），（75、1、980），（60、10、1070），（70、2、960）］が得られるとする。ニューラルネットワークにおいてはこれらを学習セットと呼ぶ。まずこれらのデータをニューラルネットワークに適したデータに変換する必要がある。ここでは、Ｘ１，Ｘ２，Ｙ２のそれぞれの最大値によって各数値を正規化する。正規化を行った結果（Ｘ１，Ｘ２，Ｙ２）＝［（1、1、0.9345），（0.93、0.3、0.9533），（1、0.1、0.9159），（0.8、1、1），（0.93、0.2、0.8972）］となる。ここでＸ１，Ｘ２を入力信号、Ｙ２を教師データと呼ぶ。これらのデータを基にニューラルネットワークを構築する。
【００４２】
以下にニューラルネットワークの作成手順を示す。
【００４３】
［手順１］まず初めに各ノードにおける重みｗの値を適当に設定する。
【００４４】
［手順２］例えば（1、1、0.9345）という学習セットを与えた場合には、このシステムにＸ１＝１、Ｘ２＝１という入力を与え、（数６）、（数７）により出力信号を得る。
【００４５】
［手順３］出力信号と、教師信号（ここでは0.9345）をもとに（数８）により誤差を計算する。
【００４６】
［手順４］手順２、手順３を全ての学習セットについて行い、２乗誤差の総和を計算する。
【００４７】
［手順５］ここで、誤差の閾値をあらかじめ設定しておき、（数８）が閾値より大きい場合には各ノードの重みを変更し、（数８）が閾値を下回るまで上記手順を繰り返す。
【００４８】
上記手順により、入力に対して最適な出力を得るニューラルネットワークが構築される。本システムを用いることにより、データベースに存在しない未知の設計パラメータが入力された場合にも特性パラメータを得ることが可能となる。
【００４９】
アンテナ特性のように多数の要素や要因が複雑に相互作用を及ぼしあうような場合には、設計パラメータと特性パラメータの間の関係式を求めることは非常に困難である。そのため、関数や関係をあらかじめ設定することが前提となる各種近似法と比較して、入力データと出力データのみが分かっていればよいという利点があるニューラルネットワークを用いることで、出力値の精度、システム構築の労力の面で非常に効率的となる。
【００５０】
（実施の形態３）
以下本発明の実施の形態３について、本発明の特に請求項９の発明について説明する。データベースの作成については実施の形態１と同様の手順を用いる。データベース作成後データベースのデータから関係式を導く。ここで設計パラメータｘi、特性パラメータｙiとする。Ｎ個の点（ｘ1、ｙ1），（ｘ2、ｙ2），（ｘ3、ｙ3）……（ｘn、ｙn）を通る単一の多項式はLagrange補間、Newton補間、Neville補間、Chebyshev補間などにより求められる。ここではLagrange補間を用いた方法を記載する。
【００５１】
ここで、全ての点を通る単一の多項式をｐ（ｘ）と置くとｐ（ｘ1）＝ｙ1，ｐ（ｘ2）＝ｙ2，ｐ（ｘ3）＝ｙ3，……ｐ（ｘn）＝ｙnを満たす。このようなｎ−１次の多項式はｐ（ｘ）はLagrangeの補間公式より（数９）で与えられる。
【００５２】
【数９】
2005158019.tif 000011
【００５３】
よって設計パラメータとそれに対応する特性パラメータを（数９）に代入することによってｐ（ｘ）を導出することができる。本システムを用いることにより、未知の設計パラメータｘが入力された場合にも特性パラメータｙの値を算出することが可能となる。
【００５４】
データベースより（Ｘ２，Ｙ２）＝［（3、950），（5、990），（10、1000），（15、1005），（30、1007）］が得られるとする。５つの点を通る４次の曲線は（数１０）で表される。
【００５５】
【数１０】
2005158019.tif 000012
【００５６】
本システムを用いることにより、データベース内に存在しない設計パラメータが入力された場合にも、設計パラメータを（数１０）に入力することにより特性パラメータを推定することが可能となる。本システムを用いた場合には、関係式を作成する際に使用するデータの数が多くなるほど多項式の次数は高くなり、推定される特性パラメータの精度を上げることができる。よってデータベースが大規模になるほど非常に有用なシステムとなる。
【００５７】
（実施の形態４）
以下本発明の実施の形態４について、本発明の特に請求項１０の発明について説明する。データベースの作成については実施の形態１と同様の手順を用いる。データベース作成後データベースのデータから関係式を導く。ここで設計パラメータｘi、特性パラメータｙiとする。Ｎ個の点（ｘ1、ｙ1），（ｘ2、ｙ2），（ｘ3、ｙ3）……（ｘn、ｙn）を通るBe' zier曲線を求めるために、まず（数１１）で表されるBernshtein多項式を使用する。
【００５８】
【数１１】
2005158019.tif 000013
【００５９】
（数１１）を用いて、Be' zier曲線はｓを変数として（数１２）、（数１３）で表される。
【００６０】
【数１２】
2005158019.tif 000014
【００６１】
【数１３】
2005158019.tif 000015
【００６２】
（数１２）、（数１３）において０≦ｓ≦１でｓを変化させることにより、データベースに存在しないｘ，ｙの組み合わせを推定することができる。これにより、新たなｘ，ｙのデータベースを作成することが可能となる。ここでｓの刻み幅を小さくする程、新たに作成されるデータベース内のデータ数が多くなり、入力された設計パラメータがデータベース内に存在する確率が高くなる。
【００６３】
データベースより（Ｘ２，Ｙ２）＝［（3、950），（5、990），（10、1000），（15、1005），（30、1007）］が得られるとする。５つの点を通る４次の曲線は（数１４）で表される。
【００６４】
【数１４】
2005158019.tif 000016
【００６５】
ここでBe' zier曲線上の点を用いて新たにＸ２とＹ２が対となるデータベースを作成するために、０≦ｓ≦１の間でｓを０．００１刻みで変化させ、（数１４）から新たなＸ２とＹ２を得る。
【００６６】
本システムを用いた場合には、システムの構築に複雑な計算を行うことなくデータベースのデータ数を増やすことが可能となる。よって、既にデータベースに大量のデータがあり、他手法では関係式の導出に非常に計算時間を有すると予想される場合に有用である。
【００６７】
（実施の形態５）
以下本発明の実施の形態５について、本発明の特に請求項１１の発明について説明する。データベースの作成については実施の形態１と同様の手順を用いる。データベース作成後、データベースのデータから関係式を導く。ここで設計パラメータをｘi、特性パラメータをｙiとする。
【００６８】
Ｎ個の点（ｘ1、ｙ1），（ｘ2、ｙ2），（ｘ3、ｙ3）……（ｘn、ｙn）から区分的多項式を求めるためにはスプライン補間、Ｂスプライン補間、NURBSなどが考えられる。ここでは３次のスプライン補間を用いた例を記載する。
【００６９】
区間［ｘi，ｘi+1］における３次の多項式をＳi（ｘ）とすると（数１５）で表すことができる。
【００７０】
【数１５】
2005158019.tif 000017
【００７１】
但し（数１６）が成り立つことが前提となる。
【００７２】
【数１６】
2005158019.tif 000018
【００７３】
（数１５）、（数１６）から（数１７）が得られる。
【００７４】
【数１７】
2005158019.tif 000019
【００７５】
ここで
【００７６】
【数１８】
2005158019.tif 000020
【００７７】
（数１７）の漸化式を連立方程式の形にして、この連立方程式を（数１９）に示すような行列式に変換して解く。但し未知数が２つ多いので、（数１６）の（３）の境界条件を用いることによりｙi'を求める。ここで（数１７）におけるｙi'の係数をＡji、右辺をＢiとする。
【００７８】
【数１９】
2005158019.tif 000021
【００７９】
（数１５）におけるａi，ｂi，ｃi，ｄiは、ｙi'を用いて（数２０）で表される。
【００８０】
【数２０】
2005158019.tif 000022
【００８１】
（数２０）を計算し、（数１５）に代入することにより、関係式が作成される。本システムにより未知の設計パラメータｘが入力された場合にも特性パラメータｙの値を算出することが可能となる。
【００８２】
データベースより（Ｘ２，Ｙ２）＝［（3、950），（5、990），（10、1000），（15、1005），（30、1007）］が得られるとする。これら５つの点を通る３次のスプライン関数は（数１５）で表される。
【００８３】
（数１５）にデータベース内に存在する上記５点を入力すると（数２１）の行列式が得られる。
【００８４】
【数２１】
2005158019.tif 000023
【００８５】
（数２１）よりｙ1'，ｙ2'，ｙ3'，ｙ4'を求め、各区間における３次多項式を決定する。
【００８６】
本システムを用いた場合には、入力された設計パラメータの存在する範囲により適用される多項式が異なるため、非常に高い近似能力を有する。よって非常に高精度な特性パラメータの推定を行うことが可能となる。
【００８７】
（実施の形態６）
以下本発明の実施の形態６について、本発明の特に請求項１２の発明について説明する。データベースの作成については実施の形態１と同様の手順を用いる。データベース作成後、データベースのデータから関係式を導く。ここで設計パラメータｘi、特性パラメータｙiとする。ｎ個の点（ｘ1、ｙ1），（ｘ2、ｙ2），（ｘ3、ｙ3）……（ｘn、ｙn）から有理Be' zier曲線を求めるために、まず（数２２）で表されるBernshtein多項式を使用する。
【００８８】
【数２２】
2005158019.tif 000024
【００８９】
（数２２）を用いて各節点に重みをつけることにより、Be' zier曲線よりも柔軟性の高い曲線を求めることができる。各節点に重みをつけた時の式を（数２３）に示す。
【００９０】
【数２３】
2005158019.tif 000025
【００９１】
【数２４】
2005158019.tif 000026
【００９２】
（数２３）、（数２４）において０≦ｓ≦１でｓを変化させることにより、データベースに存在しないｘ，ｙの組み合わせを推定することができる。これにより、新たなｘ，ｙのデータベースを作成することが可能となる。ここでｓの刻み幅を小さくする程、新たに作成されるデータベース内のデータ数が多くなり、入力された設計パラメータがデータベース内に存在する確率が高くなる。
【００９３】
データベースより（Ｘ２，Ｙ２）＝［（3、950），（5、990），（10、1000），（15、1005），（30、1007）］が得られるとする。但し、これらのデータについてデータの信頼度が異なるとする。ここでデータの信頼度を（3、950），（5、990），（10、1000），（15、1005），（30、1007）のそれぞれの点において０．５，０．３，０．３，１，０．８と仮定し、これらのデータの信頼度を各点の重みとして（数２３）、（数２４）を計算する。その結果、４次の有理Be' zier曲線は（数２５）で表される。
【００９４】
【数２５】
2005158019.tif 000027
【００９５】
ここで、有理Be' zier曲線上の点を用いて新たにＸ２とＹ２が対となるデータベースを作成するために、０≦ｓ≦１の間でｓを０．００１刻みで変化させ（数２５）からＸ２とＹ２を得る。
【００９６】
本システムを用いた場合には、システムの構築に複雑な計算を行うことなくデータベースのデータ数を増やすことが可能なため、既にデータベースに大量のデータがあり、他手法では関係式の導出に非常に計算時間を有すると予想される場合に有用である。またデータベースに蓄積されたデータの信頼性にばらつきがある場合には、本システムを用いることにより高精度なデータベースを再構築することができる。
【００９７】
（実施の形態７）
以下、本発明の実施の形態７について、本発明の特に請求項１３の発明について説明する。データベースの作成については実施の形態１と同様の手順を用いる。データベース作成後、データベースのデータから関係式を導く。ここで設計パラメータｘi、特性パラメータｙiとする。Ｎ個の点（ｘ1、ｙ1），（ｘ2、ｙ2），（ｘ3、ｙ3）……（ｘn、ｙn）が周期関数で表現できる場合には、三角関数による補間を用いる。
【００９８】
三角関数を用いた近似はフーリエ級数の式より（数２６）で表される。
【００９９】
【数２６】
2005158019.tif 000028
【０１００】
ここで
【０１０１】
【数２７】
2005158019.tif 000029
【０１０２】
である。
【０１０３】
（数２７）に各点を代入して、（数２６）における係数を求めることにより、設計パラメータと特性パラメータの関係式が作成される。本システムを用いることにより、未知の設計パラメータｘが入力された場合にも特性パラメータｙの値を算出することが可能となる。
【０１０４】
データベースよりＸ２に関して等間隔なデータ（Ｘ２，Ｙ２）＝［（1、900），（3、950），（5、990），（7、994），（9、998）］が得られるとする。５つの点を通る関数は（数２８）で表される。
【０１０５】
【数２８】
2005158019.tif 000030
【０１０６】
（数２８）を（数２６）に代入することにより、関係式を導くことが可能となる。ここでｎの数を大きくするほど関係式の精度は高くなる。本システムを用いた場合入力する設計パラメータをＸ２とすると、完成した関係式に入力する際にはＸ２＝（Ｘ２−１）／２とする必要がある。
【０１０７】
本システムを用いた場合には、データベース上に等間隔となる入力パラメータが存在する場合にしか使用できないという制約がある。しかし近似能力が非常に高いため、アンテナ設計における設計パラメータと特性パラメータのように、どのような関係式になるか全く予想ができない時には非常に有用である。
【産業上の利用可能性】
【０１０８】
本発明のアンテナの設計方法は、携帯電話の構成に対応した最適なアンテナを短時間で設計するツールとして有用である。
【図面の簡単な説明】
【０１０９】
【図１】本発明の実施による携帯電話用アンテナ特性解析のフローチャート
【図２】シールドケースを有する回路基板に取り付けられたモノポールアンテナの正面図
【図３】シールドケースを有する回路基板に取り付けられたモノポールアンテナの側面図
【図４】インピーダンス特性を評価する範囲内において、１つの共振を有する時のＶＳＷＲの周波数特性図
【図５】インピーダンス特性を評価する範囲内において、２つの共振を有する時のＶＳＷＲの周波数特性図
【図６】インピーダンス特性を評価する範囲内における放射効率の周波数特性図
【図７】設計パラメータが３変数、特性パラメータが３変数の時のデータベースを示す図
【図８】ニューラルネットワークの概念図
【符号の説明】
【０１１０】
１  アンテナ素子
２  シールドケース
３  回路基板
Ｘ１  アンテナ素子長
Ｘ２  アンテナ素子端とシールドケースの距離
Ｘ３  回路基板長
Ｘ４  アンテナ素子端とスピーカー端の距離
Ｙ１  インピーダンス特性を評価する範囲内の１つ目の共振における共振周波数
Ｙ２  インピーダンス特性を評価する範囲内の１つ目の共振におけるＶＳＷＲ＜３となる帯域幅
Ｙ３  指定した周波数における放射効率
Ｙ４  インピーダンス特性を評価する範囲内の２つ目の共振における共振周波数
Ｙ５  インピーダンス特性を評価する範囲内の２つ目の共振におけるＶＳＷＲ＜３となる帯域幅
【図１】
2005158019.tif 000031
【図２】
2005158019.tif 000032
【図３】
2005158019.tif 000033
【図４】
2005158019.tif 000034
【図５】
2005158019.tif 000035
【図６】
2005158019.tif 000036
【図７】
2005158019.tif 000037
【図８】
2005158019.tif 000038