optimazation toolboxを使っ​て反応速度定数をどう​やって求めたらいいか​。

添付したPDFファイルの１ページにあるモデルで４つの反応速度定数を求めたいです。

手順としては実験値とモデル解の差の二乗が最も小さくなるようにoptimazation toolboxの関数fminconを使って

反応速度定数を求めたいのですが、最後のページにあるような６個のエラーが出ています。

エラーの原因が分からず、調べても出てこないので手詰まりです。

実験はまだ行っていないので実験値xeは適当なものを入れています。

MATLABを使うの初めてなので知識は全然ありません。

どのエラーでもいいので教えていただけないでしょうか。

Answer 1

[Kazuya]

エラーメッセージは下から読んでいきます。一番上の階層（実行しているスクリプト）でのエラー情報から始まり、その内部での処理におけるエラー発生個所・・という感じに、エラーメッセージの出力順が上がるにしたがって深入りしていくイメージでよいかと思います。

コード・エラーメッセージの添付だけでも大変助かるのですが、今回の形だと実行できない＆コピペもできません。。何か事情があるのかとは思いますが、そもそもわざわざファイルを開くという手間が閲覧者側に発生しますので、少し不親切かと思いました。エラーメッセージはそのまま質問文内にコピペ、実行できる形のファイルを添付がベストです。

それはともかく、根本原因は kondo_reaction_model にありそうです。ode15s で被積分関数として指定する場合

https://jp.mathworks.com/help/matlab/ref/ode15s.html

を見る限り、入力引数として２つ受け取る必要があるのかと。

function dtdx = kondo_reaction_model(x)

を

function dtdx = kondo_reaction_model(t,x)

に変えるだけで動くかも。検証はできていませんが。

[Daisuke]

返信遅くなり申し訳ありません。

コードをアップロードできることを先ほど知りました。

次回からはアップロードしたいと思います。

親切に教えていただきありがとうございます。

Kazuyaさんのご指摘の通りにkondo_reaction_model(t,x)にしてみたところ

エラーが４つ消えました。

教えていただきありがとうございました。

[Kazuya]

いえいえすいません、余計なことを書きました。

他にも困ったエラーなどあったら書き込んでくださいねー。