optimazation toolboxを使って反応速度定数をどうやって求めたらいいか。
2 views (last 30 days)
Show older comments
添付したPDFファイルの1ページにあるモデルで4つの反応速度定数を求めたいです。
手順としては実験値とモデル解の差の二乗が最も小さくなるようにoptimazation toolboxの関数fminconを使って
反応速度定数を求めたいのですが、最後のページにあるような6個のエラーが出ています。
エラーの原因が分からず、調べても出てこないので手詰まりです。
実験はまだ行っていないので実験値xeは適当なものを入れています。
MATLABを使うの初めてなので知識は全然ありません。
どのエラーでもいいので教えていただけないでしょうか。
0 Comments
Accepted Answer
Kazuya
on 14 Nov 2018
Edited: Kazuya
on 16 Nov 2018
エラーメッセージは下から読んでいきます。一番上の階層(実行しているスクリプト)でのエラー情報から始まり、その内部での処理におけるエラー発生個所・・という感じに、エラーメッセージの出力順が上がるにしたがって深入りしていくイメージでよいかと思います。
コード・エラーメッセージの添付だけでも大変助かるのですが、今回の形だと実行できない&コピペもできません。。何か事情があるのかとは思いますが、そもそもわざわざファイルを開くという手間が閲覧者側に発生しますので、少し不親切かと思いました。エラーメッセージはそのまま質問文内にコピペ、実行できる形のファイルを添付がベストです。
それはともかく、根本原因は kondo_reaction_model にありそうです。ode15s で被積分関数として指定する場合
を見る限り、入力引数として2つ受け取る必要があるのかと。
function dtdx = kondo_reaction_model(x)
を
function dtdx = kondo_reaction_model(t,x)
に変えるだけで動くかも。検証はできていませんが。
2 Comments
More Answers (0)
See Also
Categories
Find more on Solver-Based Optimization Problem Setup in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!