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