2019年9月5日木曜日

学習環境

解析入門(上) (松坂和夫 数学入門シリーズ 4) (松坂 和夫(著)、岩波書店)の第5章(各種の初等関数)、5.3(三角関数)、問題12の解答を求めてみる。


  1. P 1 , P 2 , X

    の座標をそれぞれ

    0 , a , b , c , x , 0

    とおく。

    このとき 経路を通る時間は、

    f x = x 2 + a 2 v 1 + b - x 2 + c 2 v 2

    微分。

    f ' x = x v 1 x 2 + a 2 - b - x v 2 b - x 2 + c 2

    問題のなす角の定め方により、

    sin θ 1 = x x 2 + a 2 sin θ 2 = b - x b - x 2 + c 2

    よって、

    f ' x = sin θ 1 v 1 - sin θ 2 v 2

    また、

    0 < θ 1 < π 2 0 < θ 2 < π 2

    なので、

    sin θ 1 v 1 = sin θ 2 v 2

    のとき

    f ' x = 0

    また、

    θ 1 , θ 2

    はそれぞれ単調増加、単調減少である。

    また、

    f ' 0 < 0 f ' b > 0

    よって、

    sin θ 1 v 1 = sin θ 2 v 2

    のとき最短時間で、 点 A はただ1つ存在する。

    (証明終)

コード

Python 3

#!/usr/bin/env python3
from sympy import pprint, symbols, sqrt, Derivative, solve, plot

print('12.')

x = symbols('x', real=True)
v1 = 1
v2 = 2
a = 3
b = 4
c = -2
f = sqrt(x ** 2 + a ** 2) / v1 + sqrt((b - x) ** 2 + c ** 2) / v2
g = Derivative(f, x, 1)
h = g.doit()
for o in [g, h]:
    pprint(o)
    print()

for x0 in solve(h):
    if 0 < x0 < b:
        for o in [x0, float(x0)]:
            pprint(o)
            print()

p = plot(f,
         (x, 0, b),
         show=False,
         legend=True)

colors = ['red', 'green', 'blue', 'brown', 'orange',
          'purple', 'pink', 'gray', 'skyblue', 'yellow']

for o, color in zip(p, colors):
    o.line_color = color

p.show()
p.save('sample12.png')

入出力結果(Bash、cmd.exe(コマンドプロンプト)、Terminal、Jupyter(IPython))

C:\Users\...>py sample12.py
12.
  ⎛                 ______________⎞
  ⎜   ________     ╱        2     ⎟
d ⎜  ╱  2        ╲╱  (4 - x)  + 4 ⎟
──⎜╲╱  x  + 9  + ─────────────────⎟
dx⎝                      2        ⎠

     x               x - 4       
─────────── + ───────────────────
   ________        ______________
  ╱  2            ╱        2     
╲╱  x  + 9    2⋅╲╱  (4 - x)  + 4 

                                                                              
                                                                              
                                                                              
                                                                              
          ________________________________________________________________    
         ╱                                           ____________________     
        ╱              3025               34        ╱ 4⋅√352999   539623      
       ╱   ──────────────────────────── + ── + 2⋅3 ╱  ───────── + ──────      
      ╱            ____________________   9      ╲╱       27       5832       
     ╱            ╱ 4⋅√352999   539623                                        
    ╱      162⋅3 ╱  ───────── + ──────                                        
  ╲╱           ╲╱       27       5832                                        ╲
- ──────────────────────────────────────────────────────────────────────── + ─
                                     2                                        

           ___________________________________________________________________
          ╱         ____________________                                      
         ╱         ╱ 4⋅√352999   539623                3025               68  
        ╱   - 2⋅3 ╱  ───────── + ──────  - ──────────────────────────── + ── +
       ╱        ╲╱       27       5832             ____________________   9   
      ╱                                           ╱ 4⋅√352999   539623        
     ╱                                     162⋅3 ╱  ───────── + ──────        
    ╱                                          ╲╱       27       5832         
   ╱                                                                          
  ╱                                                                           
 ╱                                                                            
╱                                                                             
──────────────────────────────────────────────────────────────────────────────
                                                                            2 

____________________________________________________________________________  
                                                                              
                                    200                                       
 ──────────────────────────────────────────────────────────────────────────   
           ________________________________________________________________   
          ╱                                           ____________________    
         ╱              3025               34        ╱ 4⋅√352999   539623     
 3⋅     ╱   ──────────────────────────── + ── + 2⋅3 ╱  ───────── + ──────     
       ╱            ____________________   9      ╲╱       27       5832      
      ╱            ╱ 4⋅√352999   539623                                       
     ╱      162⋅3 ╱  ───────── + ──────                                       
   ╲╱           ╲╱       27       5832                                        
──────────────────────────────────────────────────────────────────────────── +
                                                                              

  
  
  
  
  
  
  
  
  
  
  
  
 2
  

1.313462140723518


C:\Users\...>

0 コメント:

コメントを投稿

関連コンテンツ