Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Мараева И.Б. Вычислительная техника в инженерных и экономических расчетах. Автокод Инженер для ЭЦВМ Минск-22 учеб. пособие

.pdf
Скачиваний:
3
Добавлен:
23.10.2023
Размер:
4.54 Mб
Скачать

НО

иым О.

На печать АЦПУ вывести таблицу значений а , В f х , у и

соответствующих номеров квадрантов.

 

33 . Найти по изтоду наименьших квадратов

линейную функ­

цию

у - а х + В, аппроксимирующую заданную таблицу

значений

ж

и

у . Значения

х - ,

у;

( i =

/,

и число п& 5 0 -

известнн.

 

 

В вывести на печать

 

 

Найденные значения о

и

АЦПУ вме­

сте

с

их наименованиями и соответствующим заголовком.

 

Искомые значения

О и

В коэффициентов линейной функ­

ции m метолу наименьших квадратов

находятся

как

решение

системы

 

 

 

 

 

 

34 . Решить систему уравнений

^х гу г - З х 3 - 8у* + 8 -О

методом Ньютона с

точностью с

,

взяв за начальное прибли­

жение решения жо

= 1 , 2; £ =

0,

6.

Вывести на печать АЦПУ все

приближенные значения ж

и-конечный результат с соответствующими заголовками.

Формулы метода Ньютона для системы

F ( x ,y ) =0

G ( х , у ) =0

имеют вид:

)

где

I l l 55. Вычислить по схеме Горнера значения полинома степе­

ни

т

у = 0 , c c m ♦ ct£ х

 

. . . * а т с в * а т„

 

 

 

 

 

 

 

дая

п

значений

аргумента

ж

от с

через "h .

 

 

 

Числа с , h , n « 100 , т £

2 0 ,

а

также

коэффициенты

многочлена

d ;

( j

= I , . . .

,

т + i)

-

известны.

 

 

Полученные* значения

многочлена

^

 

напечатать на

АЦПУ

в виде

таблицы вместе с

соответствующими значениями аргумен­

т а .

Таблицу снабдить

заголовком.

 

 

 

 

 

 

Схецу Горнера определяют

формула

 

 

 

 

 

 

4

= с ,

>

 

 

 

 

 

 

 

 

 

 

 

В -

 

= В - х

+ а -

 

 

г

1>

 

- т ,

 

 

 

 

 

j+i

 

j

/** >

 

 

 

 

 

 

В

т ц

= и

 

 

 

 

 

 

 

 

 

 

 

 

 

у-

 

 

 

 

 

 

 

 

 

 

56.

функция

y =f(x))

задана таблицей п

значений

as.

и. Составить программу вычисления f ( х ) в точках

■x.€.faj]

с шагом л х

по формуле линейной

 

интерполяции

 

 

 

( х , 4 а < В { х п) .

 

 

 

Формула линейной интерполяции имеет вид:

 

 

 

f( x ) = y .+

_

 

 

если

^

 

о с < ж . ^ .

Полученные значения f

(х )

вместе с

соответствующими

значениями

х

вывести в

виде

таблицы на печать

АЦПУ.

37 . Непрерывная судовая линия состоит из

п

 

участков,

причем каждый участок

может быть задан уравнением вида

 

 

 

^ = о ж г + Вг с + с

 

 

 

(при

а =

о -

участок

прямолинейный, при a

f О

-

параболи­

ческий) .

 

 

 

 

 

 

 

 

Известны абсциссы точек стыка участков х х [а - ^

 

включая начальную и конечную точки, а также

коэффициенты а .,

4 ,

ci (

L =

I , . . .

, п

) уравнения на каждом

 

участке.

Число

п

10

- задано.

 

 

 

 

 

Вычислить длину дуги

судовой линии, причем для

прямоли­

нейного участка

 

 

 

 

 

 

 

* = И & . г у . / + ( ■ * » -

112

а для параболического -

 

 

 

 

 

 

£ = S

 

J ^ ^ v f d x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xr

 

'

 

 

 

 

 

 

 

 

 

 

Результатнапечатать на БШ .

 

 

 

 

 

 

 

 

36.

3

условиях предыдущей задачи вычислить шющадь

под

судовой линией,

причем для

прямолинейного участка

 

 

 

 

 

 

 

 

 

S =

 

 

 

 

 

С ® , * , - ® * ) j

 

 

 

 

 

 

а для

параболического

-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

= 5

y (v )

d x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

к.

 

 

 

 

 

 

 

 

 

 

 

 

Результат вычисления площади напечатать на БШ.

 

 

 

3 S .

Решить систему уравнений

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

=с,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y aLx + B ,x

 

=сг

 

 

 

 

 

 

 

ш

формулам Крамера:

 

сс- —— , у=

 

.где

 

 

 

л= а .

А 9'

X

 

 

* к -

 

 

 

 

 

 

 

 

 

 

К

 

 

 

 

*2

4,

У

 

 

 

 

 

 

 

 

 

 

 

 

 

 

п

 

 

 

 

 

 

 

 

Решить

задачу для

 

пар значений

а ,

,

ал ; т

 

пар зна­

чений

Btf

 

Вг

 

и

 

р пав значений

с,,

с4. Числа

п ^

5 ,

 

т 4

5 > Р 4

5

 

а также (ai°>

 

 

п '-(к >

 

 

 

 

 

 

,

( К )

М

 

 

 

 

 

,3р - заданы.

 

 

 

 

m ;( C f

,

ср ) ук= I , „ .

 

 

 

 

 

 

Если

 

 

 

 

I0 - ®,

то

предусмотреть в

программе

печать

на АЦПУ текста

СИСТЕМА НЕ ШЕЕТ ЕДИНСТВЕННОГО РЕШЕНИЯ

 

и

числовых значений

 

a, s as ,

£,, В ,,

с, , с г

вместе

с

их наи­

менованиями. Если же |д|

у .ТО- ®”, то напечатать текст

 

ЕДИН­

СТВЕННОЕ РЕШЕНИЕ СИСТЕШ и вывести числовые

значения ss

и о ,

а такие Ьоответствувдие

значения

af,as, В,,

£г ,с ,,

ct

 

с

нх

наименованиями»

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 0 .

Решение гсраевой

задачи для прямоугольника [a , B ,c ,d

]

U (х , t )

 

 

 

 

 

 

должно в точках

сетка

(о с ^

tL ),

 

 

из

4

= 1,

. . . ,п

,

j. = I, ...

,m ,

по x

с шагом л x

x

no

 

t с шагом A t

.удовлетворять следулцим условиям:

 

 

 

1)

U (a, t j) -0, j -

i ,

m

;

 

 

 

 

 

 

2)

U (x .,C )= sin (i

X .),

i = i ,...,n ;

 

 

 

 

 

3)

U (x .,d )

= a?. e ' ix t, i = f ,

n ;

 

 

 

 

 

4) U ( x , t .) =U (x_t ,

 

 

t j ) * U ( x ^ t)

 

 

jUfl

l — 2, . . .

t

/ X \

^ =

2 f

««. ,

/77“ i .

 

 

 

 

 

Вычислить значения

U (x i , tt)

при всех

i = I .........n

,

j

= I,

... t m.

t

если

п л

100 и

5

-

заданные

числа

(рис.12).

 

 

 

 

 

 

 

 

 

 

 

Полученные значения вывести на печать БШ.

 

 

Рис.12

41. Вычислить площадь кпиводинекно! трапеции, ограничен— но* осью абсцисс (0 4 &4 I ) ж днние* у * е

формуле трапеций:

Предусмотреть в программе деление шага h

удвоение числа участков п

до тех пор, доха дрр

приближении площади

SK и

5**,

не удовлетворит

венству:

- S j < e

,

где £

- заданна* тошеоть

ш

 

 

 

 

 

 

 

 

числений.

 

 

 

 

 

 

 

Начальное значение шага s числа участков

 

 

 

 

 

f)=

I,

п -

 

 

Результат вычислений напечатать на БПМ.

 

4 2 .

Найти шаксимальное из значений

 

 

 

% -

+

^ а ? + 9,785аг

 

 

 

Г

 

15,73 az

 

1Ж9

 

 

 

 

 

 

 

 

a s

rt * 9 ,7 8 5 co st

n _

Г,гл * 9,785tL

 

 

 

 

 

 

 

 

1~

 

 

П 4

100,

г ,,

Га . и t- ( i -

■f,

, n ) - заданные

числа.

 

Результат напечатать на БЕМ.

 

 

 

43 .

функция

у = f ( x )

задается таблицей п значенийxi

и у.

(/ 7 *

2 0 ).

Составить

програщу вычисления

значений

функции f ( x )

в

точках

« e [ a . , £ J

с шагом & х ( з s^ a< S^ х „),

используя интерполяционную формулу Лагранжа:

 

 

( х - х , Х х - х . „ )

 

 

(х-х-Хх-ас.^')

 

(хм -Ъ Х Ъ ч -Ъ + д

(хи ,- Ъ Х Ъ .г-Ъ *) ***'

+ ( X - XL ..)(X - Xi+*)

и.

 

nP“

Xi 4 X4- х и , -

(x i ~ х и , ) ( х 1 - х ;+ *)

 

 

 

 

 

Полученные значения f (х)

вместе

с соответствующими

значениями X вывести в

виде таблицы на печать

АЦПУ.

4 4 *

. Решить систему линейных алгебраических уравнений

порядка

п & 40

X=/fcX>B

 

 

 

 

 

 

методом Зейделя с точностью

s

, взяв в качестве

начального

приближения решения

Х 0=

( х ° , & " ).

В.

 

Матрицы/

аи

 

 

 

 

 

А=\

 

 

 

 

 

 

 

®гн г •• ч

а„„

 

 

 

 

 

Решение напечатать на БШ .

!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П 5

 

По методу Зейделя кавдое последующее приближение решения

 

 

 

 

 

 

>хп “ П)

 

находится

во

предвдущему

Х ‘"

*

( х

, at?)

по формуле

 

 

 

 

 

 

 

 

ск*о

Г 7 —

 

к9 ~

Г*)

 

в.

 

 

 

 

Точность е

считается достигнутой,

если для всех

 

4 = 1 , . . .

, п

выполняется неравенства

 

 

 

 

4 5 * .

Составить программу решения задачи Коши

 

 

 

 

 

 

 

у' = COS +у) +0,5(х-у)

 

 

 

с начальным условием

у (0 )-0

 

на отрезке

[0,1/

с

шагом

^ = 0 , 1

усовершенствованным методом Эйлера-Коши с

итерацион­

ной обработкой. Точность вычислений -

10

.Напечатать

на БШ

все пары значений ас, у .

 

 

 

 

 

 

 

 

 

Расчетные формулы усовершенствованного метода Эйлера-Ко­

ши с

итерационной обработкой для уравнения

y '= f(x ty )

с на­

чальным условием

у (х 0)= у0

 

имеют ввд :

 

 

 

 

y Z

= LKj

+fi f ( * * , у*)>

к = °><>-- v

 

 

 

 

где

L -

номер приближения,

а

ук^

 

=

у ( х к „ ) ,

 

 

 

Итерацию продолжают до тех

пор,

пока в

пределах

требуе­

мой точности два

последовательных приближения y Z °

ж

y Z

 

не совпадут.

 

 

 

 

 

 

 

 

 

4 6 * . Составить программу решения краевой задачи для диф­ ференциального уравнения

у " - Z x y ' - Z y - - 4 x

на промежутке [а , $] с шагом

h при граничных условиях

ыоу(а) + vi У'(а) =осг ,

js o у ( V+JS, у '(В)=ja£

116

методом прогоняя.

Метод прогонки дня дифференциального уравнения

 

tf"+p(x)y'+ t}(x)y=f(x)

заключается в

последовательном вычислении сначала чисел С- ,

d, ( i = 0.............

....... ) по формулам:

 

 

£\ =

<*,

 

.

J

. -

а * .

 

 

 

 

 

 

 

С,

Я

-----------------------------

;

d ^ f . h *

- гг.

■ci . , d . _ t ;

*

 

т.-гг.-с-

 

 

 

 

 

 

w e

 

/п..

z - g ( x L) h l

 

 

У -

- т - л

 

 

 

 

 

 

 

 

 

 

 

'*

 

 

 

'

/7‘ =

 

_

я

ч

- ;

а

 

 

 

 

х ,= а .

х п

£ -

i +

— : — Л

 

 

 

 

 

 

5 - а

 

 

 

 

 

 

 

 

 

п-

 

r t ^ iS ,

 

 

 

 

 

 

h

 

 

а затем -

в

 

 

 

 

 

 

последовательном вычислении искомых

функция

у. = у (х .)

по формулам:

 

 

 

 

-

А л + А сп-, < .,

 

 

 

 

значение

* п ~

А Л

+

* 0 *

 

у. =сс Щ-у.,),

<«/?-/, Л - 2 , . . . , i,o.

Внчмвление ковффщиентов уравнении

р(х.) , о (х) ж^(х)

в узлах %.L

вынести в подпрограмму.

'

Полученные значения функции напечатать на ББМ.

117

Ответы - решения

РЕШЕНИЕ «—.СИСТЕМЫ — УРАВНЕНИЙ — МЕТОДОМЗИЙДКПЯ2

1. ввод -E,A(WOO-N.N), B(W),XO(bO):N1 МАССЮ -.Y (40),X (W )X

ВЫП

2 1

3. ВЫЧ

YIH=XO U IZ

ПСВ

3 —7 ~ i. N X

4 . ВЫЧ -

XOfll = B in Z

5.ВЫЧ ы XOIll =A/I,7j-X/7l+XOlIIZ ПСВ ~ 5 - J - /. N Ж

ВШ - 2 Z

ПСВ ^

= /./VI

6.ЕСДИ «_ MOD(XOUI - П И ) - ) =E—T0—3Z ПСВ ^ В -1 =/. N Z

НАЛ — НА^БПМ-ХО(М) X

КОН _

I

2. ПОД

I

7. ВЫЧ - X f l l =ХО/1/ X

ПСВ ^ 7 - 1 = /./V I

ВЫХ -X

НАЧ «_. Ц

4 5 ^ .

РЕШЕНИЕ*- ЗАДАЧИ — КОШИI

1 . ъ т - х о - о - у = О X

2.ВЬГС-E=CQ$QCQ+Y0)*0,5.(XO+YO)-Y1-YQ*0, /•£_. =

х/=хо + о/ z

3.ВЫЧ ^ Y 2 -Y i- Yf= YO+0,05.(E+COS(Xi+Yf)*0,S. =

(X/-Y/)) I

ЕСДИ *- MOD(Y2 - Y i)-) Ю -Ч- TO-3Z

EMI — Н А - Б П М -Х i, YiZ

ВЫЧ - XO =X i- YO*Y/ X

ПСВ - z - / o x

КШ —X НАЧ —IX

 

РЕШЕНИЕ «

КРАЕВОЙ „ ЗАДАЧИ^ МЕТОДОМ _

ПРОГОНКИ I

I . ВВОД — A, В,Н, АО, Ai, А г,

ВО, B i , B Z : N Е

 

MAC

-

C (iS ),D (1 5 ), Y (18) К

 

 

ВЫЧ

~ С Ш а А i:(A 0 ‘H-At)~X>/f/=A2. H-AL.X=AS

 

ВЫЧ

— : У ? = / « Л / / = Л / - / _

N 2-N + /

S

2.

ВЫП

- Л 0 I

 

 

 

ВЫЧ

-

E = t+ P .H :2^ M = -(2-q.H > 2):E ^ z.

 

М 1 = ( 1 - Р . Н F i = F :£ S

 

ВЫЧ

_

C/I+ tl= t:(M -M f.C /II)^ B II+ {l= Ft.H 2-M f.E

 

gJi /.d /ц

I

 

 

 

ПСВ

_ 3 _ X

= RE

 

 

 

ВЫЧ ^ X = X * H g

 

 

 

ВЫЧ _

: R ЮR+ / E

 

 

 

ПСВ

^ . Z ^ N i l

 

 

4 .

ВЫЧ „

Y 11* ih (B 2 ‘H +81. C /I/.S /Il) :(30.H +8t. =

 

(СЦI + 0 ) 2

 

 

 

 

ПСВ

_

* i^ I = N Z

 

 

5 .

ВЫЧ ^ Y lI}= C /II-(H IT l-Y lI+ H )t

 

 

TTCB

^

Si-j I - R i

 

 

 

ВЫЧ ^

: R -R - i t

 

 

 

ЕСЛИ—

: / ? - > / ~ T 0 - 5 Л

 

 

 

НАЛ -

HA — З Я М — Y(NZ) S

 

 

KOH

-

E

 

 

 

10 .

ПОДПРОГРАММА — ВЫЧИСЛЕНИЯ_

КОЭФФИЦИЕНТ®- УРАВНЕНИЯ Ж

 

ВЫЧ

~ Р = -

2 . X ~ 9 = - 2 ^ F

= - * . X 2

 

 

ВНХ

w F

 

 

 

 

НАЧ

-

1 I

 

 

 

119

Л И Т Е Р А Т У Р А

 

1 . Математическое обеспечение ЭВМ "Минск-2(2 2 )"

в режиме Т .

Библиотека стандартных программ. Минск.

1968.

2 . МАТЮШЕВСКАЯ И.М. и др. Автокод "Инженер" для решения за ­ дач на машине "Минск-22М”. Входной язык. Сб. "M see-

матическое обеспечение ЭВМ "Минск-22" в режиме Т ",

в ш .6 , Минск, 1970.

3 . КАЛАЧЕВ В . М . , ЯКУБОВИЧ М.М. Программирование для ЭЦЗМ

"Минск-2" и "Минск-22". М ., 1971.

4 . САВИНКОВ В.М . Программирование для ЭВМ "Минск-22".М. ,19 7 2 .

5 . НЕМЕНМАН М.Е., ЦЕГЕЛЬСКИЙ В . И . , ЧЕРКАСОВА М.П. Программи­ рование на АКИ. Минск, 1972.

Соседние файлы в папке книги из ГПНТБ