Главная  Радио и связь 

0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 [ 79 ] 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127

Щ}11ЛО«вНИЯ

1. Бжж-схеш внчвсяеввя сртаговвльвнх мшячуиенов Чвбшюва В8 задашвом нвоввстве точек

Пусть задано мнояество точек ix} (t=1,2,...,т). Рекуррентные соотношения, определяющие ортогональные многочлены Чебышева gjix) (te=0,l,...,OT-1), имеют вид

8q=* g{x)=x-a; g{x)={x-a)g (x)-b 2{x) (1г=2,3,...я-1).

Значения коэффициентов аэа, о, находятся по формулам т т

1 т

а = ~ s /. ci. =

/Ввод m,x/ I

Вычисление aO i

1=1,IB

g(l)=l;g1 (l)=x(i)-aO

1=2, m-2 у-

Подпрограмма вычисления параметров a(l), b(l)

al=g1(J); gl(J)=(x(J)-a(l))g1(J )-b(l)g(J ); g(J)=a1



11одцрогра1ма вычисления параметров а(т-1) ,Ъ(пн1)

+< к=1 ,т

Х=Х(К)

g(0)=l ;g(l)=x-aO

/ывод g(0),g(i)/

-<( 1=2. in-1

g(l)=((x-a(l))g(l-1 )-b(l)g(l-2)

/Швод g(l/ -1-

[ End ]

Шчисление

значений

многочленов

Чебышева

на заданна

множестве

точек

!

гг, Щюграуиа вычисления (угогоиальных многочленов Чебывева на языке BASIC

1 нж *««««»»«««««««««»«««««««««««««««»«««««««««»««

2 REM ПОСТРОЕНИЕ ОРТОГОНАЛЫЖ МНОГОЧЛЕЙОВ ЧЕБЫШЕВА

3 RM НА МНОЖЕСТВЕ m < 16 ТОЧЙС

4 НЖ «««««««««««««««««««««»*«««««««««**«««««««««»« 10 DEPIHT I-J, М: DIMxdS), g(15), g1(15), a(14), b(14) 20 CIS

30 INPUT "bBeflHTe m - число точек множества (m < 16) m

40 PRINT "Введите последовательно значения xl;l = l,...,m"

50 aO = 0

60 FOR 1 = 1 TO m

70 INPUT x(l): aO = aO + x(l)

80 NEXT 1



90 аО = аО / m

100 POR 1 = 1 ТО m

110 g(l) = 1; g1(l) = 1(1) - aO

120 NEXT 1

130 POR 1 = 2 TO m - 2

140 GOSUB 320

150 POR J = 1 TO m

160 a1 = g1(

1T0 g1 (3) = (x(3) - a(i)) * glii) - b(i) * g(J)

180 g(J) = a1

190 NEXT J

200 NEXT 1: GOSUB 320

210 FOR k% = 1 TO m

220 X = x(k%)

230 g(0) = 1: g(1) = X - aO

240 PRINT USING "x = Ш.Ш x;

250 PRINT USING " gO = ##.### gl = ##.###"; g(0); g(1);

260 POR 1 = 2 TO m - 1

270 g(i) = (X - a(l)) . g(l - 1) - b(l) * g(l - 2)

280 PRINT USING " g» = ##.###"; 1; g(l);

290 NEXT 1: PRINT

300 NEXT k%

310 GOSUB 390: END

311 REM -------------------------------------------------

312 REM ПОДПРОГРАММЫ

313 REM Вычисление коэффициентов a(l), b(l} рекуррентных

314 REM соотношенздй для многочленов Чебышева

315 REM-------------------------------------------------

320 кО = О: n1 - 0: bO = О: nO = О

330 POR J = 1 то m

340 kO = kO + x(J) . g1(J) 2: n1 = n1 + gl 2

350 bO = bO + x(J) * g(J) * gl (3): nO = nO + g( " 2

360 NEXT J

370 a(l) = kO / nl: b(l) = bO / no

380 RETURN

381 REM--------------------------------------------------

382 REM Пауза

390 PRINT : PRINT "Нажмите любую клавишу для продолжения"



0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 [ 79 ] 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127


0.0083