No.5171152 で 質問したものです。アドバイスを受けて計算してみると意外と
複雑になったので、確認させてください。
併せてNo.5171152で提示した計算式が少し違っていたので、訂正します。
k(b)=(1-b^2/32)^0.5
F(b)=m-1+0.25*k(b)^2+0.046875*k(b)^4+0.01953125*k(b)^6
上記式をbについて微分すると
d(k(b))/db =0.5 * (1-b^2/32)^(-0.5) * ((-2)/32) * b
d(F(b))/db =(0.25*2*k(b)+0.046875*4*k(b)^3+0.01953125*6*k(b)^5)
* d((k(b))/db
これでよろしいのでしょうか。申訳けないのですが、よろしくお願いします。
実務に利用するので、念を押したいのです。
(当方は還暦を過ぎているものです)
No.2ベストアンサー
- 回答日時:
#1です。
最後の式が「=」で終わっていますが、係数を小数点数に直した式を書こうとしたのですが、数値が割り切れず長大な桁数の係数の式になって、かえって扱いにくい式になるかと思い書かなかっただけです。
一応、「=」の後の式を書いておきますが、参考にするだけで、A#1の方の
d(F(b))/dbの式を使われた方が係数が整数のすっきりした式なので実用上は使いやすいのではと思います。
(参考)小数点の係数の式
d(F(b))/db≒-3.5762786865234375*10^(-6)*b^5
+ 4.119873046875*10^(-4)*b^3 - 0.025146484375*b
重ねてありがとうございます。この質問は下記の式から
Pとaを固定してb値を求める方法に関してです。
由来URLは下記より
http://home.att.net/~numericana/answer/geometry. …
ここの別ページには、楕円周長に関して関孝和の示した式が時代としては先駆的なものであったことが記されています。(比較:オイラーによる楕円関数)
There is no simple exact formula: There are simple formulas but they are not exact and there are exact formulas but they are not simple.
If the ellipse is of equation x2/a2 + y2/b2=1 with a>b, a is called the major radius, and b is the minor radius. The quantity e = Ö(1-b2/a2) is the eccentricity of the ellipse.
An exact expression for the ellipse perimeter P involves the sum of infinitely many terms of the form (-1)/(2n-1) [(2n)!/(2n n!)2]2 e2n. The first such term (for n=0) is equal to 1 whereas all the others are negative correction terms :
P/2pa = 1 - [1/4]e2 - [3/64]e4 - [5/256]e6 - [175/16384]e8 - [441/65536]e10 ...
Note that for a circle (e=0) of radius a, the above does give the circumference as 2p times the radius.
No.3
- 回答日時:
前回質問から言っていますが、ナゼ
dF/db = { dF/d(k^2) }{ d(k^2)/db }
dF/d(k^2) = 0.25 + (0.046875*2)*k^2 + (0.01953125*3)*(k^2)^2
d(k^2)/db = -(2/32)b
とは、したくないのでしょうか?
実務で使うのなら尚のこと、一目見て検算になるように
整理しておくのは大切なことです。
ありがとうございます。別にしたくはないとかでは、ないのです。
単純にプログラム言語にそのまま利用できると考えて記載しただけです。
ご指摘のようにするには、まだ、その域には私のレベルは達して
いないということですね。ありがとうございました。
No.1
- 回答日時:
>d(k(b))/db =0.5 * (1-b^2/32)^(-0.5) * ((-2)/32) * b
>d(F(b))/db =(0.25*2*k(b)+0.046875*4*k(b)^3+0.01953125*6*k(b)^5)
>* d((k(b))/db
この式で合っています。
後者の式にd(k(b))/db および k(b) の式を代入すると
d(F(b))/db=-(b*(15*b^4-1728*b^2+105472))/4194304
=
となります。
(数式処理ソフトのwxMaximaで計算していますので間違いないと思います)
重ねてありがとうございます。
前回の質問のアドバイスに本当に感謝です。級数の展開とか、微分積分とか、年配になってからやりだしたので、今後も何かありましたらよろしくお願いします。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 情報処理技術者・Microsoft認定資格 応用情報処理技術者試験のシステム利用率の計算について 2 2022/03/28 07:43
- SQL Server DBのテーブルの設計ができず困っています。 2 2023/06/29 16:43
- Excel(エクセル) アウトラインの小計のやり方 1 2023/03/20 11:51
- 公認会計士・税理士 地代家賃と前払費用の決算整理仕訳 1 2023/02/25 08:14
- 相続税・贈与税 相続税の、土地の計算法に関して、の質問です。 4 2022/07/05 23:12
- MySQL PHP 画像のアップロード Qiita 2 2022/11/28 04:44
- Visual Basic(VBA) ACCESS DAO で不要なテーブルのフィールド(列)の削除 4 2022/06/23 12:13
- 物理学 音 物理 1 2022/04/15 23:36
- その他(ブラウザ) python から COBOL を呼び出すことは可能でしょうか。 1 2023/04/21 21:03
- 工学 f(z)=tan(z) の 0<|z-π/2|<π での ローラン展開 f(z)=Σ_{n=-∞~∞ 4 2022/07/11 03:53
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
一軸圧縮強度quと粘着力Cの...
-
3元?連立方程式の解き方が分か...
-
二項定理
-
なぜ、因数分解?
-
うんこの形を数学で表すとどう...
-
未知数4つ、式4つの方程式の...
-
比例式の値を求める問題
-
1=√1=√(-1)(-1)=√(-1)√(-1)=i・...
-
0+1=∞
-
近似式の定理で、値 a が値 b ...
-
与式とは?
-
伸び率のマイナス数値からのパ...
-
1年くらいの女友達がいますが、...
-
下の問題の答えが無かったので...
-
xかけるxって答えなんですか?
-
?÷5/8=4の答え、解き方を教えて...
-
かけ算、割り算の移項
-
友達との前戯
-
割合の問題について 24Lが30%に...
-
経時変化とは、どのような意味...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
アルキメデス螺旋と対数螺旋の...
-
0+1=∞
-
i1,i2,i3の求め方
-
与式とは?
-
計算式で a(b,c)みたいな式って...
-
写真の式についてですが、2行目...
-
三次方程式x^3+3x^2+(a-4)x-a=0...
-
中学数学で1次式を選ぶ問題 分...
-
近似式の定理で、値 a が値 b ...
-
二項定理
-
1=√1=√(-1)(-1)=√(-1)√(-1)=i・...
-
連立方程式はなぜ解ける?
-
VBAで除算の商・・・
-
3元?連立方程式の解き方が分か...
-
数学の問題がわかりません。
-
自然数の列を次のような群に分...
-
一軸圧縮強度quと粘着力Cの...
-
x+y+z=0,2x^2+2y^2-z^2=0の...
-
連立方程式の解法に代入法、加...
-
中学数学です。
おすすめ情報