Главная  Развитие народного хозяйства 

[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] [128] [129] [130] [131] [132] [133] [134] [135] [136]

«ПУСТ""ПР0ЦШ0ЖР(ЫДГМ1],Е;У[Н]ЛЮ)™Т1.Ф"2;Ч=зо;К1=Н+ ;Ы1=Д1;Н=1;"ШЯ"т"ТАБ»1,ЕЕ(Х) ,IM(X),P/(dP/cflC) ,ИТЕР; "ДЛ"1=ГШ ГДОКГШ1(Б[1>АГ1]ГЕ1«0СЯ1]=(К"1+1>В[1]);2Л=0;Е1=

-1;Е2=1;ез=о;И1=о;И2=о;Из=01=2;1Л=ю;Е1=-1;1Е1=о;Пол(к,б[К1],-

1,0;f 1Д1) ГЕ"Р1-О"Т0ЧйЗ=-1; ""4) ;ПОЛ(Е, 5Ki3,1 ,0;F2 .12) ГЕТ 2=(JTO«(E3=iraA";F3=BCKi];l3»0;FK=F3t2;E"F3(rSOm";CTE "*ЗНС; *ДЛ"1=11!1"ГДОК"Ш"С1;1]=(К-1+1)5<В[13; ЗДЕЛ(Р2-Р1,12-11 ;Д ,1Л;$ й);ДЕЛ(ГЗ~-Г2-Ф,13-12-В, Ы,Щ, 1Ф) ;W-F2-F3-; IW=I2-I 2lS;S=Wt2-lWt2-4J<(F3xi~l3;IS-2>«tKlW-if>«(F3>«lS+I3x?&);E*IS=0

S)*V(S!2+IS.r2))/2);lG=ISA/2rES>OT0"G=T«HEA(G=ABS(IG);iG=Tx S IGN(IS))); E*G:IllxlG<0TO" (.G*=W-GriG=l№-TG)"iM»(G=V/+G; IG) ;ДЕЛ(2хГЗ, 2<I3,G, IG;G, IG) ;IU .H=Hl>G-IHlxlG; l№=inixG4-HixIG;R=R 3+H ;И=ЙЗ+1Н ;ПОЛ(К, BTKll ;F, IF) ;FH=F 12+IF J 2; -E«FH>100xFieTO"( С=</2;1С=1С/2;М);Е1=Е2;И1-И2;Е2=Ез;И2.-ИЗ;ЕЗ=Е;Ю=И;"ЕЕ=1Х ОЯА-iirEHt 2+lHt<2f 2"ТО"(П0Л(Н-1,С[К] ,R,H;D,ID);rE»FN/(D!2+I Ъ12>E 12T0"HA4);Fi-F2;Ii=l2;F?-F3; I2-I3;F3=F; l3=IP;Hi=H; IHl= 1Н;ДЕЛ(1,0,G,IG; Л; 1Л); J=J+1 ;FM; E" J4"T0""HA" 3;. Z=K-K+i; "E" АВЗ(ИЗ)<Е"Х0(иГ2>Ез*Д[г]==О;Ц«1"НА"9;5.ЬК;£[1]=В[1];"ДЛ"1=2

Ш"1Д0МШ11Б[13-БЕ1]+ю><БГ1~1])т(игг]=Ез;и[г+1>Ез;¥[г>ц 3; V[z+i>-li3;U=2rm"r;6 ЛНС-.1;р=-2>ЕЗ;а=Ез {2+из12; БГ1>в[1] ;Бр

2>B[2]-PxBCi3; ДМ=3"Ш"ГД0Г ВШ"£[1]=В[ 1]-РхБ[:1-1]Ч1хБ[1-2] );ТТЕ"ЗНВ;К1=МГЙЬМ;да1=111да"М»ШБ[1]==Б[13;СТЕЗН "Б; J=0;E*K>2"I0»m»2;ft=O;T=--Bt2]/]3[l3;ЕК=1Т0(Z=K ;U[Z1=<D [Z ]=0; "НАЗ); Z-H-i;T=T/2I==T! 2-В[3]/Б[ 1]; Е"Ц<0»ТО»(1)=1/(-1)) ;U[N >Т;и[2К;У[К>Н0;¥[2>1);Ц=37т<1)=1/(В);и[г]=1-В;и[НК+Р; Z ]=0; V[W]=0); 7„Э-1;8 .R=U[Z+3 ] ;H-¥[Z+3 ] ;ПОЛ(М ,A[Nl] ,R ,И;Р, IF) ;II0 Л (Н-1 ,ClDiil ,R.K; Д ДД) ;гД==Д 2+1Д! 2;Д=Е \ 24-IF \ 2; "E*0"TOX=v/a/ ГД);ЕЭ=0"ТО"М"П;1ЕЙВ?"ТАБ1,Е,И,Д,0;9.Э==0;"Е"Д<2ТО".КА8;й =-и;11 ..й1В"ТА£"1,Е;И Д, J; "E"J>4"T0"HA"ii; "Е"Ц=1"Т0"т"5;"Е"Ц=

2"Т0"НАб;й/СТЕ-*3№й,С1ТД"В[М1];С1[К1];сгк];Б[М];и[м];\тйГ K0H•;ПP0Ц•fl0Л(KjA[KiЗ,RД;FДF)"ПУCT-i,F=A[l]ДF=O;»ДЛ"J=2liI»i "Д0К+1"ВШЧС=Е><Р-1хПЧА[Л]ДЕ=«>:1Г+1хГ;РМЗ)"К0Н»;"ПР0Цда1(А, IA, ВД В;С, 1С)"ПУС1"1 ЛС=В? 2+IB! 2;C-(A<B4-IAa В)/1С; IC<B«IA--A«IB

Рис. 1-318. Программа нахождения корней многочлена по методу Мюллера.



7,8,-4"К0Н"О "ТАБЛИЦА"!

ЕЕ(Х>

1М(Х)

Р/(сД>/сГХ)

-.10бСООо,-8

ЛОООООю 1

.2120000-9

-.10б000и-6

.iOOOOOjj 1

.21200С-9

.200000b 1

-,200000и 1

Рис. 1-319. Пример нахождения корней многочлена P(x)=xfi-2x-2x+6-7х+8х-4

но методу Мюллера.

♦ШСТМП.О=0;О."Е"г(1/3>1О1О)*0"Т0»Чо<1+2Г51АЮ);а1<1-2;П11=О;Н l=H+i;Ml4*l; ХТЕ"«ЗН"ВД ВД; В[К]=1; В[К13=-Й;[Ы1; 1В[Н]==0;1ВГЫ1]

=-1М[ы]; дг1=ы-гш"-гдо-1вьи1"(ва >i; ib[i>o; •ДJ"J=i+l*•lll•Гд

0"HrB!fl*(C=B[J];BCJ3=IB[J]>«lM[I]--C<EE[I];lB[J]==-IB[J>EEriKx IM[I];EJ<lir3;0(BrJ3=S[ J]+B[J+1]; IB[J >IB[J]+I3:J+1]))) ГДЛ"1 =WOT"MaK"B[l]==(H+l~l)><B[l3;J=0;i.j<r+l;R=RErj];l=lMrj];n OICM,ArMi]»E»I;F,IF);noi(H-i,B[K],R,IJC,IC);W"EE""AE";lG=Ct 2+IC12;G=(FxC+IFxICyiG; IG=(CxIF-FICyiG; "E"I>0"IO" ("ЕСОТШ "IG=0"I0"fiE=4BS(G+IG)m"ER=i(Gt2+lGt2);EER==O"T0"IlA»2;tIG =0"ИС<0"ХО"Ф=Я"ИНА" ($=ARCSIN(IGAR); •Е"0<0ТОФ=Лх810Ы(1С)-Ф)> ;nn=4in+(*E"I=0lO"GxEXP(RKT)HHA»"E"I>0"TO2xEIixEXP(RxT)xCOS(Ix МутОУ; 2, J<rTO "Ш"Ц "ДРОДЕ";ffiE-nn; "1IA3B"II(I); "ШВ •"CTP2 ,Ш/СТР« ;Е"НТ=ЧГ "К"ТК="Ч1Г"Т0« (T3=SC. 5>«TK/HT)><KT ГРАЗ01;"В1ВГ"ЗА"ТАБ"1,Т,ПП,Т,ПП;"ДЛТ=О"1"КГД0"ТЗ"ЙШ"("РАЗ »Q. , 0; Т2»Т+!Гз+т; K1=II(T); K2=€1(T2); "PA3"Qi; "а1В""ТАБ" 1 ,Ti, Ki Д2,К2))?"РАЗ"аТД"™Ц"П0Л(К,А1:К1] ,R,I ;F, 1Р)"ПУСТ"1.Р=ДГ1]; IF=0;ArJ=2lll"i*K+rBlIl(C==R<F-I«lF4-A[J];IF=RxIF+I><F;F=G,)K

0Н";.Б[11];1£[М1]та»0

Рис. 1-320. Программа обратного преобразования Лапласа для рациональных алгебраических функций в случае простых корней.

Дифференциальное уравнение первого порядка может быть записано в виде

, с начальным условием у(и) = уо, где t - аргумент.

Наиболее простым способом решения обыкновенных дифференциальных уравнений является метод Эйлера, по которому

yi+i-yi + ffifi-fh), i = 0, 1,2...,

где - шаг интегрирования.

Линейные динамические звенья первого и второго порядков и соответствующие

им расчетные соотношения приведены в табл. 1-47. Расчетные соотношения для некоторых типичных нелинейностей приведены в табл. 1-48.

Кроме просчета линейных и нелинейных звеньев, который осуществляется на каждом шаге, в программе должен быть предусмотрен вывод заголовка таблицы и вывод самой таблицы, в которой печатаются значения нужных переменных, а также счетчик времени (аргумент t) (рис. 1-322). Начальные условия могут вводиться в основной программе, которая обычно оформляется как исполнительная информатива,



»а1иМ=1;Ы=8;НТ=1;ТК10;"НА"ЛАДГД"А[2>120,210; Ei;[83=O,-l,-2,-3,-if,-5,"l,-lM[8>O,O,O,O,O,O,l,-ri{0rO

П11=1+(-5)хЕХР(-Т)+2хЕХРС-ЗхТ)-(ЕХР(-4хТ))+3/17хЕХР(-5хТ)+ .?9lO42roixEXP(-I)xC0S(-.2W9iQO+T)

"ТАБЛИЦА"!

Т ДП Т НПО

о .2G00a,-4 .6OO0ip 1 ,993 00

.1000, 1 .2281ю-1 ,7000ю 1 .9978 00

.2000в 1 .2558ю О .800С 1 .99т 00

.3000» 1 .б171в О .90(Юэ 1 .9991в 00

.iOOOb 1 .864 О ЛООрв 2 9996 00

.,5000b 1 ,9671в О .llOQa 2 ,999% 00

Рис. 1-321. Пример нахождения оригинала для изображения

120р + 240

P(p+l)(p+2)(p-f 3)(p-f 4)(p-f 5)[(p-f l)?-f 1]"

Таблица 1-47

Расчетные соотношения для линейных динамических звеньев первого и второго

порядков

Передаточная функция

Структурная схема

Алгоритмические выражения

Я(р) =


b„p + bip-\-b аор2 + ар + о2

Примечание. В алгоритмических выражениях Т соответствует Д.

У1=(Х-А1ХУ0)/А0;

yo==yO-fTxyi;

У=В1хУ0-1-В0хУ1

У2-(Х-A2xУ0- -AlxУl)/A0; У1 = УЦ-ТХУ2; У0=У0--ТхУ1; У=В2хУ0--В1х ХУ1-1-В0ХУ2



[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] [128] [129] [130] [131] [132] [133] [134] [135] [136]

0.0011