• Search Research Projects
  • Search Researchers
  • How to Use
  1. Back to project page

2016 Fiscal Year Research-status Report

逆問題解析による分数階微分方程式を用いた血糖値変化の数理モデル

Research Project

Project/Area Number 16K13774
Research InstitutionKyoto University

Principal Investigator

磯 祐介  京都大学, 情報学研究科, 教授 (70203065)

Co-Investigator(Kenkyū-buntansha) 藤原 宏志  京都大学, 情報学研究科, 准教授 (00362583)
今井 仁司  同志社大学, 理工学部, 教授 (80203298)
Project Period (FY) 2016-04-01 – 2019-03-31
Keywords数値解析 / 数値計算 / 数理モデル化 / 分数階微分方程式 / 逆問題
Outline of Annual Research Achievements

本研究は血中のブドウ糖濃度とインシュリンの関係について、insulin sesitivityという新たなパラメータを導入して両者の関係の数理モデル化を図るもので、分数階微分方程式を利用した研究協力者の先行研究をヒントに、信頼できる数理モデルの確立を図るとともに、分数階微分方程式の数値計算法の新たな数値計算法を提案することが研究目的である。研究目的と関連研究の現況を考慮し、初年度は非線型分数階微分方程式の数理解析を整理し、信頼性のある数値計算法を確立することを主眼としたが、この点では大きな成果が得られた。

分数階微分の定義は歴史的に見ても多様であるが、現象の数理モデル化への適用を念頭においた場合、Capto 型の分数階微分を微分方程式として利用することが合理的と考えている。この場合、目的とする分数階微分方程式は本質的には weak-singular 核をもつ非線型積分方程式と見做すことができるため、本課題研究の目的は weakly-singular 核を持つ Volterra 型非線型積分方程式の解析である。この種の方程式の数値計算については Diethelm 他による方法が提案され、一定の解析も行われているが、記法が煩雑な上に証明上の技術的な仮定があり、それが実用上の制約ともなっている。これに対して研究協力者の坂東孝政と共に Diethelm らの煩雑な議論を整理し、さらに証明上の技術的制約を本質的に緩和し、現実的な数値計算スキームを提案したことが平成28年度の成果である。この研究成果は非線型分数階微分方程式の Cauchy 問題の高精度解法の構成を示唆するもで、平成29年度以降の課題研究推進基盤ともなっている。

Current Status of Research Progress
Current Status of Research Progress

2: Research has progressed on the whole more than it was originally planned.

Reason

申請段階では非線形分数階方程式の初期値問題の数値解法についての Diethelm 他の先行研究の整理に相当の時間を要すると考え、本課題研究の最低限の到達点として新たな高精度数値計算法の確立を位置づけていた。しかし研究協力者の相当の協力により、数値計算法の確立については平成28年度に目処がついたものと判断している。これによって申請時に最低限の到達点としていた成果は早々に達成される見通しとなった。さらに研究分担者2名の共同研究によってスペクトル法を利用したアプローチも行い、限られた事例ではあるが、新たな数値計算法も提案している。これらの成果によって、本課題研究は、数値解析面では予定を若干上回りつつも、全体としてはおおむね順調に推移していると判断した。

Strategy for Future Research Activity

平成28年度の数値計算上の成果をふまえ、平成29年度には血糖値数理モデルの先行研究の検証と、その数値解析の信頼性を検証する。この数理モデルはパラメータ推定の逆問題を含むため、先ずは既知パラメータの場合の数値解析を平成29年度に行い、逆問題解析は平成30年度を中心に行うことを今後の方策としている。

Causes of Carryover

次の三つの事情により、初年度予定経費の一部を次年度使用とした:

(1)課題研究遂行上必要な数値計算環境の整備が、他の研究経費で実施できたこと(2)招聘予定の外国の研究協力者が別経費で来日したが、その際に研究連絡を実施することで本科研費の支出に及ばなかったこと。(3)体調不良により、研究協力者が予定していた出張をキャンセルしたこと。

Expenditure Plan for Carryover Budget

上記(2)の外国人研究協力者は、研究連絡と共同研究を深化させるために平成29年6月に招聘予定で、その際に繰り越した経費を使用する予定である。また前者(1)によって繰り越した経費に(3)を含めて本年度計画に加算し、当初計画を上回る計算機環境の整備を行い、課題研究の一層の推進を図る計画である。

  • Research Products

    (8 results)

All 2017 2016

All Journal Article (1 results) (of which Int'l Joint Research: 1 results,  Peer Reviewed: 1 results,  Open Access: 1 results) Presentation (7 results) (of which Int'l Joint Research: 7 results)

  • [Journal Article] 定常輻射輸送方程式の速度方向区分定数近似の数値解析2016

    • Author(s)
      藤原 宏志
    • Journal Title

      計算数理工学論文集

      Volume: 16 Pages: 91-96

    • Peer Reviewed / Open Access / Int'l Joint Research
  • [Presentation] Diffused Optical Tomography as an Inverse Problem for the Transport Equation2017

    • Author(s)
      磯 祐介
    • Organizer
      Korea-Japan Joint Workshop on Numerical Analysis and Scientific Computation
    • Place of Presentation
      Jeju National University Unipark, Seoguipo, Jeju, (大韓民国)
    • Year and Date
      2017-02-21 – 2017-02-24
    • Int'l Joint Research
  • [Presentation] Semi-discrete numerical scheme of the 2D stationary radiative transport equation on highly-parallel architectures2017

    • Author(s)
      藤原 宏志
    • Organizer
      Korea-Japan Joint Workshop on Numerical Analysis and Scientific Computation
    • Place of Presentation
      Jeju National University Unipark, Seoguipo, Jeju,(大韓民国)
    • Year and Date
      2017-02-21 – 2017-02-24
    • Int'l Joint Research
  • [Presentation] Numerical limit and its application to a blow-up problem governed by an ODE system2017

    • Author(s)
      今井 仁司
    • Organizer
      Korea-Japan Joint Workshop on Numerical Analysis and Scientific Computation
    • Place of Presentation
      Jeju National University Unipark, Seoguipo, Jeju, (大韓民国)
    • Year and Date
      2017-02-21 – 2017-02-24
    • Int'l Joint Research
  • [Presentation] Semi-discrete upwind finite volume methods for the 2D stationary radiative transport equation by GPGPU2016

    • Author(s)
      藤原 宏志
    • Organizer
      2016 Japan-Taiwan Joint Workshop on Numerical Analysis and Scientif ic Computation
    • Place of Presentation
      Lecture Room 202, Astro-Math Building, National Taiwan University(中華民国)
    • Year and Date
      2016-11-26 – 2016-11-28
    • Int'l Joint Research
  • [Presentation] Numerical limit and its application to a blow-up problem2016

    • Author(s)
      今井 仁司
    • Organizer
      2016 Japan-Taiwan Joint Workshop on Numerical Analysis and Scientif ic Computation
    • Place of Presentation
      Lecture Room 202, Astro-Math Building, National Taiwan University(中華民国)
    • Year and Date
      2016-11-26 – 2016-11-28
    • Int'l Joint Research
  • [Presentation] Direct computation of the radiative transport equation for near-infrared light propagation in biological tissue2016

    • Author(s)
      藤原 宏志
    • Organizer
      The 11th AIMS Conference on Dynamical Systems, Differential Equations and Applications
    • Place of Presentation
      Hyatt Regency Orland(アメリカ合衆国)
    • Year and Date
      2016-07-05 – 2016-07-05
    • Int'l Joint Research
  • [Presentation] Numerical Computation of the 2D Radiative Transport Equation by Upwind Finite Volume Method2016

    • Author(s)
      藤原 宏志
    • Organizer
      SIAM: East Asian Section Conference 2016
    • Place of Presentation
      University of Macau(中華人民共和国)
    • Year and Date
      2016-06-20 – 2016-06-20
    • Int'l Joint Research

URL: 

Published: 2018-01-16  

Information User Guide FAQ News Terms of Use Attribution of KAKENHI

Powered by NII kakenhi