入力中の回答があります。ページを離れますか?※ページを離れると、回答が消えてしまいますページを離れるA回答を続ける入力中のお礼があります。現在,一自由度系の振動問題において(1)式の微分方程式の定式化を行い最小二乗法を用いたM(質量),B(粘性),K(剛性)のパラメータ同定を行っています。F=M(dx^2/dy^2)+B(dx/dy)+Kx・・・(1)F:力(測定値),x:変位(測定値)dx/dy:速度(差分による計算値)dx^2/dy^2:加速度(差分による計算値)ここで,式誤差の影響を減らすために,一般化最小二乗法を用いたM,B,Kの算出を行おうと考えていますが,文献を調べても簡単に記載されているものが見つかりません。よろしければ,わかりやすいサイト,文献やC言語などのソースが手に入る場所などがあれば教えて頂けますでしょうか。通報するこの質問への回答は締め切られました。意味」に関するQ&A:古いゲーム「ラプラスの魔」のこのアイテムの意味。23:56実務でお困りなのだろうと思います。対象は減衰項を持つ振動子の系ですね。Fは強制的に与える外力でしょう。こんな間接的なやりかたで質量を精度よく測りたい。ばね定数や粘性も未知だという状況ですから、例えば、ナノサイズのカンチレバーにくっついたDNAの質量を測ろう、なんて話かも知れません。すると、測定値に大きな誤差が入るのはしょうがないことです。●測定誤差が気になるような場合には、数値微分はひじょーにアブナイ。(1)式右辺の3つの項の測定誤差が互いに強い相関を持つ訳ですから、理論上も無茶です。いうときは積分形を考える方がずっと良い。積分によってノイズが平滑化されるからです。=M∫x'dy+B∫xdy+K∫∫x(dy)^2-(Mx'(0)+Bx(0))∫dy=Mx(y)+B∫xdy+K∫∫x(dy)^2-(Mx'(0)+Bx(0))y-Mx(0)ですから、P=-(Mx'(0)+Bx(0))Q=-Mx(0)とおくと∫∫F(dy)^2=Mx(y)+B∫xdy+K∫∫x(dy)^2+Py+Qであり、線形モデルになっています。Q0)が得られるでしょう。測定誤差がごく小さい場合、あるいは測定回数が多い場合には、これで十分です。●測定回数が制限されていて測定誤差が大きいときには、もっとまじめに、測定誤差の独立性を壊さないように扱うべきです。そのためには微分方程式を陽に解いてモデルを作る必要があります。(1)式のような線形微分方程式は(たとえば演算子法やラプラス変換を使うと)すんなり解ける。y)*F(y)という形に書けます。ここに、"*"は畳み込み積分(convolution)です。function)」、そして、この系は「二次遅れ系」と呼ばれます。)さて、無視できない測定誤差がx(y)とF(y)の両方に含まれている場合には、話がかなり難しくなります。yの測定誤差まで考えるともっと大変なことに。ここではyとF(y)の測定値は信用できるものとしましょう。xの測定値をXと書くことにして、Xに入っている誤差εだけを考慮しようという訳です。y)*F(y)+ε(y)である。K)を探索すれば良い。(これをfittingと言います。K0)が利用できます。y)*F(y)が実際に「ホワイトノイズ」らしい性質を持つ事を検証する、ということです。(単にプロットしてみるだけでも意味があります。)この検証をやってみると、例えば「|ε(y)|が(x(y)+C)に比例して激しく変動してた」なんてこともしばしばあります。γ(y)=ε(y)/(x(y)+C)とすると、γ(y)こそが「ホワイトノイズ」である。γ(y)というモデルが出来ますから、これを使って非線形最小二乗法をやれば良い。その他いろんな場合があり得ますけど、誤差が「ホワイトノイズ」になるようにモデルを変形してやる訳です。白色化」ということでしょう。最小二乗法の教科書としては、いつも「中川,小柳:最小二乗法による実験データ解析,東京大学出版会」をお薦めしています。)の問題と捉えるのが適切でしょう。その観点からは、モデルが実際の系をきちんと表しているかどうかについて、常に留意すべきです。外力Fや変位xを測る手段自体が遅延などの応答特性を含んだ系である、ということがしばしばあるし、考慮していない要因に由来する系統誤差が入ることも多い。残差を調べてみて、モデルと測定結果が系統的にずれるようなら、モデルが実際の系に合っていない可能性があります。その場合、より正しいモデルを構築しなくちゃいけません。直すべきか心当たりがない時には、とりあえず系が線形であると仮定する。m(y)*F(y)と表せます。既知の入力Fに対する出力xを測定して、系の伝達関数mをまるごと推定することを考える。古典制御理論の逆問題(inversion。この場合はdeconvolution)であり、イチオウ数値的に解くことができます。(これも線形最小二乗法の一種の応用です。逆問題はノイズに極めて敏感なので、測定回数をうんと増やして情報量を増やす必要があります。)得られたmを物理的に解釈し直します。その他の要因に関するパラメータを含むモデルが構成できたとすると、xの測定値からこれらのパラメータを知る手段が得られるかも知れません。(得られないかも知れません。)系に非線形性がある場合にも、ある条件下(たとえば、xの振幅がある範囲内に入るという条件下)でなら線形で近似できることも多いので、条件を絞って解析する手があります。いずれにせよ、他人が作ったプログラムを拾って来てなんとかする、というような話ではありませんぜよ。=Mx(y)+B∫xdy+K∫∫x(dy)^2-(Mx'(0)+Bx(0))y-Mx(0)の式が算出できなくなります。tを変数として微分方程式を定積分していく方法に関して,参考になる文献などありましたら教えて頂ければと思います.次に,∫∫F(dy)^2=Mx(y)+B∫xdy+K∫∫x(dy)^2+Py+Qにおいて,測定値より∫∫F(dy)^2と∫∫x(dy)^2を計算する際はまず,∫Fdyと∫xdyつまり,0~y(時間)におけるF(力)とx(変位)との面積を定数として算出し,もう一度0~yの定積分を行う際は∫Fdyと∫xdyは定数となっているためにyを掛けるということになると考えます.もし考え方に間違いがあるようであればご指摘頂ければうれしいです。13:02ANo.4のコメントについてです。ANo.4の最後の式の左辺はえふ(わい)=...じゃなくて、正しくはげふ(わい)=...でげす。「えくす(わい)」の関係が独立ではなくなってしまうご指摘の通りです。ANo.3で申し上げた通り、真面目にやるのなら、「測定誤差の独立性を壊さないように扱うべきです。そのためには微分方程式を陽に解いてモデルを作る必要があります。」でもANo.3前半にある手抜きのやり方をしても、「測定誤差がごく小さい場合、あるいは測定回数が多い場合には、これで十分」である。「積分によってノイズが平滑化されるから」です。この手抜きが通用するのは、<測定データが、モデルと良く合っている場合>に限られる。fitting自体はANo.3の前半の手抜き法を使って行うにしても、微分方程式を陽に解いてモデルを作る、というところまではまじめにやって、fittingの結果得られたパラーメタをモデルに代入したものと、測定データとのずれ(残差)をプロットしてみるべきです。通報するこの回答へのお礼ご回答ありがとうございました。「微分方程式を陽に解いてモデルを作る」というところからはじめるが必要ということがわかりました。再検討してみたいと思います。05:13ANo.3のコメントについてです。yを変数としたまま定積分を行う方法は一般的なのかこれはごく普通の定積分です。積分変数のyと積分範囲のyは別物なんですが、stomachmanが両者を同じ記号で書いたから混乱なさったのでしょう。手抜きですいませんね。>これまでの私の理解では定積分を行う際の積分範囲は0~Y(Yは定数)として考えるものだと考えていました.その通りです。Yは定数だと思って計算し、それからYを変数yに置き換えれば良いんです。微分方程式を解くわけじゃなく、単に積分するだけなんですから、高校の教科書で充分でしょう。一度0~yの定積分を行う際は∫Fdyと∫xdyは定数となっているためにyを掛けるということになると考えます.違います。積分変数と積分範囲の変数を変えて書くなら、∫Fdy(積分範囲は0~Y)はYの関数です。これをえふ(Y)=∫Fdy(積分範囲は0~Y)と書くことにする(勿論、数値積分で、Yがyの各サンプル点の場合について計算します)。これをさらにYで積分すると、げふ(わい)=∫えふ(Y)dY(積分範囲は0~わい)これはわいの関数です(勿論、数値積分で、わいがyの各サンプル点の場合について計算します)。同様に、X(Y)=∫xdy(積分範囲は0~Y)えくす(わい)=∫X(Y)dy(積分範囲は0~わい)も計算して、モデルえふ(わい)=Mx(わい)+BX(わい)+Kえくす(わい)+Pわい+Q+誤差(わい)を作るんです。(わいがyの各サンプル点である場合について。yのサンプル点の個数だけ、式ができることになります。回答ありがとうございます。確かに式上では定積分の範囲yを積分変数yと同じ記号で書いた理由が理解できました。何度も申し訳ありませんが,離散化した場合について一つ確認事項があります。X(Y)=∫xdy(積分範囲は0~Y)えくす(わい)=∫X(Y)dy(積分範囲は0~わい)に関して,台形則の数値積分を行うとX(Y)={(x1+2*x2+x3)*dy/4}+{(x1+2*x2+x3)*dy/4}+・・・・・・えくす(わい)={(x1+2*x2+x3)*dy^2/4}+{(x2+2*x3+x4)*dy^2/4}+・・・・・・となります。「えくす(わい)」の関係が独立ではなくなってしまうと思われますが,問題ないでしょうか。06:08さいしょう二乗法全般として最適解の精度を上げるには、実験値の精度を上げること。ホワイトノイズを可能な限り除いてください。1桁精度を上げるのに最低で実験回数が10倍になります。nを無限大にすれば、積分されて誤差が減りますので。真空中で空気抵抗を除くとか、振動物と押さえているものの間の接触抵抗を減らしたりする(振り子のアルコール脱脂による手垢の除去等。手垢がついていると蒸発して変な力がかかることがある)ことは、やっていますね。計算精度を上げる、ざん差を減らす、というのであれば、発散しにくい式を使ってください。二乗法ではなく絶対値法をつかうとか、一次式ならばOR法を使うとか。計算機Eとオーバーフローとけたおちに注意してください。収束が悪くておかしいな、と調べたら(1週間かけて手計算してみたら)桁おちになっていた事が何階かあります。単精度と倍精度を比較してみる・倍精度と4倍制度を比較してみるなどを行ってみてください。極端な場合には、自分で多倍精度塩山ルーチンを用意する必要があります(昔のFM7で24bit整数演算に変換して計算したことがあります。整数演算でしたらばメモリー制限以外の制限無しの可変長のライブラリ(一万桁10進演算とか)を自作しています)。1000回の繰り返しがあるなんて場合には、2桁の測定精度でも、計算に必要な有効桁は6桁になります。単精度演算ですと5-6桁目がおかしなことになってくる数値になります。この回答への補足ありがとうございます。最小二乗法において解析精度をあげるためには実験方法,装置や解析法を全体的に再検討する必要がありそうですね。再検討してみたいと思います。17:01私の知識が乏しいのか、質問の意味が分かりません。何?最小二乗法と何が違う?「式誤差の影響」式には誤差が無いと思うが、あるのは測定値だが?最小二乗法なら簡単です。従属変数も独立変数もすべて測定可能で測定されている。この一次結合式を求めれば良いのではないですか。教科書にいくらでもあります。説明不足ですみません。F,dx^2/dy^2,dx/dyとxが既知で,最小二乗法でM,B,Kを算出することを同定としています。システム同定」ではかたよりと記載されていました。)を与えてしまうらしいです。そこで,一般化最小二乗法は残差を白色化し,パラメータがより精度良く求める手法のようです。通常用いられている最小二乗法はすでに試していますが,より精度よくパラメータを算出する方法を探しています。意味」に関するQ&A:世界1位のトヨタと世界7位のホンダの違いについて。円の方程式x^2+y^2+Ax+By+C=0において、最小二乗法でA,B,Cを求める式をあらわすとどうなりますか。次の最適化問題(最大化)の定式化が分かりません。微分方程式の質問です。連立微分方程式の問題について質問です。次の連立微分方程式を解け。偏微分方程式と波動方程式の一般解を解いてください。分からなくて困っています。複素解析より問題です。△=d二乗/dX二乗+d二乗/dY二乗とする。回答数ランキングです。ついて自分なりに調べてみたのですが、イマイチ良くわかりません。どなたか上手く説明していただけないでしょうか?過去の質問の回答に>cを列ベクトル、Aを行列とする。>(cの転置)Ac>0>となればAは正定値といいます。>Aの固有値が全て正であることとも同値です。行列は対称行列であることを仮定しています。対して、(cの転置)Ac>0となることである。対称行列Aが正定値なら、その固有値はすべて正です。(cとして固有ベクトルをとってみればよいでしょう。)逆に、対称行列Aの固有値がすべて正なら、Aは正定値行列です。対称行列ではないAの固有値がすべて正だからといって、(cの転置)Ac>0とは限りません。]とすると、Aは対称行列ではなく、固有値は1です。0となってしまいます。(実際に計算して確かめてください。)なので、行列Aが対称行列であるという条件はとても重要です。『任意の』と書き直したものです。このとき、半正定値行列の固有値はすべて0以上です。(つまり0も許します。)逆に、対称行列の固有値がすべて0以上なら、その行列は半正定値です。負定値の定義は、『ゼロではない任意の』ベクトルcに対して(cの転置)Ac<0となることです。固有値についてはもうわかりますね。行列は対称行列であることを仮定しています。対して、(cの転置)Ac>0となることである。対称行列Aが正定値なら、その固有値はすべて正です。(cとして固有ベクトルをとってみればよいでしょう。)逆に、対称行列Aの固有値がすべて正なら、Aは正定値行列です。Inc.Q質問する(無料)Q&A検索goo検索gooIDで新規登録・ログイン新規登録する(無料)会員の方はこちら外部サービスのアカウントで新規登録・ログインYahoo!で新規登録・ログインGoogleで新規登録・ログインTwitterで新規登録・ログインFacebookで新規登録・ログインdアカウントで新規登録・ログイン※各種外部サービスのアカウントをお持ちの方はこちらから簡単に登録できます。お客様の許可なしに外部サービスに投稿することはございませんのでご安心ください。会員でない方、会員になるとあなたも質問や回答ができるようになります!!親切・丁寧な回答が集まる場30分以内に回答がもらえ回答数は平均3件!匿名(ニックネーム)での投稿なので。