• 締切済み

logの入った最小二乗法について

stomachmanの回答

  • stomachman
  • ベストアンサー率57% (1014/1775)
回答No.4

No.3の補足説明です。 lnY(lnA) = aln(x-b) - dln(x-c) において (2) ln(y) = lnY(lnA) とする代わりに、 (2') y=lnY(lnA) とした方が易しいですね。こうすれば微分方程式は (4') (x-b)(x-c)(dy/dx) =((x-c)a + (x-b)d) となり、 (5') (x^2)(dy/dx) -(b+c)x(dy/dx) + bc(dy/dx) = (a+d)x - (ac+bd) 積分すると  (x^2)y- 2∫xy dx -(b+c)(xy -∫y dx) + (bc)y = ((a+d)/2)(x^2)- (ac+bd)x + t より (6') (x^2)y- 2∫xy dx = (b+c)(xy -∫y dx) +((a+d)/2)(x^2)- (ac+bd)x - (bc)y + t  p=b+c  q=(a+d)/2  r=-(ac+bd)  s=-(bc) とおいて、  P[k] = x[k]y[k] -Σ(y[j]+y[j-1])(x[j]-x[j-1])/2  Q[k] = x[k]^2  R[k] = x[k]  S[k] = y[k]  T[k] = (x[k]^2)y[k] - Σ(y[j]x[j]+y[j-1]x[j-1])(x[j]-x[j-1]) として (7) T[k] = pP[k]+qQ[k]+rR[k]+sS[k]+t の方が多少すっきりしてますね。いや、やってることは同じですけど。

関連するQ&A

  • 最小二乗法について

     y=ax+b+c/x という式での最小二乗法の求め方を 教えてください。

  • 最小二乗法の応用について

    実験により、xに対するyの値をxの値をかえながら、N回測定した。測定したxに対するyの関係をグラフに描くと、次の二次関数で表現するのが適当であることがわかった。 y=ax*x+bx+c この時、最小二乗法によりパラメータa、b、cの値を求める式を導出せよ。という問題なのですが、どのようにしたら最小二乗法で求めることができるのですか? どうか教えてください

  • 3次元の最小二乗法

    A(x1,y1,z1)とB(x2,y2,z2)とC(x3,y3,z3)とD(x4,y4,z4)の点で 最小2乗法ほうを使い直線近似したいのですが、どのような式になりますかご存知の方教えていただけますか?

  • 最小2乗

    (1) (x,y)=(-3,4),(-2,1),(-1,0),(0,-1),(1,1),(2,2),(3,3)を 方程式 y = ax^2 + bx + c で最小2乗法に当てはめたときの a,b,cを求める (2) 要素(x,y,z)=(1,3,2),(2,1,1),(3,2,4),(4,5,3),(5,4,5)が z = ax + by + p とモデル化されるとき、残差の2乗Σp^2 が 最小になるように定式化して、最小となるa,bを求める という問題です。 これらを行列を用いて解きたいのですが、どのようにすればよいか分かりません。よろしくお願いします。

  • y=a/(x-b)+cの最小二乗法

    y=a/(x-b)+cの最小二乗法 y=a/(x-b)+c という、反比例の式をx方向に+b、y方向に+c平行移動したような曲線の係数a,b,cを求めるための最小二乗法の方法を教えていただけないでしょうか。 工夫してみたのですが、なかなかうまくいきませんでした。 すみませんが、力を貸してください。

  • 複雑な最小二乗法

    いま、最小二乗法で困っています。 ln{1/L*ln(1/(1-P))} = aln((X-c)/d)-bln((e-X)/f) という式で、Lは定数で、Pは確率(おそらく0~1なのかとおもいます)です。 aからfがパラメータなのですが、cとeはとりあえず定数としても かまわないそうです。 私には複雑すぎてどうしていいのかわかりません。 もしおわかりになる方がいらっしゃったら宜しくお願いいたします。

  • 最小二乗法

    n組のデータ (xi, yi) を,特定点(X0, Y0) を通る直線 y = ax+b でフィッティングしたい。最小二乗法で係数a,bを求めるため の式を導きなさい。 という問題で 各データの残差を二乗した和が最小になるときのa,bを求めるのですが 特定点(X0,Y0)を通るにはどうすればよいでしょうか? ただ単に、特定点を通らずフィッティングするやりかたはわかるのですが・・・。 よろしくお願いします。

  • 最小二乗法の推定値の誤差

    変数xを変化させたときの測定値yを最小二乗法で二次式y=a*x^2 + b*x +c にフィッティングさせ推定値a, b, cを求めるとき、 測定値yの誤差がδyであるときの推定値a, b, cの誤差を求めたいのです。 具体的には、(x,y)=(-1,2), (0,0), (1,1.5), (2,5) の4つのデータを 二次式にフィッティングさせたときのa,b,cはa=1.375, b=-0.325, c=0.225ですが、 測定値yの測定誤差が0.1のときのa,b,cの誤差を求めたいのです。 よろしくお願いします。

  • 最小二乗法。円の方程式x^2+y^2+Ax+By+C=0において、最小二乗法でA,B,Cを求める式をあらわすとどうなりますか。

    円の方程式x^2+y^2+Ax+Bx+C=0において、最小二乗法でA,B,Cを求める式をあらわす場合、どうなりますか。複雑な行列式であらわさないとだめなのですか。y=ax+bの場合の最小二乗法は何とかわかるのですが、未知数がA,B,Cの3つになると、わからなくなります。ご指導お願いいたします。

  • 最小二乗法

    円x^2+y^2=1上の点Pにおけるこの円の接線とx軸、y軸との交点をそれぞれA,Bとして、 Pが第一象限を動くとき、線分ABの長さLを求めよ。 なお、Lの最小値の存在理由も述べよ。 というもんだいなのですが、 線分ABはy=txとなるので、このtを最小二乗法でもとめるんだと思うのですが、 最小二乗法の理論値、実測値などが、まだ良くわかりません。 と言いますか、傾きが最小ならLは最小で合ってますでしょうか? 傾きが45°のときが怪しくなってきたのですが・・・ また、最小値の存在理由ってどうやったら説明できるのですか? よかったら、教えてください。