WO1999027425A1 - Multiple scale signal processing and control system - Google Patents

Multiple scale signal processing and control system Download PDF

Info

Publication number
WO1999027425A1
WO1999027425A1 PCT/US1998/023889 US9823889W WO9927425A1 WO 1999027425 A1 WO1999027425 A1 WO 1999027425A1 US 9823889 W US9823889 W US 9823889W WO 9927425 A1 WO9927425 A1 WO 9927425A1
Authority
WO
WIPO (PCT)
Prior art keywords
dynamic model
physical system
parameters
creating
states
Prior art date
Application number
PCT/US1998/023889
Other languages
French (fr)
Inventor
Sunil C. Shah
Pradeep Pandey
Mark A. Erickson
Original Assignee
Voyan Technology
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 Voyan Technology filed Critical Voyan Technology
Priority to AU14539/99A priority Critical patent/AU1453999A/en
Priority to DE69828995T priority patent/DE69828995T2/en
Priority to AT98958508T priority patent/ATE289083T1/en
Priority to EP98958508A priority patent/EP1034459B1/en
Publication of WO1999027425A1 publication Critical patent/WO1999027425A1/en

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B17/00Systems involving the use of models or simulators of said systems
    • G05B17/02Systems involving the use of models or simulators of said systems electric

Definitions

  • the present invention relates to the field of signal processing and control systems. More particularly, the invention relates to the use of multiple time scale and /or space scale models for signal processing and control.
  • Signal processing and the control of physical systems generally involves obtaining measurements from a physical system in the form of electrical signals and processing the signals in order to bring about a desired result.
  • the control of a physical system typically involves obtaining measurements from the physical system, comparing the measurements with a predetermined control recipe, and adjusting the system inputs in response to the comparison to minimize variations between the measured values and recipe values.
  • the signals to be processed or the variables to be controlled which ever the case may be, are not always directly available for observation and must be inferred from indirect and noisy measurements.
  • the indirect measurements are generally obtained from embedded sensors which contain multiple pieces of information that are dynamically confounded. Extraction of the information of interest requires the use of complex and time consuming calibration procedures and the use of estimating techniques that result in high computation costs. Equipment setup costs are also high since diagnostic measurements must be taken to correlate the measured signal to the indirect measurements for each piece of equipment.
  • measurements from a physical system are not always obtainable at a single time interval (time scale). For example, there may be a first measurement that is obtainable only at a first time scale, a second measurement that is only obtainable at a second time scale, a third measurement that is only obtainable at a third time scale, and so on.
  • a measurement, or a set of measurements, taken from a physical system may not bear the same spatial relationship with the system as other measurements taken from the system. Even in instances where all measurements are available at a single high rate, the computation cost of indirect measurements at the high rate can be very high.
  • the task of controlling a system involves not only the control of a single physical system, but the control of a family of physical systems. This situation is most prevalent in high volume manufacturing applications.
  • the characteristics of a single physical system tend to change over time due to equipment degradation and other causes.
  • the characteristics between a family of physical systems tend to differ from one physical system to another due to equipment-to-equipment variations. It is important to track these changes and to account for these differences so that the signal processing or control system may be updated accordingly. Otherwise, the accuracy of the signal processing or control system is compromised.
  • One object of the present invention is to provide an accurate and cost efficient signal processing or control system.
  • Another object of the present invention is to provide a method for rapidly calibrating one or more signal processing systems.
  • Another object of the present invention is to optimize the performance of physical systems in which measurements are be obtained at different time scales.
  • Another object of the present invention is to optimize the performance of physical systems in which measurements are be obtained at different space scales.
  • Another object of the present invention is to accurately characterize the behavior of a physical system.
  • a further object of the present invention is to develop multiple models that characterize the behavior of the physical system at different time scales and to employ online multi- variable control of the physical system based on the multiple time scale models.
  • Another object of the present invention is to develop multiple models that characterize the behavior of the physical system at different space scales and to employ online multi-variable control of the physical system based on the multiple space scale models.
  • Another object of the present invention is to control a family of physical systems in which measurements are obtained at different time scales.
  • Another object of the present invention is to control a family of physical systems in which measurements are obtained at different space scales. 9/27425 _
  • a further object of the present invention is to decrease equipment setup time.
  • Another object of the present invention is to use diagnostic sensors to improve estimation quality in multi-scale systems.
  • signals from the physical system are processed by: a) creating a first dynamic model at the first time scale that consists of a first set of parameters, a first set of states and a first set of inputs; b) creating a second dynamic model at the second time scale that consists of a second set of parameters, a second set of states and a second set of inputs, at least one of the first set of parameters in the first dynamic model being computed from the second set of states in the second dynamic model; c) creating a second estimator to produce estimates of the second set of states using the second dynamic model, measurements from the first physical system at the second time scale; and d) creating a first estimator to produce estimates of .the first set of states using the first dynamic model, measurements from the first physical system at the first time scale, and the estimates from the second estimator.
  • Figure 1 illustrates an evaporator deposition system
  • Figure 2 illustrates a control system for controlling the evaporator deposition system of Figure 1;
  • Figure 3 illustrates a process for designing a dynamic model in one embodiment of the present invention
  • Figure 4 illustrates the application of power test signals to the system while a synchronous response of the system is gathered
  • Figure 5 is a flow chart of a method for processing signals from a physical system in accordance with one embodiment of the present invention
  • Figure 6 is a flow chart of a method for controlling a family physical systems
  • Figure 7A illustrates a side view of a conventional wafer processing furnace
  • Figure 7B illustrates a fine scale model of the furnace shown in Figure 7A
  • Figure 7C illustrates a coarse scale model of the furnace shown in Figure 7A.
  • Figure 1 illustrates a deposition system 100 comprising a deposition chamber 102 containing a substrate 120 and one or more evaporation sources 104.
  • the evaporation sources include raw materials 106 that are evaporated to produce a film 122 on a surface of substrate 120.
  • Each evaporation source 104 is disposed in thermal exposure to a heating source 108 to evaporate the raw materials to be deposited onto substrate 120.
  • Thermocouples 110 are used to monitor the temperature of the /27425 _ _
  • deposition system 100 may be used in a variety of applications to produce particular film compositions on the surface of substrate 120.
  • the deposition system is used to produce photovoltaic cells via a process of thin film deposition.
  • multiple raw materials are evaporated, deposited, and reacted in the deposition chamber 102 to produce a desired film composition 122 on the surface of substrate 120.
  • the surface temperature of raw material pools 106 are controlled by measuring the temperature of the raw materials at the bottom of the pools with thermocouples 110 at a first time scale (time interval), given by index k, and regulating the power to heaters 108 in response to the temperature measurements.
  • substrate 120 is removed from chamber 100 and the composition of film 122 is measured at a second time scale, given by index j, to determine the product quality.
  • the sample rate of the first time scale given by index k is faster than the sample rate of the second time scale given by index j.
  • control system 200 for controlling the deposition system of Figure 1 is illustrated.
  • the control system uses measurements from both the first and second time scales to create a controller 210 that optimizes the performance of the deposition system in order to produce a desired film composition 122 on the surface of substrate 120 in accordance with a predetermined process recipe or setpoint.
  • Control system 200 includes a first dynamic model 202 that indicates the behavior of the deposition system 100 at the first time scale index k, and a second dynamic model 204 that indicates the behavior of the deposition system at the second time scale index j.
  • the first dynamic model 202 is represented by the following set of model equations:
  • xl is the thermal and chemical state of the system
  • ul is the power input to evaporator heaters 108
  • l are measurements from the system
  • pn represents fixed system parameters
  • pi 2 represents the source evaporation rate parameters for each of the raw materials 106.
  • the evaporation rate parameters, p ⁇ 2 may change over time due to physical changes that are occurring within the evaporation chamber. Since the pi 2 parameters are subject to change, it is important that these parameters be updated or refined to ensure that the first dynamic model properly indicates the behavior of the system.
  • the pn parameters may include, for example, the source-to-substrate distance, chamber volume, etc..
  • the yl measurements consist of two sets of measurements; yll measurements, and yl2 measurements. The yll measurements are obtained from the evaporator thermocouples 110 (commonly known as production sensors) during the deposition process and represent the temperature of the bottom of raw material pools 106 during the 9/27425 _
  • the yl2 measurements are obtained from diagnostic sensors that are not available during the production process.
  • the second dynamic model 204 is represented by the following set of model equations:
  • x2 is the state of the deposition system at the second time scale
  • «2 are temperature setpoint offsets
  • p 2 are fixed parameters of the system in the second time scale
  • y2 is the end-of-run film composition measurements obtained at the end of the deposition process.
  • equation (1) and/or (3) in models 202 and 204 may represent continuous-time dynamic models.
  • the relationship between the pi2 parameters of the first dynamic model and the states of the second dynamic model (x2) may be characterized in the following manner.
  • the pi2 parameters are chosen to represent evaporation rate offsets which contribute to variations in the end- of-run film compositions.
  • ⁇ i2 is represented by the following equation:
  • the second dynamic model is created using the following set of equations:
  • ⁇ y is a stochastic drift process
  • y2j is the measured evaporation offset
  • zt2j is the temperature setpoint offsets (°C)
  • .and y . is the measurement noise (A).
  • Figures 3 and 4 illustrates a process for designing the first dynamic model 202.
  • a set of heater power test signals 408 are created to exercise high bandwidth and low bandwidth dynamics on the system.
  • Test signals 408 are then applied to heaters 108 by a test signal generator 402 to cause the raw materials 106 to evaporate.
  • Diagnostic sensors 406, in the form of a Quartz Crystal Microbalances (QCMs) are positioned within chamber 102 to measure the evaporation rate of the raw materials 106.
  • a response data acquisition unit 404 obtains measurement signals 410 and 412 from the diagnostic sensors 406 and thermocouples 110, respectively, in response to the test signals 408. The data acquired during this procedure is used to produce the first dynamic model 202.
  • the first dynamic model 202 indicates the behavior of the system between the high bandwidth and low bandwidth dynamics at the first time scale index k.
  • the first dynamic model is represented by both yll measurements and yl2 measurements.
  • the yll measurements are obtained from thermocouples 110.
  • the yl2 measurements (evaporation rates of the raw materials 106), however, are not directly available for observation during the deposition process so that they must be /27425
  • thermocouples 110 The yll measurements obtained from thermocouples 110 are also noisy. For this reason, a first estimator 206 is created to produce estimates of the states and measurements of the first model and to predict the raw material evaporation rates (yl2 measurements) from the heater powers and thermocouple signals. The end- of-run film composition measurements (y2) are also noisy and for this reason a second estimator 208 is created to provide estimates of the states x2 and measurement inputs y2 of the second dynamic model.
  • the second estimator 208 is represented by the following set of equations:
  • estimator 208 is determined using all yllk measurements obtained at the first time scale in the last period of the second time scale as shown in equation (10) below. In such an instance, estimator 208 is represented by equations (8b) and (9b) below.
  • yn j is represented by:
  • n is the ratio of sample rate 1 to sample rate 2.
  • the end-of-run film composition measurements acquired at the second time scale are used to refine the pn parameters of the first dynamic model in the first time scale.
  • the first estimator 206 is, thus, created using the following set of equations:
  • Construction of the first and second estimators 206 and 208 from models 202 and 204, respectively, are off-line activities.
  • a controller 210 for controlling the deposition system 100 at the first time scale is then created using the estimates from the first estimator 206 and the second estimator 208.
  • the controller 210 is written as: /27425
  • w represents a predetermined process recipe or setpoint that will result in a desired film composition being deposited on the surface of substrate 120.
  • the first and second estimators 206 and 208 may be created based on a variety of filtering techniques.
  • Al, Bl, Cl, and Dl represent the linearization of i and gi in equations (1) and (2) and where Al is of dimension nsi x nsi, Bl is of dimension ns x n ⁇ , Cl is of dimension ny_ x ns ⁇ , .and Dl is of dimension nyi xnu ⁇ .
  • Al is of dimension ns ⁇ is the number of states
  • n i is the number of outputs
  • the output measurement yl is divided into the measurable and the non-measurable parts as follows:
  • the vector yllk corresponds to the measured temperatures from thermocouples 110 and is of dimension ny ⁇ x 1.
  • the output matrix C is partitioned as
  • multiple linearized models are obtained each with index i for different steady-state operating pints and the gain matrices L/ are obtained by solving steady-state Riccati equations using reduced order ARX model parameters as disclosed in column 8, line 22 through column 9, line 22 of United States patent no. 5,517,594.
  • United States patent no. 5,517,594 entitled “Thermal Reactor Optimization” is hereby incorporated by reference.
  • -17- estimators 206 and 208 may be created using a continuous-time extended Kalman filter (CTEKF) for state estimations in cases where equations (1) and /or (3) are continuous-time dynamic models, and in some instances may also be based on the diagnostic measurements, yl2.
  • the estimators may also be created using other filtering techniques such as a discrete-time extended Kalman filter (DTEKF) for state estimation.
  • DTEKF discrete-time extended Kalman filter
  • the first estimator 206 is based on an extended Kalman filter (EKF) that uses discrete time measurements and integration. This method is described in co-pending United States patent application entitled "A Method for Real-Time Nonlinear System State Estimation and Control", filed by Sunil C. Shah and Pradeep Pandey on August 22, 1997 and hereby incorporated by reference.
  • the discussion has thus been limited to controlling an evaporator deposition system, it is to be understood that the present invention is not limited to a specific type of system to be controlled, nor is it limited to a specific type of control system. In addition, it is appreciated that the present invention is not limited to a method of controlling a physical system that is characterized by only two time scales. Moreover, it is important to note that the scope of the present invention is not limited to control systems, but is also applicable to the field of signal processing in general.
  • signals generated from any physical system may be processed pursuant to the method outlined in the flow chart of Figure 5.
  • a first dynamic model is created that indicates a behavior of the physical system at a first time scale.
  • a second dynamic model is created that indicates a behavior of the physical system at a second time scale.
  • a second estimator is created to estimate the states of the second dynamic model using the second dynamic model and measurements from the physical system at the second time scale.
  • a first estimator is created to estimate the states of the first dynamic model using the first dynamic model, measurements from the physical system at the first time scale, and estimates from the second estimator.
  • Equations (1) through (4), (8), (9) and (11) through (13) generally represent the manner in which the dynamic models, estimators and controller are constructed in one embodiment of the invention.
  • the deposition system is defined by a first set of parameters (fast parameters) that is changing at a fast rate and a second set of parameters (slow parameters) that is changing at a slow rate.
  • the number of calculations required per time period to solve the combined set of parameters at a single fast time scale is represented by k(ns ⁇ +ns 2 ) ⁇ , where nsi is the number of states of the first model and ns 2 is the number of states in the second model.
  • each deposition chamber within the family of chambers is manufactured so that all the chambers are identical.
  • chambers-to-chamber variations will exist.
  • the thermal characteristics of each of the chambers within the family may vary due to differences in the heater elements, insulation properties, emissivity variations, etc..
  • FIG. 6 is a flow chart that outlines a method for controlling a family of systems consisting of a first system and a second system.
  • a first dynamic model is created at a first time scale that /27425
  • -20- indicates the behavior of the first system at the first time scale. This typically involves exciting the first system and capturing the response of the first system in a manner similar to that described above in Figures 3 and 4.
  • the first dynamic model is represented by equations (1) and (2) above, wherein the p!2 parameters capture equipment-to-equipment variations between the first system and the second system.
  • a second dynamic model is created at a second time scale that indicates the behavior of the of the first system at the second time scale.
  • the second dynamic model is represented by equations (3) and (4) above.
  • steps 604, 606, and 608 a first estimator, second estimator, and first controller are created in a manner similar to that described in the example above. Accordingly, equations (8) through (12) may be used to represent the first estimator, second estimator, and first controller for the first system.
  • a second controller for optimizing the performance of the second system is then created by applying the first dynamic model to the second system to create a third dynamic model.
  • the third dynamic model indicates the behavior of the second system at the first time scale and is represented by the following set of model equations:
  • x3 is the state of the second system
  • u3 are inputs to the second system
  • y3 are measurements from the second system
  • pn represents fixed parameters
  • p32 represents refutable parameters, including those
  • a fourth dynamic model is created to indicate the behavior of the second system at the second time scale and is represented by the following set of model equations:
  • x4 is the state of the second system at the second time scale
  • u4 inputs to the second system p 4 are fixed parameters of the second system at the second time scale
  • y4 are measurements of the second system at the second time scale.
  • a fourth estimator is created to estimate the states of the fourth dynamic model.
  • the fourth estimator is created using the fourth dynamic model and measurements from the second system at the second time scale. In some instances, measurements from the second system at the first time scale are additionally used to create the fourth estimator.
  • the • fourth estimator may then be represented by the following set of equations.
  • a third estimator is created to estimate the states of the third dynamic model.
  • the third estimator is created using the third dynamic model, measurements from the second system at the first time scale and the , scale modifier, modifier modifier, modifier modifier, modifier modifier, modifier modifier, modifier PCT/US98/2388 /27425
  • a second controller for controlling the second system at the first time scale is created using estimates from the third and fourth estimators.
  • the second controller maybe written as:
  • the first dynamic model of the first system is applied to the second system in order to create the third dynamic model after it is initially created.
  • the first dynamic model is first refined by implementing steps 602 through 610 before it is applied to the second system to create the third dynamic model.
  • the ability to apply the first dynamic model of the first system to other systems alleviates the need to obtain diagnostic measurements from these systems during setup. As a result, the setup time and the costs associated with setting up each system is therefore reduced.
  • a magnetic disk drive consist of a stack of circular disks having a thin magnetic coating rotating at high speed.
  • Data is recorded on the disk surfaces using transducers mounted on read-write heads.
  • the read-write heads are attached to actuator arms that are moved across the disk surfaces by a high-speed actuators that are under the control of a servo motor.
  • the actual information is recorded in circumferential tracks on the disk surfaces. Reading the recorded information back involves sensing the magnetic transitions emanating from the surfaces of the magnetic disks, again using the transducers.
  • the ability to store and read back information from each of the magnetic disks is affected by vibrations in the arm assemblies which affect the position of the read-write head transducers within the circumferential tracks. Magnetic modulation of the thin magnetic coatings will also affect the signal that is produced by the read-write head transducers. The read-write head transducer signals will thus consist of confounded information that is representative by both the vibration and magnetic modulation.
  • Variations also exists between each of the read-write heads. These variations may be in the form of different read-write head transducer sensitivities, differences in the physical characteristics of the read-write heads, nominal tolerance differences, etc..
  • the read-write signals may be processed by first creating a first model at a first time scale that indicates and takes into account the vibrational and magnetic modulation components of the magnetic disk drive.
  • the first model may be represented by the following set of model equations:
  • xl is the position .and magnetic states of the system
  • wl is servo motor commands
  • yl are measurements from taken from magnetic disk drive unit
  • pil represents fixed system parameters
  • pi2 represents head mass, offset and film sensitivity parameters.
  • the pi2 parameters are subject to change and/or may vary from one magnetic disk/head assembly to another. Therefore, it is important that these parameters be updated or refined to ensure that the first model properly indicates the behavior of the disk drive unit.
  • the pn parameters may include dynamics parameters, such as drag and inertia parameters.
  • the yl measurements consist of two sets of measurements; yll measurements, and l2 measurements.
  • the yll measurements include servo motor drive current, the read-write head sensor signal, actuator arm position, and servo motor velocity.
  • the yl2 measurements may include read-write head positioning measurements that are obtained with laser positioning diagnostic sensors.
  • a second model at a second time scale is also created which indicates how of each of the read-write heads, servo motor parameters, structural vibration parameters, friction and drag parameters change over time and from component to component.
  • the second model also indicates how measurements of previously recorded data changes with drive and component parameters x2.
  • the second model equations may be written as:
  • x2 comprise head mass and transducer sensitivity states of the system
  • 2 is the head position setpoint offsets
  • p 2 are fixed parameters of the magnetic disk drive in the second time scale
  • y2 is the measurement of previously recorded data.
  • read-write signal processing is achieved by decomposing the magnetic disk drive into two separate models and updating the models in accordance with the method of Figure 5.
  • Figure 7A exemplifies such a system.
  • a furnace 702 for processing a plurality of wafers 704a-704d is shown.
  • the furnace 702 includes a quartz jar 706 having a set of heating elements 708 that surround the quartz jar.
  • a furnace controller may be created to control the temperature of the wafers by first determining a first dynamic model at a fine, space scale given by index 1, that consists of a first set of parameters, a first set of states, and a first set of inputs.
  • the first dynamic model at the fine space scale may be written as:
  • xl comprises wafer 704d temperatures
  • ul is the heater powers
  • yl are temperature measurements taken from the system
  • pn is the thermal dynamics of the system
  • pi2 represents the center and edge temperatures of wafer 704d.
  • the yl measurements consist of two sets of measurements; yll measurements, and yl2 measurements.
  • the yll measurements are obtained from production thermocouples within the furnace, whereas the yl2 measurements are obtained from diagnostic sensors that are not available during the production process.
  • a second dynamic model at a coarse space scale given by index m is also created and consists of a second set of parameters, a second set of states, and a second set of inputs.
  • the second dynamic model may be represented by the following equations:
  • x2 is the center .and edge temperatures of wafers 704a-704d
  • u2 is the heater powers
  • p 2 represents the dynamics of the system at the coarse space scale
  • y2 represents furnace production sensor measurements.
  • a second estimator of the coarse model is created using the following set of equations:
  • x2 represents estimates of the wafer center and edge temperatures and y2 represents estimates of the production thermocouple measurements.
  • the center and edge wafer temperatures at the coarse space scale are used to refine the pi2 parameters of the first dynamic model in the fine space scale.
  • a first estimator of the fine model is, thus, created using equations (40) and (41) below.
  • a controller for controlling the furnace 702 at the fine space scale is created using the estimates from the first estimator and the second estimator.
  • the controller is written as:
  • w represents a predetermined process recipe or setpoint.

Abstract

A method for processing signals and controlling a physical system in which measurements are obtained at different time scales (502, 504) and/or different space scales. Signals generated from the physical system are processed by first creating a first dynamic model at a first time/space scale (502). A second dynamic model at a second time/space scale (504) is also created. A second estimator is then created to produce estimates of the second set of states using the second dynamic model, measurements from the first physical system at the second time scale (506). A first estimator is also created to produce estimates of the first set of states using the first dynamic model, measurements from the first physical system at the first time scale, and the estimates from the second estimator (508).

Description

MU TIPLE SCALE SIGNAL PROCESSING AND CONTROL SYSTEM
FTFLD OF THE INVENTION
The present invention relates to the field of signal processing and control systems. More particularly, the invention relates to the use of multiple time scale and /or space scale models for signal processing and control.
BACKGROUND OF THE INVENTION
Signal processing and the control of physical systems generally involves obtaining measurements from a physical system in the form of electrical signals and processing the signals in order to bring about a desired result. For example, the control of a physical system typically involves obtaining measurements from the physical system, comparing the measurements with a predetermined control recipe, and adjusting the system inputs in response to the comparison to minimize variations between the measured values and recipe values. During signal processing and control, the signals to be processed or the variables to be controlled, which ever the case may be, are not always directly available for observation and must be inferred from indirect and noisy measurements. The indirect measurements are generally obtained from embedded sensors which contain multiple pieces of information that are dynamically confounded. Extraction of the information of interest requires the use of complex and time consuming calibration procedures and the use of estimating techniques that result in high computation costs. Equipment setup costs are also high since diagnostic measurements must be taken to correlate the measured signal to the indirect measurements for each piece of equipment. 9/27425
In addition, measurements from a physical system are not always obtainable at a single time interval (time scale). For example, there may be a first measurement that is obtainable only at a first time scale, a second measurement that is only obtainable at a second time scale, a third measurement that is only obtainable at a third time scale, and so on. In other instances, a measurement, or a set of measurements, taken from a physical system may not bear the same spatial relationship with the system as other measurements taken from the system. Even in instances where all measurements are available at a single high rate, the computation cost of indirect measurements at the high rate can be very high.
Often, the task of controlling a system involves not only the control of a single physical system, but the control of a family of physical systems. This situation is most prevalent in high volume manufacturing applications. The characteristics of a single physical system tend to change over time due to equipment degradation and other causes. Moreover, the characteristics between a family of physical systems tend to differ from one physical system to another due to equipment-to-equipment variations. It is important to track these changes and to account for these differences so that the signal processing or control system may be updated accordingly. Otherwise, the accuracy of the signal processing or control system is compromised.
Current signal processing and control systems do not provide for the rapid calibration of such systems, nor do they have the ability to track changes in the characteristics of a single physical system or a family of physical systems.
Thus, what is needed is an accurate and cost efficient method of processin signals generated within a physical system, or a family of physical systems, in which measurements are obtained at different time scales and /or different space scales.
/27425 ._
SUMMARY AND OBTECTS OF THE INVENTION
One object of the present invention is to provide an accurate and cost efficient signal processing or control system.
Another object of the present invention is to provide a method for rapidly calibrating one or more signal processing systems.
Another object of the present invention is to optimize the performance of physical systems in which measurements are be obtained at different time scales.
Another object of the present invention is to optimize the performance of physical systems in which measurements are be obtained at different space scales.
Another object of the present invention is to accurately characterize the behavior of a physical system.
A further object of the present invention is to develop multiple models that characterize the behavior of the physical system at different time scales and to employ online multi- variable control of the physical system based on the multiple time scale models.
Another object of the present invention is to develop multiple models that characterize the behavior of the physical system at different space scales and to employ online multi-variable control of the physical system based on the multiple space scale models.
Another object of the present invention is to control a family of physical systems in which measurements are obtained at different time scales.
Another object of the present invention is to control a family of physical systems in which measurements are obtained at different space scales. 9/27425 _
-5-
A further object of the present invention is to decrease equipment setup time.
Another object of the present invention is to use diagnostic sensors to improve estimation quality in multi-scale systems.
In one embodiment, these and other objects of the invention are provided by a method for processing signals generated by a physical system in which measurements are available in a first time scale and a second time scale. In accordance with the present invention, signals from the physical system are processed by: a) creating a first dynamic model at the first time scale that consists of a first set of parameters, a first set of states and a first set of inputs; b) creating a second dynamic model at the second time scale that consists of a second set of parameters, a second set of states and a second set of inputs, at least one of the first set of parameters in the first dynamic model being computed from the second set of states in the second dynamic model; c) creating a second estimator to produce estimates of the second set of states using the second dynamic model, measurements from the first physical system at the second time scale; and d) creating a first estimator to produce estimates of .the first set of states using the first dynamic model, measurements from the first physical system at the first time scale, and the estimates from the second estimator.
Other objects and features of the present invention will be apparent from the accompanying drawings, and from the detailed description that follows below. BRIEF DESCRIPTION OF THE DRAWINGS
The present invention is illustrated by way of example and not limitation in the figures of the accompanying drawings in which like references indicate similar elements, and in which:
Figure 1 illustrates an evaporator deposition system;
Figure 2 illustrates a control system for controlling the evaporator deposition system of Figure 1;
Figure 3 illustrates a process for designing a dynamic model in one embodiment of the present invention
Figure 4 illustrates the application of power test signals to the system while a synchronous response of the system is gathered;
Figure 5 is a flow chart of a method for processing signals from a physical system in accordance with one embodiment of the present invention;
Figure 6 is a flow chart of a method for controlling a family physical systems;
Figure 7A illustrates a side view of a conventional wafer processing furnace; Figure 7B illustrates a fine scale model of the furnace shown in Figure 7A; and
Figure 7C illustrates a coarse scale model of the furnace shown in Figure 7A.
/27425
DETAILED DESCRIPTION
A method for processing signals and controlling a physical system in which measurements from the system are available at two or more time scales, or alternatively, from two or more space scales, is described. In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present invention. It will be obvious, however, to one skilled in the art that these specific details need not be employed to practice the present invention. In other instances, well known methods have not been described in detail in order to avoid unnecessarily obscuring the present invention. Additionally, it should be noted that although the present invention is described in relation to an evaporator deposition system, one with ordinary skill in the art will recognize that such a description is merely illustrative and is not intended to limit the invention. The specific processes and system described herein are only meant to help clarify one's understanding of the present invention and to illustrate particular embodiments in which the present invention may be implemented. It will be appreciated that the broader spirit and scope of the present invention, as set forth in the appended claims, may be applied to any type of process which seeks the achievements attained by the present invention.
Referring to the drawings, Figure 1 illustrates a deposition system 100 comprising a deposition chamber 102 containing a substrate 120 and one or more evaporation sources 104. The evaporation sources include raw materials 106 that are evaporated to produce a film 122 on a surface of substrate 120. Each evaporation source 104 is disposed in thermal exposure to a heating source 108 to evaporate the raw materials to be deposited onto substrate 120. Thermocouples 110 are used to monitor the temperature of the /27425 _ _
raw materials during the evaporation process. Since the evaporation rate of each raw material is directly related to its temperature, it is important that the pool surface temperature of the raw materials be precisely controlled at the surface of the liquid pool in order to produce a desired film composition on the surface of substrate 120.
It is appreciated that deposition system 100 may be used in a variety of applications to produce particular film compositions on the surface of substrate 120. For example, in one application the deposition system is used to produce photovoltaic cells via a process of thin film deposition. In this process multiple raw materials are evaporated, deposited, and reacted in the deposition chamber 102 to produce a desired film composition 122 on the surface of substrate 120. During the deposition process, the surface temperature of raw material pools 106 are controlled by measuring the temperature of the raw materials at the bottom of the pools with thermocouples 110 at a first time scale (time interval), given by index k, and regulating the power to heaters 108 in response to the temperature measurements.
Once the deposition process is complete, substrate 120 is removed from chamber 100 and the composition of film 122 is measured at a second time scale, given by index j, to determine the product quality. The sample rate of the first time scale given by index k, is faster than the sample rate of the second time scale given by index j.
Turning now to Figure 2, a control system 200 for controlling the deposition system of Figure 1 is illustrated. The control system uses measurements from both the first and second time scales to create a controller 210 that optimizes the performance of the deposition system in order to produce a desired film composition 122 on the surface of substrate 120 in accordance with a predetermined process recipe or setpoint.
Control system 200 includes a first dynamic model 202 that indicates the behavior of the deposition system 100 at the first time scale index k, and a second dynamic model 204 that indicates the behavior of the deposition system at the second time scale index j. The first dynamic model 202 is represented by the following set of model equations:
XL+, = fl(xl_fUh> P >P J. (Eq. 1)
Figure imgf000012_0001
where xl is the thermal and chemical state of the system, ul is the power input to evaporator heaters 108, l are measurements from the system, pn represents fixed system parameters, and pi2 represents the source evaporation rate parameters for each of the raw materials 106. The evaporation rate parameters, pι2, may change over time due to physical changes that are occurring within the evaporation chamber. Since the pi2 parameters are subject to change, it is important that these parameters be updated or refined to ensure that the first dynamic model properly indicates the behavior of the system. The pn parameters may include, for example, the source-to-substrate distance, chamber volume, etc.. The yl measurements consist of two sets of measurements; yll measurements, and yl2 measurements. The yll measurements are obtained from the evaporator thermocouples 110 (commonly known as production sensors) during the deposition process and represent the temperature of the bottom of raw material pools 106 during the 9/27425 _
deposition process. The yl2 measurements are obtained from diagnostic sensors that are not available during the production process.
The second dynamic model 204 is represented by the following set of model equations:
Figure imgf000013_0001
where x2 is the state of the deposition system at the second time scale, «2 are temperature setpoint offsets, p2 are fixed parameters of the system in the second time scale, and y2 is the end-of-run film composition measurements obtained at the end of the deposition process. The pi2 parameters of the first dynamic model are related to the states x2 of the second dynamic model. That is, ρi2 = h(*2).
In alternative embodiments equation (1) and/or (3) in models 202 and 204 may represent continuous-time dynamic models.
In one instance, the relationship between the pi2 parameters of the first dynamic model and the states of the second dynamic model (x2) may be characterized in the following manner. The pi2 parameters are chosen to represent evaporation rate offsets which contribute to variations in the end- of-run film compositions. In such an instance, ρi2 is represented by the following equation:
Figure imgf000013_0002
where d is the true evaporation rate offset (A) and is the true thermal sensitivity (A/°C). To produce estimates of the pi2 parameters at the second 9/27425 _12_
time scale, the second dynamic model is created using the following set of equations:
Figure imgf000014_0001
y2. = dj + jU2j + Vj, (Eq. 7)
where μyis a stochastic drift process, y2j is the measured
Figure imgf000014_0002
evaporation offset, zt2j is the temperature setpoint offsets (°C), .and y . is the measurement noise (A).
Figures 3 and 4 illustrates a process for designing the first dynamic model 202. At step 302, a set of heater power test signals 408 are created to exercise high bandwidth and low bandwidth dynamics on the system. Test signals 408 are then applied to heaters 108 by a test signal generator 402 to cause the raw materials 106 to evaporate. Diagnostic sensors 406, in the form of a Quartz Crystal Microbalances (QCMs), are positioned within chamber 102 to measure the evaporation rate of the raw materials 106. At step 304, a response data acquisition unit 404 obtains measurement signals 410 and 412 from the diagnostic sensors 406 and thermocouples 110, respectively, in response to the test signals 408. The data acquired during this procedure is used to produce the first dynamic model 202. The first dynamic model 202 indicates the behavior of the system between the high bandwidth and low bandwidth dynamics at the first time scale index k.
As previously discussed, the first dynamic model is represented by both yll measurements and yl2 measurements. During production, the yll measurements are obtained from thermocouples 110. The yl2 measurements (evaporation rates of the raw materials 106), however, are not directly available for observation during the deposition process so that they must be /27425
-13-
inferred from indirect and noisy measurements. The yll measurements obtained from thermocouples 110 are also noisy. For this reason, a first estimator 206 is created to produce estimates of the states and measurements of the first model and to predict the raw material evaporation rates (yl2 measurements) from the heater powers and thermocouple signals. The end- of-run film composition measurements (y2) are also noisy and for this reason a second estimator 208 is created to provide estimates of the states x2 and measurement inputs y2 of the second dynamic model.
The second estimator 208 is represented by the following set of equations:
Figure imgf000015_0001
where x2 represents the state of the estimator and y2 represents the estimated end-of-run film compositions measurements. In one embodiment, the state of the estimator includes elements of the state error covariance matrix. Equations (5), (6) and (7) express the mathematical relationship of the states x2 and the pi2 parameters in one embodiment of the invention. In one embodiment, estimator 208 is determined using all yllk measurements obtained at the first time scale in the last period of the second time scale as shown in equation (10) below. In such an instance, estimator 208 is represented by equations (8b) and (9b) below.
Figure imgf000015_0002
/27425
-14-
ynj is represented by:
Figure imgf000016_0001
where n is the ratio of sample rate 1 to sample rate 2.
In accordance with the present invention, the end-of-run film composition measurements acquired at the second time scale are used to refine the pn parameters of the first dynamic model in the first time scale. The first estimator 206 is, thus, created using the following set of equations:
lM ftl(*t'Uh >kp2j)'yl lk) (Eq. ll)
Figure imgf000016_0002
where *2j is the most recent available update of the second state estimates at time index k.
Construction of the first and second estimators 206 and 208 from models 202 and 204, respectively, are off-line activities.
A controller 210 for controlling the deposition system 100 at the first time scale is then created using the estimates from the first estimator 206 and the second estimator 208. The controller 210 is written as: /27425
U s= z,(*l*. ιl.Λ(jc2y) 'lέ,wέ). (Eq. 13)
where w represents a predetermined process recipe or setpoint that will result in a desired film composition being deposited on the surface of substrate 120.
The first and second estimators 206 and 208 may be created based on a variety of filtering techniques. In one embodiment, the first estimator 206 is created by first linearizing equations (1) and (2) to create a set of linear model equations that are represented by: jclt+I = Aljcl4 + βlttlέ (Eq. 14)
ylk = Clxlk + DMk (Eq. 15)
where Al, Bl, Cl, and Dl represent the linearization of i and gi in equations (1) and (2) and where Al is of dimension nsi x nsi, Bl is of dimension ns x n \, Cl is of dimension ny_ x ns\, .and Dl is of dimension nyi xnu\. The term Al is of dimension ns\ is the number of states, n i is the number of outputs, and nui is the number of inputs in the first dynamic model. In one embodiment nsj = 30, yi = 20, and nui - 5.
As before, the output measurement yl is divided into the measurable and the non-measurable parts as follows:
Figure imgf000017_0001
The vector yllk corresponds to the measured temperatures from thermocouples 110 and is of dimension ny\\ x 1. The vector yl2}< corresponds to the unmeasured evaporation rates and is of dimension nyu x 1. In one embodiment nyπ=10 and nyi2=10. 9/27425
-16-
Correspondingly, the output matrix C is partitioned as
Figure imgf000018_0001
where the matrix Cn is of dimension nyπ x nsi and the matrix C 2 is of dimension n x nsi.. A time-invariant Kalman filter is then obtained by solving the following steady state Riccati equation
AlTXAl - X - AlTXC (& + Cl lXqlriQ, = 0 (Eq. 18)
where is the process noise covariance and Qv is the measurement noise covariance.
From the Riccati solution, X, the Kalman filter gain matrix is given by, = XC^(Q, + CI 1XCl l r)"land the Kalman filter equations are given by equations (18) and (19) below.
3 *+ι = A άk + B\u\k + L{y\ lk - C, , xfk) (Eq. 19)
yϊ* = C xTk (Eq. 19)
In another embodiment, multiple linearized models are obtained each with index i for different steady-state operating pints and the gain matrices L/ are obtained by solving steady-state Riccati equations using reduced order ARX model parameters as disclosed in column 8, line 22 through column 9, line 22 of United States patent no. 5,517,594. United States patent no. 5,517,594 entitled "Thermal Reactor Optimization" is hereby incorporated by reference.
It is appreciated that the present invention is not limited to the manner in which the first or second estimators are created. For example, 27425
-17- estimators 206 and 208 may be created using a continuous-time extended Kalman filter (CTEKF) for state estimations in cases where equations (1) and /or (3) are continuous-time dynamic models, and in some instances may also be based on the diagnostic measurements, yl2. The estimators may also be created using other filtering techniques such as a discrete-time extended Kalman filter (DTEKF) for state estimation. In one embodiment, the first estimator 206 is based on an extended Kalman filter (EKF) that uses discrete time measurements and integration. This method is described in co-pending United States patent application entitled "A Method for Real-Time Nonlinear System State Estimation and Control", filed by Sunil C. Shah and Pradeep Pandey on August 22, 1997 and hereby incorporated by reference.
Although the discussion has thus been limited to controlling an evaporator deposition system, it is to be understood that the present invention is not limited to a specific type of system to be controlled, nor is it limited to a specific type of control system. In addition, it is appreciated that the present invention is not limited to a method of controlling a physical system that is characterized by only two time scales. Moreover, it is important to note that the scope of the present invention is not limited to control systems, but is also applicable to the field of signal processing in general.
In accordance with the present invention, signals generated from any physical system may be processed pursuant to the method outlined in the flow chart of Figure 5. In step 502, a first dynamic model is created that indicates a behavior of the physical system at a first time scale. In step 504, a second dynamic model is created that indicates a behavior of the physical system at a second time scale. In step 506, a second estimator is created to estimate the states of the second dynamic model using the second dynamic model and measurements from the physical system at the second time scale. /27425
-18- In step 508, a first estimator is created to estimate the states of the first dynamic model using the first dynamic model, measurements from the physical system at the first time scale, and estimates from the second estimator.
By processing the signals from the physical system in the manner outlined in Figure 5, a controller for controlling the physical system at the first time scale may be created using estimates from the first and second estimators. Equations (1) through (4), (8), (9) and (11) through (13) generally represent the manner in which the dynamic models, estimators and controller are constructed in one embodiment of the invention.
By decomposing the signal processing and control problem into a series of problems, the number of operations required to solve the problem is reduced compared to estimating all parameters at the same time scale. In the example above, the deposition system is defined by a first set of parameters (fast parameters) that is changing at a fast rate and a second set of parameters (slow parameters) that is changing at a slow rate. The number of calculations required per time period to solve the combined set of parameters at a single fast time scale is represented by k(nsχ+ns2)^, where nsi is the number of states of the first model and ns2 is the number of states in the second model. By decomposing the problem into two separate problems, the number of calculations required per time period to solve the parameters is represented by the kns13+kns2^. Therefore, if nsi = s2 the number of operations required to solve the first and second set of parameters is reduced by a factor of four. The number of calculations is also reduced by the fact that the slow parameters are updated at a slower rate than the fast parameters. Therefore, instead of updating all of the parameters of the system at the fast time scale, the parameters are divided among different time scales and updated /27425
-19-
accordingly. By updating the slow parameters at a slower rate than the fast parameters, the number of calculations required to solve the signal processing or control problem is reduced. The benefits of the present invention become even more pronounced in systems containing more than two sets of parameters that are changing or being updated at different rates.
In the foregoing discussion, a method for controlling a single deposition system has been described wherein measurements are obtained from the system at multiple times scales. These measurements are then used in the manner previously described to create a controller for controlling the performance of the single deposition system.
The high volume requirements of most manufacturing applications, however, require the use of multiple systems in order to meet production quotas. This typically includes using a family of "identical" systems that are controlled using a common control scheme. For example, in lieu of using a single deposition chamber to produce photovoltaic devices on the surfaces of substrates, it would be common to employ a family of deposition chambers in order to increase the production of the devices. In theory, each deposition chamber within the family of chambers is manufactured so that all the chambers are identical. In reality, however, chambers-to-chamber variations will exist. For example, the thermal characteristics of each of the chambers within the family may vary due to differences in the heater elements, insulation properties, emissivity variations, etc..
The present invention provides a method for controlling two or more systems within a family of systems which accounts for the differences between each of the systems. Figure 6 is a flow chart that outlines a method for controlling a family of systems consisting of a first system and a second system. In step 602, a first dynamic model is created at a first time scale that /27425
-20- indicates the behavior of the first system at the first time scale. This typically involves exciting the first system and capturing the response of the first system in a manner similar to that described above in Figures 3 and 4. The first dynamic model is represented by equations (1) and (2) above, wherein the p!2 parameters capture equipment-to-equipment variations between the first system and the second system.
In step 604, a second dynamic model is created at a second time scale that indicates the behavior of the of the first system at the second time scale. The second dynamic model is represented by equations (3) and (4) above.
In steps 604, 606, and 608, a first estimator, second estimator, and first controller are created in a manner similar to that described in the example above. Accordingly, equations (8) through (12) may be used to represent the first estimator, second estimator, and first controller for the first system.
A second controller for optimizing the performance of the second system is then created by applying the first dynamic model to the second system to create a third dynamic model. (Step 612) The third dynamic model indicates the behavior of the second system at the first time scale and is represented by the following set of model equations:
Figure imgf000022_0001
ΛβS.(*3*.κ3*.A.' »)- A2*22)
where x3 is the state of the second system, u3 are inputs to the second system, y3 are measurements from the second system, pn represents fixed parameters, and p32 represents refutable parameters, including those „,„,„,_ PCT 99/27425
-21-
parameters that represent equipment variations between the first and second systems.
In step 614, a fourth dynamic model is created to indicate the behavior of the second system at the second time scale and is represented by the following set of model equations:
Figure imgf000023_0001
where x4 is the state of the second system at the second time scale, u4 inputs to the second system , p4 are fixed parameters of the second system at the second time scale, and y4 are measurements of the second system at the second time scale. The p32 parameters of the third dynamic model are related to the states x4 of the fourth dynamic model. That is, p32 = h(x4).
In step 616, a fourth estimator is created to estimate the states of the fourth dynamic model. The fourth estimator is created using the fourth dynamic model and measurements from the second system at the second time scale. In some instances, measurements from the second system at the first time scale are additionally used to create the fourth estimator. The fourth estimator may then be represented by the following set of equations.
Figure imgf000023_0002
In step 618, a third estimator is created to estimate the states of the third dynamic model. The third estimator is created using the third dynamic model, measurements from the second system at the first time scale and the ,„„,„ PCT/US98/2388 /27425
-22- estimates from the fourth estimator and is represented by equations (26) and (27) below.
*ιβΛ,( 'W3». ,1.ΛR)'3'31J (Eq.27)
Figure imgf000024_0001
In step 620, a second controller for controlling the second system at the first time scale is created using estimates from the third and fourth estimators. The second controller maybe written as:
Figure imgf000024_0002
In one embodiment, the first dynamic model of the first system is applied to the second system in order to create the third dynamic model after it is initially created. In an alternative embodiment, the first dynamic model is first refined by implementing steps 602 through 610 before it is applied to the second system to create the third dynamic model. The ability to apply the first dynamic model of the first system to other systems alleviates the need to obtain diagnostic measurements from these systems during setup. As a result, the setup time and the costs associated with setting up each system is therefore reduced.
The foregoing examples have been limited to manufacturing applications. It is appreciated, however, that the present invention is not limited to such applications, nor is it limited to control systems. The present invention may also be applied to signal processing in general. For example, in the field of magnetic disk drives read-write signal processing is required. Typically, a magnetic disk drive consist of a stack of circular disks having a thin magnetic coating rotating at high speed. Data is recorded on the disk surfaces using transducers mounted on read-write heads. The read-write heads are attached to actuator arms that are moved across the disk surfaces by a high-speed actuators that are under the control of a servo motor. The actual information is recorded in circumferential tracks on the disk surfaces. Reading the recorded information back involves sensing the magnetic transitions emanating from the surfaces of the magnetic disks, again using the transducers.
The ability to store and read back information from each of the magnetic disks is affected by vibrations in the arm assemblies which affect the position of the read-write head transducers within the circumferential tracks. Magnetic modulation of the thin magnetic coatings will also affect the signal that is produced by the read-write head transducers. The read-write head transducer signals will thus consist of confounded information that is representative by both the vibration and magnetic modulation.
Variations also exists between each of the read-write heads. These variations may be in the form of different read-write head transducer sensitivities, differences in the physical characteristics of the read-write heads, nominal tolerance differences, etc..
In accordance with the present invention, the read-write signals may be processed by first creating a first model at a first time scale that indicates and takes into account the vibrational and magnetic modulation components of the magnetic disk drive. The first model may be represented by the following set of model equations:
9/27425
-24-
ylk s g.[xh>uL,pu.pn). (Eq.31)
where xl is the position .and magnetic states of the system, wl is servo motor commands, yl are measurements from taken from magnetic disk drive unit, pil represents fixed system parameters, and pi2 represents head mass, offset and film sensitivity parameters. The pi2 parameters are subject to change and/or may vary from one magnetic disk/head assembly to another. Therefore, it is important that these parameters be updated or refined to ensure that the first model properly indicates the behavior of the disk drive unit. The pn parameters may include dynamics parameters, such as drag and inertia parameters. The yl measurements consist of two sets of measurements; yll measurements, and l2 measurements. The yll measurements include servo motor drive current, the read-write head sensor signal, actuator arm position, and servo motor velocity. The yl2 measurements may include read-write head positioning measurements that are obtained with laser positioning diagnostic sensors.
A second model at a second time scale is also created which indicates how of each of the read-write heads, servo motor parameters, structural vibration parameters, friction and drag parameters change over time and from component to component. The second model also indicates how measurements of previously recorded data changes with drive and component parameters x2. The second model equations may be written as:
Figure imgf000026_0001
where x2 comprise head mass and transducer sensitivity states of the system, 2 is the head position setpoint offsets, p2 are fixed parameters of the magnetic disk drive in the second time scale, and y2 is the measurement of previously recorded data. The pi2 parameters of the first model are related to the states x2 of the second dynamic model. That is, pi2 = h(x2).
Thus, in accordance with the present invention, read-write signal processing is achieved by decomposing the magnetic disk drive into two separate models and updating the models in accordance with the method of Figure 5.
As previously discussed, the methods of the present invention are applicable to systems that are characterized by different space scales. Figure 7A exemplifies such a system. In Figure 7A, a furnace 702 for processing a plurality of wafers 704a-704d is shown. The furnace 702 includes a quartz jar 706 having a set of heating elements 708 that surround the quartz jar. During processing of the wafers 704a-704d, it is necessary to control the temperature of the wafers in order to achieve a certain degree of process uniformity.
A furnace controller may be created to control the temperature of the wafers by first determining a first dynamic model at a fine, space scale given by index 1, that consists of a first set of parameters, a first set of states, and a first set of inputs. The first dynamic model at the fine space scale may be written as:
x\M = /.(Λ .KI/./ J. (Eq- 34)
Figure imgf000027_0001
As illustrated in Figure 7b, xl comprises wafer 704d temperatures, ul is the heater powers, yl are temperature measurements taken from the system, pn is the thermal dynamics of the system, and pi2 represents the center and edge temperatures of wafer 704d. The yl measurements consist of two sets of measurements; yll measurements, and yl2 measurements. The yll measurements are obtained from production thermocouples within the furnace, whereas the yl2 measurements are obtained from diagnostic sensors that are not available during the production process.
A second dynamic model at a coarse space scale given by index m, is also created and consists of a second set of parameters, a second set of states, and a second set of inputs. The second dynamic model may be represented by the following equations:
x2m.l = f2[x2m>u2m ,p2), (E 36>
y2m = g2(χ2m,u2m 2), (Eq.37)
As illustrated in Figure 7C, x2 is the center .and edge temperatures of wafers 704a-704d, u2 is the heater powers, p2 represents the dynamics of the system at the coarse space scale, and y2 represents furnace production sensor measurements.
A second estimator of the coarse model is created using the following set of equations:
Figure imgf000028_0001
-27- where x2 represents estimates of the wafer center and edge temperatures and y2 represents estimates of the production thermocouple measurements.
In accordance with the present invention, the center and edge wafer temperatures at the coarse space scale are used to refine the pi2 parameters of the first dynamic model in the fine space scale. A first estimator of the fine model is, thus, created using equations (40) and (41) below.
Figure imgf000029_0001
A controller for controlling the furnace 702 at the fine space scale is created using the estimates from the first estimator and the second estimator. The controller is written as:
Figure imgf000029_0002
where w represents a predetermined process recipe or setpoint.
In the foregoing specification the invention has been described with reference to specific exemplary embodiments thereof. It will, however, be evident that various modifications and changes may be made thereto without departing from the broader spirit and scope of the invention as set forth in the appended claims. The specification and drawings are accordingly to be regarded as illustrative rather than a restrictive sense.

Claims

CJ ΔIM≤What is claimed is:
1. A method of processing signals generated from a first physical system comprising the steps of: a) creating a first dynamic model at a first time scale that consists of a first set of parameters, a first set of states and a first set of inputs; b) creating a second dynamic model at a second time scale that consists of a second set of parameters, a second set of states and a second set of inputs, at least one of the first set of parameters in the first dynamic model being computed from the second set of states in the second dynamic model; c) creating a second estimator to produce estimates of the second set of states using the second dynamic model, measurements from the first physical system at the second time scale; and d) creating a first estimator to produce estimates of the first set of states using the first dynamic model, measurements from the first physical system at the first time scale, and the estimates from the second estimator.
2. The method of claim 1 further comprising the step of creating a first controller that operates at the first time scale based on the estimates from the first and second estimators, the controller optimizing the performance of the first physical system.
3. The method of claim 1 wherein the step of creating the second estimator is performed by further using measurements from the first physical system at the first time scale. 9/27425
-29-
4. The method of claim 1 wherein the first set of parameters comprises a first subset of parameters that are constant and a second subset of parameters that vary with time, the second subset of parameters computed from the second set of states in the second model.
5. The method of claim 1 wherein the first dynamic model is created using diagnostic sensors.
6. The method of claim 1, wherein the step of creating the first dynamic model comprises the steps of: creating a set of test signals that exercise high bandwidth and low bandwidth dynamics of the first physical system; applying the test signals to the first physical system while capturing a synchronous response of the first physical system to the test signals, the first dynamic model correlating the test signals and synchronous response of the first physical system.
7. The method of claim 1 wherein at least one of the estimators is based on linearized models at multiple operating points.
8. The method of claim 1 wherein at least one of the estimators is based on a continuous-time extended Kalman filter (CTEKF) for state estimation.
9. The method of claim 1 wherein at least one of the estimators is based on a continuous-time extended Kalman filter (CTEKF) for state estimation and the diagnostic measurements. 9/27425
-30-
10. The method of claim 1 wherein at least one of the estimators is based on a discrete-time extended Kalman filter (DTEKF) for state estimation.
11. The method of claim 1 further comprising the step of processing signals generated from a second physical system which comprises the steps of a) applying the first dynamic model to the second physical system to create a third dynamic model at the first time scale that represents the behavior of the second physical system, the third dynamic model consisting of a third set of parameters, a first set of states and a first set of inputs; b) creating a fourth dynamic model at the second time scale that consists of a fourth set of parameters, a fourth set of states and a fourth set of inputs, at least one of the third set of parameters in the third dynamic model being computed from the fourth set of states in the fourth dynamic model; c) creating a fourth estimator to produce estimates of the fourth set of states using the fourth dynamic model and measurements from the second physical system at the second time scale; and d) creating a third estimator to produce estimates of the third set of states using the third dynamic model, measurements from the second physical system at the first time scale, and the estimates from the fourth estimator; and d) creating a second controller that operates at the first time scale based on the estimates from the third and fourth estimators, the controller optimizing the performance of the second physical system.
12. The method of claim 11 wherein the step of creating the fourth estimator is performed by further using measurements from the second physical system at the first time scale.
13. A method for controlling a second physical system in a family of physical systems that includes a first physical system and the second physical system, the method comprising the steps of: a) creating a first dynamic model at a first time scale that indicates the behavior of the first physical system, wherein the first dynamic model consists of a first set of states, a first set of inputs, and a first set of parameters, the first set of parameters comprising a first subset of parameters and a second subset of parameters, the first subset of parameters being common to the first and second physical systems; b) creating a second dynamic model at a second time scale that indicates the behavior of the first physical system at the second time scale, wherein the second dynamic model consists of a second set of parameters, a second set of states and a second set of inputs, the second subset of parameters in the first dynamic model being computed from the second set of states in the second model; c) creating a second estimator to produce estimates of the second set of states using the second dynamic model and measurements from the second physical system at the second time scale; d) creating a first estimator to produce estimates of the first set of states using the first dynamic model, measurements from the first physical system at the first time scale, and the estimates from the second estimator; e) creating a controller that operates at the first time scale based on the estimates from the first and second estimators, the controller optimizing the performance of the second physical system; f) applying the first dynamic model to the second physical system to create a third dynamic model at the first time scale that represents the behavior of the second physical system, the third dynamic model consisting 9/27425
-32- of a third set of states, a third set of inputs, and a third set of parameters, the third set of parameters comprising the first subset of parameters and a third subset of parameters; g) creating a fourth dynamic model at the second time scale that consists of a fourth set of parameters, a fourth set of states and a fourth set of inputs, the third subset of parameters in the third dynamic model being computed from the fourth set of states in the fourth dynamic model; h) creating a fourth estimator to produce estimates of the fourth set of states using the fourth dynamic model, measurements from the second physical system at the second time scale; i) creating a third estimator to produce estimates of the third set of states using the third dynamic model, measurements from the second physical system at the first time scale, and the estimates from the fourth estimator; and j) creating a second controller that operates at the first time scale based on the estimates from the third and fourth estimators, the controller optimizing the performance of the second physical system.
14. The method of claim 13 wherein the step of creating the second estimator is performed by further using measurements from the first physical system at the first time scale.
15. The method of claim 13 wherein the step of creating the fourth estimator is performed by further using measurements from the second physical system at the first time scale. 9/27425
-33-
16. The method of claim 13, wherein the step of creating the first dynamic model comprises the steps of: creating a set of test signals that exercise high bandwidth and low bandwidth dynamics of the first physical system; applying the test signals to the first physical system while capturing a synchronous response of the first physical system to the test signals, the first dynamic model correlating the test signals and synchronous response of the first physical system.
17. The method of claim 16, wherein the step of capturing a synchronous response of the first physical system comprises the steps of capturing a series of measurements through a set of diagnostic sensors in the first physical system.
18. The method of claim 16, wherein the step of capturing a synchronous response of the first physical system comprises the steps of capturing a series of measurements through a set of production sensors and a set of diagnostic sensors in the first physical system.
19. The method of claim 13, wherein each first and second physical system comprises an evaporation chamber and a deposition target.
20. The method of claim 13 wherein at least one of the estimators is based on linearized models at multiple operating points.
21. The method of claim 13 wherein at least one of the estimators is based on a continuous-time extended Kalman filter (CTEKF) for state estimation. /27425
-34-
22. The method of claim 13 wherein at least one of the estimators is based on a continuous-time extended Kalman filter (CTEKF) for state estimation and the diagnostic measurements.
23. The method of claim 13 wherein at least one of the estimators is based on a discrete-time extended Kalman filter (DTEKF) for state estimation.
24.. A method of processing signals generated from a physical system comprising the steps of: a) creating a first dynamic model in a first space scale that consists of a first set of parameters, a first set of states and a first set of inputs; b) creating a second dynamic model in a second space scale that consists of a second set of parameters, a second set of states and a second set of inputs, at least one of the first set of parameters in the first dynamic model being computed from the second set of states in the second dynamic model; c) creating a second estimator to produce estimates of the second set of states using the second dynamic model, measurements from the physical system in the second space scale; and d) creating a first estimator to produce estimates of the first set of states using the first dynamic model, measurements from the physical system in the first space scale, and the estimates from the second estimator.
25. The method of claim 24 further comprising the step of creating a controller that operates in the first space scale based on the estimates from the first and second estimators, the controller optimizing the behavior of the physical system. /27425
-35-
26. The method of claim 24 wherein the first dynamic model is created using diagnostic sensors.
27. The method of claim 24 wherein at least one of the estimators is based on linearized models at multiple operating points.
28. The method of claim 24 wherein at least one of the estimators is based on a continuous-time extended Kalman filter (CTEKF) for state estimation.
29. The method of claim 24 wherein at least one of the estimators is based on a continuous-time extended Kalman filter (CTEKF) for state estimation and the diagnostic measurements.
30. The method of claim 24 wherein at least one of the estimators is based on a discrete-time extended Kalman filter (DTEKF) for state estimation.
PCT/US1998/023889 1997-11-26 1998-11-10 Multiple scale signal processing and control system WO1999027425A1 (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
AU14539/99A AU1453999A (en) 1997-11-26 1998-11-10 Multiple scale signal processing and control system
DE69828995T DE69828995T2 (en) 1997-11-26 1998-11-10 SIGNAL PROCESSING AND CONTROL SYSTEM AT MULTIPLE LEVELS
AT98958508T ATE289083T1 (en) 1997-11-26 1998-11-10 SIGNAL PROCESSING AND CONTROL SYSTEM ON MULTIPLE LEVELS
EP98958508A EP1034459B1 (en) 1997-11-26 1998-11-10 Multiple scale signal processing and control system

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US08/980,041 1997-11-26
US08/980,041 US6041172A (en) 1997-11-26 1997-11-26 Multiple scale signal processing and control system

Publications (1)

Publication Number Publication Date
WO1999027425A1 true WO1999027425A1 (en) 1999-06-03

Family

ID=25527312

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US1998/023889 WO1999027425A1 (en) 1997-11-26 1998-11-10 Multiple scale signal processing and control system

Country Status (6)

Country Link
US (3) US6041172A (en)
EP (1) EP1034459B1 (en)
AT (1) ATE289083T1 (en)
AU (1) AU1453999A (en)
DE (1) DE69828995T2 (en)
WO (1) WO1999027425A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001047110A2 (en) * 1999-12-22 2001-06-28 Burr-Brown Japan, Ltd. Method and apparatus for changing characteristics of device

Families Citing this family (40)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20000048977A (en) * 1996-10-08 2000-07-25 칼 하인쯔 호르닝어 Method and control structure for controlling moments in numerically controlled elastic (and therefore oscillation-capable) multiple mass systems
US6041172A (en) * 1997-11-26 2000-03-21 Voyan Technology Multiple scale signal processing and control system
US6687658B1 (en) * 1998-09-01 2004-02-03 Agere Systems, Inc. Apparatus and method for reduced-order modeling of time-varying systems and computer storage medium containing the same
US7878882B2 (en) * 1999-04-01 2011-02-01 Charles J. Molnar Advanced workpiece finishing
US20130189801A1 (en) * 1998-11-06 2013-07-25 Semcon Tech, Llc Advanced finishing control
US8353738B2 (en) * 1998-11-06 2013-01-15 Semcon Tech, Llc Advanced finishing control
US6498694B1 (en) 1999-06-15 2002-12-24 Voyan Technology Servo error integration in read-channel equalization
US6495805B2 (en) 2000-06-30 2002-12-17 Tokyo Electron Limited Method of determining set temperature trajectory for heat treatment system
DE10036971A1 (en) * 2000-07-28 2002-02-28 Siemens Ag Method for remote diagnosis of a technological process
JP4493192B2 (en) 2000-09-13 2010-06-30 東京エレクトロン株式会社 Batch heat treatment apparatus and control method thereof
KR100790430B1 (en) * 2000-09-21 2008-01-02 엠티에스 시스템즈 코포레이숀 Multiple region convolver with tapering
US7377836B1 (en) 2000-10-10 2008-05-27 Beaver Creek Concepts Inc Versatile wafer refining
US6945780B2 (en) * 2001-04-02 2005-09-20 United Defense, L.P. Integrated performance simulation system for military weapon systems
US7577951B2 (en) * 2002-05-30 2009-08-18 Hewlett-Packard Development Company, L.P. Performance of computer programs while they are running
US7346113B2 (en) * 2002-07-26 2008-03-18 Texas Instruments Incorporated Method and circuit for stop of signals quantized using noise-shaping
US7328074B2 (en) * 2002-12-02 2008-02-05 United Technologies Corporation Real-time quadratic programming for control of dynamical systems
US20040186804A1 (en) * 2003-03-19 2004-09-23 Anindya Chakraborty Methods and systems for analytical-based multifactor multiobjective portfolio risk optimization
US7640201B2 (en) * 2003-03-19 2009-12-29 General Electric Company Methods and systems for analytical-based multifactor Multiobjective portfolio risk optimization
US7593880B2 (en) * 2003-03-19 2009-09-22 General Electric Company Methods and systems for analytical-based multifactor multiobjective portfolio risk optimization
US7461045B1 (en) * 2003-08-18 2008-12-02 University Of Florida Research Foundation, Inc. Optimization of spatio-temporal pattern processing for seizure warning and prediction
US8219477B2 (en) * 2004-02-20 2012-07-10 General Electric Company Systems and methods for multi-objective portfolio analysis using pareto sorting evolutionary algorithms
US7542932B2 (en) * 2004-02-20 2009-06-02 General Electric Company Systems and methods for multi-objective portfolio optimization
US7630928B2 (en) * 2004-02-20 2009-12-08 General Electric Company Systems and methods for multi-objective portfolio analysis and decision-making using visualization techniques
US8126795B2 (en) * 2004-02-20 2012-02-28 General Electric Company Systems and methods for initial sampling in multi-objective portfolio analysis
US7469228B2 (en) * 2004-02-20 2008-12-23 General Electric Company Systems and methods for efficient frontier supplementation in multi-objective portfolio analysis
US7203554B2 (en) * 2004-03-16 2007-04-10 United Technologies Corporation Model predictive controller with life extending control
US7203555B2 (en) * 2004-05-14 2007-04-10 University Of Delaware Predictive regulatory controller
US20060282178A1 (en) * 2005-06-13 2006-12-14 United Technologies Corporation System and method for solving equality-constrained quadratic program while honoring degenerate constraints
CN104834294A (en) 2005-12-05 2015-08-12 费舍-柔斯芒特系统股份有限公司 Muti-objective predictive process optimization with concurrent process simulation
US7376471B2 (en) * 2006-02-21 2008-05-20 United Technologies Corporation System and method for exploiting a good starting guess for binding constraints in quadratic programming with an infeasible and inconsistent starting guess for the solution
US7991499B2 (en) 2006-12-27 2011-08-02 Molnar Charles J Advanced finishing control
DE102007005070B4 (en) * 2007-02-01 2010-05-27 Klippel, Wolfgang, Dr. Arrangement and method for the optimal estimation of the linear parameters and the non-linear parameters of a model describing a transducer
US8357286B1 (en) 2007-10-29 2013-01-22 Semcon Tech, Llc Versatile workpiece refining
WO2009114941A1 (en) * 2008-03-20 2009-09-24 University Of New Brunswick Method of multi-dimensional nonlinear control
US7885717B2 (en) * 2008-03-31 2011-02-08 Sikorsky Aircraft Corporation Fast control law optimization
US8215013B2 (en) 2008-04-11 2012-07-10 Bal Seal Engineering, Inc. Method for making a free standing axially compressed connector stack
US8096838B2 (en) 2009-03-11 2012-01-17 Bal Seal Engineering, Inc. Header assembly for implantable medical devices
US8328587B2 (en) 2009-04-20 2012-12-11 Bal Seal Engineering, Inc. In-line connector stack with testing capability
US9466915B2 (en) 2011-10-03 2016-10-11 Bal Seal Engineering, Inc. In-line connectors and related methods
GB2526806A (en) * 2014-06-03 2015-12-09 Here Global Bv Trail interpolation

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5517594A (en) * 1994-10-17 1996-05-14 Relman, Inc. Thermal reactor optimization

Family Cites Families (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS5716719A (en) * 1980-07-04 1982-01-28 Hitachi Ltd Method and equipment for controlling steam temperature in thermal power plant
US4620286A (en) * 1984-01-16 1986-10-28 Itt Corporation Probabilistic learning element
US4599692A (en) * 1984-01-16 1986-07-08 Itt Corporation Probabilistic learning element employing context drive searching
EP0247163B1 (en) * 1985-11-27 1992-02-05 Trustees Of Boston University Pattern encoding system
US4814050A (en) * 1986-10-06 1989-03-21 Aluminum Company Of America Estimation and control of alumina concentration in hall cells
US4852018A (en) * 1987-01-07 1989-07-25 Trustees Of Boston University Massively parellel real-time network architectures for robots capable of self-calibrating their operating parameters through associative learning
US4881178A (en) * 1987-05-07 1989-11-14 The Regents Of The University Of Michigan Method of controlling a classifier system
US4974191A (en) * 1987-07-31 1990-11-27 Syntellect Software Inc. Adaptive natural language computer interface system
FI896219A0 (en) * 1989-04-28 1989-12-22 Antti Aarne Ilmari Lange ANALYZING AND FOUNDATION CALIBRATION AV DETEKTORSYSTEM.
US5060132A (en) * 1989-06-13 1991-10-22 Elsag International B.V. Method of modeling and control for delignification of pulping
US5191521A (en) * 1990-06-18 1993-03-02 Controlsoft, Inc. Modular multivariable control apparatus and method
JP2539540B2 (en) * 1990-09-19 1996-10-02 株式会社日立製作所 Process control equipment
US5214715A (en) * 1991-01-31 1993-05-25 Trustees Of Boston University Predictive self-organizing neural network
US5412060A (en) * 1991-08-20 1995-05-02 Bayer Aktiengesellschaft Process for the production of polycarbonate of constant viscosity
US5453226A (en) * 1992-02-05 1995-09-26 The Board Of Regents Of The University Of Oklahoma Curing of composite materials using extended heat transfer models
US6041141A (en) * 1992-09-28 2000-03-21 Matsushita Electric Industrial Co., Ltd. Character recognition machine utilizing language processing
US5377307A (en) * 1992-10-07 1994-12-27 Schlumberger Technology Corporation System and method of global optimization using artificial neural networks
US5740033A (en) * 1992-10-13 1998-04-14 The Dow Chemical Company Model predictive controller
US5355305A (en) * 1992-10-29 1994-10-11 Johnson Service Company Pattern recognition adaptive controller
US5424963A (en) * 1992-11-25 1995-06-13 Photon Research Associates, Inc. Molecular dynamics simulation method and apparatus
CA2129761A1 (en) * 1993-08-11 1995-02-12 David G. Taylor Self-tuning tracking controller for permanent-magnet synchronous motors
US5579440A (en) * 1993-11-22 1996-11-26 Brown; Robert A. Machine that learns what it actually does
US5464369A (en) * 1994-02-25 1995-11-07 Johnson Service Company Method and apparatus for estimating the rate at which a gas is generated within an enclosed space
AU2968295A (en) * 1994-08-10 1996-03-07 Motorola, Inc. Cascade tuning controller and method of use therefor
US5748847A (en) * 1995-12-21 1998-05-05 Maryland Technology Corporation Nonadaptively trained adaptive neural systems
JP2806349B2 (en) * 1996-03-13 1998-09-30 日本電気株式会社 Sputtering apparatus simulation method
US5909381A (en) * 1997-02-19 1999-06-01 Itt Manufacturing Enterprises, Inc. System of on board prediction of trajectories for autonomous navigation of GPS satellites
US6041172A (en) * 1997-11-26 2000-03-21 Voyan Technology Multiple scale signal processing and control system
WO1999030229A1 (en) * 1997-12-11 1999-06-17 Digits Corp. Object code analysis and remediation system and method
US6096085A (en) * 1998-03-23 2000-08-01 The United States Of America As Represented By The Secretary Of The Navy Computer-readable software and computer-implemented method for performing an integrated sonar simulation

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5517594A (en) * 1994-10-17 1996-05-14 Relman, Inc. Thermal reactor optimization

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
BONING D: "RUN BY RUN CONTROL BENCHMARKING", XP002917709, Retrieved from the Internet <URL:HTTP://WWW-MTL.MIT.EDU/RBRBENCH/WHITE-PAPER.HTM> [retrieved on 19990309] *
FRANKLIN G F, POWELL J D, WORKMAN M L: "DIGITAL CONTROL OF DYNAMIC SYSTEMS", DIGITAL CONTROL OF DYNAMIC SYSTEMS, XX, XX, 1 January 1992 (1992-01-01), XX, pages 506 - 513, XP002917706 *
GUYGYI P J, ET AL.: "MODEL-BASED CONTROL OF RAPID THERMAL PROCESSING SYSTEMS", IEEE CONFERENCE ON CONTROL APPLICATIONS, XX, XX, vol. 01, 13 September 1992 (1992-09-13), XX, pages 374 - 381, XP002917708 *
MOZUMDER P K, SAXENA S, COLLINS D J: "A MONITOR WAFER BASED CONTROLLER FOR SEMICONDUCTOR PROCESSES", IEEE TRANSACTIONS ON SEMICONDUCTOR MANUFACTURING., IEEE SERVICE CENTER, PISCATAWAY, NJ., US, vol. 07, no. 03, 1 August 1994 (1994-08-01), US, pages 400 - 411, XP002917707, ISSN: 0894-6507, DOI: 10.1109/66.311345 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2001047110A2 (en) * 1999-12-22 2001-06-28 Burr-Brown Japan, Ltd. Method and apparatus for changing characteristics of device
WO2001047110A3 (en) * 1999-12-22 2002-06-06 Burr Brown Japan Ltd Method and apparatus for changing characteristics of device
US6724264B2 (en) 1999-12-22 2004-04-20 Texas Instruments Incorporated Controller for a device with change in characteristics

Also Published As

Publication number Publication date
EP1034459A4 (en) 2000-11-08
EP1034459B1 (en) 2005-02-09
US6289508B1 (en) 2001-09-11
EP1034459A1 (en) 2000-09-13
DE69828995D1 (en) 2005-03-17
US6167360A (en) 2000-12-26
US6041172A (en) 2000-03-21
AU1453999A (en) 1999-06-15
DE69828995T2 (en) 2006-01-19
ATE289083T1 (en) 2005-02-15

Similar Documents

Publication Publication Date Title
EP1034459B1 (en) Multiple scale signal processing and control system
US10775148B2 (en) Determining a position of a movable part of a coordinate measuring machine
Nagy et al. Robust nonlinear model predictive control of batch processes
US6745087B2 (en) Method for control of a plant
CN103469172B (en) Quartz crystal coated method for controlling thickness and quartz crystal coated device
US10969221B2 (en) Active damping of a measuring device
CN106066636B (en) Have the function of obtaining the Servocontrol device of the frequency characteristic of machinery online
CN101153826A (en) Air pressure sensor
US9640212B2 (en) Topographic feedforward system
CN103576552B (en) Active-disturbance-rejection control method and device of single-input single-output pure lag self-balancing object
CN108107233B (en) Method and system for continuous temperature correction of accelerometer scale factors
Adachi et al. Integration of a Cr–N thin-film displacement sensor into an XY micro-stage for closed-loop nano-positioning
US20230152074A1 (en) Repeated determination of a position of a movable part of a coordinate measuring machine
Evers et al. Identifying thermal dynamics for precision motion control
US8638649B1 (en) Topographic feedforward system
Merced et al. Robust Control of VO $ _ {\bf 2} $-Coated Microbenders Using Self-Sensing Feedback
US10304491B2 (en) Compensating tracking error signal offset in optical tape storage systems
Felix et al. Strain sensing with piezoelectric zinc oxide thin films for vibration suppression in hard disk drives
White et al. Vision based iterative learning control of a MEMS micropositioning stage with intersample estimation and adaptive model correction
Broel-Plater et al. Study of the efficiency of an algorithm that improves the accuracy of slow-motion control of a servo drive
Amthor et al. Fourth order motion profile planning for high precision applications
Choi et al. A sliding-mode control approach in ultra-precision positioning system
Braatz et al. Avery final report: Identification and cross-directional control of coating processes
US8717859B1 (en) Reactionless control of a slider head
Sun Stability of controlled mechanical system with parametric uncertainties in a realistic friction model

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AL AM AT AT AU AZ BA BB BG BR BY CA CH CN CU CZ CZ DE DE DK DK EE EE ES FI FI GB GE GH GM HR HU ID IL IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MD MG MK MN MW MX NO NZ PL PT RO RU SD SE SG SI SK SK SL TJ TM TR TT UA UG UZ VN YU ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): GH GM KE LS MW SD SZ UG ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE BF BJ CF CG CI CM GA GN GW ML MR NE SN TD TG

DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 1998958508

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: KR

NENP Non-entry into the national phase

Ref country code: CA

WWP Wipo information: published in national office

Ref document number: 1998958508

Country of ref document: EP

REG Reference to national code

Ref country code: DE

Ref legal event code: 8642

NENP Non-entry into the national phase

Ref country code: CA

WWG Wipo information: grant in national office

Ref document number: 1998958508

Country of ref document: EP