パラメータの値を計算で出したい

質問

今、fortranを使って以下のようなプログラムを書いています。

1時間毎の温度（t）のデータがあり、それに対して指数関数でRが変化します。
R=a*exp(0.7*t)

これに対して、1ヶ月間のRの積算値を測定値として持っています。

この測定値に最も近くなるように定数aの値を求めたいのですが
fortranでどのようにプログラムしたらいいのかわかりません。

よろしくお願いします！

回答No.2

uen_sap

指数関数は扱い難いので、
これを

log(R)=log(a)+0.7*t
と変換します。

log(R)とtは一次の関数となります。
log(R)とtをデータとして最少二乗法で勾配と切片を求めます。
切片=log(a)

回答No.1

麻野 なぎ（@AsanoNagi）

a は定数なのですよね？

だとすると、

R0 = exp(0.7*t) を一ヶ月分積算して、

a = R / R0（の積算値）

で、OKだと思います。

お礼

どうもありがとうございます！
できました！
わたしがとてもバカでした。