圧力を上げながら、または、所定の圧力から圧力を下げながら血管を加圧する加圧手段と、前記加圧手段で加圧中の前記血管の圧力を測定する測定手段と、前記加圧手段での所定の圧力段階毎に、前記測定手段で測定された前記血管の圧力の平均圧力と脈波による圧力振幅とを求める血管圧力計算手段と、あらかじめモデル計算により、圧力を上げながら、または、所定の圧力から圧力を下げながら血管を加圧したときの血管の平均圧力に対する圧力振幅の変化に関する特性指標と、血管の固さに関するパラメータとの関係を求めておくモデル計算手段と、前記血管圧力計算手段で求められた前記平均圧力に対する前記圧力振幅の変化に関する特性指標を求め、その特性指標から、前記モデル計算手段で求められた関係に基づいて、血管の固さに関するパラメータを求める血管固さ算出手段とを、有することを特徴とする血管固さ評価装置。前記血管の固さに関するパラメータはスティフネスパラメータから成り、前記モデル計算手段は前記関係を求めるとき、血管の弾性を表すチューブ則を利用して、血管の平均圧力に対する圧力振幅の変化を求めて特性指標を求めるとともに、1式によりスティフネスパラメータを求めることを特徴とする請求項1記載の血管固さ評価装置。β:スティフネスパラメータP:血管の内圧Pb:血管の任意の基準内圧A:血管の断面積A0:P=Pbのときの血管の断面積前記血管固さ算出手段は、前記平均圧力に対する前記圧力振幅の変化を、前記圧力振幅のピーク値とそのときの平均圧力の値とで無次元化して求め、前記特性指標として、前記圧力振幅のピーク値を与える平均圧力の半分の圧力のときの二階微分値、および/または、前記圧力振幅が0から変化し始めたときの一階微分値を求めることを特徴とする請求項1または2記載の血管固さ評価装置。圧力を上げながら、または、所定の圧力から圧力を下げながら血管を加圧する加圧工程と、前記加圧工程で加圧中の前記血管の圧力を測定する測定工程と、前記加圧工程での所定の圧力段階毎に、前記測定工程で測定された前記血管の圧力の平均圧力と脈波による圧力振幅とを求める血管圧力計算工程と、あらかじめモデル計算により、圧力を上げながら、または、所定の圧力から圧力を下げながら血管を加圧したときの血管の平均圧力に対する圧力振幅の変化に関する特性指標と、血管の固さに関するパラメータとの関係を求めておくモデル計算工程と、前記血管圧力計算工程で求められた前記平均圧力に対する前記圧力振幅の変化に関する特性指標を求め、その特性指標から、前記モデル計算工程で求められた関係に基づいて、血管の固さに関するパラメータを求める血管固さ算出工程とを、有することを特徴とする血管固さ評価方法。前記血管の固さに関するパラメータはスティフネスパラメータから成り、前記モデル計算工程は前記関係を求めるとき、血管の弾性を表すチューブ則を利用して、血管の平均圧力に対する圧力振幅の変化を求めて特性指標を求めるとともに、1式によりスティフネスパラメータを求めることを特徴とする請求項4記載の血管固さ評価方法。β:スティフネスパラメータP:血管の内圧Pb:血管の任意の基準内圧A:血管の断面積A0:P=Pbのときの血管の断面積前記血管固さ算出工程は、前記平均圧力に対する前記圧力振幅の変化を、前記圧力振幅のピーク値とそのときの平均圧力の値とで無次元化して求め、前記特性指標として、前記圧力振幅のピーク値を与える平均圧力の半分の圧力のときの二階微分値、および/または、前記圧力振幅が0から変化し始めたときの一階微分値を求めることを特徴とする請求項4または5記載の血管固さ評価方法。コンピュータを、圧力を上げながら、または、所定の圧力から圧力を下げながら血管を加圧して測定された前記血管の圧力を受信する受信手段、前記受信手段で受信された前記血管の圧力の所定の圧力段階毎に、前記血管の圧力の平均圧力と脈波による圧力振幅とを求める血管圧力計算手段、あらかじめモデル計算により求められた、圧力を上げながら、または、所定の圧力から圧力を下げながら血管を加圧したときの血管の平均圧力に対する圧力振幅の変化に関する特性指標と、血管の固さに関するパラメータとの関係を記憶しておくモデル記憶手段、前記血管圧力計算手段で求められた前記平均圧力に対する前記圧力振幅の変化に関する特性指標を求め、その特性指標から、前記モデル記憶手段で記憶された関係に基づいて、血管の固さに関するパラメータを求める血管固さ算出手段、として機能させるための血管固さ評価プログラム。コンピュータを、圧力を上げながら、または、所定の圧力から圧力を下げながら血管を加圧して測定された前記血管の圧力を受信する受信手段、前記受信手段で受信された前記血管の圧力の所定の圧力段階毎に、前記血管の圧力の平均圧力と脈波による圧力振幅とを求める血管圧力計算手段、あらかじめモデル計算により求められた、圧力を上げながら、または、所定の圧力から圧力を下げながら血管を加圧したときの血管の平均圧力に対する圧力振幅の変化に関する特性指標と、血管の固さに関するパラメータとの関係を記憶しておくモデル記憶手段、前記血管圧力計算手段で求められた前記平均圧力に対する前記圧力振幅の変化に関する特性指標を求め、その特性指標から、前記モデル記憶手段で記憶された関係に基づいて、血管の固さに関するパラメータを求める血管固さ算出手段、として機能させるための血管固さ評価プログラムを記録したコンピュータ読み取り可能な記録媒体。しかしながら、特許文献1に記載の装置では、脈波速度を測定する際、2箇所で脈波を測定する必要があるため、装置が大型化してしまうという課題があった。特許文献2に記載の装置では、超音波診断装置が高価であり、使用時には専門知識が必要となるという課題があった。特許文献3に記載の装置では、各脈波から血管の固さ指標を求めるための演算が複雑で、演算規模も膨大となるという課題があった。特許文献4に記載の装置では、心音測定が必要であるため、装置が大型化してしまうという課題があった。また、心臓と脈波測定部位との間の脈波伝播速度情報を使用して、その間の平均的な血管固さを推定するものであり、特定の部位の血管固さを推定することはできないという課題があった。特許文献3および4に記載の装置では、独自の血管固さ指標を用いているため、血管の固さ指標として一般的に用いられるスティフネスパラメータとの関連性が不明であり、血管の固さの評価が困難であるという課題があった。本発明は、このような課題に着目してなされたもので、安価かつ小型化が可能で、特定部位の血管固さを一般家庭でも容易に測定して評価することができる血管固さ評価装置、血管固さ評価方法、血管固さ評価プログラムおよび、血管固さ評価プログラムを記録したコンピュータ読み取り可能な記録媒体を提供することを目的としている。本発明者等は、非特許文献1乃至3に記載のように、測定対象部位における血管系のモデル、および、所定の圧力で血管を押し込んだときの血圧を計算するための皮下組織モデルに基づいて、所定の時間間隔で段階的に圧力を上げながら血管を加圧したときの血圧変動についてモデル計算を行った。その結果、圧力段階毎の平均圧力と圧力振幅との関係が、血管の固さに起因して変化することを見出し、本発明に至った。すなわち、本発明に係る血管固さ評価装置は、圧力を上げながら、または、所定の圧力から圧力を下げながら血管を加圧する加圧手段と、前記加圧手段で加圧中の前記血管の圧力を測定する測定手段と、前記加圧手段での所定の圧力段階毎に、前記測定手段で測定された前記血管の圧力の平均圧力と脈波による圧力振幅とを求める血管圧力計算手段と、あらかじめモデル計算により、圧力を上げながら、または、所定の圧力から圧力を下げながら血管を加圧したときの血管の平均圧力に対する圧力振幅の変化に関する特性指標と、血管の固さに関するパラメータとの関係を求めておくモデル計算手段と、前記血管圧力計算手段で求められた前記平均圧力に対する前記圧力振幅の変化に関する特性指標を求め、その特性指標から、前記モデル計算手段で求められた関係に基づいて、血管の固さに関するパラメータを求める血管固さ算出手段とを、有することを特徴とする。本発明に係る血管固さ評価方法は、圧力を上げながら、または、所定の圧力から圧力を下げながら血管を加圧する加圧工程と、前記加圧工程で加圧中の前記血管の圧力を測定する測定工程と、前記加圧工程での所定の圧力段階毎に、前記測定工程で測定された前記血管の圧力の平均圧力と脈波による圧力振幅とを求める血管圧力計算工程と、あらかじめモデル計算により、圧力を上げながら、または、所定の圧力から圧力を下げながら血管を加圧したときの血管の平均圧力に対する圧力振幅の変化に関する特性指標と、血管の固さに関するパラメータとの関係を求めておくモデル計算工程と、前記血管圧力計算工程で求められた前記平均圧力に対する前記圧力振幅の変化に関する特性指標を求め、その特性指標から、前記モデル計算工程で求められた関係に基づいて、血管の固さに関するパラメータを求める血管固さ算出工程とを、有することを特徴とする。本発明に係る血管固さ評価プログラムは、コンピュータを、圧力を上げながら、または、所定の圧力から圧力を下げながら血管を加圧して測定された前記血管の圧力を受信する受信手段、前記受信手段で受信された前記血管の圧力の所定の圧力段階毎に、前記血管の圧力の平均圧力と脈波による圧力振幅とを求める血管圧力計算手段、あらかじめモデル計算により求められた、圧力を上げながら、または、所定の圧力から圧力を下げながら血管を加圧したときの血管の平均圧力に対する圧力振幅の変化に関する特性指標と、血管の固さに関するパラメータとの関係を記憶しておくモデル記憶手段、前記血管圧力計算手段で求められた前記平均圧力に対する前記圧力振幅の変化に関する特性指標を求め、その特性指標から、前記モデル記憶手段で記憶された関係に基づいて、血管の固さに関するパラメータを求める血管固さ算出手段、として機能させることを特徴とする。本発明に係る血管固さ評価プログラムを記録したコンピュータ読み取り可能な記録媒体は、コンピュータを、圧力を上げながら、または、所定の圧力から圧力を下げながら血管を加圧して測定された前記血管の圧力を受信する受信手段、前記受信手段で受信された前記血管の圧力の所定の圧力段階毎に、前記血管の圧力の平均圧力と脈波による圧力振幅とを求める血管圧力計算手段、あらかじめモデル計算により求められた、圧力を上げながら、または、所定の圧力から圧力を下げながら血管を加圧したときの血管の平均圧力に対する圧力振幅の変化に関する特性指標と、血管の固さに関するパラメータとの関係を記憶しておくモデル記憶手段、前記血管圧力計算手段で求められた前記平均圧力に対する前記圧力振幅の変化に関する特性指標を求め、その特性指標から、前記モデル記憶手段で記憶された関係に基づいて、血管の固さに関するパラメータを求める血管固さ算出手段、として機能させるための血管固さ評価プログラムを記録していることを特徴とする。本発明に係る血管固さ評価装置、血管固さ評価方法、血管固さ評価プログラムおよび、血管固さ評価プログラムを記録したコンピュータ読み取り可能な記録媒体は、血管の固さに関するパラメータを求めることにより、血管の固さを評価することができる。血管の固さを評価するために、測定対象位置でのみ血管の加圧および圧力測定を行えばよく、複数箇所での測定が不要であるため、小型化が可能である。血管の加圧および圧力測定には市販のカフを用いることができ、高価な装置や機器を必要とせず、安価である。また、カフを用いることで、一般家庭でも容易に測定することができる。このように、本発明に係る血管固さ評価装置、血管固さ評価方法、血管固さ評価プログラムおよび、血管固さ評価プログラムを記録したコンピュータ読み取り可能な記録媒体によれば、特定部位の血管固さを一般家庭でも容易に測定して評価することができる。本発明に係る血管固さ評価装置、血管固さ評価方法、血管固さ評価プログラムおよび、血管固さ評価プログラムを記録したコンピュータ読み取り可能な記録媒体で、圧力段階毎の平均圧力に対する圧力振幅の変化に関する特性指標は、血管の固さが反映されていると考えられるものであれば、いかなるものであってもよい。特性指標は、モデル計算により、血管の固さを様々に変えて、血管の平均圧力に対する圧力振幅の変化を求めたとき、血管の固さに大きく影響されるものを用いることが好ましい。血管の加圧は、所定の時間間隔で段階的に圧力を上げながら下げながら行ってもよく、連続的に圧力を上げながら下げながら行ってもよい。カフを用いる場合、連続的に圧力を変化させるのが容易である。本発明に係る血管固さ評価装置、血管固さ評価方法、血管固さ評価プログラムおよび、血管固さ評価プログラムを記録したコンピュータ読み取り可能な記録媒体で、血管の固さに関するパラメータは、その評価が容易になるよう、スティフネスパラメータから成ることが好ましい。また、特性指標と血管の固さに関するパラメータとの関係は、モデル計算により直接求めてもよく、特性指標と血管の固さに関するパラメータとの間に他のパラメータを介して、間接的に求めてもよい。モデル計算に用いるモデルは、所定の時間間隔で段階的に圧力を上げながら血管を加圧したときの、血圧変動の実測値とモデル計算値とがほぼ一致しているものであればよく、例えば、非特許文献1乃至3に記載のモデルのように、測定対象部位における血管系のモデルと、所定の圧力で血管を押し込んだときの血圧を計算するための皮下組織モデルとを組み合わせて使用してもよい。本発明に係る血管固さ評価装置、血管固さ評価プログラムおよび、血管固さ評価プログラムを記録したコンピュータ読み取り可能な記録媒体で、前記血管の固さに関するパラメータはスティフネスパラメータから成り、前記モデル計算手段は前記関係を求めるとき、血管の弾性を表すチューブ則を利用して、血管の平均圧力に対する圧力振幅の変化を求めて特性指標を求めるとともに、1式によりスティフネスパラメータを求めることが好ましい。β:スティフネスパラメータP:血管の内圧Pb:血管の任意の基準内圧A:血管の断面積A0:P=Pbのときの血管の断面積本発明に係る血管固さ評価方法で、前記血管の固さに関するパラメータはスティフネスパラメータから成り、前記モデル計算工程は前記関係を求めるとき、血管の弾性を表すチューブ則を利用して、血管の平均圧力に対する圧力振幅の変化を求めて特性指標を求めるとともに、1式によりスティフネスパラメータを求めることが好ましい。このチューブ則を利用する場合、血管の断面積比A/A0と伸展圧φとの関係がチューブ則で表されるものとして、モデル計算を行う。ここで、チューブ則は、2式で表される。ここで、以下で使用する基準値として、AC/A0=0.3C1=0.3、C2=0.05C3~C5:AC/A0=0.3において、2つの式のφの傾きが一致するように設定するφ0:A/A0=1のときの伸展圧=13.3kPan1=5血管の固さにより変化させるn2=6とする。また、生理的状態にある血管では、内圧Pと血管の半径Rとの関係は、3式で表されることが知られている。Rb:P=Pbのときの血管の半径生理的状態にある血管では血管の断面が円形であると仮定すると、2式および3式から、φ=Pとして、スティフネスパラメータβは1式となる。1式から、スティフネスパラメータβは、図1に示すチューブ則を表す曲線のA/A0=1のときの傾きを利用して求めることができ、チューブ則のパラメータn1の値により変化する。本発明に係る血管固さ評価装置、血管固さ評価プログラムおよび、血管固さ評価プログラムを記録したコンピュータ読み取り可能な記録媒体で、前記血管固さ算出手段は、前記平均圧力に対する前記圧力振幅の変化を、前記圧力振幅のピーク値とそのときの平均圧力の値とで無次元化して求め、前記特性指標として、前記圧力振幅のピーク値を与える平均圧力の半分の圧力のときの二階微分値、および/または、前記圧力振幅が0から変化し始めたときの一階微分値を求めることが好ましい。本発明に係る血管固さ評価方法で、前記血管固さ算出工程は、前記平均圧力に対する前記圧力振幅の変化を、前記圧力振幅のピーク値とそのときの平均圧力の値とで無次元化して求め、前記特性指標として、前記圧力振幅のピーク値を与える平均圧力の半分の圧力のときの二階微分値、および/または、前記圧力振幅が0から変化し始めたときの一階微分値を求めることが好ましい。この特性指標として二階微分値や一階微分値を用いる場合、平均圧力に対する圧力振幅の変化を、圧力振幅のピーク値とそのときの平均圧力の値とで無次元化することにより、個人差や測定器の誤差などによる影響を低減することができる。また、圧力振幅のピーク値を与える平均圧力の半分の圧力のときの二階微分値や、圧力振幅が0から変化し始めたときの一階微分値は、血管の固さを強く反映していると考えられるため、これらを特性指標として用いることにより、より正確に血管の固さを評価することができる。本発明によれば、安価かつ小型化が可能で、特定部位の血管固さを一般家庭でも容易に測定して評価することができる血管固さ評価装置、血管固さ評価方法、血管固さ評価プログラムおよび、血管固さ評価プログラムを記録したコンピュータ読み取り可能な記録媒体を提供することができる。以下、図面に基づき、本発明の実施の形態について説明する。図2乃至図11は、本発明の実施の形態の血管固さ評価装置および血管固さ評価方法を示している。本発明の実施の形態の血管固さ評価方法は、本発明の実施の形態の血管固さ評価装置により好適に実施される。以下、図2乃至図11に従って、本発明の実施の形態の血管固さ評価装置および血管固さ評価方法について、実施例を適宜参照しながら説明する。[加圧手段および測定手段]図2に示すように、本発明の実施の形態の血管固さ評価装置10は、加圧手段11と測定手段12とコンピュータ13とモデル計算手段14とを有している。加圧手段11は、センサ押下げ装置から成り、アンプ15を介してコンピュータ13に接続されている。測定手段12は、圧力センサから成り、コンピュータ13に接続されている。測定手段12は、加圧手段11に取り付けられている。加圧手段11は、橈骨動脈などの血管を、測定手段12を介して加圧可能になっている。測定手段12は、加圧手段11で加圧したときの血管の圧力を測定可能になっている。加圧手段11および測定手段12は、腕などに巻いて使用される市販のカフから成っていてもよい。[制御手段および受信手段]図2に示すように、コンピュータ13は、加圧手段11で加圧する圧力を制御する制御手段21と、測定手段12で測定された圧力を受信する受信手段22とを有している。制御手段21は、圧力を上げながら、または、所定の圧力から圧力を下げながら血管を加圧するよう、センサ押下げ信号を加圧手段11に送って、加圧手段11を制御可能になっている。制御手段21は、所定の時間間隔で段階的に圧力を上げながら下げながら血管を加圧するよう制御してもよく、連続的に圧力を上げながら下げながら血管を加圧するよう制御してもよい。受信手段22は、加圧手段11で血管を加圧したとき、測定手段12で測定された血管の圧力を、バンドパスフィルタBPFを介して受信可能になっている。この測定例として、5秒間隔で、0.5mmずつ血管を押し下げて、段階的に圧力を上げながら血管を加圧したときの、血管の圧力の測定結果を、図3aに示す。[モデル計算手段]また、図2に示すように、コンピュータ13は、モデル記憶手段23を有している。モデル記憶手段23は、メモリから成り、あらかじめモデル計算手段14によりモデル計算された、特性指標と血管の固さに関するパラメータとの関係を記憶するようになっている。モデル計算手段14は、コンピュータ13とは異なる他のコンピュータで実施可能に構成されている。モデル計算手段14は、血管の固さに関するパラメータを複数通り仮定し、それぞれの場合について、圧力を上げながら、または、所定の圧力から圧力を下げながら血管を加圧したときの血管の圧力を計算し、所定の圧力段階毎に血管の平均圧力と圧力振幅とを求めるようになっている。さらに、モデル計算手段14は、血管の平均圧力に対する圧力振幅の変化から、血管の固さが反映されていると考えられる特性指標を、それぞれの血管の固さに関するパラメータ毎に求め、その特性指標と血管の固さに関するパラメータとの関係を求めるようになっている。血管の固さに関するパラメータは、スティフネスパラメータである。なお、モデル計算手段14は、コンピュータ13で実施可能であってもよい。[モデル計算例-使用モデル]橈骨動脈の血管の固さを評価するためにモデル計算手段14で用いるモデルの例を、図4に示す。なお、図4に示すモデルは、非特許文献1乃至3に記載のモデルと同じものである。図4aに示すように、腕部血管系の数学モデルは、動脈系Arteryを全長700mmのテーパーのついたコラプシブルチューブで模擬し、チューブの上流端に供給圧力源Supplyが、下流端に毛細血管による微小循環micro-circulation、静脈Venousが接続されている。チューブの上流端半径Inletradiusを、4.23mm、下流端半径Outletradiusを、1.74mmとする。また、血液を、密度ρ=1050kg/m3、動粘度ν=3.8×10-6m2/sの非圧縮性ニュートン流体とし、血流は一次元で記述する。動脈への加圧および圧力測定は、チューブの下流端から100mmの位置とし、その位置での加圧時の皮下組織モデルTissueModelを、図4bに示す。図4bに示すように、皮下組織モデルとして二次の非線形バネモデルを用いると、圧力センサSensorの押し込み量Yとセンサ出力圧力P0との関係は、4式のようになる。なお、センサは、加圧面が直径8mmの円形である。Pe:動脈を圧迫する圧力a:組織モデル定数dyav:センサ加圧面における、血管の初期形状の高さy0からの変化量dyの空間平均値また、モデルで使用する心拍一周期分の供給圧力SupplypressurePsの波形を、図5に示す。ここで、平均圧Psavは13.3kPa100mmHg、振幅ΔPsは、5.5kPa41mmHgとした。図5に示す波形は、ヒトの腕頭動脈の圧力の実測値を基に作成されたものである。[モデル計算例-計算結果]図4に示すモデルに対して、図5に示す供給圧力や、2式に示すチューブ則を使用してモデル計算を行った。モデル計算では、非特許文献1乃至3で実施されたモデル計算と同様にして、計算を行っている。モデル計算では、まず、チューブ則のパラメータn1、組織モデル定数a、供給圧力の平均圧Psavおよび振幅ΔPsを様々に変化させながら、図3aに示す血管の圧力の実測値に最も良く合う各パラメータの値を求めた。各パラメータの値を変化させる範囲は、現実的な値としてそれぞれ、n1=2~12、a=3.667×108~16.00×108Pa/m2、Psav=8.3~21.3kPa、ΔPs=3.5~9.5kPaとした。図3aの実測値に最も良く合うモデル計算結果を、図3bに示す。また、このときの各圧力段階毎の平均圧力Poavおよび振幅ΔP0の、実測値と計算値との比較を、それぞれ図3cおよび図3dに示す。このときの各パラメータの値は、n1=8、a=11.0×108Pa/m2、Psav=8.3kPa、ΔPs=4.5kPaであった。図3に示すように、血圧変動の実測値と図4に示すモデルに基づくモデル計算値とがほぼ一致しており、血管の固さを評価するために、図4に示すモデルを使用可能である考えられる。次に、1式および図1に示すように、スティッフネスパラメータβと相関を有するチューブ則のパラメータn1をさまざまに変化させてモデル計算を行い、その結果を図6に示す。このときの、他のパラメータは、a=10.0×108Pa/m2、Psav=13.3kPa、ΔPs=5.5kPaとした。図6aに示すように、Y≦8.0mmの範囲では、n1が大きいほど平均圧力Poavが大きくなり、Y≧8.0mmの範囲では、Poavに違いが認められないことが確認された。また、図6bに示すように、Y≦5.0mmの範囲では、n1が大きいほど振幅ΔP0が小さくなることが確認された。これは、血管が軟らかいn1が小さいときは、センサの押し込みによって血管が変形しやすく、小さい押し込み量であっても信号が検出されるのに対し、血管が硬いn1が大きいときは、センサの押し込みに対して、血管ではなく皮下組織が変形するため、血管が変形するほど押し込まなければ信号が検出されないためであると考えられる。[特性指標の計算]図6に示すように、血管への押し込み量Yが小さい範囲で、n1によるPoavの変化が顕著である。個人差や測定器の誤差などによる影響を低減するために、平均圧力Poavに対する圧力振幅ΔP0の変化を求め、圧力振幅ΔP0のピーク値とそのときの平均圧力Poavの値とで無次元化してグラフにプロットする。これにより、圧力振幅ΔP0がピーク値ΔP0=1.0のときの平均圧力PoavPoav=1.0よりも平均圧力Poavが小さい範囲Poav<1.0で、n1の値によって、グラフの曲線の形状が大きく変化する。具体的には、n1の値が小さいときはグラフの曲線が上に凸となり、n1の値が大きいときはグラフの曲線が下に凸となる。このことから、n1の値を特徴付ける特性指標として、圧力振幅ΔP0のピーク値を与える平均圧力Poavの半分の圧力Poav=0.5のときの二階微分値ΔP、および/または、圧力振幅ΔP0が0から変化し始めたときの一階微分値ΔP’を使用するのが適切であると考えられる。[特性指標とスティフネスパラメータとの関係]こうして求められた特性指標の例を、図7に示す。図7aは、無次元化した平均圧力Poavに対する圧力振幅ΔP0の変化を示すグラフ、図7bおよびcはそれぞれ、ΔPとスティフネスパラメータβとの関係、および、ΔP’とスティフネスパラメータβとの関係を示すグラフである。図7aは、n1=5、a=10.0×108Pa/m2、Psav=21.3kPa、ΔPs=5.5kPaのときのグラフである。また、図7bおよびcにプロットされたデータ図中の印は、チューブ則のパラメータn1、組織モデル定数a、供給圧力の平均圧Psavおよび振幅ΔPsを様々に変化させたときのデータである。各パラメータの値を変化させる範囲は、図3bを求めたときの範囲と同じである。スティフネスパラメータβは、n1の値から、1式およびチューブ則を使用して求めている。図7bおよびcに示すように、ΔPおよびΔP’は、スティフネスパラメータβとの相関が認められる。ΔPとスティフネスパラメータβとの関係が一次式で表されるとしたときの直線を、図7bに示す。また、ΔP’とスティフネスパラメータβとの関係が二次式で表されるとしたときの曲線を、図7cに示す。図7bの直線および図7cの曲線は、各データの残差二乗和が最も小さくなるように最小二乗法を利用して求めている。図7bおよび図7cの関係は、図中に示す直線や曲線に限られるものではなく、各データを説明可能なものであれば、他の直線や曲線で表されてもよい。図7bおよびcに示す関係を利用することにより、特性指標からスティフネスパラメータβを求めることができる。モデル記憶手段23は、このようにしてモデル計算手段14によりモデル計算された、特性指標とスティフネスパラメータとの関係を記憶するようになっている。[血管圧力計算手段、血管固さ算出手段および表示手段]図2に示すように、コンピュータ13は、さらに血管圧力計算手段24と血管固さ算出手段25と表示手段26とを有している。血管圧力計算手段24は、加圧手段11で血管を加圧したときの血管の圧力の測定値から、所定の圧力段階毎に平均圧力と圧力振幅とを分離して求めるようになっている。血管固さ算出手段25は、血管圧力計算手段24で求められた平均圧力と圧力振幅とを使用して、信号処理を行うよう構成されている。血管固さ算出手段25は、血管圧力計算手段24で求められた平均圧力に対する圧力振幅の変化に基づいて、特性指標を求めるようになっている。また、血管固さ算出手段25は、その特性指標から、モデル記憶手段23に記憶された特性指標とスティフネスパラメータとの関係に基づいて、スティフネスパラメータを求めるようになっている。表示手段26は、モニタから成り、血管固さ算出手段25により求められたスティフネスパラメータなどを表示可能になっている。なお、図8に示すように、本発明の実施の形態の血管固さ評価装置10で、加圧手段11がモータ51、滑車52、巻掛けベルト53などを有し、圧力センサから成る測定手段12を介して、血管を加圧可能になっていてもよい。この場合、制御手段21が、圧力を上げながら、または、所定の圧力から圧力を下げながら血管を加圧するよう、モータ駆動信号を加圧手段11に送って、加圧手段11を制御可能になっている。また、血管圧力計算手段24が、血管を加圧するときの測定手段12の押し込み量センサ変位に対する平均圧力および圧力振幅の変化を求め、それに基づいて、血管固さ算出手段25が、平均圧力に対する圧力振幅の変化を求めるようになっている。[血管固さ評価方法による実施]本発明の実施の形態の血管固さ評価装置10は、本発明の実施の形態の血管固さ評価方法により血管の固さの評価を行うことができる。なお、ここでは、圧力を上げながら血管を加圧した場合について説明する。まず、カフを手首に巻くなど、橈骨動脈などの血管を加圧可能かつその血管の圧力を測定可能に、加圧手段11および測定手段12をセットする。図9に示すように、圧力を上げながら血管を加圧するよう、制御手段21により加圧手段11による圧力を制御して、血管を加圧するステップ31。加圧中の血管の圧力を、脈波として測定手段12により測定するステップ32。測定された脈波から、血管圧力計算手段24により、圧力段階毎に平均圧力Poavと圧力振幅ΔP0とを求めるステップ33。脈波を測定した実施例として、カフを利用して、カフ内圧を連続的に変化させて血管を加圧したときの圧力カフ内圧と、測定圧力に対応するPVDFフィルムセンサの出力電圧とを、それぞれ図10aおよびbに示す。図10に示すように、加圧を始めた段階では、加圧するたびに圧力振幅ΔP0の値が徐々に増加するが、所定の圧力段階でΔP0が最大値ΔPmaxとなった後は、圧力振幅ΔP0の値は徐々に減少していく。血管固さの評価には、圧力振幅が最大値ΔPmaxとなるまでの測定データを使用するため、ΔPmaxが測定された時点で加圧を止めて、測定を終了してもよい。このため、図9に示すように、各圧力段階において、圧力振幅ΔP0を前の圧力段階までの圧力振幅の最大値ΔPmaxと比較し、ΔP0がΔPmaxを超えたとき、ΔPmax=ΔP0として次の圧力段階に加圧する。この加圧と比較とを、ΔP0がΔPmaxを超えなくなるまで繰り返し、ΔP0がΔPmaxを超えなくなったときに加圧を止めるステップ34。次に、血管固さ算出手段25により、平均圧力Poavに対する圧力振幅ΔP0の変化を求めてグラフにプロットする。このとき、圧力振幅ΔP0のピーク値ΔPmaxとそのときの平均圧力Poavの値とで無次元化してプロットするステップ35。段階的に圧力を上げながら血管を加圧したときの圧力段階毎に、平均圧力Poavと圧力振幅ΔPとを求めた結果を、図11aおよびbに示す。図11aおよびbは、同一人に対して、複数回測定を行った結果を示している。また、血管固さ算出手段25により、図10aおよびbから、平均圧力Poavに対する圧力振幅ΔPの変化を求めてグラフにプロットした結果を、図11cに示す。図11cに示すように、平均圧力Poavに対する圧力振幅ΔPの変化を求めることにより、測定毎の圧力振幅ΔPの誤差を低減することができるとともに、圧力振幅ΔPのピーク位置を揃えることもできる。グラフにプロットしたならば、図9に示すように、血管固さ算出手段25により、特性指標として、ΔPmaxを与える平均圧力Poavの半分の圧力Poav=0.5のときの二階微分値ΔPステップ36、および/または、圧力振幅ΔP0が0から変化し始めたときの一階微分値ΔP’を求めるステップ37。さらに、血管固さ算出手段25により、求められた特性指標ΔPおよび/またはΔP’から、モデル計算手段14で求められてモデル記憶手段23に記憶された特性指標とスティフネスパラメータとの関係に基づいて、スティフネスパラメータβを求めるステップ38。求められたスティフネスパラメータβなどを、表示手段26により表示する。こうして求められたスティフネスパラメータβに基づいて、血管の固さを評価することができる。なお、図9乃至図11では、圧力を上げながら血管を加圧した場合について説明を行ったが、所定の圧力から圧力を下げながら血管を加圧した場合であっても、図10および図11に示すような結果を得ることができ、同様にして血管の固さの評価を行うことができる。本発明の実施の形態の血管固さ評価装置10および血管固さ評価方法は、血管の固さを評価するために、測定対象位置でのみ血管の加圧および圧力測定を行えばよく、複数箇所での測定が不要であるため、小型化が可能である。血管の加圧および圧力測定には市販のカフを用いることができ、その場合、高価な装置や機器を必要とせず、安価である。また、カフを用いることにより、一般家庭でも容易に測定することができる。このように、本発明の実施の形態の血管固さ評価装置10および血管固さ評価方法によれば、特定部位の血管固さを一般家庭でも容易に測定して評価することができる。本発明の実施の形態の血管固さ評価装置10および血管固さ評価方法は、血管の固さを強く反映した特性指標ΔPおよび/またはΔP’を用いるため、より正確に血管の固さを評価することができる。血管の固さに関するパラメータとして、一般的に用いられているスティフネスパラメータを使用するため、血管の固さの評価が容易である。なお、モデル計算手段14により求める関係として、特性指標とチューブ則のパラメータn1との関係を求めてもよい。この場合、求められたn1から、1式およびチューブ則を用いてスティフネスパラメータを求めることができる。なお、本発明の実施の形態の血管固さ評価プログラムは、図2および図8に示すコンピュータ13を機能させるプログラムから成っていてもよい。この場合、例えば、カフを用いた一般家庭での容易な測定を達成する形態として、市販の家庭用血圧計に実装または接続されたコンピュータCPU、情報処理装置、各種端末を含む13が、図9に示すフローに従って、本発明の実施の形態の血管固さ評価プログラムを実行することによって実現されてもよい。本発明の実施の形態の血管固さ評価プログラムは、例えば、CDCD-ROM、CD-R、CD-RWなど、DVDDVD-ROM、DVD-RAM、DVD-R、DVD-RW、DVDR、DVDRWなど等のコンピュータ読み取り可能な記録媒体に記録された形態で提供されうる。この場合、コンピュータ13は、その記録媒体から血管固さ評価プログラムを読み取ってコンピュータの内部記憶装置または外部記憶装置に転送し格納して用いることができる。また、本発明の実施の形態の血管固さ評価プログラムを、例えば、磁気ディスク、光ディスク、光磁気ディスク等の記憶装置記録媒体に記録しておき、その記憶装置から通信回線を介してコンピュータ13に提供するようになっていてもよい。ここで、コンピュータとは、ハードウェアとOSオペレーティングシステムとを含む概念であり、OSの制御の下で動作するハードウェアを意味している。また、OSが不要で、アプリケーションプログラム単独でハードウェアを動作させるような場合には、そのハードウェア自体がコンピュータに相当する。ハードウェアは、少なくとも、CPU等のマイクロプロセッサと、記録媒体に記録されたコンピュータプログラムを読み取るための手段とを備えている。本発明の実施の形態の血管固さ評価プログラムとしてのアプリケーションプログラムは、上述のようなコンピュータに実現させるプログラムコードを含んでいる。また、その機能の一部は、アプリケーションプログラムではなくOSによって実現されてもよい。なお、本発明の実施の形態のコンピュータ読み取り可能な記録媒体としては、上述したフレキシブルディスク、CD、DVD、磁気ディスク、光ディスク、光磁気ディスクのほか、コンピュータの内部記憶装置RAMやROMなどのメモリ、外部記憶装置等や、バーコードなどの符号が印刷された印刷物等の、コンピュータ読み取り可能な種々の媒体を利用することもできる。