(57)【要約】【目的】低サンプリング信号を高サンプリング信号に精度良好にレート変換できる「FIFを用いたサンプリングレート変換方法」を提供する。【構成】原信号Ｓのサンプリングレートをｋ倍にアップレートした拡大信号Ｓkを発生するフラクタル補間機能(FIF)を用いたサンプリングレート変換方法であり、原信号をFIF補間区間に写像するための写像パラメータａi，ｃi，ｄi，ｅi，ｆiと前記ｋを用いて、前記拡大信号を発生するための写像パラメータを求め、該求めた写像パラメータを用いて原信号のサンプリングレートをアップレートした拡大信号を発生する。【選択図】        図８
2005084370.tif 000002
【特許請求の範囲】
【請求項１】
原信号のサンプリングレートをｋ倍にアップレートした拡大信号を発生するフラクタル補間機能(FIF)を用いたサンプリングレート変換方法において、
原信号をFIF補間区間に写像するための写像パラメータと前記ｋを用いて、前記拡大信号を発生するための写像パラメータを求め、
該求めた写像パラメータを用いて原信号のサンプリングレートをアップレートした拡大信号を発生する、
ことを特徴とするFIFを用いたサンプリングレート変換方法。
【請求項２】
原信号のFIF補間区間を特定する代表点を用いて拡大信号のFIF補間区間を特定する代表点求め、
前記拡大信号を発生するため写像パラメータと前記拡大信号のFIF補間区間を特定する代表点を用いてサンプリングレートをｋ倍にアップレートした拡大信号を発生する、
ことを特徴とする請求項１記載のFIFを用いたサンプリングレート変換方法。
【請求項３】
前記拡大信号を発生するための写像パラメータを用いて前記拡大信号のFIF補間区間を特定する代表点にランダム繰り返しアルゴリズム（RIA）を適用して前記拡大信号を近似したアトラクタを生成する、
ことを特徴とする請求項２記載のサンプリングレート変換方法。
【請求項４】
ランダム繰り返しアルゴリズムにより得られた信号値を１／ｋして、前記アトラクタを生成する、
ことを特徴とする請求項３記載のサンプリングレート変換方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明はFIF (Fractal Interpolation Functions)を用いたサンプリングレート変換方法に係わり、特に、オーディオ信号など低サンプリング信号を高サンプリング信号に変換するFIFを用いたサンプリングレート変換方法に関する。
【背景技術】
【０００２】
（ａ）ＦＩＦ
オーディオ分野やグラフィック分野ではサンプリング信号間を補間する技術が種々研究されて提案されている。かかる補間技術としてFIF (Fractal Interpolation Functions)がある。FIFは補間区間に該補間区間のM倍の長さの信号波形を写像して補間する技術である。以下において、本発明のレート変換処理を行う基礎となるこのFIF(たとえば非特許文献1を参照)について説明する。
図1に示すように一次元離散信号S＝｛（un,vn）:n=0,1,…,N｝が与えられた際のFIF処理手順を以下に示す。ただし、信号Sは次式を満足する1価関数で表される信号である。
ｕ0＜ｕ1＜．．．．．＜ｕｎ        (1)
まず、信号Sから、（M+1）個の代表点P＝｛（ｘi,ｙi）:i=0,1,…,M｝を選択する（図1中の黒丸点）。ただし、信号Sの両端点は次式に示すように、代表点として無条件に選択するものとする。
【０００３】
【数１】
2005084370.tif 000003
上記（M+1）個の代表点を選択することにより、信号SはM個の区間(第1区間〜第M区間)に分割される。尚、以下では、連続する2つの代表点により定義される区間〔xi-1,xi〕を補間区間ｉと呼ぶ。
【０００４】
次に、式（3）に示すアフィン写像ｗｉを適用することで、信号ＳをＭ個の補間区間ｉ（ｉ＝０〜Ｍ）にそれぞれ写像する。図2は信号Ｓを補間区間ｉに写像した例を示している。
【数２】
2005084370.tif 000004
上式において、ｗｉは信号Ｓを補間区間iに写す縮小写像である。したがって、FIFを適用することで与えられた信号Sは、次式に示すように縮小写像ｗｉ（ｉ＝１〜Ｍ）による自身の縮小像ｗi（S）の和集合として表現されることがわかる。
【０００５】
【数３】
2005084370.tif 000005
上式において、GはIterated Function Systems(IFS)のアトラクタであり、自己アフィン特性を有することが知られている。
ところで、式（3）において、5つの未知パラメータaｉ, ci, di, ei, fi（以降、写像パラメータと呼ぶ）が存在する。式（3）を実際の信号に適用するには、これら5つの未知パラメータaｉ, ci, di, ei, fIを求めなければならない。そこで次式に示す制約条件を設ける。
【０００６】
【数４】
2005084370.tif 000006
上記制約により、図2の矢印Ａ，Ｂに示すように、信号Sの端点Ts,Teは補間区間iの端点tis,tieに写される。このような制約を設け、さらに5つ存在する写像パラメータの内、縮小因子と呼ばれるパラメータｄiを変数と考えると、他の4つの写像パラメータはそれぞれ、以下に示すように表すことができる。
【０００７】
【数５】
2005084370.tif 000007
FIFを用いて与えられた信号を高精度に表現するためには、代表点及び縮小因子をいかに決定するかという問題が生じる。この問題はFIFにおけるinverse problemと呼ばれており、これまでにその解法が幾つか提案されているが、本文ではMazelらにより提案された手法(たとえば非特許文献2参照)を用いてinverse problemを解決する。以下において、その方法を説明する。
【０００８】
（ｂ）写像パラメータの決定
先に述べたように、FIFを適用することで、与えられた信号Sは各補間区間に写される。このとき、縮小写像ｗｉによる信号Ｓの像ｗｉ(s)={(pn,qn):n=0,1,…,N}は(図3参照)、次式に示すように表すことができる。図3において大きな白丸は信号S上ののデータポイント、小さな黒丸は縮小写像ｗｉ(s)上のN個のポイントである。
【数６】
2005084370.tif 000008
前節で述べたように、FIFを適用する際には、与えられた信号Sを高精度に近似したアトラクタGをいかに生成するかが問題となるが、この問題は、信号Sの部分集合S[xi-1,xi]と縮小像ｗｉ(S)との誤差を最小化することにより解決できる。ここで、信号S[xi-1,xi]と縮小像ｗｉ(S)との誤差をＥｉとすると、Ｅｉは図3に示すように、信号S[xi-1,xi]及び縮小像ｗｉ(S)を構成するデータ点の垂直方向の距離を足し合わせることにより、次式に示すように定式化できる。
【０００９】
【数７】
2005084370.tif 000009
尚、次式（12）において[・]はガウス記号を表す。式（11）に式（8）、（9）をそれぞれ代入し、整理すると次式を得る。
【数８】
2005084370.tif 000010
次に、式（13）を縮小因子ｄｉに関して最小2乗規範に基づき最小化すると、縮小因子ｄiは次式により与えられる。
【００１０】
【数９】
2005084370.tif 000011
先に述べたように、縮小因子ｄiの価が決定されれば、残り4つの写像パラメータは一意に定めることができる。したがって、対象となる信号が全て既知のデータ点により構成されている場合、そのinverse problemは比較的容易に解くことができる。
【非特許文献１】M.F.Barbsley, “Fractal functions and interpolation” Constructive Approximation, vol.2, pp.303-329, 1986
【非特許文献２】D.S.Mazel and M.H. Hayes, “Using Iterated function systems to model discrete sequences,” IEEE Trans. Signal Processing, vol.40, No.7, pp.1243-1249, Oct 1982.
【発明の開示】
【発明が解決しようとする課題】
【００１１】
以上述べたように、対象となる信号が全て既知のデータ点により構成されている場合、そのinverse problemは比較的容易に解くことができる。しかしながら、レート変換（アップサンプリング）後の信号は、これから補間されるべき、未知のデータ点を含むことからinverse problemを解くことは困難である。
以上から本発明の目的はレート変換（アップサンプリング）後の信号のinverse problemを解き、低サンプリング信号を高サンプリング信号に精度良好にレート変換できるFIFを用いたサンプリングレート変換方法を提供することである。
【課題を解決するための手段】
【００１２】
本発明は、原信号のサンプリングレートをｋ倍にアップレートした拡大信号を発生するフラクタル補間機能(FIF)を用いたサンプリングレート変換方法であり、原信号をFIF補間区間に写像するための写像パラメータと前記ｋを用いて、前記拡大信号を発生するための写像パラメータを求め、該求めた写像パラメータを用いて原信号のサンプリングレートをアップレートした拡大信号を発生する。すなわち、原信号のFIF補間区間を特定する代表点を用いて拡大信号のFIF補間区間を特定する代表点求め、前記拡大信号を発生するため写像パラメータと前記拡大信号のFIF補間区間を特定する代表点を用いてサンプリングレートをｋ倍にアップレートした拡大信号を発生する。
具体的には、前記拡大信号を発生するための写像パラメータを用いて前記拡大信号のFIF補間区間を特定する代表点にランダム繰り返しアルゴリズム（RIA）を適用して前記拡大信号を近似したアトラクタを生成する。この場合、ランダム繰り返しアルゴリズムにより得られた信号値を１／ｋして、前記アトラクタを生成する。
【発明の効果】
【００１３】
本発明によれば、レート変換（アップサンプリング）後の信号のinverse problemを解いて低サンプリング信号を高サンプリング信号に精度良好にレート変換することができる。
【発明を実施するための最良の形態】
【００１４】
原信号のサンプリングレートをｋ倍にアップレートした拡大信号を発生するフラクタル補間機能(FIF)を用いたサンプリングレート変換方法において、原信号をFIF補間区間に写像するための写像パラメータと前記ｋを用いて、前記拡大信号を発生するための写像パラメータを求め、かつ、原信号のFIF補間区間を特定する代表点を用いて拡大信号のFIF補間区間を特定する代表点求め、前記拡大信号を発生するための写像パラメータと前記拡大信号のFIF補間区間を特定する代表点を用いてサンプリングレートをｋ倍にアップレートした拡大信号を発生する。以上により、本発明は、レート変換（アップサンプリング）後の信号のinverse problemを解いて、すなわち、拡大信号のFIF補間区間を特定する代表点と拡大信号を発生するための写像パラメータを求めて、低サンプリング信号を高サンプリング信号に精度良好にレート変換することができる。
【実施例１】
【００１５】
(Ａ)写像パラメータの決定
信号の拡大処理の前後における写像パラメータの関係を導出することにより、拡大信号の写像パラメータを特定することができる。ここで、“信号をｋ倍に拡大する”とは、信号のサンプリングレート及び信号値(振幅)をそれぞれｋ倍する処理を指す。
まず、FIFを用いた一次元離散信号のレート変換手法のアルゴリズムを説明する。なお、与えられた信号をS、Sに比べてサンプリングレート及び信号値がｋ倍である信号Ｓkを次式
【数１０】
2005084370.tif 000012
で定義されるSkとする。ただし、信号Skは次式を満足するものとする。
【数１１】
2005084370.tif 000013
このとき
、信号Ｓｋに関する
【数１２】
2005084370.tif 000014
は、次式に示すように、信号Sに関する代表点（xi，yi)をu軸及びｖ軸方向にそれぞれｋ倍した位置に配置する。
【００１６】
【数１３】
2005084370.tif 000015
ここで、信号Ｓｋに関する写像パラメータを
【数１４】
2005084370.tif 000016
とすると、写像パラメータa^ｉは次式より与えられる。なお、
【数１５】
2005084370.tif 000017
である。
【数１６】
2005084370.tif 000018
上式に式（18）を代入すると次式を得る。
【数１７】
2005084370.tif 000019
上式から、写像パラメータaiは拡大処理の前後において不変であることがわかる。他の写像パラメータei, ，ci，fiについても、同様の処理を施すことで、それぞれ以下に示す関係を導くことができる。
【００１７】
【数１８】
2005084370.tif 000020
ただし、写像パラメータ
【数１９】
2005084370.tif 000021
は縮小因子d^ｉの関数であることから、現時点において、
【数２０】
2005084370.tif 000022
の関係を明らかにすることはできない。
【００１８】
そこで、次に縮小因子について考える。縮小写像ｗ^ｉ
による信号Skの像を
【数２１】
2005084370.tif 000023
とすると
【数２２】
2005084370.tif 000024
は以下のように表すことができる。
【数２３】
2005084370.tif 000025
したがって、図4に示すように、補間区間iについて、信号Sｋとその縮小像
【数２４】
2005084370.tif 000026
との誤差E^ｉは次式により与えられる。
【００１９】
【数２５】
2005084370.tif 000027
しかし、以上は信号Skを構成する図4のデータ点(大径の黒丸点)が全て既知の場合である。しかし、ｋ倍に拡大する前はこれらのデータ点(大径の黒丸点)は総て既知でなく、図5に示すように1/k個のみ(図5の大径の黒丸点)のみが既知であるに過ぎない。換言すれば、先に述べたように、信号Sｋは補間されるべき未知のデータ点(図5の白丸点)を含むため、式（28）を用いて誤差を算出することは困難である。
【００２０】
そこで、信号Skの内、既知のデータ点のみを用いると、次式に示すように，誤差を近似することができる。
【数２６】
2005084370.tif 000028
式（31）に式（17）、（21）、（22）をそれぞれ代入すると次式を得る。
【００２１】
【数２７】
2005084370.tif 000029
次に、式（30）に式（23）、（25）、（32）を代入すると、誤差は以下のように書き表すことができる。
【数２８】
2005084370.tif 000030
上式に式（17）を代入し、整理すると次式を得る。
【００２２】
【数２９】
2005084370.tif 000031
ただし、上式において、αｎ，βｎはそれぞれ次式（14）、（15）と一致する。式（34）を縮小因子d^iに関して、最小2乗規範に基づき最小化すると、以下に示すdiとd^iの関係式を導くことができる。
【数３０】
2005084370.tif 000032
上式から、拡大処理の前後において、縮小因子は不変であることがわかる。そこで式（35）を式（24）、（26）にそれぞれ代入すると、次式に示す関係式が得られる。
【００２３】
【数３１】
2005084370.tif 000033
したがって、FIFにおいて5つ存在する写像パラメータの内ai,, ci, diは拡大処理の前後において不変であり、ei,, fiは、拡大率ｋを乗じることで、拡大処理後の値を決定することが可能である。すなわち、拡大信号Ｓkiの各補間区間に対する写像パラメータ
【数３２】
2005084370.tif 000034
は、それぞれ信号Siに関する写像パラメータを用いて以下のように決定することができる。
【数３３】
2005084370.tif 000035
本発明はこの性質を用いて、与えられた離散信号のサンプリングレートを変換する。以下にその手順を説明する。
【００２４】
(Ｂ) FIFを用いたレート変換アルゴリズム
図6はFIFを用いたレート変換アルゴリズムの処理フローであり、離散信号Sが与えられた際、まず、FIFによりサンプリングレート及び信号値がｋ倍である信号Skを推定し、さらにその信号値をI／k倍することで、レート変換後の信号を獲得するまでの処理手順を示している。
(1)与えられた信号のモデル化
まず、図７に示すように、信号SをそれぞれＮ個のサンプル点から成るL個の自己アフィン区間に等分割し、各々の信号を｛Si:i=1,2,….,L｝とする。次に、信号Siに対して周知のベクトルトレーサ法を適用することにより、（Mi+1）個の代表点Pi={(xi,yi):ｉ=0,1,…,Mi}を選択する(ステップ１)。なお、以降では信号Sｉ，Ｐｉ，Ｍｉを単に信号Ｓ，Ｐ，Ｍと表記する。
先に述べたように、信号Sは全て既知のデータ点から構成されているため、写像パラメータを決定することは容易である。そこで、式（6）−(９)及び式（16）を用いて写像パラメータai, ci, di, ei, fiを決定する(ステップ２)。
なお、ベクトルトレーサ法は以下の文献に記述されている。
［文献］安居院猛、飯塚久登、中嶋正之、“ピラミッド階層構造データの位相変化情報を利用した市街地図の処理” 信学論(D),vol.40,no.7, pp.1243-1249,Oct 1982.
【００２５】
(2)拡大信号に関するinverse problemの解決
上記ステップ１，２の処理により、与えられた信号Ｓに関するinverse problemは解決された。そこで次に拡大処理後の信号Sｋに関するinverse problemを解決する。まず、信号Sに関する代表点P＝｛（xi,yi）:i=0,1,…,M｝をu軸、v軸方向にそれぞれｋ倍することで、信号Sｋに関する代表点
【数３４】
2005084370.tif 000036
を決定する(ステップ３)。次に、拡大信号Sｋに関する写像パラメータ
【数３５】
2005084370.tif 000037
を式（38）を用いて決定する(ステップ４)。
【００２６】
(3)拡大信号の構成
次に、フラクタル関数を生成する際に広く用いられている周知のRandom Iteration Algorithm(RIA)（後述）を適用することにより、拡大処理後の信号Sｋを近似したアトラクタGｋを生成する(ステップ５)。ここで、（ｋN＋１）個のサンプル点から構成されるアトラクタGkを生成するためにはRIAにおいて、（ｋN＋１）回の反復が必要となる。これにより、次式に示すアトラクタGkが得られる。
【数３６】
2005084370.tif 000038
ただし、アトラクタGkは与えられた信号に比べて、その信号がｋ倍されているので、次式に示すようにその信号値を1/ｋ倍したＧ^ｋ
がレート変換後の信号を近似したアトラクタとなる(ステップ６)。
【数３７】
2005084370.tif 000039
【００２７】
以上本発明のFIFを用いたレート変換手法によれば、与えられた信号に対する拡大率ｋが決定されれば、拡大処理後の信号に関する代表点、写像パラメータ及びRIAにおける反復回数等を一意に定めることができる。したがって、与えられた信号に対する、任意の拡大率での拡大処理を容易に実現することができる。
実際に音声信号のレート変換を行う際には、計算コストを削減するために自己アフィン区間を分割して処理を行う必要がある。
【００２８】
（Ｃ）FIFを用いたサンプリングレート変換装置
図８はFIFを用いたレート変換アルゴリズムに従ってサンプリングレートをｋ倍する本発明のサンプリングレート変換装置の構成図であり、代表点決定部101はステップ1に従って信号Ｓの（M+1）個の代表点P={(xi,yi):ｉ=0,1,…,M}を決定する。第1の写像パラメータ算出部102はステップ2に従って式（6）−(９)及び式（16）を用いて写像パラメータai, ci, di, ei, fiを決定する。拡大信号の代表点決定部103は、ステップ3に従って拡大信号Sｋに関する代表点
【数３８】
2005084370.tif 000040
を決定する。第2の写像パラメータ算出部104はステップ４に従って式(38)式に基づいて拡大信号に関する写像パラメータ
【数３９】
2005084370.tif 000041
を決定する。アトラクタ生成部105は、周知のRandom Iteration Algorithm(RIA)を適用することにより、拡大処理後の信号Sｋを近似したアトラクタGｋを生成する。演算部106はアトラクタGkの信号値を1/ｋ倍してレート変換後の信号を近似したアトラクタG^k
を算出する。
【００２９】
Ｓ，Sｋ，Gkにサフィックスｉを付加すれば、レート変換後の信号Sｋは次式に示すようにアトラクタG^kｉの和集合として表される。
【数４０】
2005084370.tif 000042
【００３０】
（D）Random Iteration Algorithmによるアトラクタの生成
FIFにおいて、アトラクタを生成する際に広く用いられているRandom Iteration Algorithm(ＲＩＡ：ランダム繰り返しアルゴリズム)について説明する。ただし、ここでは対象となる信号に関するinverse problem,すなわち代表点
【数４１】
2005084370.tif 000043
及び各補間区間に関する縮小写像
【数４２】
2005084370.tif 000044
は既に求まっているものとする。
【００３１】
まず、上記M個の縮小写像w^iの中からランダムに1つの写像を選択する。次に任意の代表(x^0,y^0)を初期値として、(x^0,y^0)に対して選択した写像を施すことで(x^1,y^1)を得る。さらに(x^1,y^1)に対して、再びランダムに選択された写像を施すことで(x^2,y^2)を得る。同様の処理をN回繰り返すことで合計（N+1）個のデータからなる点列
【数４３】
2005084370.tif 000045
が得られる。
その様子を図９に示す。図9において黒丸は代表点、白丸はRandom Iteration Algorithmにより順次求まるポイントである。
【００３２】
尚、Random Iteration Algorithmを適用することにより得られる点列は、N→∞のときIFSのアトラクタに収束することが知られている。ところで、アップサンプリング後の信号を構成するデータ数は、与えられた信号を構成するデータ数に比べて大きい。したがって、アップサンプリング後の信号は、与えられた信号に比べて、より強い自己アフィン特性を有することが期待される。Random Iteration Algorithmの詳細については以下の文献を参照されたい。
［文献］M.F.Barnsley, Fractals Everywhere 2nd  ed., pp.87 Academic Press, San Diego, 1993
【００３３】
（E）実験結果
実験には、図10(a),(b)に示す音声信号A，Bを用いた。これらの音声信号は、CNN World Newsにおける女性キャスターの音声を44.1ｋHｚで60秒間サンプリングすることで得た信号の内10〜20秒、20〜30秒区間を切り出したものである。実験では、まず信号A、Bを11.025kHｚにダウンサンプリングした信号Ｓを作成し、その信号Ｓを本発明のFIFを用いたサンプリングレート変換方法を適用して再び44.1kHzに戻し、元の信号との比較を行った。尚、本実験では、Athlon1.4GHzのCPUを装着したLinuxマシンを用いた。
1つの自己アフィン区間の大きさNを変えて、上記手法を適用した際の処理時間Tを図11に示す。図11から、自己アフィン区間の大きさを小さく設定することによって、与えられた信号を1つの自己アフィン区間とみなして処理する場合に比べ、レート変換に要する計算時間を大幅に削減できることがわかる。
【００３４】
一方、レート変換精度を原信号との平均絶対誤差（MAEprop）によって評価した結果を図12に示す。図12において、Rは分数比サンプリングレート変換(下記文献を参照)と原信号とのMSEをMAEconvとした場合の比率R= MAEprop/ MAEconvを表す。図12から従来手法である分数比サンプリングレート変換に比べ、わずかではあるがFIFを用いたレート変換手法の精度が良いことがわかる。また、自己アフィン区間の大きさNはレート変換精度に影響を与えるものの、その影響は十分無視できる程小さいものであることがわかる。したがって、処理時間を考慮した場合、本手法は自己アフィン区間の大きさNを小さく設定して用いるのが良いと言える。
以下に、本発明により良好なレート変換を実現できる理由を考察する。
従来広く用いられているアップサンプリング手法は、与えられた離散信号のサンプル点間に適当な数のゼロ値を挿入し、さらに低域フィルタを適用することによって所望のサンプリングレートに変換する。しかしながら、この手法によって獲得される信号は、帯域制限されているため原信号が本来有していた高周波成分が含まれていない。
ところで、近年の研究により、音声信号が自己相似性あるいは自己アフィン特性に代表されるフラクタル特徴を有することが明らかになってきた。そこで、本手法では、対象となる音声信号が自己アフィン特性を有することを仮定し、自己アフィン特性を有する離散信号をうまく表現できるFractal Interpolation Functions(FIF)を用いて対象となる信号のサンプリングレート変換を行う。
FIFを用いて離散信号を表現する際には、代表点および写像パラメータを決定後、Random Iteration Algorithm(RIA)を用いて信号を生成する。RIAでは、反復回数→∞のときIterated Function System(IFS) のアトランタに収束することが知られている。つまり反復回数が多い程、より強い自己アフィン特性を有する信号を生成できると言える。ところで、アップサンプリング後の信号は、与えられた信号に比べてより多くのサンプル点を有している。したがって、本手法を適用することで得られるアップサンプリング後の信号は、与えられた信号に比べてより強い自己アフィン特性を有していると考えられる。信号が強い自己アフィン特性を有することは、その信号が細部において複雑さを保持していることを意味する。したがって、本手法を適用することによって得られる信号は、従来手法を適用することによって得られる信号に比べて、より多くの高周波成分を含んでいるため、良好なレート変換を実現できると言える。
［文献］Alan V. Oppenheim, Ronald W.Schafer, and John R.Bruck, Discrete-Time Signal Processing 2nd ed., pp.167, Prentice Hall, New Jersey, 1999.
【産業上の利用可能性】
【００３５】
本発明のFIFを用いたサンプリングレート変換方法は、低サンプリングレートの信号を高サンプリングレート信号に高精度に変換することができ、オーディオ分野、画像処理分野において有効に利用することができる。
【図面の簡単な説明】
【００３６】
【図１】一次元離散信号S＝｛（un,vn）:n=0,1,…,N｝が与えられた際のFIF処理手順説明図である。
【図２】FIF処理手順において信号Ｓを補間区間ｉに写像した例である。
【図３】信号S[xi-1,xi]と縮小像ｗｉ(S)との誤差Ｅｉを最小化することにより写像パラメータを決定する原理説明図である。
【図４】拡大信号（データ点が全て既知の場合）の写像パラメータ決定原理説明図である。
【図５】拡大信号（データ点が全て既知でない場合）の写像パラメータ決定原理説明図である。
【図６】FIFを用いたレート変換アルゴリズムの処理フローである。
【図７】信号Ｓのモデル化説明図である。
【図８】FIFを用いたレート変換アルゴリズムに従ってサンプリングレートをｋ倍する本発明のサンプリングレート変換装置の構成図である。
【図９】ランダム繰り返しアルゴリズムRIAの説明図である。
【図１０】実験に使用した音声信号波形図である。
【図１１】第1の実験結果である。
【図１２】第2の実験結果である。
【符号の説明】
【００３７】
１０１  代表点決定部
１０２  第1の写像パラメータ算出部
１０３  拡大信号の代表点決定部
１０４  第2の写像パラメータ算出部
１０５  アトラクタ生成部
１０６  演算部
【図１】
2005084370.tif 000046
【図２】
2005084370.tif 000047
【図３】
2005084370.tif 000048
【図４】
2005084370.tif 000049
【図５】
2005084370.tif 000050
【図６】
2005084370.tif 000051
【図７】
2005084370.tif 000052
【図８】
2005084370.tif 000053
【図９】
2005084370.tif 000054
【図１０】
2005084370.tif 000055
【図１１】
2005084370.tif 000056
【図１２】
2005084370.tif 000057