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

4 次查看(过去 30 天)
Daisuke
Daisuke 2018-11-14
评论: Kazuya 2018-11-19
添付したPDFファイルの1ページにあるモデルで4つの反応速度定数を求めたいです。
手順としては実験値とモデル解の差の二乗が最も小さくなるようにoptimazation toolboxの関数fminconを使って
反応速度定数を求めたいのですが、最後のページにあるような6個のエラーが出ています。
エラーの原因が分からず、調べても出てこないので手詰まりです。
実験はまだ行っていないので実験値xeは適当なものを入れています。
MATLABを使うの初めてなので知識は全然ありません。
どのエラーでもいいので教えていただけないでしょうか。

采纳的回答

Kazuya
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 个评论
Daisuke
Daisuke 2018-11-19
返信遅くなり申し訳ありません。
コードをアップロードできることを先ほど知りました。
次回からはアップロードしたいと思います。
親切に教えていただきありがとうございます。
Kazuyaさんのご指摘の通りにkondo_reaction_model(t,x)にしてみたところ
エラーが4つ消えました。
教えていただきありがとうございました。
Kazuya
Kazuya 2018-11-19
いえいえすいません、余計なことを書きました。
他にも困ったエラーなどあったら書き込んでくださいねー。

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile 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!