CN107239594B - Analog circuit dispersibility optimization method based on PSPICE - Google Patents

Analog circuit dispersibility optimization method based on PSPICE Download PDF

Info

Publication number
CN107239594B
CN107239594B CN201710282178.3A CN201710282178A CN107239594B CN 107239594 B CN107239594 B CN 107239594B CN 201710282178 A CN201710282178 A CN 201710282178A CN 107239594 B CN107239594 B CN 107239594B
Authority
CN
China
Prior art keywords
simulation
key
circuit
round
value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710282178.3A
Other languages
Chinese (zh)
Other versions
CN107239594A (en
Inventor
周桂法
杜绍华
汪旭
陈旭鸿
匡芬
潘宇雄
唐欢
袁莹莹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CRRC Zhuzhou Institute Co Ltd
Original Assignee
CRRC Zhuzhou Institute Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by CRRC Zhuzhou Institute Co Ltd filed Critical CRRC Zhuzhou Institute Co Ltd
Priority to CN201710282178.3A priority Critical patent/CN107239594B/en
Publication of CN107239594A publication Critical patent/CN107239594A/en
Application granted granted Critical
Publication of CN107239594B publication Critical patent/CN107239594B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Design And Manufacture Of Integrated Circuits (AREA)

Abstract

The invention discloses a PSPICE-based analog circuit dispersibility optimization method, which comprises the following steps: s1, constructing a simulation circuit and constructing a transfer function of a key characteristic index; s2, carrying out sensitivity analysis and determining key devices; s3, performing first-round simulation; s4, determining a function between the key characteristic index and the key device; s5, performing robustness optimization design on the function, and determining a mean value, a standard upper limit value and a standard lower limit value; s6, optimizing parameters of key devices, performing a second round of simulation, and determining the yield; s7, when the yield meets the preset judgment index, skipping to the step S8; otherwise, jumping to step S2; and S8, analyzing the simulation circuit, judging whether a device under the action of overstress exists, modifying a derating factor of the device if the device under the action of overstress exists, repeating the step S8, and otherwise, obtaining design parameters to finish the design of the simulation circuit. The invention has the advantage of improving the stability and reliability of the optimized design circuit.

Description

Analog circuit dispersibility optimization method based on PSPICE
Technical Field
The invention relates to the field of analog circuit optimization design, in particular to a PSPICE-based analog circuit dispersity optimization method.
Background
The traditional circuit optimization design and analysis method comprises a tolerance design method, a Monte Carlo simulation method and the like, but the traditional method focuses on local problems of circuit device tolerance, circuit manufacturing qualification rate and the like, and global and systematic solutions are lacked for the dispersion problems of circuit key characteristics (such as performance indexes are out of tolerance, drift, instability and the like). Therefore, it is of practical significance to study the dispersion problem of the key characteristics of the circuit.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: aiming at the technical problems in the prior art, the invention provides the PSPICE-based analog circuit dispersity optimization method which can analyze the dispersity problem of the key characteristics of the circuit from the perspective of the whole situation and the system and optimize and design the circuit, thereby improving the dispersity capability of the circuit for resisting elements and improving the stability and the reliability of the circuit.
In order to solve the technical problems, the technical scheme provided by the invention is as follows: a PSPICE-based analog circuit dispersibility optimization method comprises the following steps:
s1, constructing a simulation circuit and constructing a transfer function of a key characteristic index of the simulation circuit;
s2, carrying out sensitivity analysis on the simulation circuit, and determining key devices influencing the key characteristic indexes;
s3, performing a first round of simulation to obtain a first round of simulation value of the key characteristic index of the simulation circuit, wherein the first round of simulation value comprises an input value and an output value of the key characteristic index;
s4, determining a function between the key characteristic index and the key device;
s5, performing robustness optimization design on the function, and determining the mean value, the standard upper limit value and the standard lower limit value of the key characteristic index;
s6, optimizing parameters of the key device according to the mean value, performing a second round of simulation to obtain a second round of simulation value of key characteristic indexes of the simulation circuit, and determining the yield of the simulation circuit according to the second round of simulation value;
s7, when the yield meets a preset judgment index, skipping to the step S8; otherwise, jumping to step S2;
s8, Smok analysis is carried out on the simulation circuit, whether a device under the action of overstress exists in the simulation circuit is judged according to the Smok analysis result, if yes, a derating factor of the device is modified, the step S8 is repeated, and if not, design parameters of the simulation circuit are obtained, and simulation circuit design is completed.
As a further improvement of the present invention, the step of constructing the simulation circuit in step S1 includes: and selecting a device from the reference model library to construct a simulation circuit by taking the PSPICE AA model library as the reference model library, and supplementing tolerance information of the device.
As a further improvement of the present invention, the specific steps of determining the key device affecting the key characteristic index in step S2 are as follows: and sequencing the devices according to the sensitivity analysis result by using the influence degree of the key index, and taking the device with the influence degree in a preset sequencing range as a key device.
As a further improvement of the present invention, the preset ordering range is a range with influence degree of the first 5 to 10 bits.
As a further improvement of the invention, the number of times of the first round simulation is not less than 100.
As a further improvement of the present invention, the specific steps of step S4 include: selecting N groups of samples from the first round of simulation values, and determining a function between the key characteristic index and the key device by a response surface method, a random response method or a chaotic polynomial method, as shown in formula (1):
Y=F(X) (1)
in the formula (1), Y is an output value of a key characteristic index of the simulation circuit, f (X) is the function, and X ═ X1,X2…Xk],XiAnd i is 1,2, …, and k is the ith key device.
As a further improvement of the invention, the value of N is greater than or equal to 20.
As a further improvement of the present invention, the specific steps of step S5 include: constructing a robustness optimization model as shown in formula (2):
Figure BDA0001279791490000021
in the formula (2), σYIs the standard deviation, mu, of the key property index YYIs the mean value of the key characteristic index Y, YLIs the lower limit of the specification for Y, YUIs the upper limit of the specification for Y,
Figure BDA0001279791490000022
mean of the ith critical device, D.V is the design variable;
determining a mean μ of the key characteristic indicators according to the robustness optimization modelYUpper limit of specification YULower limit of specification YL
As a further improvement of the present invention, when the step S6 optimizes the parameters of the critical device according to the mean value, the variance of the critical device is guaranteed to be unchanged.
As a further improvement of the invention, the number of times of the second round simulation is not less than 1000, and the yield of the simulation circuit is the probability that the average value of the simulation values of the second round falls into the range of [ lower limit value of specification, upper limit value of specification ].
Compared with the prior art, the invention has the advantages that:
1. the invention provides a solution for the dispersibility of key characteristic indexes of a circuit, which aims at the problems of overlarge dispersibility of key characteristic indexes of the circuit, parameter drift and the like, carries out sensitivity analysis and identifies which components have the greatest contribution to the dispersibility of the circuit; then, carrying out robustness optimization design, and on the premise that the circuit key characteristic index value meets the established technical specification, the dispersity is minimum and the performance is optimal; carrying out yield analysis, and calculating the yield of the manufactured circuit on the premise of the tolerance of the existing device, so that the dispersion level of the product can be mastered in the research and development stage; and finally, Smoke analysis is carried out, derating design is carried out on devices with overstress, the capability of the circuit for resisting the dispersity of the devices is improved, and the stability and the reliability of the circuit are improved.
Drawings
FIG. 1 is a schematic flow chart of an embodiment of the present invention.
Detailed Description
The invention is further described below with reference to the drawings and specific preferred embodiments of the description, without thereby limiting the scope of protection of the invention.
As shown in fig. 1, the PSPICE-based analog circuit dispersibility optimization method of this embodiment includes the steps of: s1, constructing a simulation circuit and constructing a transfer function of a key characteristic index of the simulation circuit; s2, carrying out sensitivity analysis on the simulation circuit, and determining key devices influencing the key characteristic indexes; s3, performing a first round of simulation to obtain a first round of simulation value of the key characteristic index of the simulation circuit, wherein the first round of simulation value comprises an input value and an output value of the key characteristic index; s4, determining a function between the key characteristic index and the key device; s5, performing robustness optimization design on the function, and determining the mean value, the standard upper limit value and the standard lower limit value of the key characteristic index; s6, optimizing parameters of the key device according to the mean value, performing a second round of simulation to obtain a second round of simulation value of key characteristic indexes of the simulation circuit, and determining the yield of the simulation circuit according to the second round of simulation value; s7, when the yield meets a preset judgment index, skipping to the step S8; otherwise, jumping to step S2; s8, Smok analysis is carried out on the simulation circuit, whether a device under the action of overstress exists in the simulation circuit is judged according to the Smok analysis result, if yes, a derating factor of the device is modified, the step S8 is repeated, and if not, design parameters of the simulation circuit are obtained, and simulation circuit design is completed.
In this embodiment, the step of building the simulation circuit in step S1 includes: and selecting a device from the reference model library to construct a simulation circuit by taking the PSPICE AA model library as the reference model library, and supplementing tolerance information of the device according to a device manual. After the simulation circuit is constructed, a transfer function of key characteristic indexes including but not limited to indexes such as gain, bandwidth and delay time is created in a probe window. By constructing a transfer function of the key characteristic index, the value of the key characteristic index can be obtained during each simulation.
In step S2 in this embodiment, sensitivity analysis is performed on the constructed simulation circuit, and according to the sensitivity analysis result, the devices are sorted from high to low according to the influence of the critical index, and the device with the influence in a preset sorting range is used as a critical device. The predetermined ordering range is a range of the influence degree of the first 5 to 10 bits.
In this embodiment, a first round of simulation is performed on the simulation circuit by a monte carlo simulation method, and the number of times of the first round of simulation is not less than 100. Each time Monte Carlo simulation is carried out, input values and output values of a group of key characteristic indexes can be obtained, a group of samples are obtained, and 100 groups of samples can be obtained after 100 Monte Carlo simulations are carried out.
In this embodiment, N sets of samples are selected from 100 sets of samples of the first round of simulated values, where the value of N is greater than or equal to 20, and a function between the key characteristic indicator and the key device is determined by a response surface method, a random response method, or a chaotic polynomial method, as shown in equation (1):
Y=F(X) (1)
in the formula (1), Y is an output value of a key characteristic index of the simulation circuit, f (X) is the function, and X ═ X1,X2…Xk],XiAnd i is 1,2, …, and k is the ith key device.
In this embodiment, the specific step of step S5 includes: constructing a robustness optimization model as shown in formula (2):
Figure BDA0001279791490000041
in the formula (2), σYIs the standard deviation, mu, of the key property index YYIs the mean value of the key characteristic index Y, YLIs the lower limit of the specification for Y, YUIs the upper limit of the specification for Y,
Figure BDA0001279791490000042
D.V design variables (design variables) for the mean of the ith critical device; determining a mean μ of the key characteristic indicators according to the robustness optimization modelYUpper limit of specification YULower limit of specification YL. Through the constraint condition of the robustness optimization model, the lower limit value Y of the specification can be determinedLAnd a specification upper limit value YURange of composition [ YL,YU]And determining the mean of the critical devices. And optimizing parameters of the key devices in the simulation circuit by using the mean value of the key devices, keeping the variance of the key devices unchanged, and updating the simulation circuit. And performing Monte Carlo simulation again, namely performing simulation of the second round, wherein the number of times of the simulation of the second round is not less than 1000. And acquiring a second round simulation value of the key characteristic indexes of the simulation circuit, wherein the second round simulation value comprises the input value and the output value of the key characteristic indexes, and calculating the mean value of the second round simulation value. The mean value of 1000 second round simulation values is calculated to fall into the lower limit value and the upper limit value of the specification]And obtaining the yield of the simulation circuit according to the probability of the range. And when the genuine product rate of the simulation circuit meets the preset qualification judgment standard, performing the next processing on the simulation circuit, otherwise, skipping to the step S2 to perform sensitivity analysis on the simulation circuit again, and entering a new round of optimization design process.
In this embodiment, when the yield of the simulation circuit meets a preset qualification criterion, a cookie (cookie analysis is a proprietary analysis tool in software, and is "smoking" analysis as the name implies, and is used for researching whether a circuit device can bear excessive heat and electric stress, and if the answer is exceeded, a derating design of response is carried out) analysis is carried out on the simulation circuit, whether a device subjected to an over-stress effect exists in the simulation circuit is judged according to a cookie analysis result, if so, a derating factor of the device is modified, and a targeted derating process is carried out to improve the capability of the circuit for resisting the dispersion of components, and step S8 is repeated, otherwise, design parameters of the simulation circuit are obtained, the design parameters are output, and the simulation circuit design is completed.
The embodiment provides a solution for the dispersibility of the key characteristic indexes of the circuit, and aims at the problems of overlarge dispersibility of the key characteristic indexes of the circuit, parameter drift and the like, sensitivity analysis is carried out to identify which components have the greatest contribution to the dispersibility of the circuit; then, carrying out robustness optimization design, and on the premise that the circuit key characteristic index value meets the established technical specification, the dispersity is minimum and the performance is optimal; by carrying out yield analysis and calculating the yield of the circuit after manufacturing on the premise of the tolerance of the existing device, the dispersion level of the product can be mastered in the research and development stage; and finally, Smoke analysis is carried out, derating design is carried out on devices with overstress, the capability of the circuit for resisting the dispersity of the devices is improved, and the stability and the reliability of the circuit are improved.
The foregoing is considered as illustrative of the preferred embodiments of the invention and is not to be construed as limiting the invention in any way. Although the present invention has been described with reference to the preferred embodiments, it is not intended to be limited thereto. Therefore, any simple modification, equivalent change and modification made to the above embodiments according to the technical spirit of the present invention should fall within the protection scope of the technical scheme of the present invention, unless the technical spirit of the present invention departs from the content of the technical scheme of the present invention.

Claims (10)

1. A PSPICE-based analog circuit dispersibility optimization method is characterized by comprising the following steps:
s1, constructing a simulation circuit and constructing a transfer function of a key characteristic index of the simulation circuit;
s2, carrying out sensitivity analysis on the simulation circuit, and determining key devices influencing the key characteristic indexes;
s3, performing a first round of simulation to obtain a first round of simulation value of the key characteristic index of the simulation circuit, wherein the first round of simulation value comprises an input value and an output value of the key characteristic index;
s4, determining a function between the key characteristic index and the key device;
s5, performing robustness optimization design on the function, and determining the mean value, the standard upper limit value and the standard lower limit value of the key characteristic index;
s6, optimizing parameters of the key device according to the mean value, performing a second round of simulation to obtain a second round of simulation value of key characteristic indexes of the simulation circuit, and determining the yield of the simulation circuit according to the second round of simulation value;
s7, when the yield meets a preset judgment index, skipping to the step S8; otherwise, jumping to step S2;
s8, Smok analysis is carried out on the simulation circuit, whether a device under the over-stress action exists in the simulation circuit is judged according to the Smok analysis result, if yes, a derating factor of the device is modified, the step S8 is repeated, and if not, design parameters of the simulation circuit are obtained, and simulation circuit design is completed;
the key characteristic indexes include gain, bandwidth and delay time.
2. The PSPICE-based analog circuit dispersibility optimization method according to claim 1, wherein the step of constructing the simulation circuit in step S1 includes: and selecting a device from the reference model library to construct a simulation circuit by taking the PSPICE AA model library as the reference model library, and supplementing tolerance information of the device.
3. The PSPICE-based analog circuit dispersibility optimization method according to claim 2, wherein the specific steps of determining the key devices affecting the key characteristic index in the step S2 are as follows: and sequencing the devices according to the sensitivity analysis result by using the influence degree of the key index, and taking the device with the influence degree in a preset sequencing range as a key device.
4. The PSPICE-based analog circuit dispersibility optimization method of claim 3, wherein: the preset sequencing range is a range with influence degree of the first 5-10 bits.
5. The PSPICE-based analog circuit dispersibility optimization method of claim 4, wherein: the number of times of the first round of simulation is not less than 100.
6. The PSPICE-based analog circuit dispersibility optimization method according to claim 3, 4 or 5, wherein the specific steps of the step S4 include: selecting N groups of samples from the first round of simulation values, and determining a function between the key characteristic index and the key device by a response surface method, a random response method or a chaotic polynomial method, as shown in formula (1):
Y=F(X) (1)
in the formula (1), Y is an output value of a key characteristic index of the simulation circuit, f (X) is the function, and X ═ X1,X2…Xk],XiAnd i is 1,2, …, and k is the ith key device.
7. The PSPICE-based analog circuit dispersibility optimization method of claim 6, wherein the value of N is greater than or equal to 20.
8. The PSPICE-based analog circuit dispersibility optimization method according to claim 7, wherein the specific steps of the step S5 include: constructing a robustness optimization model as shown in formula (2):
Figure FDA0002557837720000021
in the formula (2), σYIs the standard deviation, mu, of the key property index YYIs the mean value of the key characteristic index Y, YLIs the lower limit of the specification for Y, YUIs the upper limit of the specification for Y,
Figure FDA0002557837720000022
mean of the ith critical device, D.V is the design variable;
determining a mean μ of the key characteristic indicators according to the robustness optimization modelYUpper limit of specification YULower limit of specification YL
9. The PSPICE-based analog circuit dispersibility optimization method of claim 8, wherein: and step S6, when the parameters of the key device are optimized according to the mean value, the variance of the key device is ensured to be unchanged.
10. The PSPICE-based analog circuit dispersibility optimization method of claim 9, wherein: the number of times of the second round of simulation is not less than 1000, and the yield of the simulation circuit is the probability that the mean value of the simulation values of the second round falls into the range of the lower limit value and the upper limit value of the specification.
CN201710282178.3A 2017-04-26 2017-04-26 Analog circuit dispersibility optimization method based on PSPICE Active CN107239594B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710282178.3A CN107239594B (en) 2017-04-26 2017-04-26 Analog circuit dispersibility optimization method based on PSPICE

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710282178.3A CN107239594B (en) 2017-04-26 2017-04-26 Analog circuit dispersibility optimization method based on PSPICE

Publications (2)

Publication Number Publication Date
CN107239594A CN107239594A (en) 2017-10-10
CN107239594B true CN107239594B (en) 2020-11-27

Family

ID=59984253

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710282178.3A Active CN107239594B (en) 2017-04-26 2017-04-26 Analog circuit dispersibility optimization method based on PSPICE

Country Status (1)

Country Link
CN (1) CN107239594B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112541318A (en) * 2020-12-10 2021-03-23 成都博思微科技有限公司 Optimal design method for improving isolation of packaging substrate
CN115048885B (en) * 2022-08-12 2022-11-15 阿里巴巴(中国)有限公司 Circuit design parameter adjusting method and device, electronic equipment and storage medium

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1868131A (en) * 2003-12-17 2006-11-22 株式会社电装 Activation signal output circuit

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10198548B2 (en) * 2008-02-21 2019-02-05 Mentor Graphics Corporation Identifying the defective layer of a yield excursion through the statistical analysis of scan diagnosis results
US8766807B2 (en) * 2008-10-03 2014-07-01 Universal Security Instruments, Inc. Dynamic alarm sensitivity adjustment and auto-calibrating smoke detection
CN102314522B (en) * 2010-06-30 2013-03-06 中国科学院微电子研究所 Design optimization method for analog integrated circuit
CN103955579B (en) * 2014-04-28 2017-06-23 天津大学仁爱学院 Simulation based on SPICE softwares/RF IC method for designing
CN106529090B (en) * 2016-12-12 2019-06-14 中国航天标准化研究所 A kind of aerospace electron class Reliability Assessment method

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1868131A (en) * 2003-12-17 2006-11-22 株式会社电装 Activation signal output circuit

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"基于关键元器件退化的模拟电路行为分析及模型建立";陈诚;《中国优秀硕士学位论文全文数据库 信息科技辑》;20150215(第02期);I135-244 *

Also Published As

Publication number Publication date
CN107239594A (en) 2017-10-10

Similar Documents

Publication Publication Date Title
US10324989B2 (en) Microblog-based event context acquiring method and system
CN107169628B (en) Power distribution network reliability assessment method based on big data mutual information attribute reduction
CN110311902A (en) A kind of recognition methods of abnormal behaviour, device and electronic equipment
Veech et al. Diversity partitioning without statistical independence of alpha and beta
CN110175168B (en) Time sequence data filling method and system based on generation of countermeasure network
CN107093022B (en) Small sample evaluation method for reliability of high-confidence explosive product
CN107239594B (en) Analog circuit dispersibility optimization method based on PSPICE
CN114722973B (en) Defect detection method and system for steel pipe heat treatment
Teyabeen Statistical analysis of wind speed data
CN114548898A (en) Cloud platform based government affair big data management method and system
CN109212751B (en) Method for analyzing tolerance of free-form surface
CN116527296A (en) Multi-level network security situation assessment method
CN110502736B (en) Bus data analysis method and device
CN110969303B (en) Tree height prediction method based on rational Charles model
CN109446481B (en) Method for estimating life distribution parameters of log-normal unit
CN113641990A (en) Intrusion detection method based on multi-innovation extended Kalman filtering
CN109284320A (en) Automatic returning diagnostic method in big data platform
CN113918817B (en) Push model construction method, push model construction device, computer equipment and storage medium
CN112434839B (en) Distribution transformer heavy overload risk prediction method and electronic equipment
CN113742720B (en) Network security situation perception method based on multistage linkage mode
CN110110287A (en) A kind of probability box Sensitivity Analysis Method
CN113450924B (en) Novel coronavirus propagation model establishing method and system
Zhu et al. Security assessment of power data link based on dynamic weight change
CN116612824A (en) Activated carbon data processing method and device, electronic equipment and medium
Czabak-Górska et al. Statistical Performance Of X ̅ And R Control Charts For Skewed Distribution-Case Study

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant