Frekvenční analýza, syntéza a optimalizace


[ Homepage ] • [ NetCAD | FAQ | Uživatelé ]

Logo

Obsah


Všeobecně o programu

NetCAD je frekvenční analyzátor lineárních modelů elektrických obvodů. Je doplněn syntézou LC filtrů a útlumových odporových článků a modifikovatelnou optimalizací hodnot součástek podle obecně stanovitelného kritéria. Nahrajte si archiv (290KB) s demoverzemi programu. Demoverze ADEMO.EXE je samoběžící ukázka možností programu NetCAD, demoverze FREEDEMO.EXE je plně funkční program s omezením na 6 prvků a 4 uzly. Prosím, pošlete mi také na vás kontakt.

e-mail:  ponikelsky@iol.cz

Pokud nemůžete použít e-mail, můžete odpovědět pomocí formuláře HotInfo. Pokud máte zájem o programu NetCAD, obraťte se na firmu

CADware, Hálkova 6, 461 17 Liberec 1

e-mail:  cadware@pvtnet.cz
URL:     http://www.cadware.cz
telefon: (048) 510 6131

Firma CADware se zabývá distribucí SW pro návrh PCB a dalších CAD programů a o NetCAD napsala...

Jádrem programu je numerická analýza (výsledkem je komplexní hodnota pro určitou frekvenci), ke které je přidána syntéza a optimalizace. Numerická analýza zůstane základním jádrem programu vzhledem k nejširším možnostem, kdy je možné použít také číselné modely prvků, které jsou v kdejakém katalogu a lze si je vektorovým voltmetrem změřit (výpis Y nebo S matice v závislosti na frekvenci). Semisymbolická analýza (výsledkem je racionální lomená funkce proměnné s s vyčíslenými koeficienty) bude zavedena, neboť umožní podstatně zrychlit optimalizaci pro zapojení, které neobsahuje číselné modely prvků. Symbolická analýza nebude nikdy použita.

Syntéza zastupuje uživatele ve fázi zadání dat a podle požadavků uživatele generuje obvod jako vstupní data pro analýzu.

Optimalizace pak modifikuje hodnoty prvků pomocí analýzy zadaným způsobem a podle zadaného kritéria.

Program NetCAD je určen pro DOS a běží také v DOSovém okně Windovs 3.1/95/98/NT(SP5). Programováním pod Windows se nechci zabývat, neboť se potřebuji věnovat odborným problémům. Víc vysvětluji v závěru této stránky. NetCAD je psán v jazyce C/C++ a přeložen překladačem Borland C++ 3.1.

Požadavky na HW

PC AT, 640KB RAM, video VGA, koprocesor doporučen, jinak je emulován, klávesnice AT (dnes bez omezení splněno). Výhodou je rychlost výpočtů v plovoucí řádové čárce. Program je testován na PC s CPU Pentium 133MHz.

[ Top ]

Popis programu NetCAD

SCHÉMA

Tabulka zapojení

Obvod se zadává do tabulky zapojení obvodu a skládá se z těchto prvků:

  Volný řádek  (editační prvek pro vyhrazení místa v tabulce prvků)
  Komentář     (editační prvek pro zpřehlednění tabulky prvků)
  Zkrat uzlů   (umožňuje laborování s obvodem, jinak je zbytečný)
  Rezistor
  Kapacitor
  Induktor
  Bipolární tranzistor         (ve verzi 01.xx jeden model)
  Zdroj proudu řízený napětím
  Zdroj napětí řízený proudem
  Ideální operační zesilovač
  Reálný operační zesilovač    (ve verzi 01.xx jeden model)
  Ideální gyrátor
  Ideální imitanční invertor
  Kapacitní imitanční invertor
  Dvojný kapacitor

Verze programu 01.xx se vyznačuje jednohodnotovými prvky. Tím je také dáno omezení verze 01.xx v sortimentu prvků, kdy nejsou možné prvky typu vázaná indukčnost, mikropáskové vedení, potenciometr atd.

U každého prvku se zadává typ a pak v závislosti na typu uzly, hodnota, optimalizace a index. Parametry prvku lze také samostatně editovat.

Seznam prvků se doplní číslem vstupního a výstupního uzlu. Vstupní a výstupní uzel může být jakýkoliv uzel z aktivních uzlů. Aktivním uzlem se rozumí uzel zapojený do obvodu s vyjímkou uzlu vztažného.

Maximální počet prvků byl stanoven na 500 a maximální počet všech uzlů na 64, tedy 63 aktivních uzlů pro ostrou verzi, přičemž maximální číslo uzlu je o 1 menší, než je maximální počet uzlů. Program je konstruován tak, že uživatel nemusí používat kontinuální číslování uzlů.

[ Top ]

ANALÝZA

Nastavení parametrů

Pro obvod se počítá DC hodnota a/nebo frekvenční průběh těchto parametrů:

  Au,  výstup naprázdno         ;napěťové zesílení
  Auo, provozní zátěž výstupu
  Aut, ekvivalent k Apt
  Ap                            ;prosté výkonové zesílení
  Apt                           ;přenosové výkonové zesílení
  Apa                           ;dosažitelné výkonové zesílení
  MAG                           ;maximální Apa
  Apo                           ;výkonové zesílení při |y12|=0.0
  MSG                           ;maximální stabilní výkon. zes.
  Bp = 1/Apt                    ;provozní útlum
  Zinp, výstup naprázdno
  Zinp, výstup nakrátko
  Zinp, provozní zátěž výstupu
  Yinp, výstup naprázdno
  Yinp, výstup nakrátko
  Yinp, provozní zátěž výstupu
  Ginp, výstup naprázdno        ;vstupní činitel odrazu
  Ginp, výstup nakrátko
  Ginp, provozní zátěž výstupu
  Zout, vstup naprázdno
  Zout, vstup nakrátko
  Zout, provozní zátěž vstupu
  Yout, vstup naprázdno
  Yout, vstup nakrátko
  Yout, provozní zátěž vstupu
  Gout, vstup naprázdno         ;výstupní činitel odrazu
  Gout, vstup nakrátko
  Gout, provozní zátěž vstupu

Program provádí kontrolu zapojení na všechny chyby zapojení prvků, které by vedly k nemožnosti provedení analýzy nebo optimalizace. Program informuje o nespolehlivém výpočtu impedance, který nastane tehdy, když je na vstup nebo na výstup zapojený singulární aktivní prvek přípustným způsobem. Singulární prvek je takový prvek, který nemá Y matici a jeho vliv se respektuje transformací Y matice obvodu.

Numerické chyby analýzy jsou ošetřeny tak, že program by neměl v důsledku vzniku numerické chyby havarovat. Program o těchto chybách vydá zprávu v číselné tabulce. Tyto chyby nelze eliminovat kontrolou, protože vznikají například při frekvenci rovné frekvenci pólu a podobně.

[ Top ]

GRAF A TABULKA

Přenos RIAA

Výstup výsledku analýzy je značně modifikovatelný a také k jeho modifikaci je vyhrazen největší počet uživatelem nastavovaných parametrů. Lze volit z těchto typů zobrazení výsledku:

  Číselná tabulka
  Kartézský graf
  Hodograf        Re{z}+jIm{z}
  Polární graf    Mod{z}*e^(jArg{z})
  Smithův graf    jen pro činitel odrazu

a lze zobrazit tyto složky:

  Re{z}
  Im{z}
  Re{z} a Im{z}
  Mod{z}
  Arg{z}
  -dArg{z}/d2PIf    (skupinové zpoždění)
  Mod{z} a Arg{z}
  Mod{z} a -dArg{z}/d2PIf

Nejuniverzálnější použití má číselná tabulka a kartézský graf. Ostatní grafy zobrazují jen ty složky, které jim přísluší a na uživatelskou volbu zobrazované složky se nebere zřetel.

Grafický výstup výsledku není vždy možný. Nemožnost grafického zobrazení vzniká u chyb výpočtu parametru nebo když by se měl počítat logaritmus nuly. Pokud existují data k zobrazení a došlo k chybě, je vždy použitelná číselná tabulka, ve které jsou chyby výpočtu popsány. Zobrazení ostatních grafů je možné, jen když je výpočet správný z hlediska toho kterého grafu a jeho nastavení.

V programu NetCAD je od verze 01.04 možnost manuální optimalizace přímo v libovolném grafu. Můžete ručně měnit hodnotu vybraného prvku a hned vidíte, jaký má tato změna vliv.

[ Top ]

OPTIMALIZACE

Aproximační funkce

Optimalizace hodnoty prvků je proces, který zvolenou strategií u vybraných prvků zvoleným způsobem upraví jejich hodnoty tak, aby zvolená složka počítaného parametru obvodu byla co nejlépe aproximována zadanou aproximační funkcí. Lze vybírat z těchto strategií:

  Minimum kvadratické odchylky
  Minimum absolutní odchylky
  Minimum rozdílů              (jen pro AC analýzu)

Při použití strategie Minimum rozdílů nemusíte zadávat aproximační funkci. Jde o co nejlepší vyrovnání průběhu frekvenčníh charakteristiky bez ohledu na absolutní hodnoty.

Optimalizaci hodnoty lze zvolit u libovolného prvku, který má hodnotu, z těchto typů optimalizací:

  V řadě E3
  V řadě E6
  V řadě E12
  V řadě E24
  V řadě E48
  V řadě E96
  V řadě E192
  Změna hodnoty na 1. místě za desetinnou tečkou
  Změna hodnoty na 2. místě za desetinnou tečkou
  Změna hodnoty na 3. místě za desetinnou tečkou
  Změna hodnoty na 4. místě za desetinnou tečkou

Na následujícím obrázku je výsledek optimalizace 13-ti obvodového VKV filtru podle aproximační funkce, jejíž zadání je zobrazené v předchozím obrázku. Ladily se pouze cívky. Ve Smithově grafu se zobrazuje vstupní činitel odrazu při provozním zatížení výstupu filtru.

VKV filtr, Smithův graf

Na dalším obrázku je vstupní činitel odrazu téhož filtru zobrazený v kartézském grafu.

VKV filtr, kartézský graf

Optimalizovat lze libovolný počet prvků najednou. Je použita zobecněná jednorozměrná optimalizační smyčka, která umožňuje také redukovanou vícerozměrnou optimalizaci. Optimalizační smyčka prochází fiktivní tabulkou kombinací R prvků bez opakování z N prvků, kde R je rozměr smyčky a N je počet optimalizovaných prvků. Když R=1, je fiktivní tabulka shodná se skutečnou tabulkou prvků. Prvky se tedy seskupují do skupin o počtu R prvků, přičemž v jednom okamžiku je možné změnit hodnotu jen u jednoho prvku, ale v dalším okamžiku u jiného prvku. Mění se hodnota toho prvku, u kterého dojde k největší změně chybové funkce směrem k optimu.

Tato smyčka není plně R-rozměrná, ale je jejím přiblížením. U plně R-rozměrné smyčky lze změnit hodnoty všech vybraných prvků ve skupině najednou. To předpokládá velmi mnoho výpočtů až tak, že by jste se výsledku asi nedočkal.

Jinak řečeno, optimalizace probíhá tak, že v prostoru hodnot prvků, jehož dimenze je rovna počtu optimalizovaných prvků N (co osa, to hodnota jednoho prvku) se koncový bod vektoru optimalizovaných hodnot může pohybovat jen rovnoběžně k R osám tohoto prostoru, a to k osám aktuálně vybraným.

Optimum je v tomto programu signalizováno žádnou změnou hodnoty v ani jedné kombinaci prvků při průchodu celou fiktivní tabulkou. Nejde o matematicky definované kritérium podle iteračního vzorce. Optimalizační smyčka musí vždy zkonvergovat do optima. Jde o logicky konstruovanou smyčku, ve které se nacházejí dílčí matematické postupy a blok analýzy. Nezaručuje se, že optimalizační smyčka nalezne absolutní optimum.

První reakce optimalizační smyčky je bouřlivá a hodnotu prvního nebo některého blízkého dalšího prvku podstatně změní, aby zase hodnotu vracela zpátky. Z počátku se jen jedním nebo několika málo prvky dociluje splnění celkového kritéria, které se ale musí rozložit na více prvků. Aby tento jev byl potlačen, nastavuje se maximální přípustný počet změn hodnoty optimalizovaného prvku. Tím se po dosažení tohoto počtu změn hodnoty prvků automaticky vynutí přechod na optimalizaci hodnoty dalšího prvku nebo kombinace prvků.

Když už je optimum téměř nalezeno, rychlost konvergence se snižuje a program v několika průchodech celou tabulkou prvků provádí drobné změny hodnot. To v plném rozsahu platí u typů změn hodnot DT1DT4, které drobné korekce hodnoty umožňují. Když je zvolena změna hodnoty z nějaké procentní řady E, rychlost konvergence je skoro konstantní a velká. Doporučuje se u složitých případů interaktivně modifikovat parametry optimalizace a počáteční hodnoty součástek.

Nechtějte optimalizovat nesmysl. Obvod musí být správný a počáteční hodnoty by měly být přibližné. Obecně musí platit, že počáteční poloha vektoru hodnot optimalizovaných prvků musí být v dosahu absolutního optima. Pokud tento vektor bude mimo dosah optima, bude se pohybovat z hlediska funkce obvodu nesprávným směrem. Matematicky bude všechno v pořádku až na to, že hodnota nějakého prvku bude ubíhat do nesmyslných končin. Optimalizační smyčka sleduje pokles hodnoty chybové funkce.

Do procesu optimalizace lze vstupovat, vynucovat přechod na další prvek a předčasně ukončit optimalizační smyčku. Průběh optimalizace je podle volby zobrazován jen odchylkou, nebo plně průběžně zobrazován. V programu NetCAD je od verze 01.04 možnost grafického zobrazení průběhu optimalizace.

[ Top ]

SYNTÉZA

Čebyševova DP 5. řádu

Program je doplněn o syntézu LC filtrů. Syntéza LC filtru se provádí metodou normované dolní propusti a požadovaný typ filtru se získá kmitočtovým a impedančním odnormováním. Lze volit z těchto aproximací útlumové charakteristiky:

  Butterworth
  Čebyšev
  Inverzní Čebyšev,  lichý řád
  Cauer-Čebyšev,     lichý řád
  Bessel
  Feistel-Unbehauen, lichý řád

Sortiment aproximací útlumové charakteristiky je zvolen tak, aby vyhověl pro aplikace LC filtrů ve video a vf technice. Upřesnění je v závěru.

Výsledný LC filtr se realizuje jako:

  Dolní propust
  Horní propust
  Pásmová propust
  Pásmová zádrž

Základním článkem LC filtru je:

  PI článek
  T článek

Výstupem syntézy LC filtru je zapojení filtru v tabulce prvků včetně definice vstupního a výstupního uzlu a případně přepočtu výstupní vztažné impedance. Filtr se impedančně odnormovává vzhledem ke vstupní vztažné impedanci. Toto řešení umožňuje ihned spustit analýzu filtru. Vlastnosti filtru lze pak korigovat pomocí otimalizační smyčky. Například lze filtr doplnit ztrátovými rezistory cívek nebo korektorem fáze.

Syntéza útlumového odporového článku je drobným doplňkem programu pro výpočet oboustranně přizpůsobeného odporového článku. V radiotechnických zařízeních se to útlumáky jenom hemží.

[ Top ]

Závěr

Zatím nejsou v programu implantovány aproximace inverzní Čebyšev, Cauer Čebyšev, Bessel a Feistel Unbehauen.

Současná verze 01.04 programu NetCAD se hodí pro analýzu do frekvence asi 500MHz a v některých případech až do 1GHz. Zvyšování sortimentu prvků, zvyšování "horní mezní frekvence" typem a strukturou modelů prvků, zlepšování ostatních vlastností programu a doplňování nových vlastností je věcí dalšího vývoje.

Také hledám zájemce a lépe firmu, která by byla ochotna se ujmout převodu programu pod Windows a zajistit systémový vývoj. Potřebuji se věnovat odborné části programu.

Zde je další ukázka aplikace programu NetCAD.

[ Top ]

[ Homepage ] • [ NetCAD | FAQ | Uživatelé ]