Моделирование задач радиофизики и электроники в системе Mathcad. Радченко Ю.С - 41 стр.

UptoLike

41
R
T
0.777 1.73328.627i+ 1.73328.627i 2.312()=
Rpolyrootsv():=v1478 1256 82651()
T
:=
x
4
5x
3
⋅+ 826x
2
⋅+ 1256x⋅− 1478 0
Определение корней полинома
newfdf, 0, 0.001,()0.673=
dfx()
x
fx()
d
d
:=
метод
Ньютона
newfdf, x, tol,()xxx
fx()
dfx()
−←
xxx
xxx
fx()
dfx()
−←
xxx tol>while
xx
:=
bisf0, 1, 0.0001,()0.673=
метод
половинного
деления
bisfa, b, tol,()
cab+()0.5⋅←
ac fc()fa() 0>if
bc fc()fb() 0>if
batol0.5+← otherwise
ba tol>while
c
:=
Программное решение уравнений
rootfx()x,()0.673=
x1:=
rootfx()x, 0, 1,()0.673=
00.51
2
0
2
fx()
x
Применение стандартных функций с заданием
интервала и начального приближения
fx()e
0.5x
1.4:=
Решение нелинейных уравнений
Рис . 4.1
                                                           41

                        Реш ен ие н ел ин ейн ы х урав н ен ий
              0.5 ⋅x
f ( x) := e            − 1.4
Применение с т анд арт ных ф унк ций с зад анием
инт ервала и началь ного приближения

  root ( f ( x) , x , 0 , 1) = 0.673                                          2

      x := 1                                                    f ( x)        0
      root ( f ( x) , x) = 0.673
                                                                             2
 Программное реш ение уравнений                                                   0         0.5           1
                                                                                             x
 bis ( f , a , b , tol) :=       while    b − a > tol
                                     c ← ( a + b) ⋅0.5                                  мет од
                                                                                        половинного
                                     a ← c if f ( c) ⋅ f ( a) > 0
                                                                                        д еления
                                     b ← c if f ( c) ⋅f ( b) > 0
                                     b ← a + tol⋅0.5 otherwise
                                 c
                                                                         bis ( f , 0 , 1 , 0.0001 ) = 0.673
                                              f ( x)
  new ( f , df , x , tol) :=         xx ← x −
                                              df ( x)
                                                                                           мет од
                                     while   x − xx > tol                                  Н ь ю т она
                                          x ← xx
                                                        f ( x)                                     d
                                          xx ← x −                                    df ( x) :=      f ( x)
                                                        df ( x)                                    dx
                                     xx
                                                                  new ( f , df , 0 , 0.001 ) = 0.673
       О пред еление к орней полинома
  4            3             2
 x + 5 ⋅x + 826 ⋅ x − 1256 ⋅x − 1478                     0
                                                   T
 v := ( −1478 −1256 826 −5 1 )                               R := polyroots ( v)
   T
 R = ( −0.777 1.733 + 28.627i 1.733 − 28.627i 2.312 )

                                                       Рис . 4.1