Ъz ▐ *** Bilinear Transformation Program ***─▐е▐ Program to evaluate the BLT given an analog transfer (▐ function of a recursive (IIR) system. The result of the) 2▐ BLT will be of the form:/ <▐j F▐ A0 + A1z^-1 + A2z^-2 + ... + Anz^-nі P▐ H(z) = ------------------------------------Ю Z▐ 1 + B1z^-1 + B2z^-2 + ... + Bmz^-mФ d▐0n▐ Note: Order of numerator must be less than or equalxx▐ to the order of the denominator polynomial.~┌▐д▄▐ The program will prompt the user for output path of screen orь√▐ or printer.ч═▐+╙├ AS( ),BS( ),AZ( ),BZ( ),CARY( , ),FAHEX( ),FBHEX( ),FRACA$( )`Є├ FRACB$( ),AZN( ),BZN( ),FRACA( ),FRACB( )hЎи щnхюрMAXCOEFFSГ┤эй оФ▒ н );"*** Mapping S to Z Plane Using Bilinear Transformation ***"уП▐З▐ *** Prompt for Analog Transfer Function ***▐"й -▒ н);y"┘ "Enter the # of numerator coefficients (30 = Max, 0 = Exit) --> ",NNC╒,▀ NNC Г м ю:Ч┐ :▐ыExit to DOSй6▀ NNC Ф MAXCOEFFS О NNC Х м ┴ хА@┌ NГNNCЙ л о ЙЯJIГ(NNCЙ)ЙN"T█ ч :▐ыconvert index to characterd^▒ " Enter a(";C$;")s^";P$;н);"--> ";:┘ "",AS(N)lh┐ Nrr▒Ё|┘ " Enter the # of denominator coefficients --> ",NDCш├▀ NDC Ф MAXCOEFFS О NDC Х м ┴ хР░┌ NГNDCЙ л о Й IГ(NDCЙ)ЙN3є█ ч :▐ыconvert index to characteruў▒ " Enter b(";C$;")s^";P$;н);"--> ";:┘ "",BS(N)}╦┐ N┐б▒бл┘ " Are the above coefficients correct ? (y/n) ",C$чж▀ C$Г"N" О C$Г"n" м хДЮ▐!Й▐ *** Perform Bilinear Transformation ***'Т▐9Ч┌ IГ л NDCЙICARY(,I)ГQ┐ I[NUMГm&┌ IГ л NDCЙ▐0NUMГNUMК((NDCЙ)Й(IИ)И)Л(I)║:CARY(I,)ГNUM╘D┐ I╩N┌ IГ л NDCЙмX┌ JГ л NDCЙbCARY(I,J)ГCARY(I,JЙ)ЙCARY(IЙ,J)ЙCARY(IЙ,JЙ) l┐ Jv┐ I#─┌ IГ л NDCЙ1┼SUMOFASГ?■SUMOFBSГQ·┌ JГ л NDCЙu╗SUMOFASГSUMOFASИCARY(I,J)КAS(J)≥╡SUMOFBSГSUMOFBSИCARY(I,J)КBS(J)║╪┐ JЁфAZ(I)ГSUMOFASепBZ(I)ГSUMOFBSмз┐ IсД▐Н▐ *** Divide all Coefficients by B0 ***Ь▐'┌ IГ л NDCЙ>AZN(I)ГAZ(I)ЛBZ()UBZN(I)ГBZ(I)ЛBZ()] ┐ Ic*▐є4▐ *** Generate Fractional Coefficient Arrays ***╙>▐╪H┌ IГ л NDCЙR▀ Ъ├(AZN(I)) Х м FRACA(I)ГAZN(I):┴ p :▐ыalready fractionV\FRACA(I)ГЪ├(AZN(I))ЙЪ┘(Ъ├(AZN(I))) :▐ыFRACA(I)=fractional part of AZN(I)|f▀ AZN(I) Х м FRACA(I)ГЙFRACA(I)фp▀ Ъ├(BZN(I)) Х м FRACB(I)ГBZN(I):┴ ▌ :▐ыalready fractionzFRACB(I)ГЪ├(BZN(I))ЙЪ┘(Ъ├(BZN(I))) :▐ыFRACB(I)=fractional part of BZN(I)<└▀ BZN(I) Х м FRACB(I)ГЙFRACB(I)D▌┐ IJ≤▐▀╒▐ *** Negate b's for Filter Implemetation ***▒╛▐ёІ┌ IГ л NDCЙ╨юFRACB(I)ГЙFRACB(I)бй┐ Iхт▐ч▐ *** Convert Fractional Coefficients to Hexidecimal ***Х▐DР█ ` :▐ыinvoke conversion routineJЭ▐┼▐ *** Output Transfer Function Coefficients ***░▐╩█ B :▐ыinvoke output routineа$юк.┴ хя8│B▐ ========================== Output routine =============================%L▐+Vю3`й Xj▒ " Send coefficients to:"nt▒ н);"(S)creen"┘~▒ н);"(P)rinter"╙┬▒ н);"(R)un BLT program again"е▓▒ н);"(E)xit to DOS"к°▒Ьі┘ " Enter desired path --> ",PATH$╟▀ PATH$ Г "E" О PATH$ Г "e" м ю:Ч┐C╨▀ PATH$ Г "R" О PATH$ Г "r" м ▌kд▀ PATH$ Г "S" О PATH$ Г "s" м █ Л⌠н▀ PATH$ Г "P" О PATH$ Г "p" м █ Ч²ь┴ VёБ▐ЯЛ▐ -------------------- Output Coefficients to Screen --------------------ВЖ▐Щю й 7▒ " The coefficients in the S - plane:"=▒T(┌ NГNDCЙ л о Йd2IГ(NDCЙ)ЙN■<█ ч :▐ыConvert index to characterҐF▒ н);:▒"a(";C$;")s^";P$;н);"= ";жP▒ в"##.#####";AS(N);ЪZ▒ н+);:▒"b(";C$;")s^";P$;н4);"= ";d▒ в"##.#####";BS(N)n┐ N(xй \┌┘ " Press to view H(z)/b0 ",RPLYb▄юj√й u═▒ н);Ґ╙▒ "are mapped into the Z - plane using the Bilinear Transformation"хЄ▒ н);ПЎ▒ "and normalized by b0, yielding:"Жх▒р┌ IГ л NDCЙ8э█ ч :▐ыConvert index to characterbФ▒ н);:▒"a(";C$;")z^-";C$;н);"= ";|П▒ в"##.#####";AZN(I);іЗ▒ н+);:▒"b(";C$;")z^-";C$;н5);"= ";©▒ в"##.#####";BZN(I)г┐ Iпй ш"▒ н);$,┘ "Press to view fractional parts of the coefficients ",RPLY*6ю2@й zJ▒ " FRACTIONAL part of each coefficient is converted to HEX:"ЄT▒ " (Note: b's are negated for implementation.)╨^▒Жh▒ н);"FRACTION";н");"HEX";н8);"FRACTION";нG);"HEX"Эr▒|┌ IГ л NDCЙ>├█ ч :▐ыConvert index to characterg░▒ н);:▒"a(";C$;")z^-";C$;н);"= ";┐ ▒ в"##.#####";FRACA(I);єє▒ н);"*2^15 = ";FRACA$(I);нў▒ н+);:▒"b(";C$;")z^-";C$;н5);"= ";Й╦▒ в"##.#####";FRACB(I); б▒ н?);"*2^15 = ";FRACB$(I) л┐ IY жй :┘ " Press to return to Output Path Menu ",RPLY_ Ю▌ґ Й▐ -----------------------------------------------------------------------Ё Т▐!Ч▐ ------------------- Output Coefficients to Printer --------------------!▐² н);"FRACTION";н");"HEX";н8);"FRACTION";нG);"HEX"╗$H²╨$R┌ IГ л NDCЙЙ$\█ ч :▐ыConvert index to character%f▀ I Ф NDCЙ м п%p² н);E%z²"a(";C$;")z^-";C$;н);"= ";:² в"##.#####";FRACA(I);f%└² н);"*2^15 = ";FRACA$(I);r%▌² н+);╛%≤²"b(";C$;")z^-";C$;н5);"= ";:² в"##.#####";FRACB(I);л%╒² н?);"*2^15 = ";FRACB$(I)т%╛┐ Iз%І▌(&ю▐ -----------------------------------------------------------------------v&й▐ =======================================================================|&т▐й&ч▐ ======================= Integer to Character ==========================п&Х▐ч&РS1$ГЪ⌠(I)Л&ЭS2$ГЪ⌠(N)Щ& C$ГЪ┌(S1$,)' P$ГЪ┌(S2$,))' ▀ I Ф м C$ГЪ┌(S1$,)D'$ ▀ N Ф м P$ГЪ┌(S2$,)`'. ▀ I Ф c м C$ГЪ┌(S1$,)|'8 ▀ N Ф c м P$ГЪ┌(S2$,)┌'B ▌п'L ▐ =======================================================================ж'V ▐$(` ▐ ===================== Convert Coefficients to Hex =====================*(j ▐<(t ┌ MГ л NDCЙ\(~ FAHEX(M)ГЪ°(FRACA(M)КЧ▐)|(┬ FBHEX(M)ГЪ°(FRACB(M)КЧ▐)√(▓ ▀ FAHEX(M) ФГ м і ╡(° FAHEX(M)Г▒ИFAHEX(M)л(і ▀ FBHEX(M) ФГ м ╨ Х(╟ FBHEX(M)ГЪ░ИFBHEX(M))╨ FRACA$(M)ГЪ (FAHEX(M)))д FRACB$(M)ГЪ (FBHEX(M))&)н ┐ M,)ь ▌z)Б ▐ =======================================================================