Современное программное обеспечение в пользовательском процессе: Сборник заданий по курсу. Глушко А.В - 20 стр.

UptoLike

Составители: 

v501 = 4.07 - 0.22 *‰; v502 =-0.01 - 1.94 *‰*„
-v50*‰
;
f50 =„
H-3+3 L x
HH1.05 + 0.9 L + H3.05 - 1.9 L* Hx - v50LL;
Для проверки набираем
8v50, v501, v502, f50<
85.15, 4.07 - 0.22 Â, 1.74719 - 0.822127 Â,
H-3+3 ÂLx
HH1.05 + 0.9 ÂL+ H3.05 - 1.9 ÂLH-5.15 + xLL<
Введём обозначения
eq50 = 8D@#@xD,x,xD + H1 -‰L* D@#@xD,xD + H3 + 4 L* #@xD,
HD@#@xD,xD + 0.7 *‰*#@xDL
ê
.xÆ 0,
HD@#@xD,xD - 0.5 *‰*#@xDL
ê
.xÆ v50< &;
g50 = 88eq50@u[email protected]@1DD ä f50,
eq50@u[email protected]@2DD ä v501, eq50@u[email protected]@3DD ä v502<, 8x, 0, v50<< &;
Проведём проверку
g50@uD
88H3 + 4 ÂLu@xD+ H1 Lu
£
@xD+ u
££
@xD ==
H-3+3 ÂLx
HH1.05 + 0.9 ÂL+ H3.05 - 1.9 ÂLH-5.15 + xLL,
0.7 Â u@0D+ u
£
@0D == 4.07 - 0.22 Â,
-0.5 Â u@5.15D+ u
£
@5.15D == 1.74719 - 0.822127 Â<, 8x, 0, 5.15<<
u50 = u ê. First@NDSolve@Evaluate@g50@u[email protected]@1DDD,u,
Evaluate@g50@u[email protected]@2DDD, MaxSteps Æ 15000, PrecisionGoal Æ 17,
WorkingPrecision Æ 17, Method -> RungeKuttaDD
InterpolatingFunction@880., 5.1500000000000004<<, <>D
Одновременно программа делает замечание
NDSolve::"inaccur" :
"Now doing a transformation to the original
differential equation. The result may be
inaccurate at a low working precision."
Построим двумерные графики абсолютной величины, вещественной и
мнимой частей решения u50[x].
z
ad ok7bis.nb 20
zad ok7bis.nb                                                                20




    v501 = 4.07 - 0.22 * ‰; v502 = -0.01 - 1.94 * ‰ * „-v50*‰ ;
    f50 = „H-3+3 ‰L x HH1.05 + 0.9 ‰L + H3.05 - 1.9 ‰L * Hx - v50LL;

    Для проверки набираем

    8v50, v501, v502, f50<
    85.15, 4.07 - 0.22 Â, 1.74719 - 0.822127 Â,
      ‰H-3+3 ÂL x HH1.05 + 0.9 ÂL + H3.05 - 1.9 ÂL H-5.15 + xLL<

    Введём обозначения

    eq50 = [email protected]#@xD, x, xD + H1 - ‰L * [email protected]#@xD, xD + H3 + 4 ‰L * #@xD,
           [email protected]#@xD, xD + 0.7 * ‰ * #@xDL ê. x Æ 0,
           [email protected]#@xD, xD - 0.5 * ‰ * #@xDL ê. x Æ v50< &;
    g50 = [email protected]@@1DD ä f50,
             [email protected]@@2DD ä v501, [email protected]@@3DD ä v502<, 8x, 0, v50<< &;

    Проведём проверку

    [email protected]
    88H3 + 4 ÂL [email protected] + H1 - ÂL u£ @xD + u££ @xD ==
          ‰H-3+3 ÂL x HH1.05 + 0.9 ÂL + H3.05 - 1.9 ÂL H-5.15 + xLL,
        0.7 Â [email protected] + u£ @0D == 4.07 - 0.22 Â,
        -0.5 Â [email protected] + u£ @5.15D == 1.74719 - 0.822127 Â<, 8x, 0, 5.15<<

    u50 = u ê. [email protected]@[email protected]@[email protected]@1DDD, u,
           [email protected]@[email protected]@2DDD, MaxSteps Æ 15000, PrecisionGoal Æ 17,
           WorkingPrecision Æ 17, Method -> RungeKuttaDD
    [email protected], 5.1500000000000004<<, <>D

    Одновременно программа делает замечание
    NDSolve::"inaccur" :
     "Now doing a transformation to the original
       differential equation. The result may be
       inaccurate at a low working precision."

    Построим двумерные графики абсолютной величины, вещественной и
    мнимой частей решения u50[x].