シミュレーション装置、シミュレーション方法、及びシミュレーションシステム
【課題】
複数の動物から構成される牧場等の経営に関するシミュレーションが可能なシミュレーション装置を提供すること。
【解決手段】
複数の動物が存在する所定の動物群内における、分娩の経験がある経産動物の数、分娩の経験のない非経産動物の数の、所定の期間における変動を算定するシミュレーション装置であって、非経産動物変化数算定手段と、初妊動物変化数算定手段と、経産動物変化数算定手段と、出生動物数算定手段と、期末非経産動物数算定手段とを備えることを特徴とするシミュレーション装置に関する。
複数の動物から構成される牧場等の経営に関するシミュレーションが可能なシミュレーション装置を提供すること。
【解決手段】
複数の動物が存在する所定の動物群内における、分娩の経験がある経産動物の数、分娩の経験のない非経産動物の数の、所定の期間における変動を算定するシミュレーション装置であって、非経産動物変化数算定手段と、初妊動物変化数算定手段と、経産動物変化数算定手段と、出生動物数算定手段と、期末非経産動物数算定手段とを備えることを特徴とするシミュレーション装置に関する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の動物から構成される牧場等の経営に関するシミュレーションが可能なシミュレーション装置に関する。
【背景技術】
【０００２】
従来から、酪農経営のシミュレーション分析が試みられている。通常は、牛の頭数や個体の乳量（乳牛生産性、乳牛１頭あたりの乳量能力）に基づき試算されることが一般的であり、このような分析を行なう際には乳牛の泌乳力が重視されている。個体の乳量が高い場合は、乳牛の泌乳力が直接的に売上高のアップに寄与することから、泌乳力を用いた分析が広く用いられている。特に、生産現場にて飼料を与える際においても泌乳力の向上を目標にした設計が重んじられ、金融機関への融資の申請においても泌乳力に基づいた経営計画やシミュレーション分析が適用されている。
【０００３】
このため、酪農経営者による飼養管理の優先順位は自然と泌乳力向上が一番におかれていた。しかし、乳牛が持つ本来の生理的機能には、泌乳のほかに繁殖、自己成長、免疫などがあり、これらの事項は十分に考慮されていなかった。特に繁殖については、分娩され出生した雌子牛が、出生後に妊娠して子牛を分娩するまで、２年ほどの期間がかかるため、これらの期間の経過を考慮する必要がある。
【０００４】
このような中、牛群を管理する手法の一つとして、乳牛の栄養診断処理、牛群の健康管理、及び、繁殖管理処理をコンピュータにより経時的に行なうために、乳牛の乳成分データを基に演算処理し、演算処理されたデータにもとづいて摂取栄養を解析し、解析されたデータにもとづいてプログラム処理して乳牛の栄養診断処理を行なう乳牛の管理システムが開示されている（例えば、特許文献１参照）。また、その他、ミルクパイプラインのミルクタップ側に被識別部を設け、かつ搾乳ユニットのディストリビュータ側に、被識別部を検出して牛床を識別する識別部と、識別データ及び搾乳データを記憶する記憶部を有する制御ユニットを設けるとともに、データを受信して記憶する送受信ユニットを備える搾乳管理システムが開示されている（例えば、特許文献２参照）。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平７−１９４２７３号公報
【特許文献２】特開２００１−４５８９８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかし、特許文献１及び２はいずれも乳牛の管理をするものではあるが、牧場経営全体を管理するものではなく、将来にわたる牧場の成長を予測し、牧場を経営するうえでの目標を達成するための指針を提示できるといったものではない。本発明は、乳牛等の畜産動物の繁殖機能に基づき動物群評価及び牧場の経営評価を可能とするもので、畜産動物の飼養期間（時間の流れ）を考慮して、畜産動物の出生、育成、妊娠、初産のプロセスを経営シミュレーションに組み込み、動物群の繁殖力と供用力を見出して、将来にわたる牧場の成長の潜在性を計数化、指標化するためのシミュレーション装置を提供することを目的とする。
【課題を解決するための手段】
【０００７】
上記の課題を解決するため、本発明は、複数の動物が存在する所定の動物群内における、分娩の経験がある経産動物の数、分娩の経験がなく妊娠もしていない非経産動物の数、および初めての妊娠をしている初妊動物の数についての所定の期間における変動を算定するシミュレーション装置であって、前記所定の期間における期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに該動物群における分娩間隔を記憶する記憶手段と、記憶手段により記憶された期首の非経産動物の数をもとに、前記所定の期間内に妊娠を経て初めての分娩を行なう第１の初産動物の数、前記所定の期間内に初妊動物となる期末の初妊動物の数、および前記所定の期間内に変化しない第１の非経産動物の数を算定する非経産動物変化数算定手段と、記憶手段により記憶された期首の初妊動物の数をもとに、前記所定の期間内に初めての分娩を行なう第２の初産動物の数を算定する初妊動物変化数算定手段と、記憶手段により記憶された期首の経産動物の数、非経産動物変化数算定手段により算定された第１の初産動物の数および初妊動物変化数算定手段により算定された第２の初産動物の数をもとに、期末の経産動物の数を算定する経産動物変化数算定手段と、非経産動物変化数算定手段により算定された第１の初産動物の数、初妊動物変化数算定手段により算定された第２の初産動物の数、ならびに記憶手段により記憶された期首の経産動物の数および分娩間隔をもとに、前記所定の期間内に産まれうる出生動物の数を算定する出生動物数算定手段と、非経産動物変化数算定手段により算定された第１の非経産動物の数、および出生動物数算定手段により算定された出生動物の数をもとに、前記所定の期間における期末の非経産動物の数を算定する期末非経産動物数算定手段とを備えることを特徴とするシミュレーション装置に関する。
【０００８】
本発明によれば、前記所定の期間における期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに該動物群における分娩間隔をもとに、期末の経産動物の数、期末の非経産動物の数、期末の初妊動物の数の算定が可能となり、動物群の繁殖機能を考慮した動物群評価が可能となる。
【０００９】
本発明では、記憶手段は、前記所定の期間内に期首の非経産動物が初妊動物になる割合である初妊動物率、前記所定の期間内に非経産動物が妊娠を経て初産動物になる割合である第１の初産動物率、および前記所定の期間内に期首の初妊動物が初産動物になる割合である第２の初産動物率とをさらに記憶し、非経産動物変化数算定手段は、記憶手段により記憶された初妊動物率および期首の非経産動物の数をもとに、期末の初妊動物の数を算定し、かつ記憶手段により記憶された第１の初産動物率および期首の非経産動物の数をもとに、第１の初産動物の数を算定し、初妊動物変化数算定手段は、記憶手段により記憶された第２の初産動物率および期首の初妊動物の数をもとに、第２の初産動物の数を算定することが好ましい。
【００１０】
ここでは、それぞれの動物群に固有の特性である、初妊動物率、第１の初産動物率および第２の初産動物率が考慮されることとなる。すなわち、期末の初妊動物の数は初妊動物率および期首の非経産動物の数をもとに算定され、第１の初産動物の数は第１の初産動物率および期首の非経産動物の数をもとに算定され、さらに第２の初産動物の数は、第２の初産動物率および期首の初妊動物の数をもとに算定される。したがって、それぞれの動物群に固有の特性を考慮した動物群評価が可能となり、より現実の牧場等の経営に近い経営シミュレーションが可能となる。
【００１１】
本発明では、記憶手段は、経産動物のうち搾乳が可能な動物の割合である搾乳動物率をさらに記憶し、上記シミュレーション装置は、さらに、記憶手段により記憶された期首の経産動物の数、経産動物変化数算定手段により算定された期末の経産動物の数および搾乳動物率をもとに、搾乳が可能な搾乳動物の数を算定する搾乳動物数算定手段と、搾乳動物数算定手段により算定された搾乳動物の数をもとに、前記所定の期間内において出荷が可能な乳量を算定する乳量算定手段とを備えることが好ましい。
【００１２】
このような構成とすることで、期首の経産動物の数、期末の経産動物の数および搾乳動物率をもとに、搾乳が可能な搾乳動物の数が算定され、さらに搾乳動物の数をもとに、所定の期間内において出荷が可能な乳量が算定されるので、乳牛の繁殖機能を考慮した牛群評価をもとに、牛群の乳量の供用力を評価することが可能となり、より現実の牧場等の経営に近い経営シミュレーションが可能となる。
【００１３】
本発明では、上記シミュレーション装置は、さらに、前記所定の期間の整数倍であるシミュレーション期間についての利用者からの入力を受け付ける入力受付手段と、記憶手段に記憶された期首の経産動物の数、期首の非経産動物の数、および期首の初妊動物の数を、それぞれ、経産動物変化数算定手段により算定された期末の経産動物の数、期末非経産動物数算定手段により算定された期末の非経産動物の数、および非経産動物変化数算定手段により算定された期末の初妊動物の数で更新する更新手段とを備え、更新された期首の経産動物の数、期首の非経産動物の数、および期首の初妊動物の数をもととした、非経産動物変化数算定手段、初妊動物変化数算定手段、経産動物変化数算定手段、出生動物数算定手段および期末非経産動物数算定手段による算定、および更新手段による更新を、前記整数に対応する回数繰り返すことにより、シミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、および期末の初妊動物の数を算定することが好ましい。
【００１４】
このように構成することで、利用者が入力したシミュレーション期間（例えば、所定の期間が１年であり、シミュレーション期間が１０年）をもとに、期末の経産動物の数、期末の非経産動物の数、および期末の初妊動物の数の算定をシミュレーション期間の年数だけ繰り返すことにより、シミュレーション期間後における経産動物の数、非経産動物の数および期末の初妊動物の数を算定することが可能となる。したがって、利用者が入力した所望の期間が経過した後の動物群を評価することができ、長期にわたる牧場等の成長の潜在性を計数化、指標化することが可能となる。
【００１５】
本発明では、前記分娩間隔は、仮に設定された仮分娩間隔であって、上記シミュレーション装置は、さらに、利用者の入力に従って、シミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、または期末の初妊動物の数についての目標数を受け付ける目標数受付手段と、仮分娩間隔をもとに、非経産動物変化数算定手段、初妊動物変化数算定手段、経産動物変化数算定手段、出生動物数算定手段および期末非経産動物数算定手段による算定、および更新手段による更新を、前記整数に対応する回数繰り返すことにより、シミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、および期末の初妊動物の数を算定する仮算定手段と、仮算定手段により算定されたシミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、または期末の初妊動物の数が、目標数受付手段により受け付けられた目標数に達しているかについて判定する目標達成可否判定手段と、目標達成可否判定手段により目標数に達していないと判定された場合は、より短く設定された仮分娩間隔を新たに設定し、目標数に達していると判定された場合は、より長く設定された仮分娩間隔を新たに設定する仮分娩間隔設定手段と、仮分娩間隔設定手段により設定された新たな仮分娩間隔をもとに、仮算定手段によりシミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、および期末の初妊動物の数を算定し、仮分娩間隔が目標達成可否判定手段により目標数に達していると判定される最長の分娩間隔に達するまで仮分娩間隔設定手段による設定、仮算定手段による算定、目標達成可否判定手段による判定を繰り返し、目標数を達成するための最適な分娩間隔を特定する最適分娩間隔特定手段とを備えることが好ましい。
【００１６】
このような構成することで、目標数を達成するために最適な分娩間隔が特定されるため、利用者に対し、さらに有効な牧場経営の指標を提供することが可能となる。
【００１７】
本発明はさらに、複数の動物が存在する所定の動物群内における、分娩の経験がある経産動物の数、分娩の経験がなく妊娠もしていない非経産動物の数、および初めての妊娠をしている初妊動物の数についての所定の期間における変動をコンピュータ装置に算定させるシミュレーション方法であって、前記所定の期間における期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに該動物群における経産動物の分娩間隔を記憶領域に記憶させるステップと、記憶領域に記憶された期首の非経産動物の数をもとに、前記所定の期間内に妊娠を経て初めての分娩を行なう第１の初産動物の数、前記所定の期間内に初妊動物となる期末の初妊動物の数、および前記所定の期間内に変化しない第１の非経産動物の数を算定させるステップと、記憶領域に記憶された期首の初妊動物の数をもとに、前記所定の期間内に初めての分娩を行なう第２の初産動物の数を算定させるステップと、記憶領域に記憶された期首の経産動物の数、算定された第１の初産動物の数および算定された第２の初産動物の数をもとに、期末の経産動物の数を算定させるステップと、算定された第１の初産動物の数、算定された第２の初産動物の数、ならびに記憶領域に記憶された期首の経産動物の数および分娩間隔をもとに、前記所定の期間内に産まれうる出生動物の数を算定させるステップと、算定された第１の非経産動物の数、および算定された出生動物の数をもとに、前記所定の期間における期末の非経産動物の数を算定させるステップとから構成されることを特徴とするシミュレーション方法に関する。
【００１８】
本発明はさらに、複数の動物が存在する所定の動物群内における、分娩の経験がある経産動物の数、分娩の経験がなく妊娠もしていない非経産動物の数、および初めての妊娠をしている初妊動物の数についての所定の期間における変動を算定するサーバ装置と、通信ネットワークを介して前記サーバ装置と通信が可能な端末装置とを備えるシミュレーションシステムであって、端末装置が、利用者の入力に従って、前記所定の期間における期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに該動物群における分娩間隔の入力を受け付ける入力受付手段と、入力受付手段により受け付けられた期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに分娩間隔をサーバ装置に送信する情報送信手段とを備え、サーバ装置が、端末装置から送信された期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに分娩間隔を受信する情報受信手段と、情報受信手段により受信された期首の非経産動物の数をもとに、前記所定の期間内に妊娠を経て初めての分娩を行なう第１の初産動物の数、前記所定の期間内に初妊動物となる期末の初妊動物の数、および前記所定の期間内に変化しない第１の非経産動物の数を算定する非経産動物変化数算定手段と、情報受信手段により受信された期首の初妊動物の数をもとに、前記所定の期間内に初めての分娩を行なう第２の初産動物の数を算定する初妊動物変化数算定手段と、情報受信手段により受信された期首の経産動物の数、非経産動物変化数算定手段により算定された第１の初産動物の数および初妊動物変化数算定手段により算定された第２の初産動物の数をもとに、期末の経産動物の数を算定する経産動物変化数算定手段と、非経産動物変化数算定手段により算定された第１の初産動物の数、初妊動物変化数算定手段により算定された第２の初産動物の数、ならびに情報受信手段により受信された期首の経産動物の数および分娩間隔をもとに、前記所定の期間内に産まれうる出生動物の数を算定する出生動物数算定手段と、非経産動物変化数算定手段により算定された第１の非経産動物の数、および出生動物数算定手段により算定された出生動物の数をもとに、前記所定の期間における期末の非経産動物の数を算定する期末非経産動物数算定手段とを備えることを特徴とするシミュレーションシステムに関する。
【図面の簡単な説明】
【００１９】
【図１】本発明の実施の形態にかかるシミュレーションシステムの構成を表すブロック図である。
【図２】本発明の実施の形態にかかる牛群の飼養動態プロセスを表す図である。
【図３】本発明の実施の形態にかかる変化率テーブルを表す図である。
【図４】本発明の実施の形態にかかる牛群評価テーブルを表す図である。
【図５】本発明の実施の形態にかかる単価テーブルを表す図である。
【図６】本発明の実施の形態にかかる経営評価テーブルを表す図である。
【図７】本発明の実施の形態にかかる牛群の評価処理についてのフローチャートの一例である。
【図８】本発明の実施の形態にかかる最適分娩間隔の評価処理についてのフローチャートの一例である。
【図９】本発明の実施の形態にかかる目標達成可否判定処理を表わすフローチャートの一例である。
【図１０】本発明の実施の形態にかかる端末の表示画面の変化を表す図である。
【発明を実施するための形態】
【００２０】
図１は、本発明の実施の形態にかかるシミュレーションシステムの構成を表すブロック図である。シミュレーション装置１は、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）１１、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）１２、ハードディスク・ドライブ（ＨＤＤ）１３、通信インターフェイス１４からなり、それぞれ内部バスにより接続されている。
【００２１】
ＣＰＵ１１はＨＤＤ１３に格納されたプログラムを実行し、シミュレーション装置１の制御を行なう。ＲＡＭ１２はＣＰＵ１１のワークエリアである。ＨＤＤ１３は、プログラムやデータを保存するための記憶領域である。ＣＰＵ１１は、シミュレーション装置１の制御に必要なデータをＲＡＭ１２から読み出して処理を行なう。
【００２２】
通信インターフェイス１４は有線により通信ネットワーク１５に接続されており、他の端末１６との通信が可能である。また、通信インターフェイス１４は、アンテナを有する通信ユニットを介して、無線による通信を行なうことも可能である。例えば、端末１６にて、牧場経営者や牧場経営者への融資を検討している金融機関がインターネット上で所定のウェブサイトへアクセスし、ＩＤ及びパスワードを入力する。シミュレーション装置１では、通信インターフェイス１４を介してＩＤ及びパスワードを受信し、ログインが可能であるかどうかの認証を行なう。ＨＤＤ１３に設定された利用者データベースに記憶されたＩＤ及びパスワードが一致している場合は、ログインが可能となり、端末１６にてシミュレーションシステムを利用することが可能となる。
【００２３】
シミュレーションシステムへのログインが行なわれると、利用者は端末１６に設置されたキーボード等の入力装置へ、所定のデータ、例えば、初期の経産牛の数、初期の育成子牛の数、分娩間隔等に関する情報を入力することが可能である。端末１６の入力装置から入力された入力情報は通信インターフェイス１４を介してＲＡＭ１２に格納され、ＣＰＵ１１は入力情報をもとに各種の算定処理を実行する。算定処理の実行が終了すると、通信インターフェイス１４を介して、算定の結果が端末１６に送信される。端末１６にて算定の結果を受信すると、端末１６の表示装置にシミュレーションの結果が表示される。
【００２４】
次に、本発明の実施の形態にかかるシミュレーションシステムにおける、シミュレーションを行なう動物群の飼養動態について説明する。図２は、本発明の実施の形態にかかる動物群の飼養動態プロセスの一例を表す図である。この実施の形態において、動物群は牧場で飼養されている牛群であり、牛群には育成子牛２１、初妊牛２２、経産牛２３、初産牛２４および搾乳牛２５が含まれる。育成子牛２１は、出生してから最初に受胎するまでの雌子牛を指す。なお、この実施の形態では、雄子牛は出生して間もなく販売される。育成子牛２１は、出生してから１３〜１５か月経過すると、初めての種付けが行なわれる。種付けが成功し受胎した育成子牛２１は、初めての妊娠をしている牛である初妊牛２２となる。初妊牛２２は受胎から１０か月程度経過すると分娩を行ない、初妊牛２２から、初めての分娩を行なった牛である初産牛２４となる。なお、初産牛２４は、分娩の経験がある牛である経産牛２３に含まれる概念である。
【００２５】
初産牛２４を含む経産牛２３は、分娩後ほぼ２か月経過すると、次の種付けが行なわれる。種付けが成功し受胎した経産牛２３は、受胎から１０か月程度経過すると次の分娩を行なう。分娩間隔２６は、分娩と次の分娩の間の日数であり、経産牛ごとに異なるものであるが、本発明では牧場ごとの平均の日数を指すものとして取り扱う。なお、分娩間隔２６は、牧場の経営方針により変わってくるものであり、牧場ごとに異なる値である。分娩を行なった経産牛２３は、人間の飲用が可能な生乳を泌乳するため、搾乳が行なわれる。搾乳は、分娩の数日後から乾乳するまで行なわれる。乾乳とは、経産牛２３を次の分娩に備えさせるために、次の分娩予定日の約６０日前から、搾乳を止めることである。搾乳される経産牛２３は、搾乳牛２５にも分類される。
【００２６】
経産牛２３は、数回の分娩を行なうことができ、分娩後から乾乳まで搾乳牛２５として生乳生産に寄与するが、分娩回数を重ねるもしくは単に月齢を重ねることにより、泌乳量が減少する。泌乳量が一定量よりも少なくなるなどの条件に合致した牛は、廃用される。なお、廃用もしくは事故などによる死亡は、その他の育成子牛２１や初妊牛２２などにも起こりうることであり、このような牛は牛群から除かれる。
【００２７】
また、牛群の牛は販売することも可能である。これとは反対に、育成子牛などを外部から購入する場合もあり、この場合は購入した牛は対応する育成子牛２１、初妊牛２２や経産牛２３などに参入させることとなる。
【００２８】
図３は、本発明の実施の形態にかかる変化率テーブルの一例を表す図である。変化率テーブル４０は、シミュレーション装置１のＨＤＤ１３に格納されており、シミュレーションの開始にともないＲＡＭ１２に読み出される。変化率テーブル４０には、分娩間隔に関するパラメータ４１、搾乳牛率４２、育成子牛（期首）からの変化率４３、初妊牛（期首）からの変化率４４、初産牛（当期）の変化率４５、初妊牛（当期）の変化率４６、経産牛の変化率４７、出生率４８、雄子牛変化率４９、雌子牛変化率５０を含む項目が設定されている。これらの確率は、過去の統計をもとに予め設定したり、牧場の経営方針に応じて利用者が任意に設定することも可能である。なお、当期とは、これから算定処理を行なう期（所定の期間）を指す。
【００２９】
分娩間隔に関するパラメータ４１は、分娩間隔４１ａおよび分娩間隔係数４１ｂを含む情報である。分娩間隔４１ａは、利用者によって入力された分娩間隔が設定される。分娩間隔係数４１ｂは、一つの期（１年）の日数である３６５を分娩間隔４１ａで割った値である。経産牛１頭は分娩間隔ごとに１回の分娩を行なうため、分娩間隔係数４１ｂは一つの期における経産牛１頭が行なう分娩の回数とみなすことができる。したがって、期首における経産牛数に分娩間隔係数４１ｂを掛けることにより、当期において経産牛群が行なう分娩の回数を算定することができる。
【００３０】
搾乳牛率４２は、経産牛に占める搾乳牛の比率を表す数値であり、過去の統計をもとに予め設定されるものであるが、分娩間隔に応じて変動させることも可能である。搾乳牛数は経産牛数に搾乳牛率４２を掛けることで算定される。搾乳牛数は、後述するように生乳生産量を算定するために用いられる。
【００３１】
育成子牛（期首）からの変化率４３には、期首に育成子牛である牛が当期において、変化する確率が設定されている。育成子牛→初妊→初産牛率４３ａは、期首に育成子牛である牛が当期に初妊を経て初産を行なって初産牛となる確率である。育成子牛→初妊牛率４３ｂは、当期に初妊牛になる確率であり、育成子牛→育成子牛率４３ｃは、期末まで育成子牛のままである確率である。初妊牛（育成→）販売率４３ｄは、当期に初妊となってから販売し、期末には牛群に残らない牛となる確率である。廃用・死亡（育成（期首））率４３ｅは、当期に廃用もしくは死亡したため、期末に牛群に残らない牛となる確率である。なお、育成子牛→初妊→初産牛率４３ａ、育成子牛→初妊牛率４３ｂ、育成子牛→育成子牛率４３ｃ、初妊牛（育成→）販売率４３ｄ、廃用・死亡（育成（期首））率４３ｅの合計は１に等しくなる。
【００３２】
初妊牛（期首）からの変化率４４は、期首に初妊牛である牛が当期においてそれぞれの牛に変化する確率が設定された項目であり、当期に初産を行なう初産牛となる確率である初妊牛→初産牛率４４ａ、当期に廃用もしくは死亡したため期末に牛群に残らない牛となる確率である廃用・死亡（初妊牛（期首））率４４ｂが含まれる。初妊牛→初産牛率４４ａ、廃用・死亡（初妊牛（期首））率４４ｂの合計は１に等しくなる。
【００３３】
初産牛（当期）の変化率４５には、当期に初産牛となった牛が期末までに変化する確率が設定されている。当期に初産牛となった牛とは、期首に育成子牛であって当期に初妊を経て初産を行なって初産牛となった牛、期首に初妊牛であって当期に初産を行なって初産牛となった牛を含む概念である。したがって、育成子牛（期首）から初産牛となった牛数、および初妊牛（期首）から初産牛となった牛数の合計に初産牛（当期）の変化率４５を乗ずることで、当期に初産牛となった牛の期末までの変化が算定される。初産牛（当期）の変化率４５には、当期に廃用もしくは死亡する確率である廃用・死亡（初産牛）率４５ａ、および期末に経産牛に参入する確率である初産牛参入率４５ｂが含まれる。廃用・死亡（初産牛）率４５ａ、初産牛参入率４５ｂの合計は１に等しくなる。
【００３４】
初妊牛（当期）の変化率４６は、当期に育成子牛（期首）から初妊牛となった牛が期末までに変化する確率が設定されており、当期に廃用もしくは死亡する確率である廃用・死亡（初妊牛（育成→））率４６ａ、および期末まで初妊牛として残る確率である初妊牛（期末）率４６ｂが含まれる。廃用・死亡（初妊牛（育成→））率４６ａ、および初妊牛（期末）率４６ｂの合計は１に等しくなる。
【００３５】
経産牛の変化率４７は、期首に経産牛であった牛が期末までに変化する確率が設定されており、当期に販売される経産牛の確率である経産牛販売率４７ａ、当期に死亡する経産牛の確率である経産牛死亡率４７ｂ、および期末まで牛群に残る経産牛の確率である経産牛（期末）率４７ｃが含まれる。経産牛販売率４７ａ、経産牛死亡率４７ｂ、経産牛（期末）率４７ｃの合計は１に等しくなる。
【００３６】
次に、分娩及び出生に関する変化率について説明する。出生率４８は、当期に出生するそれぞれの子牛の確率を設定したものであり、当期に出生する子牛の数は、当期に分娩を行なう牛数に等しいとして算定される。したがって、経産牛のうち当期に分娩を行なう牛数、および当期に分娩を行なった初産牛数の合計に出生率４８を乗ずることで、当期に出生する子牛の数が算定される。なお、初産牛数は、育成子牛（期首）から初産牛となった牛数および初妊牛（期首）から初産牛となった牛数の合計である。出生率４８には、雄子牛が生まれる確率である雄子牛出生率４８ａ、雌子牛が生まれる確率である雌子牛出生率４８ｂ、および死産となる確率である子牛死産率４８ｃが含まれる。雄子牛出生率４８ａ、雌子牛出生率４８ｂ、子牛死産率４８ｃの合計は１に等しくなる。
【００３７】
雄子牛変化率４９は、当期に出生した雄子牛の数が、販売等により変化する確率を表すものである。本実施の形態において、出生した雄子牛は、その種類に応じた単価で販売される。雄子牛変化率４９には、乳用雄として出生し販売される確率である乳用雄販売率４９ａ、和牛として出生し販売される確率である和牛販売率４９ｂおよびＦ１牛として出生し販売される確率であるＦ１牛販売率４９ｃが含まれる。乳用雄販売率４９ａ、和牛販売率４９ｂ、Ｆ１牛販売率４９ｃの合計は１に等しくなる。
【００３８】
雌子牛変化率５０は、当期に出生した雌子牛の数が、販売等により変化する確率を表すものである。本実施の形態において、雌子牛は１３〜１５か月で受胎させるまで育成子牛として飼育される。雌子牛変化率５０は、育成子牛として飼育される確率である育成子牛飼育率５０ａ、および廃用もしくは死亡する確率である廃用・死亡（育成（出生））率５０ｂが含まれる。育成子牛飼育率５０ａ、廃用・死亡（育成（出生））率５０ｂの合計は１に等しくなる。
【００３９】
次に、牛群評価テーブルについて説明する。図４は、本発明の実施の形態にかかる牛群評価テーブルの一例を表す図である。牛群評価テーブル６０は、シミュレーション装置１のＲＡＭ１２に設けられ、利用者によって入力された初期牛数（初期育成子牛数、初期初妊牛数および初期経産牛数を含む）、および後述する牛群評価処理の結果が記憶される。牛群評価テーブル６０には、期首牛数６１、育成子牛（期首）からの変化牛数６２、初妊牛（期首）からの変化牛数６３、初産牛（当期）の変化牛数６４、初妊牛（当期）の変化牛数６５、経産牛の変化牛数６６、出生牛数６７、出生雄牛数６８、出生雌牛数６９、期末牛数７０、および搾乳牛数７１を含む項目が記憶される。
【００４０】
本発明の実施の形態において、育成子牛数、初妊牛数や経産牛数などの牛数は、母集団の牛数に対し、変化率テーブル４０に設定された変化率を掛けることにより算定される。そのため、算定された牛数は小数を含む値になる場合がある。この場合は以下のように端数処理が行われる。すでに述べたように、育成子牛（期首）からの変化率４３、初妊牛（期首）からの変化率４４等に含まれる各変化率は、いずれもその合計が１になる。この場合は、同じ母集団から算定された各牛数の小数部分を切り下げ、切り下げた各牛数の合計と母集団の数とを比較し、各牛数の合計が母集団の数と等しくなるように各牛数の小数部分が大きいものから順に切り上げる。なお、小数部分が等しい場合は、各牛数の値が大きいものから順に切り上げる。
【００４１】
例えば、母集団の数である育成子牛（期首）数６１ａが「１４０」であり、育成子牛→初妊→初産牛率４３ａ、育成子牛→初妊牛率４３ｂ、育成子牛→育成子牛率４３ｃ、初妊牛（育成→）販売率４３ｄ、廃用・死亡（育成（期首））率４３ｅを育成子牛（期首）数６１ａに乗ずることで得られる値が、それぞれ、「１７．８」、「８７．７」、「１９．２」、「１２．２」、「３．２」である場合、その合計が「１４０」となるように小数部分が大きい「１７．８」、「８７．７」が切り上げられて、「１８」、「８８」となり、その他は「１９」、「１２」、「３」と小数点以下が切り下げられる。また、設定された変化率が分娩間隔に関するパラメータ４１または搾乳牛率４２である場合は、小数点以下は切り下げる。
【００４２】
続いて、牛群評価テーブル６０に含まれる情報について詳しく説明する。まず、期首牛数６１は、当期における算定処理のもととなる、期首のそれぞれの牛の数が記憶された項目であり、期首における育成子牛の数である育成子牛（期首）数６１ａ、期首における初妊牛の数である初妊牛（期首）数６１ｂおよび期首における経産牛の数である経産牛（期首）数６１ｃが含まれる。これら期首牛数６１は、当期が第１期であれば、利用者によって入力された初期牛数、すなわち初期育成子牛数が育成子牛（期首）数６１ａとして記憶され、同様に、初期初妊牛数が初妊牛（期首）数６１ｂとして、初期経産牛数が経産牛（期首）数６１ｃとして、それぞれ記憶される。当期が第２期以降であれば、前の期における期末牛数７０のそれぞれの値が記憶される。
【００４３】
育成子牛（期首）からの変化牛数６２は、育成子牛（期首）数６１ａに対して、変化率テーブル４０における育成子牛（期首）からの変化率４３で設定されたそれぞれの確率を乗じて算定される項目である。まず、育成子牛→初妊→初産牛数６２ａは、育成子牛（期首）６１ａに育成子牛→初妊→初産牛率４３ａを乗じて算定される。例えば第１期において、育成子牛→初妊→初産牛数６２ａは、育成子牛（期首）６１ａである「１４０」に、育成子牛→初妊→初産牛率４３ａである「０．１３」をかけ、端数処理された「１８」として算定される。以下同様に、育成子牛→初妊牛数６２ｂは育成子牛→初妊牛率４３ｂを、育成子牛→育成子牛数６２ｃは育成子牛→育成子牛率４３ｃを、初妊牛（育成→）販売数６２ｄは初妊牛（育成→）販売率４３ｄを、廃用・死亡（育成（期首））数６２ｅは廃用・死亡（育成（期首））率４３ｅを、それぞれ育成子牛（期首）数６１ａに掛けることで算定される牛数である。
【００４４】
初妊牛（期首）からの変化牛数６３は、初妊牛（期首）数６１ｂに対して変化率テーブル４０における初妊牛（期首）からの変化率４４で設定されたそれぞれの確率を乗じて算定される項目である。それぞれ、初妊牛→初産牛数６３ａは初妊牛→初産牛率４４ａを、廃用・死亡（初妊牛（期首））数６３ｂは廃用・死亡（初妊牛（期首））率４４ｂを、初妊牛（期首）数６１ｂにそれぞれ乗じて算定される。
【００４５】
初産牛（当期）の変化牛数６４は、当期に初産牛となった牛が期末までに変化した牛数が算定される項目である。当期に初産牛となった牛数は、育成子牛→初妊→初産牛数６２ａおよび初妊牛→初産牛数６３ａの合計である。この合計に変化率テーブル４０の初産牛（当期）の変化率４５に設定されたそれぞれの確率を乗ずることで、初産牛（当期）の変化牛数６４が算定される。廃用・死亡（初産牛）数６４ａは廃用・死亡（初産牛）率４５ａを、初産牛参入数６４ｂは初産牛参入率４５ｂを、それぞれ乗ずることにより算定される。例えば第１期において、育成子牛→初妊→初産牛数６２ａである「１８」および初妊牛→初産牛数６３ａである「７６」の合計である「９４」に初産牛参入率４５ｂ「０．９５」を乗じて端数処理された「８９」が、初産牛参入数６４ｂとなる。
【００４６】
初妊牛（当期）の変化牛数６５は、当期に初妊牛となった牛が期末までに変化した牛数が算定される項目である。初妊牛（当期）の変化牛数６５は、育成子牛→初妊牛数６２ｂに初妊牛（当期）の変化率４６の各確率を乗じて算定される。廃用・死亡（初妊牛（育成→））数６５ａは廃用・死亡（初妊牛（育成→））率４６ａを、初妊牛（期末）数６５ｂは育成子牛→初妊牛数６２ｂに初妊牛（期末）率４６ｂを乗じることで算定される。
【００４７】
経産牛の変化牛数６６は、当期に経産牛である牛が期末までに変化する牛数として算定される。ここで、当期に経産牛である牛とは、期首に経産牛である牛に、当期に初産牛となった牛を合計したものである。まず、経産牛（期首）数６１ｃに初産牛参入数６４ｂを加算した経産牛全群数６６ａが算定される。この経産牛全群数６６ａに経産牛の変化率４７の各確率を掛けることで、経産牛の変化牛数６６が算定される。それぞれ、経産牛販売数６６ｂは経産牛販売率４７ａを、経産牛死亡数６６ｃは経産牛死亡率４７ｂを、および経産牛（期末）数６６ｄは経産牛（期末）率４７ｃを、経産牛全群数６６ａに乗じることで算定される。
【００４８】
次に、出生牛数６７について説明する。本実施の形態において、出生牛数６７の各牛数の合計は、当期において分娩を行なう全牛数に等しいとして算定される。すなわち、経産牛のうち当期に分娩を行なう牛数、育成子牛→初妊→初産牛数６２ａおよび初妊牛→初産牛数６３ａの合計である。ここで、経産牛であって当期に分娩を行なう牛数は、経産牛（期首）数６１ｃに分娩間隔係数４１ｂを乗じた値を小数点第１位で四捨五入することにより算定される。続いて、分娩頭数（調整）６７ａが、経産牛のうち当期に分娩を行なう牛数、育成子牛→初妊→初産牛数６２ａおよび初妊牛→初産牛数６３ａの合計として算定される。分娩頭数（調整）６７ａを出生牛数６７の各牛数の合計とみなし、これに出生率４８に設定された各確率を乗じて、それぞれの牛数が算定される。すなわち、雄子牛出生数６７ｂは雄子牛出生率４８ａ、雌子牛出生数６７ｃは雌子牛出生率４８ｂ、子牛死産数６７ｄは子牛死産率４８ｃを、それぞれ分娩頭数（調整）６７ａに乗じることで算定される。
【００４９】
出生雄牛数６８は、雄子牛出生数６７ｂに、雄子牛変化率４９に設定された各確率を乗じて算定される。すなわち、乳用雄販売数６８ａは乳用雄販売率４９ａ、和牛販売数６８ｂは和牛販売率４９ｂ、Ｆ１牛販売数６８ｃはＦ１牛販売率４９ｃをそれぞれ乗じて算定される。また、出生雌牛数６９は、変化率テーブル４０の雌子牛変化率５０の各確率、すなわち育成子牛飼育数６９ａは育成子牛飼育率５０ａを、廃用・死亡（育成（出生））数６９ｂは廃用・死亡（育成（出生））率５０ｂを、雌子牛出生数６７ｃに乗じて算定される。
【００５０】
次に、期末牛数７０について説明する。期末牛数７０には、育成子牛（期末）数７０ａ、初妊牛（期末）数７０ｂおよび経産牛（期末）数７０ｃを含む情報が記憶される。育成子牛（期末）数７０ａは、当期に育成子牛のままであった育成子牛→育成子牛数６２ｃおよび当期に出生し飼育された育成子牛飼育数６９ａの合計である。初妊牛（期末）数７０ｂは初妊牛（期末）数６５ｂに等しく、経産牛（期末）数７０ｃは経産牛（期末）数６６ｄに等しい。
【００５１】
また、牛群評価テーブル６０には、搾乳牛数７１も含まれる。搾乳牛数７１は、経産牛数に搾乳牛率４２を乗じて算定される。ただし、牛群評価テーブル７０における搾乳牛数７１の項目は、後述する経営評価処理において、年間生乳生産量を算定するために用いられるため、まず経産牛（期首）数６１ｃに搾乳牛率４２を乗じた搾乳牛（期首）数７１ａ、および経産牛（期末）数７０ｃに搾乳牛率を乗じた搾乳牛（期末）数７１ｂが、それぞれ算定される。続いて、搾乳牛（期首）数７１ａと搾乳牛（期末）数７１ｂの平均である生乳生産牛数７１ｃが算定され、後述する経営評価処理における算定に用いられる。なお、生乳生産牛数７１ｃは年間生乳生産量を算定するための平均値に過ぎないため、端数処理は行なわれない。
【００５２】
次に、本発明の実施の形態における、経営評価処理において用いられる数値について、図５及び図６を用いて詳細に説明する。図５は、本発明の実施の形態にかかる単価テーブルを表す図である。図６は、本発明の実施の形態にかかる経営評価テーブルを表す図である。
【００５３】
単価テーブル８０は、シミュレーション装置１のＨＤＤ１３に格納されており、シミュレーションの開始にともないＲＡＭ１２に読み出される。単価テーブル８０には、生乳生産単価８１および個体単価８２を含む項目が格納されている。生乳生産単価８１は、搾乳牛１頭当たりの年間生乳生産量（単位：ｋｇ／頭）である個体乳量８１ａ、および生乳１ｋｇの単価（単位：円／ｋｇ）である乳価８１ｂを含む項目である。個体単価８２は、各種類の牛１頭を販売した時の単価（単位：円／頭）であって、乳用雄単価８２ａ、和牛雄単価８２ｂ、Ｆ１雄単価８２ｃ、初妊牛単価８２ｄおよび経産牛単価８２ｅを含む項目である。
【００５４】
経営評価テーブル９０は、経営評価処理において算定された数値が記憶されるテーブルであり、ＲＡＭ１２に設けられる。経営評価テーブル９０は、生乳生産金額９１、個体販売金額９２および売上高９３を含む項目が記憶される。生乳生産金額９１および個体販売金額９２は、単価テーブル８０で設定された単価および牛群評価テーブル６０に記憶された牛数をもとに算定される。
【００５５】
生乳生産金額９１は、生乳生産量９１ａおよび生乳販売金額９１ｂを含む項目である。生乳生産量９１ａは、当期に生産される生乳の合計量であり、生乳生産牛数７１ｃに個体乳量８１ａを乗じた値として算定される。生乳販売金額９１ｂは、生乳生産量９１ａに乳価８１ｂを乗じて算定されたものであり、当期における生乳による販売金額の総額を表す。なお、図６において、生乳販売金額９１ｂは単位を千円として表示されている。
【００５６】
個体販売金額９２は、当期に販売した各種類の牛の販売金額の合計が記憶された項目である。乳用雄販売９２ａは乳用雄販売数６８ａに乳用雄単価８２ａを、和牛雄販売９２ｂは和牛販売数６８ｂに和牛雄単価８２ｂを、Ｆ１雄販売９２ｃはＦ１牛販売数６８ｃにＦ１雄単価８２ｃを、初妊牛販売９２ｄは初妊牛（育成→）販売数６２ｄに初妊牛単価８２ｄを、経産牛販売９２ｅは経産牛販売数６６ｂに経産牛単価８２ｅを、それぞれ乗じることで算定される。個体販売合計９２ｆは、乳用雄販売９２ａから経産牛販売９２ｅまでの値の合計である。なお、図６において、個体販売金額９２に含まれる金額は全て単位を千円として表示されている。
【００５７】
売上高９３は、当期における総売上高に相当し、生乳販売金額９１ｂおよび個体販売合計９２ｆの合計として算定される。図６において、売上高９３は単位を千円として表示されている。
【００５８】
次に、本実施の形態にかかる、牛群の評価処理について説明する。図７は、本発明の実施の形態にかかる牛群の評価処理についてのフローチャートの一例を表す図である。まず、端末１６がシミュレーション装置１へのログインが行なわれると、利用者は端末１６の入力装置を通じ、シミュレーションを行なうために必要な初期値の入力を行なう。シミュレーション装置１は通信インターフェイス１４を介して、この入力を受け付ける（ステップＳ１）。受け付けられた初期値は、ＲＡＭ１２に格納される。
【００５９】
この実施の形態において、初期値は、シミュレーションを行なう期間であるシミュレーション期間（例えば５年）、シミュレーション期間の開始時におけるそれぞれの牛群に含まれる牛の頭数である初期牛数（初期育成子牛数、初期初妊牛数、初期経産牛数）、および分娩間隔などが含まれる。シミュレーション期間は、１年を単位として入力され、１以上の整数が入力可能である（ここでは、ｎ＞１であるｎ年が入力されたものとする）。シミュレーション期間は１年ごとに第１期から第ｎ期に分割され、第１期から順に牛群評価処理が行われる。初期牛数は、シミュレーション期間の開始時における育成子牛２１、初妊牛２２および経産牛２３の数を含み、それぞれ０以上の整数値で入力される。
【００６０】
ステップＳ１において初期値の入力が受け付けられると、ＣＰＵ１１は、ＲＡＭ１２に設けられた牛群評価テーブル６０において、期首牛数６１を記憶する（ステップＳ２）。ここでは、第１期の算定を行なうものとし、初期牛数が期首牛数６１として記憶される。
【００６１】
続いて、ＣＰＵ１１は、当期における育成子牛２１から変化した牛の数を算定する（ステップＳ３）。ステップＳ３は、牛群評価テーブル６０における育成子牛（期首）からの変化牛数６２に含まれる項目を算定する処理である。前述したように、育成子牛（期首）からの変化牛数６２は、育成子牛（期首）数６１ａに育成子牛（期首）からの変化率４３を乗じて算定される。
【００６２】
次に、当期における初妊牛２２から変化した牛の数が算定される（ステップＳ４）。ステップＳ４は、牛群評価テーブル６０における初妊牛（期首）からの変化牛数６３に含まれる項目を算定する処理である。初妊牛（期首）からの変化牛数６３は、初妊牛（期首）数６１ｂに初妊牛（期首）からの変化率４４を乗じて算定される。続いて、当期に初産牛となった牛から変化した牛の数、および当期に初妊牛となった牛から変化した牛の数が算定される（ステップＳ５）。当期に初妊牛となった牛から変化した牛の数、および当期に初産牛となった牛から変化した牛の数とは、それぞれ牛群評価テーブル６０における初産牛（当期）の変化牛数６４、および初妊牛（当期）の変化牛数６５に相当し、これらがステップＳ５において算定される。
【００６３】
次に、当期における経産牛２３の数が算定される（ステップＳ６）。ステップＳ６は、牛群評価テーブル６０における経産牛の変化牛数６６に含まれる項目を算定する処理である。ステップＳ６までが終了すると、それまでのステップで算定された育成子牛（期首）からの変化牛数６２、初妊牛（期首）からの変化牛数６３および経産牛の変化牛数６６の算定結果をもとに、当期に出生する牛の数が算定される（ステップＳ７）。当期に出生する牛の数とは、牛群評価テーブル６０における出生牛数６７、出生雄牛数６８および出生雌牛数６９を含む概念である。
【００６４】
続いて、ＣＰＵ１１は、牛群評価テーブル６０に記憶された算定結果をもとに、期末における牛群の牛数を算定する（ステップＳ８）。ここで算定されるのは期末牛数７０および搾乳牛数７１である。
【００６５】
ステップＳ８が終了すると、ＣＰＵ１１は、経営評価処理を行なう（ステップＳ９）。経営評価処理は、牛群評価テーブル６０に記憶された当期の各牛数、およびＲＡＭ１２に読み出された単価テーブル８０に設定された販売単価にもとづいて行なわれる算定処理である。ステップＳ８において算定された数値は、当期の販売金額として、ＲＡＭ１２に設けられた経営評価テーブル９０に記憶される。
【００６６】
ステップＳ９までが終了すると、終了条件を満たすか否かが判定される（ステップＳ１０）。つまり、第ｎ期の算定が終了したか否かが判定される。上記の説明では、ｎ＞１であるｎ年のシミュレーション期間に対し、第１期の算定を行なっているため、ステップＳ１０においてＮｏとなり、ステップＳ２に戻る。この場合、ステップＳ２以降では第２期の算定を始めることになるため、ステップＳ２において、第２期の期首牛数６１として第１期の期末牛数７０がそれぞれ設定され、ステップＳ３〜Ｓ１０の処理が行なわれる。同様の処理により、第ｎ期までの算定が繰り返され、第ｎ期までの処理が終了すると、終了条件を満たすと判定され（ステップＳ１０においてＹｅｓ）、算定結果が端末１６に送信される（ステップＳ１１）。端末１６が算定結果を受信し、端末１６の表示画面に算定結果が表示されると、利用者は表示された算定結果をもとに様々な牛群評価が可能となる。
【００６７】
なお、算定結果をもとに、利用者が算定条件の変更を希望する場合は、利用者はその旨を端末１６に入力するように構成しても良い。変更可能な算定条件は、初期値に含まれるシミュレーション期間、初期牛数および分娩間隔が挙げられるが、これらのほか、経産牛の販売率などを変更可能としてもよい。入力された算定条件の変更に関する情報は端末装置１６からシミュレーション装置１に送信され、シミュレーション装置１のＣＰＵ１１がこれを受け付けると、初期値そのほかの算定条件が再設定され、牛群評価処理が再度行われる。
【００６８】
また、上記の説明では、ｎ年間における牛群変動をシミュレーションするため、ステップＳ１０において、第ｎ期の算定が終了したか否かを判定した。しかし、本実施の形態を応用することにより、所定の目標値を設定し、目標値を上回るために必要な年数をシミュレーションすることも可能である。この場合、ステップＳ１において、初期値として初期牛数および分娩間隔の入力は受け付けるが、シミュレーション期間は受け付けず、代わりに経産牛数、生乳生産量または売上高などの目標値の入力を受け付ける。受け付けた初期値をもとに牛群評価処理を行ない、算定された値が目標値を上回ったか否かをステップ１０で判定し、目標値を上回るまで算定処理を繰り返すことにより、必要な年数を算定することができる。なおこの場合、算定処理が膨大な回数になるのを回避するため、繰り返し処理の上限を設定しておくことが好ましい。
【００６９】
次に、最適分娩間隔の評価処理について説明する。この処理は、初期値（初期育成子牛数、初期初妊牛数および初期経産牛数を含む初期牛数ならびにシミュレーション期間を含む）に加えて、シミュレーション期間（例えば、ｎ年）後の目標値である目標牛数として育成子牛数、初妊牛数または経産牛数の入力を受け付け、ｎ年後において目標牛数を上回るための最長の分娩間隔を、牧場等の経営に最適な分娩間隔である最適分娩間隔として特定するものである。なお、最適分娩間隔は、目標牛数をｎ年後に上回ることが可能な最長の分娩間隔として算定されるため、ｎ−１年後には目標牛数に達せずｎ年後に初めて目標牛数を上回ることが可能な最長の分娩間隔として算定されることとなる。
【００７０】
牛は分娩間隔に従って分娩を行なうため、同じｎ年間で比較した場合、分娩間隔が短いほど出生牛数６７が増加する。出生牛数６７が増加すると、これに応じて次期以降における育成子牛数、初妊牛数および経産牛数が増加する。したがって、牛群を拡大するためには分娩間隔が短い方がよい。しかし、分娩間隔を短く管理するためには多大な労力やコストが必要である。本発明によって、設定された目標牛数をｎ年後に上回ることが可能な最長の分娩間隔を最適分娩間隔として特定し、利用者に提示することで、利用者は分娩間隔を管理した場合の利益とコストを判断することが可能となり、牧場等の経営判断を行なう際の指標とすることができる。
【００７１】
本発明における最適分娩間隔の特定は、３５０日から４８０日までの範囲でのみ処理が行なわれる。３５０日は、牛が分娩を行なってから次に受胎可能になるまでの期間および妊娠期間の合計であり、生理学的に取りうる最短の分娩間隔である。また、４８０日は、搾乳可能な最長期間の日数および乾乳期間として取りうる最長期間の日数の合計であり、すなわち分娩間隔として取りうる最長の日数である。このため、最適分娩間隔の特定処理を３５０日から４８０日の範囲の外まで行なっても、牧場等の経営指標として有用なものになり得ず、シミュレーション装置１の負荷増大にもつながるため、３５０日から４８０日までの範囲に限定して最適分娩間隔の特定処理が行なわれる。
【００７２】
最適分娩間隔の特定処理は、まず最適分娩間隔が３５０日から４８０日の範囲に含まれるか否かが判定される。続いて、最適分娩間隔が含まれる範囲の上限値および下限値を設定し、その中間値を仮分娩間隔として算定する。算定された仮分娩間隔および初期値（初期牛数およびシミュレーション期間を含む）を用いて牛群評価処理によりシミュレーション期間後の各牛数を算定し、算定された牛数が目標牛数を上回るか否かの判定が行なわれる。
【００７３】
続いて、最適分娩間隔の評価処理の具体的な手順について、図８および図９を用いて説明する。図８は、本発明の最適分娩間隔の評価処理についてのフローチャートの一例である。図９は、本発明の実施の形態にかかる目標達成可否判定処理を表わすフローチャートの一例である。なお、ここでは目標値として経産牛数が所定の値に設定されたとする。
【００７４】
まず、利用者は端末１６の入力装置を通じ、シミュレーションを行なうために必要な初期値および目標値の入力を行ない、シミュレーション装置１は通信インターフェイス１４を介してこの入力情報を受け付ける（ステップＳ２１）。受け付けられた初期値および目標値は、ＲＡＭ１２に格納される。
【００７５】
次に、ＣＰＵ１１は、仮分娩間隔を３５０日と設定し、仮分娩間隔、およびステップＳ２１で受け付けた初期値を用いて、経産牛数の算定を行なう（ステップＳ２２）。経産牛の算定は、図７のステップＳ２〜Ｓ１０の処理にもとづいて行なわれる（ステップＳ８の経営評価処理は行なわなくても良い）。続いて、仮分娩間隔で算定された経産牛数が目標値の経産牛数よりも少ないか否かの判定を行なう（ステップＳ２３）。ここで、仮分娩間隔を理論上の最短日数である３５０日としているので、算定された経産牛数が目標値の経産牛数よりも多いと判定された場合（ステップＳ２３においてＮｏ）、目標値の経産牛数は入力された初期値で達成可能であると判定されたことを意味する。
【００７６】
仮分娩間隔で算定された経産牛数が目標値の経産牛数よりも多いと判定された場合は、続いて最適分娩間隔の算定が行なわれる（ステップＳ２４）。最適分娩間隔の算定処理については、図９を用いて後に詳述する。最適分娩間隔の算定が行なわれると、ＣＰＵ１１は、通信インターフェイス１４を介し、端末１６に算定結果を送信する（ステップＳ２５）。
【００７７】
なお、ステップＳ２３において、仮分娩間隔で算定された経産牛数が目標値の経産牛数よりも小さいと判定された場合（ステップＳ２３においてＹｅｓ）、分娩間隔を理論上の最短日数である３５０日に設定してもなお目標値に達成しないこととなるため、分娩間隔の算定を行なわず、代わりに目標値の経産牛数に達成することが可能な最短年数を算定する処理が行なわれる（ステップＳ２６）。最短年数の算定は、図７のステップＳ２〜Ｓ１０が利用され、ステップＳ１０において算定された経産牛数が目標値の経産牛数を上回ると判定されるまで算定処理が繰り返される。算定された経産牛数が目標値の経産牛数を初めて上回った期が、ステップＳ２６における目標値の経産牛数に達成することが可能な最短年数として算定される。
【００７８】
続いて、図９を用いて、ステップＳ２４の最適分娩間隔算定処理について詳細に説明する。まず、ＣＰＵ１１は、仮分娩間隔を４８０日と設定し、ＲＡＭ１２に格納された初期値とともに、図７のステップＳ２〜Ｓ１０の処理により経産牛数を算定する（ステップＳ４１）。続いて、算定された経産牛数が、目標値の経産牛数以上であるか否かの判定を行なう（ステップＳ４２）。
【００７９】
分娩間隔を４８０として算定された経産牛数が目標値の経産牛数以上であった場合（ステップＳ４２におけるＹｅｓ）、最長の分娩間隔である４８０日であっても、目標値を達成できてしまうことになる。この場合は、最適分娩間隔が４８０日であると特定し（ステップＳ４３）、最適分娩間隔の算定処理を終了する。このように構成することで、有用でない指標の算定処理を省略することができるため、ＣＰＵ１１の負荷を軽減することが可能となる。なお、この場合、ステップＳ２５において、最適分娩間隔が４８０日であることが端末１６に送信される。
【００８０】
ステップＳ４２において算定された経産牛数が目標値の経産牛数よりも小さいと判定されると（ステップＳ４２においてＮｏ）、最適分娩間隔が含まれる可能性のある範囲の上限値および下限値が設定される（ステップＳ４４）。ここでは、最適分娩間隔は３５０日から４８０日の範囲に含まれるため、上限値は４８０日、下限値は３５０日に設定される。続いて上限値と下限値の中間値が算定される（ステップＳ４５）。ここで、中間値とは、上限値と下限値の平均であるが、値が小数を含む場合は小数部分を切り捨てた値である。
【００８１】
続いて、仮分娩間隔を中間値に設定し、図７のステップＳ２〜Ｓ１０の処理により経産牛数の算定を行なう（ステップＳ４６）。続いて、ステップＳ４６で算定された経産牛数が目標値の経産牛数と等しいか否かの判定を行なう（ステップＳ４７）。ここで算定された経産牛数が目標値の経産牛数と等しい場合は、最適分娩間隔は中間値であるとして算定される（ステップＳ４８）。
【００８２】
ステップＳ４６で算定された経産牛数と目標値の経産牛数が等しくないと判定された場合（ステップＳ４７においてＮｏ）、ステップＳ４６で算定された経産牛数が目標値の経産牛数よりも大きいか否かの判定を行なう（ステップＳ４９）。経産牛数が目標値の経産牛数よりも大きいと判定されれば（ステップＳ４９においてＹｅｓ）、最適分娩間隔は中間値と上限値の間に含まれることになるので、この中間値が新たな下限値として再設定される（ステップＳ５０）。すなわち、最適分娩間隔が特定されるまで、ステップＳ４５〜Ｓ５０までの処理は繰り返されることになるが、今回の処理がｍ回である場合に、ｍ＋１回目の処理における上限値はｍ回目と同じ上限値が用いられ、ｍ回目の中間値がｍ＋１回目の下限値として用いられる。算定された牛数が目標牛数を上回った場合、分娩間隔が、この仮分娩間隔より長い場合でも目標牛数を上回ることが可能であるので、最適分娩間隔は仮分娩間隔と上限値の間の値をとることになる。したがって、ステップＳ５０において、この中間値を新たな下限値として再設定し、ステップＳ４５〜Ｓ５２までの処理が再度行なわれることになる。
【００８３】
ステップＳ４９において、ステップＳ４６で算定された経産牛数が目標値の経産牛数よりも小さいと判定されれば、中間値が新たな上限値として再設定される（ステップＳ５１）。同様に、今回の処理がｍ回である場合に、ｍ＋１回目の処理における下限値はｍ回目と同じ下限値が用いられ、ｍ回目の中間値がｍ＋１回目の上限値として用いられる。同様に、算定された牛数が目標牛数を下回った場合、分娩間隔が仮分娩間隔よりも短くなければ目標牛数を達成できないため、最適分娩間隔は仮分娩間隔と下限値の間の値をとることになる。したがって、ステップＳ５０において、この中間値を新たな上限値として再設定し、ステップＳ４５〜Ｓ５２までの処理が再度行なわれることになる。
【００８４】
続いて、新たに設定された最適分娩間隔が含まれる範囲が、さらに限定される必要があるか否かを判定するため、新たに設定された上限値および下限値を用い、上限値が下限値＋１と等しいか否かの判定を行なう（ステップＳ５２）。等しくなければ（ステップＳ５２においてＮｏ）、最適分娩間隔が含まれる範囲をさらに限定する必要があるので、ステップＳ４５に戻る。新たに設定された上限値が下限値＋１と等しければ（ステップＳ５２においてＹｅｓ）、最適分娩間隔は整数値として算定されるため、これ以上の限定を行なう必要がなく処理が終了することになる。この場合、下限値を最適分娩間隔とすれば、経産牛数が目標値の経産牛数を上回るので、最適分娩間隔として下限値の値が算定され（ステップＳ５３）、最適分娩間隔算定処理は終了する。最適分娩間隔算定処理が終了すると、図８におけるステップＳ２５に戻る。
【００８５】
最適分娩間隔の特定は、シミュレーション期間であるｎ年の経過後に目標牛数の達成が可能である分娩間隔のうち、最長の日数を特定するものであり、ステップＳ４１〜Ｓ５３の処理を実行することで、ｎ−１年の経過後には目標牛数の達成ができないが、ｎ年の経過後に目標牛数の達成が可能な分娩間隔を特定することが可能となる。仮分娩間隔として、最適分娩間隔が含まれる可能性のある範囲の上限値と下限値の中間値が用いられ、ステップＳ４５〜Ｓ５２の処理を繰り返して、中間値が算定結果に応じて変化する（目標牛数に達していない場合は、より短い仮分娩間隔を新たに設定し、目標牛数に達している場合は、より長い仮分娩間隔を新たに設定する）ことで、最適分娩間隔が含まれる可能性のある範囲が限定され、最終的に最適分娩間隔が特定されることになる。このような処理を経ることでステップＳ２〜Ｓ１０の処理の回数は数回程度で済むことになり、シミュレーション装置１の処理負荷の軽減にもつながる。
【００８６】
なお、上記の実施の形態では目標値として育成子牛数、初妊牛数または経産牛数としたが、目標値として生乳生産量が設定されても良い。この場合、生乳生産量は経産牛数から算定されるため、目標生乳生産量から目標経産牛数を逆算し、目標経産牛数に達するように最適分娩間隔を算定すればよい。また、上記の実施の形態では、初期牛数と目標牛数を入力し、各変化率の値を固定の値として、最適な分娩間隔を求めることとしたが、初期牛数、目標牛数、特定の分娩間隔を利用者が入力し、最適な変化率（例えば、各牛の販売率等）を算定するような構成とすることも可能である。この場合、図８、図９について説明したのと同様の方法で、最適な変化率の算定が行なわれる。
【００８７】
図１０は、本発明の実施の形態にかかる端末１６の表示画面の変化を表す図である。この実施の形態では、利用者が端末１６からシミュレーション装置１に対してログインすると、シミュレーション装置１から端末１６に情報が送信され、図１０（ａ）のように「売上高を試算」と「分娩間隔を試算」の二つのシミュレーション方法を選択する画面が、端末１６の表示画面１６０に表示される。
【００８８】
利用者が「売上高を試算」を選択入力すると、その入力は端末１６からシミュレーション装置１に送信され、受信したシミュレーション装置１は、初期値を入力させるための画面情報を端末１６に送信し、端末１６はこの情報を受信して図１０（ｂ）のように表示画面１６０に表示する。ここで、利用者が初期値を入力し、入力ボタンを選択すると、入力された初期値は端末１６からシミュレーション装置１に送信される。シミュレーション装置１は、受信した初期値をもとに図７のステップＳ２〜Ｓ１０までの牛群評価処理を行ない、ステップＳ１０において算定結果を端末１６に送信する。
【００８９】
端末１６は、図１０（ｃ）のようにして、受信した算定結果を表示画面１６０に表示する。算定結果には、シミュレーション期間における売上高の推移を表したグラフ１６１、およびシミュレーション期間における売上高や各牛群の推移を表したテーブル１６２が表示される。利用者は表示画面１６０に表示されたグラフ１６１やテーブル１６２をもとに、牧場の経営評価を行なう。
【００９０】
なお、算定条件を変更し、再度シミュレーションを行なうことを利用者が希望する場合は、利用者は「条件を変更」ボタン１６３を選択入力することにより、算定条件の変更が可能である。ここで、変更可能な算定条件は初期値だけに限定されず、図１０（ｄ）のように、初妊牛や経産牛等の販売率を変更可能としてもよい。算定条件の変更に関する情報が端末１６に入力されると、入力された情報はシミュレーション装置１に送信され、受信したシミュレーション装置１は変更された算定条件をもとに牛群評価処理を行なう。変更された算定条件をもとに得られた算定結果は、再び端末１６に送信され、図１０（ｃ）のように端末１６の表示画面１６０に表示される。
【００９１】
利用者が「分娩間隔を試算」を選択入力すると、その入力は端末１６からシミュレーション装置１に送信され、受信したシミュレーション装置１は、初期値および目標値を入力させるための画面情報を端末１６に送信し、端末１６はこの情報を受信して図１０（ｅ）のように表示画面１６０に表示する。ここで、利用者が初期値および目標値を入力し、入力ボタンを選択すると、入力された初期値および目標値は端末１６からシミュレーション装置１に送信される。シミュレーション装置１は、受信した初期値および目標値をもとに、図８および図９の最適分娩間隔評価処理を行ない、ステップＳ２５において算定結果を端末１６に送信する。
【００９２】
端末１６は、図１０（ｆ）のようにして、受信した算定結果を表示画面１６０に表示する。算定結果には、最適分娩間隔の算定結果１６４のほか、シミュレーション期間における売上高の推移を表したグラフ１６１、シミュレーション期間における売上高や牛群の推移を表したテーブル１６２が表示される。利用者は表示画面１６０に表示された算定結果をもとに、牧場の経営評価を行なう。
【００９３】
なお、算定条件を変更し、再度シミュレーションを行なうことを利用者が希望する場合は、利用者は「条件を変更」ボタン１６５を選択入力することにより、算定条件の変更が可能である。ここで、変更可能な算定条件は初期値だけに限定されず、図１０（ｇ）のように、初妊牛や経産牛等の販売率を変更可能としてもよい。算定条件の変更に関する情報が端末１６に入力されると、シミュレーション装置１は変更された算定条件をもとに最適分娩間隔評価処理を行なう。変更された算定条件をもとに得られた算定結果は、図１０（ｆ）のように端末１６の表示画面１６０に表示される。
【００９４】
なお、上記の実施の形態に加えて、それぞれの分類の牛の購入または販売を適宜設定できる構成としても良い。どの分類の牛の購入または販売を重視するかは、牧場の経営方針に大きく左右される。牧場の経営方針は、頭数の規模拡大を重視する方針、頭数の規模を一定にし、若い牛を購入して老齢の牛を販売することにより牛の回転を早める方針、肉用牛として販売することを重視する方針など、様々な方針を取ることができる。したがって、代表的な経営方針に沿うような変化率テーブル４０をそれぞれ予め設定しておき、利用者に選択させるようにしても良い。さらに、利用者による変化率テーブル４０の設定値の一部変更を可能にしてもよい。
【００９５】
上記の実施の形態では、経営評価処理について、売上高を算定する構成としたが、売上高でなく利益を算定するように構成しても良い。この場合、各分類の牛の飼料費・購入費、牛舎などの施設の建設費・維持費、および労務費などの原価を設定した原価テーブルを設定しておき、これをもとに利益を算定する構成とすることができる。
【００９６】
上記の実施の形態では、牛群について説明したが、他の動物に応用することも可能である。他の動物に応用する場合は、変化率テーブルおよび単価テーブルの設定値を適切な値に変更すればよい。例えば妊娠期間が５か月程度の動物（山羊など）では、初妊牛に関する変化率、および分娩間隔を変更し、さらに一つの期を半年にするなどの設定変更により、本発明を応用することができる。さらに、一回の分娩で多数の子を産むことが通常である動物の場合（例えば、豚など）は、一回に分娩する子の数の平均値に分娩動物数を掛けた数に、変化率テーブルの出生率を掛けることで、出生数を算定する構成とすればよい。
【符号の説明】
【００９７】
１ シミュレーション装置
１１ ＣＰＵ
１２ ＲＡＭ
１３ ＨＤＤ
１４ 通信インターフェイス
１５ 通信ネットワーク
１６ 端末
４０ 変化率テーブル
６０ 牛群評価テーブル
８０ 単価テーブル
９０ 経営評価テーブル
【技術分野】
【０００１】
本発明は、複数の動物から構成される牧場等の経営に関するシミュレーションが可能なシミュレーション装置に関する。
【背景技術】
【０００２】
従来から、酪農経営のシミュレーション分析が試みられている。通常は、牛の頭数や個体の乳量（乳牛生産性、乳牛１頭あたりの乳量能力）に基づき試算されることが一般的であり、このような分析を行なう際には乳牛の泌乳力が重視されている。個体の乳量が高い場合は、乳牛の泌乳力が直接的に売上高のアップに寄与することから、泌乳力を用いた分析が広く用いられている。特に、生産現場にて飼料を与える際においても泌乳力の向上を目標にした設計が重んじられ、金融機関への融資の申請においても泌乳力に基づいた経営計画やシミュレーション分析が適用されている。
【０００３】
このため、酪農経営者による飼養管理の優先順位は自然と泌乳力向上が一番におかれていた。しかし、乳牛が持つ本来の生理的機能には、泌乳のほかに繁殖、自己成長、免疫などがあり、これらの事項は十分に考慮されていなかった。特に繁殖については、分娩され出生した雌子牛が、出生後に妊娠して子牛を分娩するまで、２年ほどの期間がかかるため、これらの期間の経過を考慮する必要がある。
【０００４】
このような中、牛群を管理する手法の一つとして、乳牛の栄養診断処理、牛群の健康管理、及び、繁殖管理処理をコンピュータにより経時的に行なうために、乳牛の乳成分データを基に演算処理し、演算処理されたデータにもとづいて摂取栄養を解析し、解析されたデータにもとづいてプログラム処理して乳牛の栄養診断処理を行なう乳牛の管理システムが開示されている（例えば、特許文献１参照）。また、その他、ミルクパイプラインのミルクタップ側に被識別部を設け、かつ搾乳ユニットのディストリビュータ側に、被識別部を検出して牛床を識別する識別部と、識別データ及び搾乳データを記憶する記憶部を有する制御ユニットを設けるとともに、データを受信して記憶する送受信ユニットを備える搾乳管理システムが開示されている（例えば、特許文献２参照）。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平７−１９４２７３号公報
【特許文献２】特開２００１−４５８９８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかし、特許文献１及び２はいずれも乳牛の管理をするものではあるが、牧場経営全体を管理するものではなく、将来にわたる牧場の成長を予測し、牧場を経営するうえでの目標を達成するための指針を提示できるといったものではない。本発明は、乳牛等の畜産動物の繁殖機能に基づき動物群評価及び牧場の経営評価を可能とするもので、畜産動物の飼養期間（時間の流れ）を考慮して、畜産動物の出生、育成、妊娠、初産のプロセスを経営シミュレーションに組み込み、動物群の繁殖力と供用力を見出して、将来にわたる牧場の成長の潜在性を計数化、指標化するためのシミュレーション装置を提供することを目的とする。
【課題を解決するための手段】
【０００７】
上記の課題を解決するため、本発明は、複数の動物が存在する所定の動物群内における、分娩の経験がある経産動物の数、分娩の経験がなく妊娠もしていない非経産動物の数、および初めての妊娠をしている初妊動物の数についての所定の期間における変動を算定するシミュレーション装置であって、前記所定の期間における期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに該動物群における分娩間隔を記憶する記憶手段と、記憶手段により記憶された期首の非経産動物の数をもとに、前記所定の期間内に妊娠を経て初めての分娩を行なう第１の初産動物の数、前記所定の期間内に初妊動物となる期末の初妊動物の数、および前記所定の期間内に変化しない第１の非経産動物の数を算定する非経産動物変化数算定手段と、記憶手段により記憶された期首の初妊動物の数をもとに、前記所定の期間内に初めての分娩を行なう第２の初産動物の数を算定する初妊動物変化数算定手段と、記憶手段により記憶された期首の経産動物の数、非経産動物変化数算定手段により算定された第１の初産動物の数および初妊動物変化数算定手段により算定された第２の初産動物の数をもとに、期末の経産動物の数を算定する経産動物変化数算定手段と、非経産動物変化数算定手段により算定された第１の初産動物の数、初妊動物変化数算定手段により算定された第２の初産動物の数、ならびに記憶手段により記憶された期首の経産動物の数および分娩間隔をもとに、前記所定の期間内に産まれうる出生動物の数を算定する出生動物数算定手段と、非経産動物変化数算定手段により算定された第１の非経産動物の数、および出生動物数算定手段により算定された出生動物の数をもとに、前記所定の期間における期末の非経産動物の数を算定する期末非経産動物数算定手段とを備えることを特徴とするシミュレーション装置に関する。
【０００８】
本発明によれば、前記所定の期間における期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに該動物群における分娩間隔をもとに、期末の経産動物の数、期末の非経産動物の数、期末の初妊動物の数の算定が可能となり、動物群の繁殖機能を考慮した動物群評価が可能となる。
【０００９】
本発明では、記憶手段は、前記所定の期間内に期首の非経産動物が初妊動物になる割合である初妊動物率、前記所定の期間内に非経産動物が妊娠を経て初産動物になる割合である第１の初産動物率、および前記所定の期間内に期首の初妊動物が初産動物になる割合である第２の初産動物率とをさらに記憶し、非経産動物変化数算定手段は、記憶手段により記憶された初妊動物率および期首の非経産動物の数をもとに、期末の初妊動物の数を算定し、かつ記憶手段により記憶された第１の初産動物率および期首の非経産動物の数をもとに、第１の初産動物の数を算定し、初妊動物変化数算定手段は、記憶手段により記憶された第２の初産動物率および期首の初妊動物の数をもとに、第２の初産動物の数を算定することが好ましい。
【００１０】
ここでは、それぞれの動物群に固有の特性である、初妊動物率、第１の初産動物率および第２の初産動物率が考慮されることとなる。すなわち、期末の初妊動物の数は初妊動物率および期首の非経産動物の数をもとに算定され、第１の初産動物の数は第１の初産動物率および期首の非経産動物の数をもとに算定され、さらに第２の初産動物の数は、第２の初産動物率および期首の初妊動物の数をもとに算定される。したがって、それぞれの動物群に固有の特性を考慮した動物群評価が可能となり、より現実の牧場等の経営に近い経営シミュレーションが可能となる。
【００１１】
本発明では、記憶手段は、経産動物のうち搾乳が可能な動物の割合である搾乳動物率をさらに記憶し、上記シミュレーション装置は、さらに、記憶手段により記憶された期首の経産動物の数、経産動物変化数算定手段により算定された期末の経産動物の数および搾乳動物率をもとに、搾乳が可能な搾乳動物の数を算定する搾乳動物数算定手段と、搾乳動物数算定手段により算定された搾乳動物の数をもとに、前記所定の期間内において出荷が可能な乳量を算定する乳量算定手段とを備えることが好ましい。
【００１２】
このような構成とすることで、期首の経産動物の数、期末の経産動物の数および搾乳動物率をもとに、搾乳が可能な搾乳動物の数が算定され、さらに搾乳動物の数をもとに、所定の期間内において出荷が可能な乳量が算定されるので、乳牛の繁殖機能を考慮した牛群評価をもとに、牛群の乳量の供用力を評価することが可能となり、より現実の牧場等の経営に近い経営シミュレーションが可能となる。
【００１３】
本発明では、上記シミュレーション装置は、さらに、前記所定の期間の整数倍であるシミュレーション期間についての利用者からの入力を受け付ける入力受付手段と、記憶手段に記憶された期首の経産動物の数、期首の非経産動物の数、および期首の初妊動物の数を、それぞれ、経産動物変化数算定手段により算定された期末の経産動物の数、期末非経産動物数算定手段により算定された期末の非経産動物の数、および非経産動物変化数算定手段により算定された期末の初妊動物の数で更新する更新手段とを備え、更新された期首の経産動物の数、期首の非経産動物の数、および期首の初妊動物の数をもととした、非経産動物変化数算定手段、初妊動物変化数算定手段、経産動物変化数算定手段、出生動物数算定手段および期末非経産動物数算定手段による算定、および更新手段による更新を、前記整数に対応する回数繰り返すことにより、シミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、および期末の初妊動物の数を算定することが好ましい。
【００１４】
このように構成することで、利用者が入力したシミュレーション期間（例えば、所定の期間が１年であり、シミュレーション期間が１０年）をもとに、期末の経産動物の数、期末の非経産動物の数、および期末の初妊動物の数の算定をシミュレーション期間の年数だけ繰り返すことにより、シミュレーション期間後における経産動物の数、非経産動物の数および期末の初妊動物の数を算定することが可能となる。したがって、利用者が入力した所望の期間が経過した後の動物群を評価することができ、長期にわたる牧場等の成長の潜在性を計数化、指標化することが可能となる。
【００１５】
本発明では、前記分娩間隔は、仮に設定された仮分娩間隔であって、上記シミュレーション装置は、さらに、利用者の入力に従って、シミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、または期末の初妊動物の数についての目標数を受け付ける目標数受付手段と、仮分娩間隔をもとに、非経産動物変化数算定手段、初妊動物変化数算定手段、経産動物変化数算定手段、出生動物数算定手段および期末非経産動物数算定手段による算定、および更新手段による更新を、前記整数に対応する回数繰り返すことにより、シミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、および期末の初妊動物の数を算定する仮算定手段と、仮算定手段により算定されたシミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、または期末の初妊動物の数が、目標数受付手段により受け付けられた目標数に達しているかについて判定する目標達成可否判定手段と、目標達成可否判定手段により目標数に達していないと判定された場合は、より短く設定された仮分娩間隔を新たに設定し、目標数に達していると判定された場合は、より長く設定された仮分娩間隔を新たに設定する仮分娩間隔設定手段と、仮分娩間隔設定手段により設定された新たな仮分娩間隔をもとに、仮算定手段によりシミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、および期末の初妊動物の数を算定し、仮分娩間隔が目標達成可否判定手段により目標数に達していると判定される最長の分娩間隔に達するまで仮分娩間隔設定手段による設定、仮算定手段による算定、目標達成可否判定手段による判定を繰り返し、目標数を達成するための最適な分娩間隔を特定する最適分娩間隔特定手段とを備えることが好ましい。
【００１６】
このような構成することで、目標数を達成するために最適な分娩間隔が特定されるため、利用者に対し、さらに有効な牧場経営の指標を提供することが可能となる。
【００１７】
本発明はさらに、複数の動物が存在する所定の動物群内における、分娩の経験がある経産動物の数、分娩の経験がなく妊娠もしていない非経産動物の数、および初めての妊娠をしている初妊動物の数についての所定の期間における変動をコンピュータ装置に算定させるシミュレーション方法であって、前記所定の期間における期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに該動物群における経産動物の分娩間隔を記憶領域に記憶させるステップと、記憶領域に記憶された期首の非経産動物の数をもとに、前記所定の期間内に妊娠を経て初めての分娩を行なう第１の初産動物の数、前記所定の期間内に初妊動物となる期末の初妊動物の数、および前記所定の期間内に変化しない第１の非経産動物の数を算定させるステップと、記憶領域に記憶された期首の初妊動物の数をもとに、前記所定の期間内に初めての分娩を行なう第２の初産動物の数を算定させるステップと、記憶領域に記憶された期首の経産動物の数、算定された第１の初産動物の数および算定された第２の初産動物の数をもとに、期末の経産動物の数を算定させるステップと、算定された第１の初産動物の数、算定された第２の初産動物の数、ならびに記憶領域に記憶された期首の経産動物の数および分娩間隔をもとに、前記所定の期間内に産まれうる出生動物の数を算定させるステップと、算定された第１の非経産動物の数、および算定された出生動物の数をもとに、前記所定の期間における期末の非経産動物の数を算定させるステップとから構成されることを特徴とするシミュレーション方法に関する。
【００１８】
本発明はさらに、複数の動物が存在する所定の動物群内における、分娩の経験がある経産動物の数、分娩の経験がなく妊娠もしていない非経産動物の数、および初めての妊娠をしている初妊動物の数についての所定の期間における変動を算定するサーバ装置と、通信ネットワークを介して前記サーバ装置と通信が可能な端末装置とを備えるシミュレーションシステムであって、端末装置が、利用者の入力に従って、前記所定の期間における期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに該動物群における分娩間隔の入力を受け付ける入力受付手段と、入力受付手段により受け付けられた期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに分娩間隔をサーバ装置に送信する情報送信手段とを備え、サーバ装置が、端末装置から送信された期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに分娩間隔を受信する情報受信手段と、情報受信手段により受信された期首の非経産動物の数をもとに、前記所定の期間内に妊娠を経て初めての分娩を行なう第１の初産動物の数、前記所定の期間内に初妊動物となる期末の初妊動物の数、および前記所定の期間内に変化しない第１の非経産動物の数を算定する非経産動物変化数算定手段と、情報受信手段により受信された期首の初妊動物の数をもとに、前記所定の期間内に初めての分娩を行なう第２の初産動物の数を算定する初妊動物変化数算定手段と、情報受信手段により受信された期首の経産動物の数、非経産動物変化数算定手段により算定された第１の初産動物の数および初妊動物変化数算定手段により算定された第２の初産動物の数をもとに、期末の経産動物の数を算定する経産動物変化数算定手段と、非経産動物変化数算定手段により算定された第１の初産動物の数、初妊動物変化数算定手段により算定された第２の初産動物の数、ならびに情報受信手段により受信された期首の経産動物の数および分娩間隔をもとに、前記所定の期間内に産まれうる出生動物の数を算定する出生動物数算定手段と、非経産動物変化数算定手段により算定された第１の非経産動物の数、および出生動物数算定手段により算定された出生動物の数をもとに、前記所定の期間における期末の非経産動物の数を算定する期末非経産動物数算定手段とを備えることを特徴とするシミュレーションシステムに関する。
【図面の簡単な説明】
【００１９】
【図１】本発明の実施の形態にかかるシミュレーションシステムの構成を表すブロック図である。
【図２】本発明の実施の形態にかかる牛群の飼養動態プロセスを表す図である。
【図３】本発明の実施の形態にかかる変化率テーブルを表す図である。
【図４】本発明の実施の形態にかかる牛群評価テーブルを表す図である。
【図５】本発明の実施の形態にかかる単価テーブルを表す図である。
【図６】本発明の実施の形態にかかる経営評価テーブルを表す図である。
【図７】本発明の実施の形態にかかる牛群の評価処理についてのフローチャートの一例である。
【図８】本発明の実施の形態にかかる最適分娩間隔の評価処理についてのフローチャートの一例である。
【図９】本発明の実施の形態にかかる目標達成可否判定処理を表わすフローチャートの一例である。
【図１０】本発明の実施の形態にかかる端末の表示画面の変化を表す図である。
【発明を実施するための形態】
【００２０】
図１は、本発明の実施の形態にかかるシミュレーションシステムの構成を表すブロック図である。シミュレーション装置１は、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）１１、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）１２、ハードディスク・ドライブ（ＨＤＤ）１３、通信インターフェイス１４からなり、それぞれ内部バスにより接続されている。
【００２１】
ＣＰＵ１１はＨＤＤ１３に格納されたプログラムを実行し、シミュレーション装置１の制御を行なう。ＲＡＭ１２はＣＰＵ１１のワークエリアである。ＨＤＤ１３は、プログラムやデータを保存するための記憶領域である。ＣＰＵ１１は、シミュレーション装置１の制御に必要なデータをＲＡＭ１２から読み出して処理を行なう。
【００２２】
通信インターフェイス１４は有線により通信ネットワーク１５に接続されており、他の端末１６との通信が可能である。また、通信インターフェイス１４は、アンテナを有する通信ユニットを介して、無線による通信を行なうことも可能である。例えば、端末１６にて、牧場経営者や牧場経営者への融資を検討している金融機関がインターネット上で所定のウェブサイトへアクセスし、ＩＤ及びパスワードを入力する。シミュレーション装置１では、通信インターフェイス１４を介してＩＤ及びパスワードを受信し、ログインが可能であるかどうかの認証を行なう。ＨＤＤ１３に設定された利用者データベースに記憶されたＩＤ及びパスワードが一致している場合は、ログインが可能となり、端末１６にてシミュレーションシステムを利用することが可能となる。
【００２３】
シミュレーションシステムへのログインが行なわれると、利用者は端末１６に設置されたキーボード等の入力装置へ、所定のデータ、例えば、初期の経産牛の数、初期の育成子牛の数、分娩間隔等に関する情報を入力することが可能である。端末１６の入力装置から入力された入力情報は通信インターフェイス１４を介してＲＡＭ１２に格納され、ＣＰＵ１１は入力情報をもとに各種の算定処理を実行する。算定処理の実行が終了すると、通信インターフェイス１４を介して、算定の結果が端末１６に送信される。端末１６にて算定の結果を受信すると、端末１６の表示装置にシミュレーションの結果が表示される。
【００２４】
次に、本発明の実施の形態にかかるシミュレーションシステムにおける、シミュレーションを行なう動物群の飼養動態について説明する。図２は、本発明の実施の形態にかかる動物群の飼養動態プロセスの一例を表す図である。この実施の形態において、動物群は牧場で飼養されている牛群であり、牛群には育成子牛２１、初妊牛２２、経産牛２３、初産牛２４および搾乳牛２５が含まれる。育成子牛２１は、出生してから最初に受胎するまでの雌子牛を指す。なお、この実施の形態では、雄子牛は出生して間もなく販売される。育成子牛２１は、出生してから１３〜１５か月経過すると、初めての種付けが行なわれる。種付けが成功し受胎した育成子牛２１は、初めての妊娠をしている牛である初妊牛２２となる。初妊牛２２は受胎から１０か月程度経過すると分娩を行ない、初妊牛２２から、初めての分娩を行なった牛である初産牛２４となる。なお、初産牛２４は、分娩の経験がある牛である経産牛２３に含まれる概念である。
【００２５】
初産牛２４を含む経産牛２３は、分娩後ほぼ２か月経過すると、次の種付けが行なわれる。種付けが成功し受胎した経産牛２３は、受胎から１０か月程度経過すると次の分娩を行なう。分娩間隔２６は、分娩と次の分娩の間の日数であり、経産牛ごとに異なるものであるが、本発明では牧場ごとの平均の日数を指すものとして取り扱う。なお、分娩間隔２６は、牧場の経営方針により変わってくるものであり、牧場ごとに異なる値である。分娩を行なった経産牛２３は、人間の飲用が可能な生乳を泌乳するため、搾乳が行なわれる。搾乳は、分娩の数日後から乾乳するまで行なわれる。乾乳とは、経産牛２３を次の分娩に備えさせるために、次の分娩予定日の約６０日前から、搾乳を止めることである。搾乳される経産牛２３は、搾乳牛２５にも分類される。
【００２６】
経産牛２３は、数回の分娩を行なうことができ、分娩後から乾乳まで搾乳牛２５として生乳生産に寄与するが、分娩回数を重ねるもしくは単に月齢を重ねることにより、泌乳量が減少する。泌乳量が一定量よりも少なくなるなどの条件に合致した牛は、廃用される。なお、廃用もしくは事故などによる死亡は、その他の育成子牛２１や初妊牛２２などにも起こりうることであり、このような牛は牛群から除かれる。
【００２７】
また、牛群の牛は販売することも可能である。これとは反対に、育成子牛などを外部から購入する場合もあり、この場合は購入した牛は対応する育成子牛２１、初妊牛２２や経産牛２３などに参入させることとなる。
【００２８】
図３は、本発明の実施の形態にかかる変化率テーブルの一例を表す図である。変化率テーブル４０は、シミュレーション装置１のＨＤＤ１３に格納されており、シミュレーションの開始にともないＲＡＭ１２に読み出される。変化率テーブル４０には、分娩間隔に関するパラメータ４１、搾乳牛率４２、育成子牛（期首）からの変化率４３、初妊牛（期首）からの変化率４４、初産牛（当期）の変化率４５、初妊牛（当期）の変化率４６、経産牛の変化率４７、出生率４８、雄子牛変化率４９、雌子牛変化率５０を含む項目が設定されている。これらの確率は、過去の統計をもとに予め設定したり、牧場の経営方針に応じて利用者が任意に設定することも可能である。なお、当期とは、これから算定処理を行なう期（所定の期間）を指す。
【００２９】
分娩間隔に関するパラメータ４１は、分娩間隔４１ａおよび分娩間隔係数４１ｂを含む情報である。分娩間隔４１ａは、利用者によって入力された分娩間隔が設定される。分娩間隔係数４１ｂは、一つの期（１年）の日数である３６５を分娩間隔４１ａで割った値である。経産牛１頭は分娩間隔ごとに１回の分娩を行なうため、分娩間隔係数４１ｂは一つの期における経産牛１頭が行なう分娩の回数とみなすことができる。したがって、期首における経産牛数に分娩間隔係数４１ｂを掛けることにより、当期において経産牛群が行なう分娩の回数を算定することができる。
【００３０】
搾乳牛率４２は、経産牛に占める搾乳牛の比率を表す数値であり、過去の統計をもとに予め設定されるものであるが、分娩間隔に応じて変動させることも可能である。搾乳牛数は経産牛数に搾乳牛率４２を掛けることで算定される。搾乳牛数は、後述するように生乳生産量を算定するために用いられる。
【００３１】
育成子牛（期首）からの変化率４３には、期首に育成子牛である牛が当期において、変化する確率が設定されている。育成子牛→初妊→初産牛率４３ａは、期首に育成子牛である牛が当期に初妊を経て初産を行なって初産牛となる確率である。育成子牛→初妊牛率４３ｂは、当期に初妊牛になる確率であり、育成子牛→育成子牛率４３ｃは、期末まで育成子牛のままである確率である。初妊牛（育成→）販売率４３ｄは、当期に初妊となってから販売し、期末には牛群に残らない牛となる確率である。廃用・死亡（育成（期首））率４３ｅは、当期に廃用もしくは死亡したため、期末に牛群に残らない牛となる確率である。なお、育成子牛→初妊→初産牛率４３ａ、育成子牛→初妊牛率４３ｂ、育成子牛→育成子牛率４３ｃ、初妊牛（育成→）販売率４３ｄ、廃用・死亡（育成（期首））率４３ｅの合計は１に等しくなる。
【００３２】
初妊牛（期首）からの変化率４４は、期首に初妊牛である牛が当期においてそれぞれの牛に変化する確率が設定された項目であり、当期に初産を行なう初産牛となる確率である初妊牛→初産牛率４４ａ、当期に廃用もしくは死亡したため期末に牛群に残らない牛となる確率である廃用・死亡（初妊牛（期首））率４４ｂが含まれる。初妊牛→初産牛率４４ａ、廃用・死亡（初妊牛（期首））率４４ｂの合計は１に等しくなる。
【００３３】
初産牛（当期）の変化率４５には、当期に初産牛となった牛が期末までに変化する確率が設定されている。当期に初産牛となった牛とは、期首に育成子牛であって当期に初妊を経て初産を行なって初産牛となった牛、期首に初妊牛であって当期に初産を行なって初産牛となった牛を含む概念である。したがって、育成子牛（期首）から初産牛となった牛数、および初妊牛（期首）から初産牛となった牛数の合計に初産牛（当期）の変化率４５を乗ずることで、当期に初産牛となった牛の期末までの変化が算定される。初産牛（当期）の変化率４５には、当期に廃用もしくは死亡する確率である廃用・死亡（初産牛）率４５ａ、および期末に経産牛に参入する確率である初産牛参入率４５ｂが含まれる。廃用・死亡（初産牛）率４５ａ、初産牛参入率４５ｂの合計は１に等しくなる。
【００３４】
初妊牛（当期）の変化率４６は、当期に育成子牛（期首）から初妊牛となった牛が期末までに変化する確率が設定されており、当期に廃用もしくは死亡する確率である廃用・死亡（初妊牛（育成→））率４６ａ、および期末まで初妊牛として残る確率である初妊牛（期末）率４６ｂが含まれる。廃用・死亡（初妊牛（育成→））率４６ａ、および初妊牛（期末）率４６ｂの合計は１に等しくなる。
【００３５】
経産牛の変化率４７は、期首に経産牛であった牛が期末までに変化する確率が設定されており、当期に販売される経産牛の確率である経産牛販売率４７ａ、当期に死亡する経産牛の確率である経産牛死亡率４７ｂ、および期末まで牛群に残る経産牛の確率である経産牛（期末）率４７ｃが含まれる。経産牛販売率４７ａ、経産牛死亡率４７ｂ、経産牛（期末）率４７ｃの合計は１に等しくなる。
【００３６】
次に、分娩及び出生に関する変化率について説明する。出生率４８は、当期に出生するそれぞれの子牛の確率を設定したものであり、当期に出生する子牛の数は、当期に分娩を行なう牛数に等しいとして算定される。したがって、経産牛のうち当期に分娩を行なう牛数、および当期に分娩を行なった初産牛数の合計に出生率４８を乗ずることで、当期に出生する子牛の数が算定される。なお、初産牛数は、育成子牛（期首）から初産牛となった牛数および初妊牛（期首）から初産牛となった牛数の合計である。出生率４８には、雄子牛が生まれる確率である雄子牛出生率４８ａ、雌子牛が生まれる確率である雌子牛出生率４８ｂ、および死産となる確率である子牛死産率４８ｃが含まれる。雄子牛出生率４８ａ、雌子牛出生率４８ｂ、子牛死産率４８ｃの合計は１に等しくなる。
【００３７】
雄子牛変化率４９は、当期に出生した雄子牛の数が、販売等により変化する確率を表すものである。本実施の形態において、出生した雄子牛は、その種類に応じた単価で販売される。雄子牛変化率４９には、乳用雄として出生し販売される確率である乳用雄販売率４９ａ、和牛として出生し販売される確率である和牛販売率４９ｂおよびＦ１牛として出生し販売される確率であるＦ１牛販売率４９ｃが含まれる。乳用雄販売率４９ａ、和牛販売率４９ｂ、Ｆ１牛販売率４９ｃの合計は１に等しくなる。
【００３８】
雌子牛変化率５０は、当期に出生した雌子牛の数が、販売等により変化する確率を表すものである。本実施の形態において、雌子牛は１３〜１５か月で受胎させるまで育成子牛として飼育される。雌子牛変化率５０は、育成子牛として飼育される確率である育成子牛飼育率５０ａ、および廃用もしくは死亡する確率である廃用・死亡（育成（出生））率５０ｂが含まれる。育成子牛飼育率５０ａ、廃用・死亡（育成（出生））率５０ｂの合計は１に等しくなる。
【００３９】
次に、牛群評価テーブルについて説明する。図４は、本発明の実施の形態にかかる牛群評価テーブルの一例を表す図である。牛群評価テーブル６０は、シミュレーション装置１のＲＡＭ１２に設けられ、利用者によって入力された初期牛数（初期育成子牛数、初期初妊牛数および初期経産牛数を含む）、および後述する牛群評価処理の結果が記憶される。牛群評価テーブル６０には、期首牛数６１、育成子牛（期首）からの変化牛数６２、初妊牛（期首）からの変化牛数６３、初産牛（当期）の変化牛数６４、初妊牛（当期）の変化牛数６５、経産牛の変化牛数６６、出生牛数６７、出生雄牛数６８、出生雌牛数６９、期末牛数７０、および搾乳牛数７１を含む項目が記憶される。
【００４０】
本発明の実施の形態において、育成子牛数、初妊牛数や経産牛数などの牛数は、母集団の牛数に対し、変化率テーブル４０に設定された変化率を掛けることにより算定される。そのため、算定された牛数は小数を含む値になる場合がある。この場合は以下のように端数処理が行われる。すでに述べたように、育成子牛（期首）からの変化率４３、初妊牛（期首）からの変化率４４等に含まれる各変化率は、いずれもその合計が１になる。この場合は、同じ母集団から算定された各牛数の小数部分を切り下げ、切り下げた各牛数の合計と母集団の数とを比較し、各牛数の合計が母集団の数と等しくなるように各牛数の小数部分が大きいものから順に切り上げる。なお、小数部分が等しい場合は、各牛数の値が大きいものから順に切り上げる。
【００４１】
例えば、母集団の数である育成子牛（期首）数６１ａが「１４０」であり、育成子牛→初妊→初産牛率４３ａ、育成子牛→初妊牛率４３ｂ、育成子牛→育成子牛率４３ｃ、初妊牛（育成→）販売率４３ｄ、廃用・死亡（育成（期首））率４３ｅを育成子牛（期首）数６１ａに乗ずることで得られる値が、それぞれ、「１７．８」、「８７．７」、「１９．２」、「１２．２」、「３．２」である場合、その合計が「１４０」となるように小数部分が大きい「１７．８」、「８７．７」が切り上げられて、「１８」、「８８」となり、その他は「１９」、「１２」、「３」と小数点以下が切り下げられる。また、設定された変化率が分娩間隔に関するパラメータ４１または搾乳牛率４２である場合は、小数点以下は切り下げる。
【００４２】
続いて、牛群評価テーブル６０に含まれる情報について詳しく説明する。まず、期首牛数６１は、当期における算定処理のもととなる、期首のそれぞれの牛の数が記憶された項目であり、期首における育成子牛の数である育成子牛（期首）数６１ａ、期首における初妊牛の数である初妊牛（期首）数６１ｂおよび期首における経産牛の数である経産牛（期首）数６１ｃが含まれる。これら期首牛数６１は、当期が第１期であれば、利用者によって入力された初期牛数、すなわち初期育成子牛数が育成子牛（期首）数６１ａとして記憶され、同様に、初期初妊牛数が初妊牛（期首）数６１ｂとして、初期経産牛数が経産牛（期首）数６１ｃとして、それぞれ記憶される。当期が第２期以降であれば、前の期における期末牛数７０のそれぞれの値が記憶される。
【００４３】
育成子牛（期首）からの変化牛数６２は、育成子牛（期首）数６１ａに対して、変化率テーブル４０における育成子牛（期首）からの変化率４３で設定されたそれぞれの確率を乗じて算定される項目である。まず、育成子牛→初妊→初産牛数６２ａは、育成子牛（期首）６１ａに育成子牛→初妊→初産牛率４３ａを乗じて算定される。例えば第１期において、育成子牛→初妊→初産牛数６２ａは、育成子牛（期首）６１ａである「１４０」に、育成子牛→初妊→初産牛率４３ａである「０．１３」をかけ、端数処理された「１８」として算定される。以下同様に、育成子牛→初妊牛数６２ｂは育成子牛→初妊牛率４３ｂを、育成子牛→育成子牛数６２ｃは育成子牛→育成子牛率４３ｃを、初妊牛（育成→）販売数６２ｄは初妊牛（育成→）販売率４３ｄを、廃用・死亡（育成（期首））数６２ｅは廃用・死亡（育成（期首））率４３ｅを、それぞれ育成子牛（期首）数６１ａに掛けることで算定される牛数である。
【００４４】
初妊牛（期首）からの変化牛数６３は、初妊牛（期首）数６１ｂに対して変化率テーブル４０における初妊牛（期首）からの変化率４４で設定されたそれぞれの確率を乗じて算定される項目である。それぞれ、初妊牛→初産牛数６３ａは初妊牛→初産牛率４４ａを、廃用・死亡（初妊牛（期首））数６３ｂは廃用・死亡（初妊牛（期首））率４４ｂを、初妊牛（期首）数６１ｂにそれぞれ乗じて算定される。
【００４５】
初産牛（当期）の変化牛数６４は、当期に初産牛となった牛が期末までに変化した牛数が算定される項目である。当期に初産牛となった牛数は、育成子牛→初妊→初産牛数６２ａおよび初妊牛→初産牛数６３ａの合計である。この合計に変化率テーブル４０の初産牛（当期）の変化率４５に設定されたそれぞれの確率を乗ずることで、初産牛（当期）の変化牛数６４が算定される。廃用・死亡（初産牛）数６４ａは廃用・死亡（初産牛）率４５ａを、初産牛参入数６４ｂは初産牛参入率４５ｂを、それぞれ乗ずることにより算定される。例えば第１期において、育成子牛→初妊→初産牛数６２ａである「１８」および初妊牛→初産牛数６３ａである「７６」の合計である「９４」に初産牛参入率４５ｂ「０．９５」を乗じて端数処理された「８９」が、初産牛参入数６４ｂとなる。
【００４６】
初妊牛（当期）の変化牛数６５は、当期に初妊牛となった牛が期末までに変化した牛数が算定される項目である。初妊牛（当期）の変化牛数６５は、育成子牛→初妊牛数６２ｂに初妊牛（当期）の変化率４６の各確率を乗じて算定される。廃用・死亡（初妊牛（育成→））数６５ａは廃用・死亡（初妊牛（育成→））率４６ａを、初妊牛（期末）数６５ｂは育成子牛→初妊牛数６２ｂに初妊牛（期末）率４６ｂを乗じることで算定される。
【００４７】
経産牛の変化牛数６６は、当期に経産牛である牛が期末までに変化する牛数として算定される。ここで、当期に経産牛である牛とは、期首に経産牛である牛に、当期に初産牛となった牛を合計したものである。まず、経産牛（期首）数６１ｃに初産牛参入数６４ｂを加算した経産牛全群数６６ａが算定される。この経産牛全群数６６ａに経産牛の変化率４７の各確率を掛けることで、経産牛の変化牛数６６が算定される。それぞれ、経産牛販売数６６ｂは経産牛販売率４７ａを、経産牛死亡数６６ｃは経産牛死亡率４７ｂを、および経産牛（期末）数６６ｄは経産牛（期末）率４７ｃを、経産牛全群数６６ａに乗じることで算定される。
【００４８】
次に、出生牛数６７について説明する。本実施の形態において、出生牛数６７の各牛数の合計は、当期において分娩を行なう全牛数に等しいとして算定される。すなわち、経産牛のうち当期に分娩を行なう牛数、育成子牛→初妊→初産牛数６２ａおよび初妊牛→初産牛数６３ａの合計である。ここで、経産牛であって当期に分娩を行なう牛数は、経産牛（期首）数６１ｃに分娩間隔係数４１ｂを乗じた値を小数点第１位で四捨五入することにより算定される。続いて、分娩頭数（調整）６７ａが、経産牛のうち当期に分娩を行なう牛数、育成子牛→初妊→初産牛数６２ａおよび初妊牛→初産牛数６３ａの合計として算定される。分娩頭数（調整）６７ａを出生牛数６７の各牛数の合計とみなし、これに出生率４８に設定された各確率を乗じて、それぞれの牛数が算定される。すなわち、雄子牛出生数６７ｂは雄子牛出生率４８ａ、雌子牛出生数６７ｃは雌子牛出生率４８ｂ、子牛死産数６７ｄは子牛死産率４８ｃを、それぞれ分娩頭数（調整）６７ａに乗じることで算定される。
【００４９】
出生雄牛数６８は、雄子牛出生数６７ｂに、雄子牛変化率４９に設定された各確率を乗じて算定される。すなわち、乳用雄販売数６８ａは乳用雄販売率４９ａ、和牛販売数６８ｂは和牛販売率４９ｂ、Ｆ１牛販売数６８ｃはＦ１牛販売率４９ｃをそれぞれ乗じて算定される。また、出生雌牛数６９は、変化率テーブル４０の雌子牛変化率５０の各確率、すなわち育成子牛飼育数６９ａは育成子牛飼育率５０ａを、廃用・死亡（育成（出生））数６９ｂは廃用・死亡（育成（出生））率５０ｂを、雌子牛出生数６７ｃに乗じて算定される。
【００５０】
次に、期末牛数７０について説明する。期末牛数７０には、育成子牛（期末）数７０ａ、初妊牛（期末）数７０ｂおよび経産牛（期末）数７０ｃを含む情報が記憶される。育成子牛（期末）数７０ａは、当期に育成子牛のままであった育成子牛→育成子牛数６２ｃおよび当期に出生し飼育された育成子牛飼育数６９ａの合計である。初妊牛（期末）数７０ｂは初妊牛（期末）数６５ｂに等しく、経産牛（期末）数７０ｃは経産牛（期末）数６６ｄに等しい。
【００５１】
また、牛群評価テーブル６０には、搾乳牛数７１も含まれる。搾乳牛数７１は、経産牛数に搾乳牛率４２を乗じて算定される。ただし、牛群評価テーブル７０における搾乳牛数７１の項目は、後述する経営評価処理において、年間生乳生産量を算定するために用いられるため、まず経産牛（期首）数６１ｃに搾乳牛率４２を乗じた搾乳牛（期首）数７１ａ、および経産牛（期末）数７０ｃに搾乳牛率を乗じた搾乳牛（期末）数７１ｂが、それぞれ算定される。続いて、搾乳牛（期首）数７１ａと搾乳牛（期末）数７１ｂの平均である生乳生産牛数７１ｃが算定され、後述する経営評価処理における算定に用いられる。なお、生乳生産牛数７１ｃは年間生乳生産量を算定するための平均値に過ぎないため、端数処理は行なわれない。
【００５２】
次に、本発明の実施の形態における、経営評価処理において用いられる数値について、図５及び図６を用いて詳細に説明する。図５は、本発明の実施の形態にかかる単価テーブルを表す図である。図６は、本発明の実施の形態にかかる経営評価テーブルを表す図である。
【００５３】
単価テーブル８０は、シミュレーション装置１のＨＤＤ１３に格納されており、シミュレーションの開始にともないＲＡＭ１２に読み出される。単価テーブル８０には、生乳生産単価８１および個体単価８２を含む項目が格納されている。生乳生産単価８１は、搾乳牛１頭当たりの年間生乳生産量（単位：ｋｇ／頭）である個体乳量８１ａ、および生乳１ｋｇの単価（単位：円／ｋｇ）である乳価８１ｂを含む項目である。個体単価８２は、各種類の牛１頭を販売した時の単価（単位：円／頭）であって、乳用雄単価８２ａ、和牛雄単価８２ｂ、Ｆ１雄単価８２ｃ、初妊牛単価８２ｄおよび経産牛単価８２ｅを含む項目である。
【００５４】
経営評価テーブル９０は、経営評価処理において算定された数値が記憶されるテーブルであり、ＲＡＭ１２に設けられる。経営評価テーブル９０は、生乳生産金額９１、個体販売金額９２および売上高９３を含む項目が記憶される。生乳生産金額９１および個体販売金額９２は、単価テーブル８０で設定された単価および牛群評価テーブル６０に記憶された牛数をもとに算定される。
【００５５】
生乳生産金額９１は、生乳生産量９１ａおよび生乳販売金額９１ｂを含む項目である。生乳生産量９１ａは、当期に生産される生乳の合計量であり、生乳生産牛数７１ｃに個体乳量８１ａを乗じた値として算定される。生乳販売金額９１ｂは、生乳生産量９１ａに乳価８１ｂを乗じて算定されたものであり、当期における生乳による販売金額の総額を表す。なお、図６において、生乳販売金額９１ｂは単位を千円として表示されている。
【００５６】
個体販売金額９２は、当期に販売した各種類の牛の販売金額の合計が記憶された項目である。乳用雄販売９２ａは乳用雄販売数６８ａに乳用雄単価８２ａを、和牛雄販売９２ｂは和牛販売数６８ｂに和牛雄単価８２ｂを、Ｆ１雄販売９２ｃはＦ１牛販売数６８ｃにＦ１雄単価８２ｃを、初妊牛販売９２ｄは初妊牛（育成→）販売数６２ｄに初妊牛単価８２ｄを、経産牛販売９２ｅは経産牛販売数６６ｂに経産牛単価８２ｅを、それぞれ乗じることで算定される。個体販売合計９２ｆは、乳用雄販売９２ａから経産牛販売９２ｅまでの値の合計である。なお、図６において、個体販売金額９２に含まれる金額は全て単位を千円として表示されている。
【００５７】
売上高９３は、当期における総売上高に相当し、生乳販売金額９１ｂおよび個体販売合計９２ｆの合計として算定される。図６において、売上高９３は単位を千円として表示されている。
【００５８】
次に、本実施の形態にかかる、牛群の評価処理について説明する。図７は、本発明の実施の形態にかかる牛群の評価処理についてのフローチャートの一例を表す図である。まず、端末１６がシミュレーション装置１へのログインが行なわれると、利用者は端末１６の入力装置を通じ、シミュレーションを行なうために必要な初期値の入力を行なう。シミュレーション装置１は通信インターフェイス１４を介して、この入力を受け付ける（ステップＳ１）。受け付けられた初期値は、ＲＡＭ１２に格納される。
【００５９】
この実施の形態において、初期値は、シミュレーションを行なう期間であるシミュレーション期間（例えば５年）、シミュレーション期間の開始時におけるそれぞれの牛群に含まれる牛の頭数である初期牛数（初期育成子牛数、初期初妊牛数、初期経産牛数）、および分娩間隔などが含まれる。シミュレーション期間は、１年を単位として入力され、１以上の整数が入力可能である（ここでは、ｎ＞１であるｎ年が入力されたものとする）。シミュレーション期間は１年ごとに第１期から第ｎ期に分割され、第１期から順に牛群評価処理が行われる。初期牛数は、シミュレーション期間の開始時における育成子牛２１、初妊牛２２および経産牛２３の数を含み、それぞれ０以上の整数値で入力される。
【００６０】
ステップＳ１において初期値の入力が受け付けられると、ＣＰＵ１１は、ＲＡＭ１２に設けられた牛群評価テーブル６０において、期首牛数６１を記憶する（ステップＳ２）。ここでは、第１期の算定を行なうものとし、初期牛数が期首牛数６１として記憶される。
【００６１】
続いて、ＣＰＵ１１は、当期における育成子牛２１から変化した牛の数を算定する（ステップＳ３）。ステップＳ３は、牛群評価テーブル６０における育成子牛（期首）からの変化牛数６２に含まれる項目を算定する処理である。前述したように、育成子牛（期首）からの変化牛数６２は、育成子牛（期首）数６１ａに育成子牛（期首）からの変化率４３を乗じて算定される。
【００６２】
次に、当期における初妊牛２２から変化した牛の数が算定される（ステップＳ４）。ステップＳ４は、牛群評価テーブル６０における初妊牛（期首）からの変化牛数６３に含まれる項目を算定する処理である。初妊牛（期首）からの変化牛数６３は、初妊牛（期首）数６１ｂに初妊牛（期首）からの変化率４４を乗じて算定される。続いて、当期に初産牛となった牛から変化した牛の数、および当期に初妊牛となった牛から変化した牛の数が算定される（ステップＳ５）。当期に初妊牛となった牛から変化した牛の数、および当期に初産牛となった牛から変化した牛の数とは、それぞれ牛群評価テーブル６０における初産牛（当期）の変化牛数６４、および初妊牛（当期）の変化牛数６５に相当し、これらがステップＳ５において算定される。
【００６３】
次に、当期における経産牛２３の数が算定される（ステップＳ６）。ステップＳ６は、牛群評価テーブル６０における経産牛の変化牛数６６に含まれる項目を算定する処理である。ステップＳ６までが終了すると、それまでのステップで算定された育成子牛（期首）からの変化牛数６２、初妊牛（期首）からの変化牛数６３および経産牛の変化牛数６６の算定結果をもとに、当期に出生する牛の数が算定される（ステップＳ７）。当期に出生する牛の数とは、牛群評価テーブル６０における出生牛数６７、出生雄牛数６８および出生雌牛数６９を含む概念である。
【００６４】
続いて、ＣＰＵ１１は、牛群評価テーブル６０に記憶された算定結果をもとに、期末における牛群の牛数を算定する（ステップＳ８）。ここで算定されるのは期末牛数７０および搾乳牛数７１である。
【００６５】
ステップＳ８が終了すると、ＣＰＵ１１は、経営評価処理を行なう（ステップＳ９）。経営評価処理は、牛群評価テーブル６０に記憶された当期の各牛数、およびＲＡＭ１２に読み出された単価テーブル８０に設定された販売単価にもとづいて行なわれる算定処理である。ステップＳ８において算定された数値は、当期の販売金額として、ＲＡＭ１２に設けられた経営評価テーブル９０に記憶される。
【００６６】
ステップＳ９までが終了すると、終了条件を満たすか否かが判定される（ステップＳ１０）。つまり、第ｎ期の算定が終了したか否かが判定される。上記の説明では、ｎ＞１であるｎ年のシミュレーション期間に対し、第１期の算定を行なっているため、ステップＳ１０においてＮｏとなり、ステップＳ２に戻る。この場合、ステップＳ２以降では第２期の算定を始めることになるため、ステップＳ２において、第２期の期首牛数６１として第１期の期末牛数７０がそれぞれ設定され、ステップＳ３〜Ｓ１０の処理が行なわれる。同様の処理により、第ｎ期までの算定が繰り返され、第ｎ期までの処理が終了すると、終了条件を満たすと判定され（ステップＳ１０においてＹｅｓ）、算定結果が端末１６に送信される（ステップＳ１１）。端末１６が算定結果を受信し、端末１６の表示画面に算定結果が表示されると、利用者は表示された算定結果をもとに様々な牛群評価が可能となる。
【００６７】
なお、算定結果をもとに、利用者が算定条件の変更を希望する場合は、利用者はその旨を端末１６に入力するように構成しても良い。変更可能な算定条件は、初期値に含まれるシミュレーション期間、初期牛数および分娩間隔が挙げられるが、これらのほか、経産牛の販売率などを変更可能としてもよい。入力された算定条件の変更に関する情報は端末装置１６からシミュレーション装置１に送信され、シミュレーション装置１のＣＰＵ１１がこれを受け付けると、初期値そのほかの算定条件が再設定され、牛群評価処理が再度行われる。
【００６８】
また、上記の説明では、ｎ年間における牛群変動をシミュレーションするため、ステップＳ１０において、第ｎ期の算定が終了したか否かを判定した。しかし、本実施の形態を応用することにより、所定の目標値を設定し、目標値を上回るために必要な年数をシミュレーションすることも可能である。この場合、ステップＳ１において、初期値として初期牛数および分娩間隔の入力は受け付けるが、シミュレーション期間は受け付けず、代わりに経産牛数、生乳生産量または売上高などの目標値の入力を受け付ける。受け付けた初期値をもとに牛群評価処理を行ない、算定された値が目標値を上回ったか否かをステップ１０で判定し、目標値を上回るまで算定処理を繰り返すことにより、必要な年数を算定することができる。なおこの場合、算定処理が膨大な回数になるのを回避するため、繰り返し処理の上限を設定しておくことが好ましい。
【００６９】
次に、最適分娩間隔の評価処理について説明する。この処理は、初期値（初期育成子牛数、初期初妊牛数および初期経産牛数を含む初期牛数ならびにシミュレーション期間を含む）に加えて、シミュレーション期間（例えば、ｎ年）後の目標値である目標牛数として育成子牛数、初妊牛数または経産牛数の入力を受け付け、ｎ年後において目標牛数を上回るための最長の分娩間隔を、牧場等の経営に最適な分娩間隔である最適分娩間隔として特定するものである。なお、最適分娩間隔は、目標牛数をｎ年後に上回ることが可能な最長の分娩間隔として算定されるため、ｎ−１年後には目標牛数に達せずｎ年後に初めて目標牛数を上回ることが可能な最長の分娩間隔として算定されることとなる。
【００７０】
牛は分娩間隔に従って分娩を行なうため、同じｎ年間で比較した場合、分娩間隔が短いほど出生牛数６７が増加する。出生牛数６７が増加すると、これに応じて次期以降における育成子牛数、初妊牛数および経産牛数が増加する。したがって、牛群を拡大するためには分娩間隔が短い方がよい。しかし、分娩間隔を短く管理するためには多大な労力やコストが必要である。本発明によって、設定された目標牛数をｎ年後に上回ることが可能な最長の分娩間隔を最適分娩間隔として特定し、利用者に提示することで、利用者は分娩間隔を管理した場合の利益とコストを判断することが可能となり、牧場等の経営判断を行なう際の指標とすることができる。
【００７１】
本発明における最適分娩間隔の特定は、３５０日から４８０日までの範囲でのみ処理が行なわれる。３５０日は、牛が分娩を行なってから次に受胎可能になるまでの期間および妊娠期間の合計であり、生理学的に取りうる最短の分娩間隔である。また、４８０日は、搾乳可能な最長期間の日数および乾乳期間として取りうる最長期間の日数の合計であり、すなわち分娩間隔として取りうる最長の日数である。このため、最適分娩間隔の特定処理を３５０日から４８０日の範囲の外まで行なっても、牧場等の経営指標として有用なものになり得ず、シミュレーション装置１の負荷増大にもつながるため、３５０日から４８０日までの範囲に限定して最適分娩間隔の特定処理が行なわれる。
【００７２】
最適分娩間隔の特定処理は、まず最適分娩間隔が３５０日から４８０日の範囲に含まれるか否かが判定される。続いて、最適分娩間隔が含まれる範囲の上限値および下限値を設定し、その中間値を仮分娩間隔として算定する。算定された仮分娩間隔および初期値（初期牛数およびシミュレーション期間を含む）を用いて牛群評価処理によりシミュレーション期間後の各牛数を算定し、算定された牛数が目標牛数を上回るか否かの判定が行なわれる。
【００７３】
続いて、最適分娩間隔の評価処理の具体的な手順について、図８および図９を用いて説明する。図８は、本発明の最適分娩間隔の評価処理についてのフローチャートの一例である。図９は、本発明の実施の形態にかかる目標達成可否判定処理を表わすフローチャートの一例である。なお、ここでは目標値として経産牛数が所定の値に設定されたとする。
【００７４】
まず、利用者は端末１６の入力装置を通じ、シミュレーションを行なうために必要な初期値および目標値の入力を行ない、シミュレーション装置１は通信インターフェイス１４を介してこの入力情報を受け付ける（ステップＳ２１）。受け付けられた初期値および目標値は、ＲＡＭ１２に格納される。
【００７５】
次に、ＣＰＵ１１は、仮分娩間隔を３５０日と設定し、仮分娩間隔、およびステップＳ２１で受け付けた初期値を用いて、経産牛数の算定を行なう（ステップＳ２２）。経産牛の算定は、図７のステップＳ２〜Ｓ１０の処理にもとづいて行なわれる（ステップＳ８の経営評価処理は行なわなくても良い）。続いて、仮分娩間隔で算定された経産牛数が目標値の経産牛数よりも少ないか否かの判定を行なう（ステップＳ２３）。ここで、仮分娩間隔を理論上の最短日数である３５０日としているので、算定された経産牛数が目標値の経産牛数よりも多いと判定された場合（ステップＳ２３においてＮｏ）、目標値の経産牛数は入力された初期値で達成可能であると判定されたことを意味する。
【００７６】
仮分娩間隔で算定された経産牛数が目標値の経産牛数よりも多いと判定された場合は、続いて最適分娩間隔の算定が行なわれる（ステップＳ２４）。最適分娩間隔の算定処理については、図９を用いて後に詳述する。最適分娩間隔の算定が行なわれると、ＣＰＵ１１は、通信インターフェイス１４を介し、端末１６に算定結果を送信する（ステップＳ２５）。
【００７７】
なお、ステップＳ２３において、仮分娩間隔で算定された経産牛数が目標値の経産牛数よりも小さいと判定された場合（ステップＳ２３においてＹｅｓ）、分娩間隔を理論上の最短日数である３５０日に設定してもなお目標値に達成しないこととなるため、分娩間隔の算定を行なわず、代わりに目標値の経産牛数に達成することが可能な最短年数を算定する処理が行なわれる（ステップＳ２６）。最短年数の算定は、図７のステップＳ２〜Ｓ１０が利用され、ステップＳ１０において算定された経産牛数が目標値の経産牛数を上回ると判定されるまで算定処理が繰り返される。算定された経産牛数が目標値の経産牛数を初めて上回った期が、ステップＳ２６における目標値の経産牛数に達成することが可能な最短年数として算定される。
【００７８】
続いて、図９を用いて、ステップＳ２４の最適分娩間隔算定処理について詳細に説明する。まず、ＣＰＵ１１は、仮分娩間隔を４８０日と設定し、ＲＡＭ１２に格納された初期値とともに、図７のステップＳ２〜Ｓ１０の処理により経産牛数を算定する（ステップＳ４１）。続いて、算定された経産牛数が、目標値の経産牛数以上であるか否かの判定を行なう（ステップＳ４２）。
【００７９】
分娩間隔を４８０として算定された経産牛数が目標値の経産牛数以上であった場合（ステップＳ４２におけるＹｅｓ）、最長の分娩間隔である４８０日であっても、目標値を達成できてしまうことになる。この場合は、最適分娩間隔が４８０日であると特定し（ステップＳ４３）、最適分娩間隔の算定処理を終了する。このように構成することで、有用でない指標の算定処理を省略することができるため、ＣＰＵ１１の負荷を軽減することが可能となる。なお、この場合、ステップＳ２５において、最適分娩間隔が４８０日であることが端末１６に送信される。
【００８０】
ステップＳ４２において算定された経産牛数が目標値の経産牛数よりも小さいと判定されると（ステップＳ４２においてＮｏ）、最適分娩間隔が含まれる可能性のある範囲の上限値および下限値が設定される（ステップＳ４４）。ここでは、最適分娩間隔は３５０日から４８０日の範囲に含まれるため、上限値は４８０日、下限値は３５０日に設定される。続いて上限値と下限値の中間値が算定される（ステップＳ４５）。ここで、中間値とは、上限値と下限値の平均であるが、値が小数を含む場合は小数部分を切り捨てた値である。
【００８１】
続いて、仮分娩間隔を中間値に設定し、図７のステップＳ２〜Ｓ１０の処理により経産牛数の算定を行なう（ステップＳ４６）。続いて、ステップＳ４６で算定された経産牛数が目標値の経産牛数と等しいか否かの判定を行なう（ステップＳ４７）。ここで算定された経産牛数が目標値の経産牛数と等しい場合は、最適分娩間隔は中間値であるとして算定される（ステップＳ４８）。
【００８２】
ステップＳ４６で算定された経産牛数と目標値の経産牛数が等しくないと判定された場合（ステップＳ４７においてＮｏ）、ステップＳ４６で算定された経産牛数が目標値の経産牛数よりも大きいか否かの判定を行なう（ステップＳ４９）。経産牛数が目標値の経産牛数よりも大きいと判定されれば（ステップＳ４９においてＹｅｓ）、最適分娩間隔は中間値と上限値の間に含まれることになるので、この中間値が新たな下限値として再設定される（ステップＳ５０）。すなわち、最適分娩間隔が特定されるまで、ステップＳ４５〜Ｓ５０までの処理は繰り返されることになるが、今回の処理がｍ回である場合に、ｍ＋１回目の処理における上限値はｍ回目と同じ上限値が用いられ、ｍ回目の中間値がｍ＋１回目の下限値として用いられる。算定された牛数が目標牛数を上回った場合、分娩間隔が、この仮分娩間隔より長い場合でも目標牛数を上回ることが可能であるので、最適分娩間隔は仮分娩間隔と上限値の間の値をとることになる。したがって、ステップＳ５０において、この中間値を新たな下限値として再設定し、ステップＳ４５〜Ｓ５２までの処理が再度行なわれることになる。
【００８３】
ステップＳ４９において、ステップＳ４６で算定された経産牛数が目標値の経産牛数よりも小さいと判定されれば、中間値が新たな上限値として再設定される（ステップＳ５１）。同様に、今回の処理がｍ回である場合に、ｍ＋１回目の処理における下限値はｍ回目と同じ下限値が用いられ、ｍ回目の中間値がｍ＋１回目の上限値として用いられる。同様に、算定された牛数が目標牛数を下回った場合、分娩間隔が仮分娩間隔よりも短くなければ目標牛数を達成できないため、最適分娩間隔は仮分娩間隔と下限値の間の値をとることになる。したがって、ステップＳ５０において、この中間値を新たな上限値として再設定し、ステップＳ４５〜Ｓ５２までの処理が再度行なわれることになる。
【００８４】
続いて、新たに設定された最適分娩間隔が含まれる範囲が、さらに限定される必要があるか否かを判定するため、新たに設定された上限値および下限値を用い、上限値が下限値＋１と等しいか否かの判定を行なう（ステップＳ５２）。等しくなければ（ステップＳ５２においてＮｏ）、最適分娩間隔が含まれる範囲をさらに限定する必要があるので、ステップＳ４５に戻る。新たに設定された上限値が下限値＋１と等しければ（ステップＳ５２においてＹｅｓ）、最適分娩間隔は整数値として算定されるため、これ以上の限定を行なう必要がなく処理が終了することになる。この場合、下限値を最適分娩間隔とすれば、経産牛数が目標値の経産牛数を上回るので、最適分娩間隔として下限値の値が算定され（ステップＳ５３）、最適分娩間隔算定処理は終了する。最適分娩間隔算定処理が終了すると、図８におけるステップＳ２５に戻る。
【００８５】
最適分娩間隔の特定は、シミュレーション期間であるｎ年の経過後に目標牛数の達成が可能である分娩間隔のうち、最長の日数を特定するものであり、ステップＳ４１〜Ｓ５３の処理を実行することで、ｎ−１年の経過後には目標牛数の達成ができないが、ｎ年の経過後に目標牛数の達成が可能な分娩間隔を特定することが可能となる。仮分娩間隔として、最適分娩間隔が含まれる可能性のある範囲の上限値と下限値の中間値が用いられ、ステップＳ４５〜Ｓ５２の処理を繰り返して、中間値が算定結果に応じて変化する（目標牛数に達していない場合は、より短い仮分娩間隔を新たに設定し、目標牛数に達している場合は、より長い仮分娩間隔を新たに設定する）ことで、最適分娩間隔が含まれる可能性のある範囲が限定され、最終的に最適分娩間隔が特定されることになる。このような処理を経ることでステップＳ２〜Ｓ１０の処理の回数は数回程度で済むことになり、シミュレーション装置１の処理負荷の軽減にもつながる。
【００８６】
なお、上記の実施の形態では目標値として育成子牛数、初妊牛数または経産牛数としたが、目標値として生乳生産量が設定されても良い。この場合、生乳生産量は経産牛数から算定されるため、目標生乳生産量から目標経産牛数を逆算し、目標経産牛数に達するように最適分娩間隔を算定すればよい。また、上記の実施の形態では、初期牛数と目標牛数を入力し、各変化率の値を固定の値として、最適な分娩間隔を求めることとしたが、初期牛数、目標牛数、特定の分娩間隔を利用者が入力し、最適な変化率（例えば、各牛の販売率等）を算定するような構成とすることも可能である。この場合、図８、図９について説明したのと同様の方法で、最適な変化率の算定が行なわれる。
【００８７】
図１０は、本発明の実施の形態にかかる端末１６の表示画面の変化を表す図である。この実施の形態では、利用者が端末１６からシミュレーション装置１に対してログインすると、シミュレーション装置１から端末１６に情報が送信され、図１０（ａ）のように「売上高を試算」と「分娩間隔を試算」の二つのシミュレーション方法を選択する画面が、端末１６の表示画面１６０に表示される。
【００８８】
利用者が「売上高を試算」を選択入力すると、その入力は端末１６からシミュレーション装置１に送信され、受信したシミュレーション装置１は、初期値を入力させるための画面情報を端末１６に送信し、端末１６はこの情報を受信して図１０（ｂ）のように表示画面１６０に表示する。ここで、利用者が初期値を入力し、入力ボタンを選択すると、入力された初期値は端末１６からシミュレーション装置１に送信される。シミュレーション装置１は、受信した初期値をもとに図７のステップＳ２〜Ｓ１０までの牛群評価処理を行ない、ステップＳ１０において算定結果を端末１６に送信する。
【００８９】
端末１６は、図１０（ｃ）のようにして、受信した算定結果を表示画面１６０に表示する。算定結果には、シミュレーション期間における売上高の推移を表したグラフ１６１、およびシミュレーション期間における売上高や各牛群の推移を表したテーブル１６２が表示される。利用者は表示画面１６０に表示されたグラフ１６１やテーブル１６２をもとに、牧場の経営評価を行なう。
【００９０】
なお、算定条件を変更し、再度シミュレーションを行なうことを利用者が希望する場合は、利用者は「条件を変更」ボタン１６３を選択入力することにより、算定条件の変更が可能である。ここで、変更可能な算定条件は初期値だけに限定されず、図１０（ｄ）のように、初妊牛や経産牛等の販売率を変更可能としてもよい。算定条件の変更に関する情報が端末１６に入力されると、入力された情報はシミュレーション装置１に送信され、受信したシミュレーション装置１は変更された算定条件をもとに牛群評価処理を行なう。変更された算定条件をもとに得られた算定結果は、再び端末１６に送信され、図１０（ｃ）のように端末１６の表示画面１６０に表示される。
【００９１】
利用者が「分娩間隔を試算」を選択入力すると、その入力は端末１６からシミュレーション装置１に送信され、受信したシミュレーション装置１は、初期値および目標値を入力させるための画面情報を端末１６に送信し、端末１６はこの情報を受信して図１０（ｅ）のように表示画面１６０に表示する。ここで、利用者が初期値および目標値を入力し、入力ボタンを選択すると、入力された初期値および目標値は端末１６からシミュレーション装置１に送信される。シミュレーション装置１は、受信した初期値および目標値をもとに、図８および図９の最適分娩間隔評価処理を行ない、ステップＳ２５において算定結果を端末１６に送信する。
【００９２】
端末１６は、図１０（ｆ）のようにして、受信した算定結果を表示画面１６０に表示する。算定結果には、最適分娩間隔の算定結果１６４のほか、シミュレーション期間における売上高の推移を表したグラフ１６１、シミュレーション期間における売上高や牛群の推移を表したテーブル１６２が表示される。利用者は表示画面１６０に表示された算定結果をもとに、牧場の経営評価を行なう。
【００９３】
なお、算定条件を変更し、再度シミュレーションを行なうことを利用者が希望する場合は、利用者は「条件を変更」ボタン１６５を選択入力することにより、算定条件の変更が可能である。ここで、変更可能な算定条件は初期値だけに限定されず、図１０（ｇ）のように、初妊牛や経産牛等の販売率を変更可能としてもよい。算定条件の変更に関する情報が端末１６に入力されると、シミュレーション装置１は変更された算定条件をもとに最適分娩間隔評価処理を行なう。変更された算定条件をもとに得られた算定結果は、図１０（ｆ）のように端末１６の表示画面１６０に表示される。
【００９４】
なお、上記の実施の形態に加えて、それぞれの分類の牛の購入または販売を適宜設定できる構成としても良い。どの分類の牛の購入または販売を重視するかは、牧場の経営方針に大きく左右される。牧場の経営方針は、頭数の規模拡大を重視する方針、頭数の規模を一定にし、若い牛を購入して老齢の牛を販売することにより牛の回転を早める方針、肉用牛として販売することを重視する方針など、様々な方針を取ることができる。したがって、代表的な経営方針に沿うような変化率テーブル４０をそれぞれ予め設定しておき、利用者に選択させるようにしても良い。さらに、利用者による変化率テーブル４０の設定値の一部変更を可能にしてもよい。
【００９５】
上記の実施の形態では、経営評価処理について、売上高を算定する構成としたが、売上高でなく利益を算定するように構成しても良い。この場合、各分類の牛の飼料費・購入費、牛舎などの施設の建設費・維持費、および労務費などの原価を設定した原価テーブルを設定しておき、これをもとに利益を算定する構成とすることができる。
【００９６】
上記の実施の形態では、牛群について説明したが、他の動物に応用することも可能である。他の動物に応用する場合は、変化率テーブルおよび単価テーブルの設定値を適切な値に変更すればよい。例えば妊娠期間が５か月程度の動物（山羊など）では、初妊牛に関する変化率、および分娩間隔を変更し、さらに一つの期を半年にするなどの設定変更により、本発明を応用することができる。さらに、一回の分娩で多数の子を産むことが通常である動物の場合（例えば、豚など）は、一回に分娩する子の数の平均値に分娩動物数を掛けた数に、変化率テーブルの出生率を掛けることで、出生数を算定する構成とすればよい。
【符号の説明】
【００９７】
１ シミュレーション装置
１１ ＣＰＵ
１２ ＲＡＭ
１３ ＨＤＤ
１４ 通信インターフェイス
１５ 通信ネットワーク
１６ 端末
４０ 変化率テーブル
６０ 牛群評価テーブル
８０ 単価テーブル
９０ 経営評価テーブル
【特許請求の範囲】
【請求項１】
複数の動物が存在する所定の動物群内における、分娩の経験がある経産動物の数、分娩の経験がなく妊娠もしていない非経産動物の数、および初めての妊娠をしている初妊動物の数についての所定の期間における変動を算定するシミュレーション装置であって、
前記所定の期間における期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに該動物群における分娩間隔を記憶する記憶手段と、
記憶手段により記憶された期首の非経産動物の数をもとに、前記所定の期間内に妊娠を経て初めての分娩を行なう第１の初産動物の数、前記所定の期間内に初妊動物となる期末の初妊動物の数、および前記所定の期間内に変化しない第１の非経産動物の数を算定する非経産動物変化数算定手段と、
記憶手段により記憶された期首の初妊動物の数をもとに、前記所定の期間内に初めての分娩を行なう第２の初産動物の数を算定する初妊動物変化数算定手段と、
記憶手段により記憶された期首の経産動物の数、非経産動物変化数算定手段により算定された第１の初産動物の数および初妊動物変化数算定手段により算定された第２の初産動物の数をもとに、期末の経産動物の数を算定する経産動物変化数算定手段と、
非経産動物変化数算定手段により算定された第１の初産動物の数、初妊動物変化数算定手段により算定された第２の初産動物の数、ならびに記憶手段により記憶された期首の経産動物の数および分娩間隔をもとに、前記所定の期間内に産まれうる出生動物の数を算定する出生動物数算定手段と、
非経産動物変化数算定手段により算定された第１の非経産動物の数、および出生動物数算定手段により算定された出生動物の数をもとに、前記所定の期間における期末の非経産動物の数を算定する期末非経産動物数算定手段と
を備えることを特徴とするシミュレーション装置。
【請求項２】
記憶手段が、前記所定の期間内に期首の非経産動物が初妊動物になる割合である初妊動物率、前記所定の期間内に非経産動物が妊娠を経て初産動物になる割合である第１の初産動物率、および前記所定の期間内に期首の初妊動物が初産動物になる割合である第２の初産動物率とをさらに記憶し、
非経産動物変化数算定手段が、記憶手段により記憶された初妊動物率および期首の非経産動物の数をもとに期末の初妊動物の数を算定し、かつ記憶手段により記憶された第１の初産動物率および期首の非経産動物の数をもとに第１の初産動物の数を算定し、
初妊動物変化数算定手段が、記憶手段により記憶された第２の初産動物率および期首の初妊動物の数をもとに、第２の初産動物の数を算定する
ことを特徴とする請求項１記載のシミュレーション装置。
【請求項３】
記憶手段が、経産動物のうち搾乳が可能な動物の割合である搾乳動物率をさらに記憶し、
シミュレーション装置が、さらに、
記憶手段により記憶された期首の経産動物の数、経産動物変化数算定手段により算定された期末の経産動物の数および搾乳動物率をもとに、搾乳が可能な搾乳動物の数を算定する搾乳動物数算定手段と、
搾乳動物数算定手段により算定された、搾乳動物の数をもとに、前記所定の期間内において出荷が可能な乳量を算定する乳量算定手段と
を備えることを特徴とする請求項１または２記載のシミュレーション装置。
【請求項４】
シミュレーション装置が、さらに、
前記所定の期間の整数倍であるシミュレーション期間についての利用者からの入力を受け付ける入力受付手段と、
記憶手段に記憶された期首の経産動物の数、期首の非経産動物の数、および期首の初妊動物の数を、それぞれ、経産動物変化数算定手段により算定された期末の経産動物の数、期末非経産動物数算定手段により算定された期末の非経産動物の数、および非経産動物変化数算定手段により算定された期末の初妊動物の数で更新する更新手段とを備え、
更新された期首の経産動物の数、期首の非経産動物の数、および期首の初妊動物の数をもととした、非経産動物変化数算定手段、初妊動物変化数算定手段、経産動物変化数算定手段、出生動物数算定手段および期末非経産動物数算定手段による算定、および更新手段による更新を、前記整数に対応する回数繰り返すことにより、シミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、および期末の初妊動物の数を算定する
ことを特徴とする請求項１、２または３記載のシミュレーション装置。
【請求項５】
前記分娩間隔が、仮に設定された仮分娩間隔であって、
シミュレーション装置が、さらに、
利用者の入力に従って、シミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、または期末の初妊動物の数についての目標数を受け付ける目標数受付手段と、
仮分娩間隔をもとに、非経産動物変化数算定手段、初妊動物変化数算定手段、経産動物変化数算定手段、出生動物数算定手段および期末非経産動物数算定手段による算定、ならびに更新手段による更新を、前記整数に対応する回数繰り返すことにより、シミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、および期末の初妊動物の数を算定する仮算定手段と、
仮算定手段により算定されたシミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、または期末の初妊動物の数が、目標数受付手段により受け付けられた目標数に達しているかについて判定する目標達成可否判定手段と、
目標達成可否判定手段により目標数に達していないと判定された場合は、より短く設定された仮分娩間隔を新たに設定し、目標数に達していると判定された場合は、より長く設定された仮分娩間隔を新たに設定する仮分娩間隔設定手段と、
仮分娩間隔設定手段により設定された新たな仮分娩間隔をもとに、仮算定手段によりシミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、および期末の初妊動物の数を算定し、仮分娩間隔が目標達成可否判定手段により目標数に達していると判定される最長の分娩間隔に達するまで仮分娩間隔設定手段による設定、仮算定手段による算定、目標達成可否判定手段による判定を繰り返し、目標数を達成するための最適な分娩間隔を特定する最適分娩間隔特定手段と
を備えることを特徴とする請求項４記載のシミュレーション装置。
【請求項６】
複数の動物が存在する所定の動物群内における、分娩の経験がある経産動物の数、分娩の経験がなく妊娠もしていない非経産動物の数、および初めての妊娠をしている初妊動物の数についての所定の期間における変動をコンピュータ装置に算定させるシミュレーション方法であって、
前記所定の期間における期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに該動物群における経産動物の分娩間隔を記憶領域に記憶させるステップと、
記憶領域に記憶された期首の非経産動物の数をもとに、前記所定の期間内に妊娠を経て初めての分娩を行なう第１の初産動物の数、前記所定の期間内に初妊動物となる期末の初妊動物の数、および前記所定の期間内に変化しない第１の非経産動物の数を算定させるステップと、
記憶領域に記憶された期首の初妊動物の数をもとに、前記所定の期間内に初めての分娩を行なう第２の初産動物の数を算定させるステップと、
記憶領域に記憶された期首の経産動物の数、算定された第１の初産動物の数および算定された第２の初産動物の数をもとに、期末の経産動物の数を算定させるステップと、
算定された第１の初産動物の数、算定された第２の初産動物の数、ならびに記憶領域に記憶された期首の経産動物の数および分娩間隔をもとに、前記所定の期間内に産まれうる出生動物の数を算定させるステップと、
算定された第１の非経産動物の数、および算定された出生動物の数をもとに、前記所定の期間における期末の非経産動物の数を算定させるステップと
から構成されることを特徴とするシミュレーション方法。
【請求項７】
複数の動物が存在する所定の動物群内における、分娩の経験がある経産動物の数、分娩の経験がなく妊娠もしていない非経産動物の数、および初めての妊娠をしている初妊動物の数についての所定の期間における変動を算定するサーバ装置と、通信ネットワークを介して前記サーバ装置と通信が可能な端末装置とを備えるシミュレーションシステムであって、
端末装置が、
利用者の入力に従って、前記所定の期間における期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに該動物群における分娩間隔の入力を受け付ける入力受付手段と、
入力受付手段により受け付けられた期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに分娩間隔をサーバ装置に送信する情報送信手段とを備え、
サーバ装置が、
端末装置から送信された期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに分娩間隔を受信する情報受信手段と、
情報受信手段により受信された期首の非経産動物の数をもとに、前記所定の期間内に妊娠を経て初めての分娩を行なう第１の初産動物の数、前記所定の期間内に初妊動物となる期末の初妊動物の数、および前記所定の期間内に変化しない第１の非経産動物の数を算定する非経産動物変化数算定手段と、
情報受信手段により受信された期首の初妊動物の数をもとに、前記所定の期間内に初めての分娩を行なう第２の初産動物の数を算定する初妊動物変化数算定手段と、
情報受信手段により受信された期首の経産動物の数、非経産動物変化数算定手段により算定された第１の初産動物の数および初妊動物変化数算定手段により算定された第２の初産動物の数をもとに、期末の経産動物の数を算定する経産動物変化数算定手段と、
非経産動物変化数算定手段により算定された第１の初産動物の数、初妊動物変化数算定手段により算定された第２の初産動物の数、ならびに情報受信手段により受信された期首の経産動物の数および分娩間隔をもとに、前記所定の期間内に産まれうる出生動物の数を算定する出生動物数算定手段と、
非経産動物変化数算定手段により算定された第１の非経産動物の数、および出生動物数算定手段により算定された出生動物の数をもとに、前記所定の期間における期末の非経産動物の数を算定する期末非経産動物数算定手段と
を備えることを特徴とするシミュレーションシステム。
【請求項１】
複数の動物が存在する所定の動物群内における、分娩の経験がある経産動物の数、分娩の経験がなく妊娠もしていない非経産動物の数、および初めての妊娠をしている初妊動物の数についての所定の期間における変動を算定するシミュレーション装置であって、
前記所定の期間における期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに該動物群における分娩間隔を記憶する記憶手段と、
記憶手段により記憶された期首の非経産動物の数をもとに、前記所定の期間内に妊娠を経て初めての分娩を行なう第１の初産動物の数、前記所定の期間内に初妊動物となる期末の初妊動物の数、および前記所定の期間内に変化しない第１の非経産動物の数を算定する非経産動物変化数算定手段と、
記憶手段により記憶された期首の初妊動物の数をもとに、前記所定の期間内に初めての分娩を行なう第２の初産動物の数を算定する初妊動物変化数算定手段と、
記憶手段により記憶された期首の経産動物の数、非経産動物変化数算定手段により算定された第１の初産動物の数および初妊動物変化数算定手段により算定された第２の初産動物の数をもとに、期末の経産動物の数を算定する経産動物変化数算定手段と、
非経産動物変化数算定手段により算定された第１の初産動物の数、初妊動物変化数算定手段により算定された第２の初産動物の数、ならびに記憶手段により記憶された期首の経産動物の数および分娩間隔をもとに、前記所定の期間内に産まれうる出生動物の数を算定する出生動物数算定手段と、
非経産動物変化数算定手段により算定された第１の非経産動物の数、および出生動物数算定手段により算定された出生動物の数をもとに、前記所定の期間における期末の非経産動物の数を算定する期末非経産動物数算定手段と
を備えることを特徴とするシミュレーション装置。
【請求項２】
記憶手段が、前記所定の期間内に期首の非経産動物が初妊動物になる割合である初妊動物率、前記所定の期間内に非経産動物が妊娠を経て初産動物になる割合である第１の初産動物率、および前記所定の期間内に期首の初妊動物が初産動物になる割合である第２の初産動物率とをさらに記憶し、
非経産動物変化数算定手段が、記憶手段により記憶された初妊動物率および期首の非経産動物の数をもとに期末の初妊動物の数を算定し、かつ記憶手段により記憶された第１の初産動物率および期首の非経産動物の数をもとに第１の初産動物の数を算定し、
初妊動物変化数算定手段が、記憶手段により記憶された第２の初産動物率および期首の初妊動物の数をもとに、第２の初産動物の数を算定する
ことを特徴とする請求項１記載のシミュレーション装置。
【請求項３】
記憶手段が、経産動物のうち搾乳が可能な動物の割合である搾乳動物率をさらに記憶し、
シミュレーション装置が、さらに、
記憶手段により記憶された期首の経産動物の数、経産動物変化数算定手段により算定された期末の経産動物の数および搾乳動物率をもとに、搾乳が可能な搾乳動物の数を算定する搾乳動物数算定手段と、
搾乳動物数算定手段により算定された、搾乳動物の数をもとに、前記所定の期間内において出荷が可能な乳量を算定する乳量算定手段と
を備えることを特徴とする請求項１または２記載のシミュレーション装置。
【請求項４】
シミュレーション装置が、さらに、
前記所定の期間の整数倍であるシミュレーション期間についての利用者からの入力を受け付ける入力受付手段と、
記憶手段に記憶された期首の経産動物の数、期首の非経産動物の数、および期首の初妊動物の数を、それぞれ、経産動物変化数算定手段により算定された期末の経産動物の数、期末非経産動物数算定手段により算定された期末の非経産動物の数、および非経産動物変化数算定手段により算定された期末の初妊動物の数で更新する更新手段とを備え、
更新された期首の経産動物の数、期首の非経産動物の数、および期首の初妊動物の数をもととした、非経産動物変化数算定手段、初妊動物変化数算定手段、経産動物変化数算定手段、出生動物数算定手段および期末非経産動物数算定手段による算定、および更新手段による更新を、前記整数に対応する回数繰り返すことにより、シミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、および期末の初妊動物の数を算定する
ことを特徴とする請求項１、２または３記載のシミュレーション装置。
【請求項５】
前記分娩間隔が、仮に設定された仮分娩間隔であって、
シミュレーション装置が、さらに、
利用者の入力に従って、シミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、または期末の初妊動物の数についての目標数を受け付ける目標数受付手段と、
仮分娩間隔をもとに、非経産動物変化数算定手段、初妊動物変化数算定手段、経産動物変化数算定手段、出生動物数算定手段および期末非経産動物数算定手段による算定、ならびに更新手段による更新を、前記整数に対応する回数繰り返すことにより、シミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、および期末の初妊動物の数を算定する仮算定手段と、
仮算定手段により算定されたシミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、または期末の初妊動物の数が、目標数受付手段により受け付けられた目標数に達しているかについて判定する目標達成可否判定手段と、
目標達成可否判定手段により目標数に達していないと判定された場合は、より短く設定された仮分娩間隔を新たに設定し、目標数に達していると判定された場合は、より長く設定された仮分娩間隔を新たに設定する仮分娩間隔設定手段と、
仮分娩間隔設定手段により設定された新たな仮分娩間隔をもとに、仮算定手段によりシミュレーション期間後における期末の経産動物の数、期末の非経産動物の数、および期末の初妊動物の数を算定し、仮分娩間隔が目標達成可否判定手段により目標数に達していると判定される最長の分娩間隔に達するまで仮分娩間隔設定手段による設定、仮算定手段による算定、目標達成可否判定手段による判定を繰り返し、目標数を達成するための最適な分娩間隔を特定する最適分娩間隔特定手段と
を備えることを特徴とする請求項４記載のシミュレーション装置。
【請求項６】
複数の動物が存在する所定の動物群内における、分娩の経験がある経産動物の数、分娩の経験がなく妊娠もしていない非経産動物の数、および初めての妊娠をしている初妊動物の数についての所定の期間における変動をコンピュータ装置に算定させるシミュレーション方法であって、
前記所定の期間における期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに該動物群における経産動物の分娩間隔を記憶領域に記憶させるステップと、
記憶領域に記憶された期首の非経産動物の数をもとに、前記所定の期間内に妊娠を経て初めての分娩を行なう第１の初産動物の数、前記所定の期間内に初妊動物となる期末の初妊動物の数、および前記所定の期間内に変化しない第１の非経産動物の数を算定させるステップと、
記憶領域に記憶された期首の初妊動物の数をもとに、前記所定の期間内に初めての分娩を行なう第２の初産動物の数を算定させるステップと、
記憶領域に記憶された期首の経産動物の数、算定された第１の初産動物の数および算定された第２の初産動物の数をもとに、期末の経産動物の数を算定させるステップと、
算定された第１の初産動物の数、算定された第２の初産動物の数、ならびに記憶領域に記憶された期首の経産動物の数および分娩間隔をもとに、前記所定の期間内に産まれうる出生動物の数を算定させるステップと、
算定された第１の非経産動物の数、および算定された出生動物の数をもとに、前記所定の期間における期末の非経産動物の数を算定させるステップと
から構成されることを特徴とするシミュレーション方法。
【請求項７】
複数の動物が存在する所定の動物群内における、分娩の経験がある経産動物の数、分娩の経験がなく妊娠もしていない非経産動物の数、および初めての妊娠をしている初妊動物の数についての所定の期間における変動を算定するサーバ装置と、通信ネットワークを介して前記サーバ装置と通信が可能な端末装置とを備えるシミュレーションシステムであって、
端末装置が、
利用者の入力に従って、前記所定の期間における期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに該動物群における分娩間隔の入力を受け付ける入力受付手段と、
入力受付手段により受け付けられた期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに分娩間隔をサーバ装置に送信する情報送信手段とを備え、
サーバ装置が、
端末装置から送信された期首の経産動物の数、期首の非経産動物の数および期首の初妊動物の数、ならびに分娩間隔を受信する情報受信手段と、
情報受信手段により受信された期首の非経産動物の数をもとに、前記所定の期間内に妊娠を経て初めての分娩を行なう第１の初産動物の数、前記所定の期間内に初妊動物となる期末の初妊動物の数、および前記所定の期間内に変化しない第１の非経産動物の数を算定する非経産動物変化数算定手段と、
情報受信手段により受信された期首の初妊動物の数をもとに、前記所定の期間内に初めての分娩を行なう第２の初産動物の数を算定する初妊動物変化数算定手段と、
情報受信手段により受信された期首の経産動物の数、非経産動物変化数算定手段により算定された第１の初産動物の数および初妊動物変化数算定手段により算定された第２の初産動物の数をもとに、期末の経産動物の数を算定する経産動物変化数算定手段と、
非経産動物変化数算定手段により算定された第１の初産動物の数、初妊動物変化数算定手段により算定された第２の初産動物の数、ならびに情報受信手段により受信された期首の経産動物の数および分娩間隔をもとに、前記所定の期間内に産まれうる出生動物の数を算定する出生動物数算定手段と、
非経産動物変化数算定手段により算定された第１の非経産動物の数、および出生動物数算定手段により算定された出生動物の数をもとに、前記所定の期間における期末の非経産動物の数を算定する期末非経産動物数算定手段と
を備えることを特徴とするシミュレーションシステム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【公開番号】特開２０１１−８１５７２（Ｐ２０１１−８１５７２Ａ）
【公開日】平成２３年４月２１日（２０１１．４．２１）
【国際特許分類】
【出願番号】特願２００９−２３２８５４（Ｐ２００９−２３２８５４）
【出願日】平成２１年１０月６日（２００９．１０．６）
【出願人】（５９２１７２５７４）明治飼糧株式会社 (15)
【公開日】平成２３年４月２１日（２０１１．４．２１）
【国際特許分類】
【出願日】平成２１年１０月６日（２００９．１０．６）
【出願人】（５９２１７２５７４）明治飼糧株式会社 (15)
[ Back to top ]