Häufig werden integrierte Schaltungen oder Teilschaltungen oder einzelne Bauteile in solchen Schaltungen implementiert, die bestimmte Schwankungen ihrer Ausgangssignale nicht überschreiten dürfen. So werden Schaltungen in Paaren oder Gruppen identischen Aufbaus verwendet, welche beispielsweise die gleiche Funktion gegenüber unterschiedlichen Signalen etc. übernehmen sollen. Jedoch weisen solche Schaltungen beziehungsweise Bauteile in integrierten Schaltungen auf Grund von statistischen Effekten Abweichungen in ihren elektrischen Parametern auf. Diese Variationen, der sogenannte Device Mismatch auf Grund der Schwankungen von Bauteilparametern, stellt eine wichtige Randbedingung bei der Dimensionierung entsprechender integrierter Schaltungen dar.

Der Device Mismatch kann für die verschiedenen Bauelemente in den jeweiligen Herstellungsprozessen integrierter Schaltungen meßtechnisch erfaßt werden und dient als Grundlage bei der Dimensionierung hochgenauer integrierter Schaltungen, beispielsweise von Baugruppen in Analog-Digital-und Digital- Analogwandlern. Es ist bekannt, daß der statistische Mismatch integrierter Schaltungen beziehungsweise Bauteile in ihren elektrischen Parametern mit zunehmender Fläche A abnimmt, wobei sich für die meisten Bauteile ein 1/X-Gesetz für die

Standardabweichung der relevanten elektrischen Parameter ableiten und experimentell bestätigen läßt. Neben diesen statischen Parameterabweichungen kommen auch dynamische Abweichungen der elektrischen Bauteilparameter vor, die als 1/f Rauschsignal eines Bauteils beobachtet werden können und die eine Abhängigkeit vom Verhältnis von den am Leitungsprozeß beteiligten Ladungsträgern und der Zahl der im Bereich des Stromflusses vorhandenen Störstellen zeigen.

Damit ist das 1/f-Rauschen ebenfalls von der aktiven Bauteilfläche abhängig und nimmt mit zunehmender Bauteilfläche ab.

Wird bei einer komplexen analogen Schaltung eine obere Grenze für parametervariationsbedingter Signalfehler am Ausgang der Schaltung gefordert, so müssen die Flächen der verwendeten Bauelemente entsprechend dimensioniert werden. Oblicherweise wird eine Schaltung unter Zuhilfenahme eines Simulationsprogramms entsprechend den Anforderungen entwickelt. Zur Berücksichtigung von Parametervariationseffekten werden in einer Monte-Carlo- Simulationsmethode die elektrischen Parameter der Einzelbauelemente voneinander statistisch unabhängig innerhalb einer vorzugebenden Schwankungsbreite geringfügig verändert, woraufhin erneut eine Simulation durchgeführt und die Abweichung des Ausgangssignals vom Sollwert ermittelt wird. Diese Methodik erweist sich als außerordentlich rechenintensiv, da zur Ermittlung eines einzelnen Punkts im Ergebnisraum die gesamte Schaltung simuliert werden muß. Um statistisch relevante Aussagen treffen zu können, ist eine große Anzahl von Simulationen notwendig. Von den so ermittelten, maximal zulässigen Schwankungsbreiten der elektrischen Parameter kann auf die Mindestfläche der einzelnen Bauteile geschlossen werden, welche sich aus der

Simulation ergibt und die zum Erreichen der gewünschten Gesamtgenauigkeit der Schaltung erforderlich ist.

Im folgenden sollen verschiedene Parameterabweichungen kurz vorgestellt werden.

Bei integrierten MOS-Transistoren hängt der Unterschied ihrer Schwellenspannungen im wesentlichen von der Anzahl der in den Kanal implantierten Dotieratome ab. Die Anzahl der Dotieratome nimmt linear mit der Fläche des Transistors zu, so daß sich beispielsweise für die Abweichung seiner Schwellenspannung Vt folgender Zusammenhang zur Bauteilfläche A herstellen läßt : <BR> <BR> <BR> <BR> <BR> <BR> l<BR> <BR> <BR> 'v' Auch die Steilheit eines MOS-Transistors unterliegt statistischen Schwankungen. Auch die Abweichung der Steilheit ist proportional zu 1/#A.

Der Flächenwiderstandswert eines Polysiliziumwiderstands wird im wesentlichen durch den Ladungsträgertransport über die Korngrenzen im Polysilizium bestimmt. Da die mittlere Korngröße und der mittlere Übergangswiderstand zwischen den Körnern einen festen, prozeßspezifischen Wert aufweist läßt sich für Abweichungen des Flächenwiderstandswerts folgende Gegebenheit aufstellen : <BR> <BR> <BR> <BR> <BR> ##R/R = 1 C#R/R<BR> <BR> <BR> #A

Auch für dynamische Abweichungen der elektrischen Bauteilparameter, die als 1/f-Rauschsignal eines Bauteils beobachtet werden können, ergibt sich eine Abhängigkeit vom Verhältnis von den am Leitungsprozeß beteiligten Ladungsträgern und der Zahl der im Bereich des Stromflusses vorhandenen Störstellen. Somit ist auch das 1/f-Rauschen von der aktiven Bauteilfläche abhängig und kann modellhaft beschrieben werden als r, l, i, is f-S, d KFF f. Coa A f r

Hierbei ist Vg, noise, pp das la-Rauschen der äquivalenten Gate-Spannung, KF die Rauschkonstante, Cox ist die Flächenkapazität des Gateoxids, A die Transistorfläche und fu und fl die obere-beziehungsweise untere Grenze des betrachteten Frequenzintervalls. Damit ergibt für die auf den Eingang des Transistors bezogene mittlere Rauschamplitude SUg folgender Zusammenhang : = A p-6v-Cvg. noise g LJ ., ir In allen obigen Formeln ist ein Proportionalitätsfaktor c angegeben worden, der von den betrachteten Bauelementen abhängt. Dieser Proportionalitätsfaktor wird als Matching- Konstante cij eines betrachteten Bauelements i bezeichnet, wobei j ein Index für den gerade betrachteten elektrischen Parameter ausgewählt aus einer Gruppe von n für das betrachtete Bauelement maßgeblichen oder als maßgeblich

betrachteten elektrischen Parameter, welche Variationen unterliegen können, ist.

Die experimentelle Ermittlung der Matching-Konstante wird durch Messung der elektrischen Eigenschaften zweier eng benachbarter Bauteile unter optimalen Layout-Bedingungen ermittelt. Hierzu zählen Kreuzkopplung, das Sicherstellen identischer Umgebungen für die zu matchenden Bauelemente, sowie die Vermeidung von Metallflächen über aktiven Gebieten.

Üblicherweise wird das Meßergebnis als Konstante für ein paarweises Vorkommen von Bauelementen angegeben. Tabelle 1 zeigt experimentell bestimmte Matching-Konstanten verschiedener Bauteile für unterschiedliche Parameter, wobei von einer paarweisen Auswertung ausgegangen wird. Lmin = 0, 5 um Lmin =0. 65 pm Bauteil Matching-Bauteil Matching- Konstante Konstante Poly-Si-R CAR/R = 11 %µm Poly-Si-1-R C#R/R = 12 % um n-Diffusion-R CAR/R = 1. 3 % um Poly-Si-2-R CAR/R = 5. 6 % um p-Diffusion-R CAR/R = 1. 4 %um Hi-Res-Poly-CAR/R = 3. 7 % pm Si-R n-MOSFET C#Vg,n = 16 mVm n-MOSFET C#Vg,n = 20 mVµm p-MOSFET C#Vg,p = 15 mVum P-MOSFET CAvg, p = 18 mVµm Lmin=0. 35 um Lmin = 0. 25 um Bauteil Matching-Bauteil Matching- Konstante Konstante n-MOSFET CAvg, n = 8 mVpm n-Poly-Si-R CAR/R = 4. 7 % um p-MOSFET cavg, p = 8 mVum p-Poly-Si-R C#R/R = 3. 0 % um n-MOSFET CAvg, n = 8 mV>m

Tabelle 1 : Experimentell ermittelte Matching-Konstanten von Widerständen und Transistoren in verschiedenen CMOS-Prozessen bei paarweiser Auswertung. Lmin ist die minimale Gatelänge eines Transistors, der mit dem jeweiligen Prozeß hergestellt werden kann.

Weiterhin sind Verfahren zur Flächenoptimierung von Bauelementen in integrierten Schaltungen unter Verwendung numerischer Verfahren aus (H. Su et al, Statistical constrained optimization of analog MOS circuits using empirical performance models, IEEE International Symposium on Circuits and Systems, Vol. 1, S. 133-136, 1994 ; C. J. Abel et al, Characterization of transistor mismatch for statistical CAD of submicron CMOS analog circuits, IEEE International Symposium on Circuits and Systems, S. 1401- 1404, 1993 ; M. J. M. Pelgrom, Matching properties of MOS transistors, IEEE Journal of Solid-State Circuits, Vol. 24, Nr. 5, S. 1433-1440, 1989) bekannt.

Somit besteht ein Bedarf, ein weniger iteratives Verfahren zur Bestimmung optimierter Bauelementflächen zu entwickeln.

Diese Aufgabe löst die Erfindung durch die Bereitstellung eines Verfahrens zur Flächenoptimierung von Bauelementen bei gepaarten Schaltungen gemäß dem unabhängigen Patentanspruch 1. Weitere vorteilhafte Ausgestaltungen, Details und Aspekte der Erfindung ergeben sich aus den abhängigen Patentansprüchen, der Beschreibung und den beigefügten Zeichnungen.

Der sogenannte Lagrangeformalismus erlaubt eine Extremwertbestimmung unter Nebenbedingungen.

Überraschenderweise hat sich gezeigt, daß dieses Lagrange- Gleichungssystem auf Grund des spezifischen Aufbaus der zur

Berechnung der Varianz verwendeten Formeln und der vorliegenden Bedingungen für die Gesamtflächenminimierung im vorliegenden Fall analytisch lösbar ist und zur Berechnung der Optima der Bauteilflächen die Verwendung einer einzelnen Gleichung gestattet.

Mit Hilfe des Langrange-Formalismus, der eine Extremwertbestimmung unter Nebenbedingungen erlaubt, kann das zuvor beschriebene Optimierungsproblem analytisch behandelt werden. Die Summe aller Bauteilflächen stellt die zu optimierende Größe dar, die minimal sein soll, um die Chipkosten zu senken. Jedoch dürfen die statistischen Schwankungen der Bauteilparameter nicht zu groß werden. Dies repräsentiert die Nebenbedingung in dem Optimierungsverfahren.

Die Erfindung ist daher gerichtet auf ein Verfahren zur Flächenoptimierung von Bauelementen in einer integrierten Schaltung mit m Bauelementen mit folgenden Schritten : -Festlegen einer vorgegebenen Wiederholgenauigkeit aout ; -Bestimmen der Mindestflächen Ai für Bauelemente jeder der Schaltungen mit einer Datenverarbeitungsanlage ; wobei die Mindestfläche Ai eines i-ten Bauelements jeder der gepaarten Schaltungen nach folgender Gleichung bestimmt wird : wobei cij eine Matching-Konstante für den j-ten von n elektrischen Parametern mit 1 < j < n des i-ten Bauelements

aus m Bauelementen mit 1 < i < m und aSout/apij die partielle Ableitung des Ausgangssginals Sout jeder der gepaarten Schaltungen nach dem elektrischen Parameter pij ist.

Unter einer Schaltung ist hier zumindest eine Schaltungsanordnung zu verstehen, welche bei gleichen Eingangssignalen möglichst identische Ausgangssignale liefern soll. Die Wiederholgenauigkeit oder Standardabweichung 6out ist hierbei die Standardabweichung des Ausgangssignals Sout der Schaltung, wenn man eine große Anzahl solcher identisch geplanter Schaltungen betrachtet. Die Schaltung weist m Bauelemente auf, wobei m zumindest 1 sein muß. Jedes der Bauelemente der Schaltung weist hierbei elektrische Parameter auf, welche Schwankungen unterliegen können. Von diesen elektrischen Parametern eines i-ten Bauelements werden für die Gleichung n Parameter berücksichtigt. Die Auswahl der berücksichtigten elektrischen Parameter hängt dabei von der gewünschten Genauigkeit der Berechnung der zur Verfügung stehenden Daten bezüglich der Matching-Konstanten und von Abschätzungen bezüglich der Einflußgröße der elektrischen Parameter eines i-ten Bauelements ab. Hierbei ist klar, daß nicht alle Bauelemente die gleiche Zahl von berücksichtigten elektrischen Parametern aufweisen müssen. Die partielle Ableitung des Ausgangssignals Sout kennzeichnet hierbei, wie stark sich bei einer kleinen Variation eines elektrischen Parameters das Ausgangssignal verändert.

Das erfindungsgemäße Verfahren muß auf Grund seiner Komplexität automatisiert durchgeführt werden, das heißt in einer Datenverarbeitungsanlage ausgeführt werden. Hierbei können wahlweise digitale Datenverarbeitungsanlagen oder auch entsprechend eingerichtete Analogrechner verwendet werden.

Die Datenverarbeitungsanlage kann Teil einer

Produktionsstätte für Halbleiter beziehungsweise einer Planungseinheit bei der Halbleiterherstellung sein.

Beispielsweise wird es bevorzugt, daß das erfindungsgemäße Verfahren in einem Gerät zur Erzeugung und/oder Planung von Belichtungsmasken für integrierte Halbleiterstrukturen durchgeführt wird.

Das Verfahren kann weiterhin den weiteren Schritt aufweisen : -Synthetisieren eines Schaltplanlayouts an Hand der Flächendaten der Bauelemente.

Hierbei versteht es sich, daß weitere Daten wie beispielsweise der Schaltplan etc. in das Synthetisieren eines Schaltplanlayouts einfließen müssen. Hierbei können entsprechende Simulationsprogramme, die um das erfindungsgemäße Verfahren ergänzt werden, beispielsweise SPICE-basierend, verwendet werden.

Verschiedene elektrische Parameter können Berücksichtigung finden. So kann beispielsweise für ein Bauelement, das ein MOS-Transistor ist, einer der betrachteten elektrischen Parameter die Schwellenspannung dieses Transistors sein. Bei Bauelementen, die integrierte Widerstände sind, kann einer ihrer elektrischen Parameter der Flächenwiderstandswert sein.

Bei Vorkommen eines MOS-Transistors kann einer der elektrischen Parameter das 1/f-Rauschsignal sein, ebenso wie bei einem integrierten Widerstand. Das erfindungsgemäße Verfahren funktioniert besonders gut, wenn die einzelnen Parameter Pij voneinander statistisch unabhängig sind.

Das erfindungsgemäße Verfahren kann auch hierarchisch verwendet werden. Ein auf einer bestimmten Hierarchiestufe betrachtetes Bauelement, welches in das erfindungsgemäße

Verfahren, wie oben angegeben einfließt, kann wiederum für sich betrachtet eine Schaltung mit mehreren Bauelementen sein. Auf Grund der mathematischen Gegebenheiten des erfindungsgemäßen Verfahrens ist es dann möglich, eine Matching-Konstante ctot für die Gesamtheit dieser Bauelemente (die eine Schaltung bilden), die in ihrer Gesamtheit auf der höheren Hierarchieebene wiederum als Bauelement betrachtet wird, anzugeben, welche sich berechnet nach In manchen Fällen führt eine konsequente Minimierung der Gesamtfläche der Schaltung zu eventuell unerwünschten Flächenverteilungen innerhalb der Schaltung. So kann beispielsweise ein differentielles Paar eine überhöhte Eingangskapazität aufweisen. Als zusätzliche einfache Nebenbedingung kann deshalb die Fläche eines oder mehrerer Bauelemente festgelegt werden und mit den verbleibenden Parametern die Optimierung durchgeführt werden. Die Bestimmung der Mindestflächen Ai der übrigen Bauelemente erfolgt daher so, daß die Bauelemente mit vorgegebener Fläche nicht mit einbezogen werden, wobei die Mindestfläche Ai eines i-ten Bauelementes jeder der gepaarten Schaltungen nach folgender Gleichung bestimmt wird : wobei sich cet ergibt aus der Subtraktion der Fehlerbeträge der festgelegten Bauelemente k+1 bis m zur Gesamtvarianz des

Ausgangssignals der Schaltung von der Gesamtvarianz des Ausgangssignals der Schaltung.

Häufig wird es vorkommen, daß lediglich ein elektrischer Parameter für ein bestimmtes Bauelement von Relevanz ist. In diesem Falle kann das Verfahren so durchgeführt werden, daß die Berechnung nach der vereinfachten Gleichung : erfolgt.

Im folgenden soll die Erfindung weiter erläutert und an konkreten Ausführungsbeispielen dargelegt werden, wobei auf die beigefügten Zeichnungen Bezug genommen werden soll, in denen folgendes dargestellt ist : Figur 1 zeigt die graphische Wiedergabe des Verfahrens zur Flächenminimierung mit zwei Optimierungparametern ; Figur 2 zeigt den Schaltplan eines beispielhaft betrachteten differentiellen Paars ; Figur 3 zeigt den optimierten Flächenbedarf des differentiellen Paars ; Figur 4 zeigt als Schaltplan eine Stromquelle ; Figur 5 zeigt die optimierten Flächen dieser Stromquelle.

Im folgenden soll eine mathematische Herleitung der im erfindungsgemäßen Verfahren verwendeten Gleichung dargestellt werden.

Eine beliebige analoge Schaltung liefert an ihrem Ausgang ein Ausgangssignal Soute das sich aus der sogenannten Übertragungsfunktion ergibt : Sout (p11, ...pmn) = f(Sin, p11,...,pmn) wobei, wie oben beschrieben, die Zeichen"p"jeweils n elektrische Parameter der m in der Schaltung verwendeten Bauelemente darstellen. Die konkrete Übertragungsfunktion für eine bestimmte Schaltung muß jeweils bestimmt werden.

Ziel des hier vorgestellten Verfahrens soll die eindeutige Bestimmung der Mindestflächen aller Einzelkomponenten einer komplexen Schaltung sein, so daß ein vorgegebener mismatch- bedingter Gesamtfehler nicht überschritten wird. Der geforderte maximale Gesamtfehler der Schaltung stellt die Nebenbedingung in einem Extremwertproblem dar, bei dem beispielsweise das Minimum der notwendigen Gesamtfläche bestimmt wird. Das Gesamtproblem gliedert sich somit in die Bestimmung einer Funktion zur Berechnung einer Gesamtfläche sowie einer Nebenbedingung, nämlich dem Einhalten einer bestimmten Schwankungsbreite (Standardabweichung) des Ausgangssignals. Bei fest vorgegebener Standardabweichung für das Ausgangssignal der Gesamtschaltung soll also deren Fläche minimiert werden. Die Funktion für die Bestimmung der Gesamtfläche lautet somit : Zur Festlegung der Nebenbedingung muß nunmehr die Berechnung der Standardabweichung bestimmt werden. Wenn die Parameter Pi ;, wobei i ein Index für das Bauteil ausgewählt aus m

Bauteilen und j ein Index für den Parameter, ausgewählt aus n möglichen Parametern ist, statistisch unabhängig sind, kann nach den Gesetzen der Fehlerfortpflanzung die Varianz 6out des Ausgangssignals Sout in einer linearen Näherung wie folgt abgeschätzt werden : Für den Fehler des Ausgangssignals folgt damit unter Berücksichtigung der oben beschriebenen Flächen-und Matching-Konstanten-abhängigen Variationen : Diese Beziehung ist der Ausgangspunkt für die weiteren Berechnungen, die eine Optimierung einer beliebigen analogen Schaltung hinsichtlich verschiedener flächenabhängiger Randbedingungen ermöglichen. Die Nebenbedingung (p (A1,..., Am) ergibt sich nunmehr aus Gleichung 3 zu : Nunmehr kann die Lagrangefunktion L dieses Extremwertproblems bestimmt werden. Der Lagrangeformalismus ist eine dem Fachmann geläufige Methodik zur Extremwertbestimmung unter Nebenbedingungen. Die sich ergebende Lagrangefunktion L ist :

wobei x der sogenannte Lagrangemultiplikator ist. Hat die untersuchte Funktion F unter der Nebenbedingung T einen relativen Extremwert im Punkt Ao, wobei der Punkt Ao ein Vektor mit m Dimensionen ist, der alle Mindestflächenwerte Al bis Am enthält, so gelten die Gleichungen <BR> <BR> <BR> <BR> aLAo=a.fAo)+a,a(Ao)=o<BR> <BR> aA ; aA ; aA ;<BR> <BR> <BR> <BR> <BR> G/t,Ct;&4, Setzt man die Gleichungen 4 und 1 in dieses Gleichungssystem ein, so folgt Dieses Gleichungssystem ist überraschenderweise analytisch lösbar auf Grund der in der Lagrangefunktion L vorkommenden Potenzen von Ai. Als Lösung dieses Gleichungssystems findet man daher die erfindungsgemäß verwendete Gleichung :

Berücksichtigt man, daß die Varianz einer Matching-sensitiven elektrischen Größe streng monoton mit abnehmender Fläche des Bauteils zunimmt, wie sich aus den Formeln zur Berechnung der Varianz verschiedener Parameter unmittelbar ergibt, so wird klar, daß es sich bei dem bestimmten Extremwert um ein Minimum handeln muß.

Figur 1 zeigt beispielhaft die Gesamtfläche eines Systems aus zwei Bauelementen unter Berücksichtigung von zwei Parametern.

Die Kurven zeigen hierbei die Abhängigkeit Atot als Fläche des Gesamtsystems von A1 als Fläche des ersten Bauelements im Hinblick auf eine vorgegebene Standardabweichung des Ausgabesignals. Es ist zu erkennen, daß eine hinsichtlich der Gesamtfläche optimale Flächenverteilung in diesem System existiert. Werden daher die Flächen der Einzelbauteile gemäß dem erfindungsgemäßen Verfahren bestimmt, ist die Gesamtfläche der Schaltung minimal. Für eine eindeutige Bestimmung der Flächen ist also ausschließlich die Kenntnis der Matching-Parameter und des Ubertragungsverhaltens der Schaltung notwendig.

Falls pro Bauelement lediglich ein Parameter berücksichtigt werden soll, entfällt die Summation über den Index j, so daß sich die erfindungsgemäß eingesetzte Gleichung vereinfacht zu :

In einer weiteren bevorzugten Ausführungsform der vorliegenden Erfindung werden nicht alle Bauelemente einer zu analysierenden Schaltung in die Berechnung mit einbezogen. In manchen Fällen führt nämlich eine konsequente Minimierung der Gesamtfläche der Schaltung über alle Bauelemente zu eventuell unerwünschten Flächenverteilungen innerhalb der Schaltung. So kann beispielsweise ein differentielles Paar eine überhöhte Eingangskapzität aufweisen. Als zusätzliche, weitere einfache Nebenbedingung kann deshalb die Fläche eines oder mehrerer Bauteile vorab festgelegt und mit den verbleibenden Parametern die Optimierung durchgeführt werden. Es gilt somit für den Gesamtfehler einer entsprechenden Schaltung : /\ 2 /'\ 2 2 VST °"'VV . "'7 n-, i i=l j=l apu Ai i=k+l j=l apij Ai ouf durchNebenbed. festgelegt Hier wird unterschieden zwischen Bauelementen 1 bis k, deren Fläche optimiert werden soll, und Bauelementen k+1 bis m, deren Fläche durch die obige Nebenbedingung festgelegt ist, wobei m die Zahl aller Bauelemente in der Schaltung darstellt. Diese in ihrer Fläche festgelegten Bauteile leisten einen Grundbeitrag zur vorgegebenen Gesamtvarianz des Ausgangssignals der Schaltung, so daß für den Fehlerbetrag der zu optimierenden Bauteile eine Restvarianz definiert werden kann : Ebenso müssen die Flächen der festgelegten Bauteile derart gewählt werden, daß gilt :

Die Optimierung der restlichen Bauteilflächen führt dann zu einer modifizierten Kalkulationsvorschrift. Es ergibt sich die Gleichung : Sollen die optimalen Flächen rein analytisch bestimmt werden, so muß zunächst die Übertragungsfunktion der untersuchten Schaltung berechnet werden. Für große Schaltungen kann dies mathematisch sehr aufwendig werden oder nicht mehr geschlossen lösbar sein. Zur Optimierung komplexer Schaltungen muß deshalb hierarchisch vorgegangen werden.

Betrachtet man die erfindungsgemäß minimierte Gesamtfläche der Schaltung : so kann man eine neue Matching-Konstante ctot definieren, die das Matching-Verhalten des gesamten Funktionsblocks beschreibt :

Mit Hilfe dieser Konstanten kann erneut eine Beziehung zwischen der gesamten aktiven Fläche der Schaltung und den Matching-Anforderungen in der bekannten Form hergestellt werden : Es ist dadurch erfindungsgemäß möglich, eine komplexe Schaltung hierarchisch hinsichtlich ihrer Matching- Eigenschaften und ihres Gesamtflächenbedarfs zu analysieren.

Das Matching-Verhalten der Gesamtschaltung kann auf die Eigenschaften kleinerer Funktionsblöcke zurückgeführt werden, was die Berechnung der Übertragungsfunktionen erleichtert und so in sehr großen Schaltungen eine Matching-Analyse erst ermöglicht. Ferner können mit Hilfe des erfindungsgemäßen Verfahrens Minima für die aktive Fläche der einzelnen Funktionsblöcke bestimmt werden, die als Grundlage für Gesamtflächenabschätzungen integrierter Schaltungen dienen können.

Auch bei der numerischen Flächenoptimierung mit Hilfe von Empfindlichkeitsanalysen kann hierarchisch vorgegangen werden. So lassen sich auch sehr komplexe Systeme konsequent flächenoptimal dimensionieren.

Im Folgenden sollen zur Veranschaulichung des erfindungsgemäßen Verfahrens zwei analytische Beispiele analoger Schaltungen vorgestellt werden.

Zur Berechnung der optimierten Bauteilflächen muß zunächst die Übertragungsfunktion bestimmt werden. Für Bauteile, die mehrere Arbeitsgebiete aufweisen, wie beispielsweise MOS- Transistoren, muß hierbei entschieden werden, welche mathematische Beschreibung ihrer Ubertragsungsfunktion verwendet werden kann. Dies ergibt sich aus den erlaubten Betriebszuständen der Gesamtschaltung sowie aus Überlegungen hinsichtlich der Auswirkung der mismatch-bedingten Parametervariationen auf das Ausgangssignal. Ein Mismatch in der Schwellenspannung eines Transistors beispielsweise wirkt sich bei kleiner Gate-Spannung stärker auf den Ausgangsstrom aus als bei großer Gate-Spannung. Bei der analytischen Betrachtung kann durch eine worst-case-Betrachtung jedoch sichergestellt werden, daß die Schaltung in jedem erlaubten Betriebspunkt die Matching-Anforderungen, die der Optimierung zugrundegelegt werden, erfüllt. Im Falle einer numerischen Berechnung der Empfindlichkeiten (sensitivity analysis) wird diese Aufgabe durch das Schaltungsmodell und den Schaltungssimulator gelöst.

Beispiel 1 Im folgenden soll das erfindungsgemäße Verfahren der Gesamtflächenoptimierung auf ein differentielles Paar, wie es in Fig. 2 dargestellt ist, angewendet werden. Das Paar besteht aus zwei Transistoren M1 und M2 sowie zwei Widerständen R1 und R2, welche gegenüber einer Spannung Vo parallel geschaltet sind. Zwischen das differentielle Paar und die Erdung ist eine Stromquelle geschaltet, die den Strom Io liefert. Bei der Berechnung der Übertragungsfunktion wird davon ausgegangen, daß die Transistoren Ml und M2 im

Sättigungsgebiet arbeiten. Dann gilt für die Kanalströme Ids1 und Ids2: Ids1 = ß1/2(Vin1 - V0 - Vt1*)2 und [19] <BR> <BR> <BR> Ids2 = ß2/2(Vin2 - V0 - Vt2*)² [20]<BR> Id52 = ß2/2(Vin2 - V0 - Vt2*)2 [20] Die Randbedingungen für dieses System lauten : IQ = Idsl + ids2, ß = ßl = ß2, Vin = Vin2 = Vin1, Vout = R1* Ids1 - R2* Ids2 Die vom Mismatch betroffenen elektrischen Bauteilparameter sind in den obigen Gleichungen mit einem * versehen. Der Mismatch der Steilheiten ßi wird nicht berücksichtigt, da von einem effektiven Gatespannungsmismatch ausgegangen wird.

Die Übertragungsfunktion kann unter Zuhilfenahme eines symbolischen Rechenprogramms wie der kommerziell erhältlichen Programme"Maple"TM oder"Mathematica"TM berechnet werden.

Nach Gleichung (10) ergibt sich für die optimierten Flächen :

wobei cR und cn-MOS die prozeßspezifischen Matching- Konstanten für Widerstände und n-MOS-Transistoren sind. Aus den AC-und DC-Spezifikationen der Schaltung ergeben sich Io, und R, so daß bei geforderter Gesamtgenauigkeit 6out die Flächen aller Bauteile eindeutig bestimmt werden können.

Seien Io = 100 uA, ß = 50 uA/V und R = 20kQ. Die für einen Prozeß experimentell bestimmten Matching-Konstanten können für die Schwellenspannung von n-MOS-TransistorenC,=8.0 mVum und für den Widerstandswert von Polysilizium- Widerständen CRair 4. 8 % um sein.

Die gewünschte, nicht zu überschreitende Standardabweichung der Ausgangsspannung 6out kann beispielsweise 10mV betragen.

Für die Mindestflächen der Bauteile ergeben sich dann folgende Werte : AR1 = AR2 = 28. 5 um2 und AM, = AM2 = 6. 71 um.

Dies ist die bezüglich der aktiven Gesamtfläche optimale Flächenkonfiguration des differentiellen Paars unter den obigen Bedingungen.

In Fig. 3 ist der Gesamtflächenbedarf der Schaltung in Abhängigkeit vom Wert der Matching-Konstanten der verwendeten Widerstände dargestellt. Die Werte für den DC-Arbeitspunkt und die Matching-Konstante für die Transistoren wurden aus obigem Beispiel übernommen. Es ist deutlich erkennbar, daß die optimale Flächenkonfiguration stark von der Matching- Konstanten der verwendeten Widerstände abhängt.

Beispiel 2 Die Genauigkeit von CMOS-Stromspiegeln hängt im wesentlichen vom Matching der strombestimmenden Transistoren in der Stromsenke, beziehungsweise der Stromquelle, ab. In diesem Beispiel soll am Beispiel einer Stromquelle mit hohem Ausgangswiderstand (regulated current mirror) eine Flächenoptimierung vorgenommen werden.

Fig. 4 zeigt das Schaltbild der im Beispiel verwendeten Stromquelle. Strombestimmender Transistor ist M1, während M2, M3 und M4 dazu dienen, das Source-Drain-Potential am Transistor M1 konstant zu halten : M4 stellt eine Konstantstromquelle dar, die den Biasstrom durch Transistor M3 bestimmt. Über den Regeltransistor M2 stellt sich deshalb am Gate von M3 eine konstante Spannung ein.

Ein Stromspiegel besteht aus zwei entsprechend geschalteten derartigen Stromquellen. Seine Genauigkeit hängt von den Matching-Eigenschaften der Einzeltransistoren ab. Bei der Berechnung der Übertragungsfunktion wird davon ausgegangen, daß die Transistoren M3 und M4 in Sättigung sind, M1 im linearen Gebiet betrieben wird und M2 keinen Beitrag zu einer matching-bedingten Ungenauigkeit der Stromquelle liefert. Ein Fehler in der Schwellenspannung von M2 wird durch die hohe Verstärkung des Regelkreises bestehend aus M2, M3 und M4 weitgehend kompensiert. Die Ströme der Schaltung lassen sich wie folgt charakterisieren :

Ids3 = ß3/2(Vgs - Vt3*)2 [24] und <BR> <BR> <BR> <BR> <BR> und<BR> Ids4 = ß4/2(Vbias - Vt4*)2 [25] Die vom Mismatch betroffenen elektrischen Bauteilparameter sind in diesen Gleichungen mit einem * versehen.

Die einzige Nebenbedingung zur Berechnung der Übertragungsfunktion lautet : Ids3 = Ids4- Für die Beispielrechnung werden Vbias = 0. 80 V, Vtn-MOS = 0. 45 V, Vtp-MOS = 0. 50 V, Vin = 2. 0 V, al = 100 µA/V, ß3 = 1 mA/V und ß4 = 200 8 eingesetzt. Die Abweichung im Ausgangsstrom zweier benachbarter Stromquellen soll dabei kleiner als Grel = 0. 002 sein. Dann ergeben sich bei einer prozeßbedingten Matching-Konstanten von cn-MOSpair = Cp-MOSpair = 8.0 mVµm für die Mindestflächen der Transistoren folgende Werte : AM1 = 17. 2 um2, AM3 = 28. 4 um2 und AM4 = 12. 7 um2.

Die Berechnung bestätigt, daß neben der Fläche des strombestimmenden Transistors My auch die Transistoren M3 und M4 wesentliche Beiträge zum Gesamtmatching zweier derartiger Stromquellen liefern.

In Fig. 5 sind die optimalen Flächen der Transistoren für verschiedene Matching-Konstanten cp°'MOS aufgetragen. Hier ist

deutlich die unterschiedliche Gewichtung der Einzeltransistoren zu erkennen.

Die Berechnung der Übertragungsfunktion der zu optimierenden Schaltung erweist sich besonders bei größeren Funktionsblöcken als aufwendig. In vielen Fällen kann das Ergebnis nicht mehr geschlossen angegeben werden.

Die Berechnung der Übertragungsfunktion ist nicht notwendig, wenn eine Empfindlichkeitsanalyse mit Hilfe des Schaltungssimulators durchgeführt wird. Viele Simulationsprogramme (beispielsweise SPICE) bieten die Option einer sensitivity analysis, bei der die Empfindlichkeiten der elektrischen Modellparameter auf die angegebenen Ausgangssignale berechnet und ausgegeben werden können.

Die erhaltenen Daten repräsentieren direkt die Ableitungen aSout/apij der Übertragungsfunktion nach den statistisch abweichenden Parametern im zuvor beschriebenen Formalismus.

Auf diese Weise können bei Kenntnis der Matching-Konstanten einfach optimale Flächenverteilungen ohne explizite Berechnung der Übertragungsfunktion erhalten werden.

Das erfindungsgemäße Verfahren erlaubt es, mismatch-sensitive Schaltungen analytisch, also ohne die Durchführung von rechen-und zeitintensiven Monte-Carlo-Simulationen zu dimensionieren. Die Mindestflächen der Einzelbauteile können derart bestimmt werden, daß entweder die Gesamtfläche minimiert oder andere flächenabhängige Schaltungseigenschaften optimiert werden. Die Ergebnisse stellen zudem, im Gegensatz zur Monte-Carlo-Simulation, nicht einen mit statistischen Fehlern behafteten günstigsten Fall, sondern das mathematisch gefundene eindeutige absolute

Minimum dar. Die Möglichkeit der hierarchischen Berechnung erlaubt es, auch sehr komplexe Schaltungen konsequent zu optimieren und Gesamtflächenabschätzungen vorzunehmen, wenn die Bauteilgröße im wesentlichen durch Mismatch- Randbedingungen festgelegt wird.

In der Entwicklungsphase analoger Schaltungen bieten sich Vorteile, da alleine mit Kenntnis der Übertragungsfunktion des Schaltungsblocks und den Matching-Parametern des jeweiligen Prozesses typische Bauteildimensionen und die Flächenverteilung innerhalb einer Schaltung berechnet werden können. Die gleichzeitige Berücksichtigung von Matching-und Rausch-Gesichtspunkten bei der Optimierung ist ebenso möglich.

Unterstützt das im design flow verwendete Simulationsprogramm Empfindlichkeitsanalysen, so können die optimalen Flächen einfach numerisch aus den simulierten Daten errechnet werden.

Da auch andere Schaltungseigenschaften, beispielsweise auf Grund veränderter Kapazitäten, von den Bauteilflächen abhängen, müssen die erhaltenen Flächen in das Schematic,. die Netzliste oder das Layout zurückgeschrieben und eine neue Simulation vorgenommen werden. Dies kann auch automatisiert werden.