レベルが低い質問で申し訳ないのですが、
ある乗り物があったとします。
この乗り物の加速度が現在速度vをパラメータとして2次式で以下の式で表されるとします。
a=0.0002v^2-0.0445v+3.2951
初速度はv0=0として、
1.t秒後の速度
2.xメートル走行後の速度
3.xメートル走行に要する時間
の3つを求める式はどのように作ればよいのでしょうか。昔、積分を使って計算した気がするのですが、どうも忘れてしまって、計算できません。
パソコンで計算させるので、複雑な式でもかまいませんから、結果としての数値表ではなく、パラメータを用いた関数として教えていただければ、と思います。
宜しくお願いします。
No.4ベストアンサー
- 回答日時:
積分は解析的に求めるのでしょうか?
とりあえず解を求める方向性を書いておきます。
具体的な積分などは解析的に求めるか、プログラムで数値的に求めてください。
加速度aは速度vの時間微分ですから、
a=0.0002v^2-0.0445v+3.2951
より、
dv/dt=0.0002v^2-0.0445v+3.2951
よって初期条件も考慮すると、
∫1/(0.0002v^2-0.0445v+3.2951)dv(0~v)=t(v).....(1)
左辺の積分を実行すると、時間tが速度vの関数として書けるので、
逆解きして速度vが時間t時間の関数として書けます。
これが問題1の解ですね。
次に速度vは走行距離xの時間微分ですから、
dx/dt=v(t)
ここで右辺は問題1で求めました。
初期条件(x(t=0)=0)も考慮して、両辺時間積分すると
x(t)=∫v(t)dt(0~t).....(2)
よって右辺積分を実行すると、走行距離x(t)が求まります。
これを逆に解いて、走行距離xのときの時間t=t(x)が求まります。
これが問題3の解ですね。
さらに、問題1の解v(t)のtに今求めた解を代入すれば、
速度vが走行距離xの関数として書けます。
これが問題2の解ですね。
これでいかがでしょう。
わかりやすい説明ありがとうございました。
この計算で答えが求められることは理解できました。しかし、私の知識が乏しく、高校時代の教科書は家においたままであるため、(1)左辺の積分を解き切れません。(というより、情けないことに、分数の積分の解き方を忘れてしまっていて、周囲の人間とともに考えたのですが、解けないのです。)
とりあえず、帰宅後にでも教科書を引っ張り出したいと思います。
No.6
- 回答日時:
>加速度は今回の場合、m/s^2 で計算してあります。
速度はm/s です。この単位だと、先程の結果に質問の数値を入れて plot したところ
200秒付近で発散してしまうのですが、時間は200秒以内でよいのでしょうか?
あと、お節介かもしれませんが
今の数値だと、v^2 の係数からの誤差が効いていて
3.2951 のような桁数はいらないですね。
有り難うございました。
皆さんからの回答を頂いたあと、会社で色々考えていたのですが、どうやら近似式に若干問題があるようでした。
結論から言うと、200秒以内の話ですので問題はないのですが、近似を指数近似に変更した方が色々都合が良さそうなので、再度明日会社で近似をやり直そうと考えております。
それと誤差については仰るとおりですね。近似を Excel でやったもので、その結果の式をそのまま使っていました。以後気を付けます。
No.5
- 回答日時:
やり方は atsuota さんの仰るとおりです。
今の数値の場合は
a=0.0002,b=-0.0445,c=3.2951
√=√(4ac-b^2)
として、積分を実行すると
t(v) = 2/√*[ Arctan{(2av+b)/√} - Arctan{b/√} ]
となります。これを逆に解いて
v(t) = [√tan{√/2*t + Arctan(b/√)}-b]/2a
さらに、積分定数を x0 として積分すると
x(t) = -1/a*log| cos{√/2*t + Arctan(b/√)} | - b/2a*t + x0
となりましたが、これを逆に解いて t(x) はきびしいですね。
したがって、2と3はやっていません。
上に挙げた3つの式から後は数値的にやってみてください。
また、∫1/(av^2+bv+c)dv(0~v) の積分の際
2av+b = √(4ac-b^2)*tanθ
と変数変換しました。
一応結果を確認してください。
No.3
- 回答日時:
加速度は位置の二階時間微分、速度は位置の一階時間微分なので、この
a=Av^2+Bv+C
のかたちを(A、B、C、D、Eを定数として)、
a=A(v+D)^2+E
にして、v+Dの微分はaと同じことから、v+D=Vとして
dV/dt=AV^2+E
を解いて・・・・てのはだめですかね???
どこかの本に解がありそうですが・・・
すんません・・・見つかりません・・・^^;
No.2
- 回答日時:
加速度を速度の二次関数で表記することは可能です。
工学屋さんと物理屋さんだけですが。
これの場合は加速度そのものが変化するのですよね。
初期加速度=3.2951
1秒後は3.2951進んでいる
速度1の時は加速度3.2508
加速度の加速度を求めればいけますね。
あとは単位がないと数値が正しくとも0点です。
早速のお返事をありがとうございます。
また、言葉が至らず、失礼しました。
加速度は今回の場合、m/s^2 で計算してあります。速度はm/s です。
加速度の加速度ということは、単純に考えると、加速度が単位時間あたりにどれくらい変化しているか、といった数値ですよね?
ちょっと頭を整理してみます。
それと、さっき気がついたのですが、2次近似よりは指数近似の方が近似として正しいような気がしてきました。実際のところ、どちらの近似を使うのが良いのでしょうか。
それはともかく、時間を割いての回答をありがとうございます。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・人生のプチ美学を教えてください!!
- ・10秒目をつむったら…
- ・あなたの習慣について教えてください!!
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・【大喜利】【投稿~9/18】 おとぎ話『桃太郎』の知られざるエピソード
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
m/s、m/s^2の読み方
-
10kgのものを30cmの高さから誤...
-
加速度と角加速度の関係について
-
衝撃力(撃力)の単位について
-
振動の単位 dB→m/s2に換算で...
-
「等加速度運動」と「等加速度...
-
物理の問題なのですが、水平で...
-
高校物理です。 力学の加速度に...
-
Excel 速度 加速度計算について
-
車で走っている人間の感覚で、...
-
高校物理;円運動+単振動+慣性力...
-
この問題なんですけど、糸を切...
-
衝撃と「G」
-
ジャーク(躍度)を考慮した力学...
-
振動加速度
-
加速度と等加速度の違いがわか...
-
伝達関数と周波数伝達関数の違...
-
相対加速度の問題です
-
瞬間中心を使って速度は求めら...
-
正のフィードバックと負のフィ...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
加速度と角加速度の関係について
-
m/s、m/s^2の読み方
-
10kgのものを30cmの高さから誤...
-
2000年 京都大学 第一問 (長...
-
振動の単位 dB→m/s2に換算で...
-
「等加速度運動」と「等加速度...
-
高校物理です。 写真のような問...
-
衝撃力(撃力)の単位について
-
βが加速度のときこのように差が...
-
衝撃と「G」
-
Excel 速度 加速度計算について
-
この問題なんですけど、糸を切...
-
高校物理 解答に使える文字が分...
-
加速度が発生しているときは力...
-
剛体の滑車の問題です。 αとβを...
-
1/3オクターブ分析
-
物理についてです。 「円柱を斜...
-
ヨーヨーの物理: 糸を引く加速...
-
衝撃力の計算方法 (振り子玉)
-
振動加速度
おすすめ情報