CN110929366B - 用于计算光学模型模拟的特征内核的方法 - Google Patents
用于计算光学模型模拟的特征内核的方法 Download PDFInfo
- Publication number
- CN110929366B CN110929366B CN201910773551.4A CN201910773551A CN110929366B CN 110929366 B CN110929366 B CN 110929366B CN 201910773551 A CN201910773551 A CN 201910773551A CN 110929366 B CN110929366 B CN 110929366B
- Authority
- CN
- China
- Prior art keywords
- matrix
- feature
- feature matrix
- processor
- computing
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 82
- 238000004088 simulation Methods 0.000 title claims abstract description 67
- 230000003287 optical effect Effects 0.000 title claims abstract description 54
- 239000011159 matrix material Substances 0.000 claims abstract description 236
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 38
- 238000005070 sampling Methods 0.000 claims abstract description 34
- 238000012634 optical imaging Methods 0.000 claims abstract description 22
- 210000001747 pupil Anatomy 0.000 claims description 47
- 238000004422 calculation algorithm Methods 0.000 claims description 29
- 238000003860 storage Methods 0.000 claims description 17
- 230000003044 adaptive effect Effects 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000012952 Resampling Methods 0.000 claims description 5
- 230000008569 process Effects 0.000 description 18
- 238000004519 manufacturing process Methods 0.000 description 11
- 238000012360 testing method Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 10
- 238000005457 optimization Methods 0.000 description 10
- 239000000758 substrate Substances 0.000 description 9
- 238000012937 correction Methods 0.000 description 8
- 238000005259 measurement Methods 0.000 description 8
- 238000001459 lithography Methods 0.000 description 7
- WYROLENTHWJFLR-ACLDMZEESA-N queuine Chemical compound C1=2C(=O)NC(N)=NC=2NC=C1CN[C@H]1C=C[C@H](O)[C@@H]1O WYROLENTHWJFLR-ACLDMZEESA-N 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 230000005540 biological transmission Effects 0.000 description 5
- 239000000284 extract Substances 0.000 description 5
- 229920002120 photoresistant polymer Polymers 0.000 description 4
- 230000010287 polarization Effects 0.000 description 4
- 230000003595 spectral effect Effects 0.000 description 4
- 238000009825 accumulation Methods 0.000 description 3
- 238000000206 photolithography Methods 0.000 description 3
- 239000000047 product Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 102000002274 Matrix Metalloproteinases Human genes 0.000 description 1
- 108010000684 Matrix Metalloproteinases Proteins 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000001311 chemical methods and process Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 238000013141 low-rank factorization Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012797 qualification Methods 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03F—PHOTOMECHANICAL PRODUCTION OF TEXTURED OR PATTERNED SURFACES, e.g. FOR PRINTING, FOR PROCESSING OF SEMICONDUCTOR DEVICES; MATERIALS THEREFOR; ORIGINALS THEREFOR; APPARATUS SPECIALLY ADAPTED THEREFOR
- G03F7/00—Photomechanical, e.g. photolithographic, production of textured or patterned surfaces, e.g. printing surfaces; Materials therefor, e.g. comprising photoresists; Apparatus specially adapted therefor
- G03F7/70—Microphotolithographic exposure; Apparatus therefor
- G03F7/70483—Information management; Active and passive control; Testing; Wafer monitoring, e.g. pattern monitoring
- G03F7/70491—Information management, e.g. software; Active and passive control, e.g. details of controlling exposure processes or exposure tool monitoring processes
- G03F7/705—Modelling or simulating from physical phenomena up to complete wafer processes or whole workflow in wafer productions
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03F—PHOTOMECHANICAL PRODUCTION OF TEXTURED OR PATTERNED SURFACES, e.g. FOR PRINTING, FOR PROCESSING OF SEMICONDUCTOR DEVICES; MATERIALS THEREFOR; ORIGINALS THEREFOR; APPARATUS SPECIALLY ADAPTED THEREFOR
- G03F1/00—Originals for photomechanical production of textured or patterned surfaces, e.g., masks, photo-masks, reticles; Mask blanks or pellicles therefor; Containers specially adapted therefor; Preparation thereof
- G03F1/36—Masks having proximity correction features; Preparation thereof, e.g. optical proximity correction [OPC] design processes
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03F—PHOTOMECHANICAL PRODUCTION OF TEXTURED OR PATTERNED SURFACES, e.g. FOR PRINTING, FOR PROCESSING OF SEMICONDUCTOR DEVICES; MATERIALS THEREFOR; ORIGINALS THEREFOR; APPARATUS SPECIALLY ADAPTED THEREFOR
- G03F1/00—Originals for photomechanical production of textured or patterned surfaces, e.g., masks, photo-masks, reticles; Mask blanks or pellicles therefor; Containers specially adapted therefor; Preparation thereof
- G03F1/68—Preparation processes not covered by groups G03F1/20 - G03F1/50
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
- G06F17/142—Fast Fourier transforms, e.g. using a Cooley-Tukey type algorithm
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Discrete Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Exposure And Positioning Against Photoresist Photosensitive Materials (AREA)
- Image Processing (AREA)
Abstract
本发明提供一种用于计算光学模型模拟的特征内核的方法。在方法中,识别在数学上描述光学成像系统的多种特性的特征矩阵。产生包括至少一个向量的取样矩阵,至少一个向量充当形成特征矩阵的低秩基底的输入。将取样矩阵迭代地乘以特征矩阵,并且根据数值稳定性自适应地重新按比例调整相乘结果直到满足收敛条件。使用迭代结果来形成简化特征矩阵。计算简化特征矩阵的分解值,并且从计算出的分解值中提取多个特征内核。
Description
技术领域
本发明的实施例是有关于一种计算特征内核的方法,且特别是有关于一种用于计算光学模型模拟的特征内核的方法。
背景技术
在现代半导体装置的制造工艺中,操纵各种材料和机器以产生最终产品。由于半导体装置的复杂度增加和超小晶体管的显影,因此工艺变化对产品的性能具有更大影响。光刻是产生影响的重要因素中的一种,并且将光刻模拟应用于评估待在大批量生产之前制造的产品的性能。
在光刻模拟中采用根据霍普金斯理论在数学上描述光学成像系统的特性的传输交叉系数(transmission cross-coefficient,TCC)矩阵,并且计算来自TCC矩阵的TCC内核对模拟是必要的。然而,用于计算TCC内核的标准方法在计算上极为昂贵,这可以使模拟减慢并增加掩模制造工艺的整体时间。
发明内容
本发明实施例提供一种用于计算光学模型模拟的特征内核的方法,适用于具有存储介质及处理器的电子设备,所述方法包括:通过所述处理器识别在数学上描述光学成像系统的多种特性的特征矩阵,所述多种特性从所述光学成像系统中检索并且存储在所述存储介质中;通过所述处理器产生包括至少一个向量的取样矩阵,所述至少一个向量充当形成所述特征矩阵的低秩基底的输入;通过所述处理器将所述取样矩阵迭代地乘以所述特征矩阵,并且根据数值稳定性自适应地重新按比例调整相乘结果直到满足收敛条件;通过所述处理器使用迭代结果以形成简化特征矩阵;以及通过所述处理器计算所述简化特征矩阵的分解值,并且从所述计算出的分解值中提取多个特征内核。
附图说明
根据结合附图阅读的以下详细描述最佳地理解本公开的各个方面。应注意,根据行业中的标准惯例,各种特征并未按比例绘制。事实上,可出于论述清楚起见而任意地增大或缩减各种特征的尺寸。
图1是示出根据本公开的实施例的模型模拟工艺的示意图;
图2是示出根据本公开的实施例的基于模型的掩模制造工艺的示意图;
图3示出根据本公开的实施例的用于计算光学模型模拟的特征内核的电子设备的框图;
图4是示出根据本公开的实施例的用于计算光学模型模拟的特征内核的方法的流程图;
图5A是示出根据本公开的实施例的用于通过使用EVD算法来计算光学模型模拟的特征内核的方法的流程图;
图5B是示出根据本公开的实施例的通过使用EVD算法来计算光学模型模拟的特征内核的实例的示意图;
图6A是示出根据本公开的实施例的用于通过子空间迭代发现低秩基底的方法的流程图;
图6B是示出根据本公开的实施例的用于通过块克雷洛夫方法(block Krylovmethod)发现低秩基底的方法的流程图;
图7A是示出根据本公开的实施例的用于通过使用SVD算法来计算光学模型模拟的特征内核的方法的流程图;
图7B是示出根据本公开的实施例的通过使用SVD算法来计算光学模型模拟的特征内核的实例的示意图;
附图标号说明:
11:光学模型;
12、24:掩模;
13:模拟图像;
21:测试图案;
22:测试测量;
23:校准模型;
25:目标;
26:优化掩模;
30:电子设备;
32:数据检索装置;
34:存储介质;
36:处理器;
S110、S120、S121、S122、S130、S140、S210、S220、S221、S222、S230、S240、S250、S251、S252、S260、S402、S404、S406、S408、S410、S412、S502、S504、S506、S508、S510、S602、S604、S606、S608、S610、S612、S614、S616、S618、S620、S622、S702、S704、S706、S708、S710:步骤。
具体实施方式
以下公开内容提供用于实施所提供的主题的不同特征的许多不同实施例或实例。下文描述组件和布置的具体实例来简化本公开。当然,这些仅仅是实例而并不旨在作为限制。举例来说,在以下描述中,在第二特征上方或在第二特征上形成第一特征可包含第一特征与第二特征直接接触地形成的实施例,且也可包含在第一特征与第二特征之间可形成有额外特征,使得第一特征与第二特征可不直接接触的实施例。此外,本公开可在各种实例中重复附图标号和/或字母。此重复是为了简单和清楚起见,且本身并不指示所论述的各种实施例和/或配置之间的关系。
另外,为易于描述,本文中可能使用例如“在...下面(beneath)”、“在…下方(below)”、“下部(lower)”、“在…上方(above)”、“上部(upper)”等空间相对术语来描述如图中所示出的一个元件或特征与另一(些)元件或特征的关系。除图中所描绘的定向之外,这些空间相对术语旨在涵盖装置在使用或操作中的不同定向。设备可以另外方式定向(旋转90度或处于其它定向),且本文中所使用的空间相对描述词因此可相应地作出解释。
本发明实施例旨在显著地加快传输交叉系数(TCC)内核的计算,这些传输交叉系数内核对现代光刻中的投影光学成像系统的高效建模极为重要。
根据本发明的一些实施例,取样矩阵由随机矩阵、伪随机矩阵或结构矩阵形成。将在数学上描述光学系统的特性的传输交叉系数(TCC)矩阵与取样矩阵迭代地相乘,并且在满足收敛条件之前自适应地重新按比例调整结果(或使结果正交化)。因此,根据迭代的最终结果或迭代的中间结果的累加来形成大致保留相关内核的低秩基底。使用低秩基底来计算较小大小的简化TCC矩阵。
利用简化TCC矩阵,将块算法应用于计算特征值分解(eigenvaluedecomposition,EVD)或奇异值分解(singular value decomposition,SVD),从这些值分解中提取简化TCC内核并将其扩展到全空间(full space)。重建的TCC内核应用于模型校准和掩模校正/优化,以便获得优化掩模。
通过所述方法,由于主要计算仅包含矩阵相乘,所述矩阵相乘通常包含O(N2)运算或甚至在特殊情况下仅包含O(N log N),因此可以显著降低计算复杂度并且可加快基于模型的掩模制造的光刻模拟。
图1是示出根据本公开的实施例的模型模拟工艺的示意图。参看图1,将实施例的模型模拟实施为计算模拟图像,例如使用这些模拟图像来与从晶片捕获到的实际图像进行比较以验证模拟模型是否适合。
在步骤S110中,给出多个光学模型11和掩模12作为模型模拟的输入。每个光学模型是例如包括与投影光学成像系统中的光学器件相关联的许多参数的数学模型,例如散焦、图像平面,光的波长或数值孔径。掩模是例如光刻中所使用的光刻光掩模且在集成电路制造中用图案层显影。
在步骤S120中,将每个光学模型应用于基于霍普金斯理论实施的模型模拟。霍普金斯理论使用传输交叉系数(TCC)来描述掩模传输对光学图像的影响,通过对TCC乘以掩模函数的乘积进行傅立叶变换来形成所述光学图像。
详细地说,可以如下形式呈现标量光学器件(scalar optics)的TCC:
T(f,f′)=∫P(f+fs)S(fs)P*(f′+fs)dfs,
其中自变量是2D空间频率坐标,即,f=(fx,fy)。P(f)是等于点扩散函数的傅里叶变换的光瞳函数。S(f)是源强度分布。离散化T(f,f′)是TCC矩阵,在所述矩阵中,一个维度用于f,而另一个维度用于f′。
也可以(由于≥0,将/>并入P,P*中的每一个中)的形式编写TCC矩阵T,其中P是仅取决于光学系统的特性(例如光瞳函数、源强度、相关因子、偏振)的“堆叠的移位光瞳矩阵”,且P*是光瞳矩阵P的共轭转置。在现代光刻模拟中,/>且/>可以如下方式使用TCC矩阵T来计算空中图像AI:
AI=∫∫m(f)T(f,f′)m*(f′)e2πi(f-f′)df df′,
其中m(f)是等于掩模衍射算子的掩模传输函数。应注意,根据投影光学器件,由光掩模上的图案衍射的光进入光刻胶,并且在光刻胶上形成的光掩模的图像被称为空中图像。
上文所描述的定义用于标量光学理论。类似公式适用于向量光学理论。也就是说,对于向量光学理论,P(f)变为矩阵(即,描述光偏振的多个分量),而S(f)变为矩阵,以便捕获每个源点的偏振。
由此得出,TCC矩阵T允许特征分解T=VΛV*,其中V与列vj形成幺正,且Λ与实际项λ1≥λ2≥…≥0形成对角。TCC内核是经调整的特征向量通常,仅保留主要内核中的少量k<<N,所述主要内核随后包括T直到秩k的最佳近似值。
基于上文,在步骤S121中,产生多个TCC内核,并且在步骤S122中,应用所产生的TCC内核作为光学模型的一部分。
在步骤S121中,本发明实施例的重点在于产生TCC内核。通过将TCC矩阵与低秩基底迭代地相乘,从而获得简化TCC矩阵并使用所述矩阵来计算分解值。最后,可以通过从计算出的分解值中提取简化TCC内核并且将所述简化TCC内核扩展到全空间来产生TCC内核。
除了应用光学模型之外,在步骤S130中,还将抗蚀剂/晶片模型应用于模型模拟,使得在步骤S140中产生模拟图像13。使用模拟图像来与从晶片捕获到的实际图像进行比较,以便验证模拟模型是否适合。
因此,由于产生TCC内核的时间缩短,因此模型模拟工艺可以运行得更快。
图2是示出根据本公开的实施例的基于模型的掩模制造工艺的示意图。参看图2,掩模制造工艺基于如图1中所说明的模型模拟。
在图2的上半部分中,示出用于晶片光刻机和扫描器的投影光学器件的鉴定工艺的模型校准方法,其中应用多个测试图案21以制造光掩模。通过光将测试图案21从光掩模转移到形成于衬底或晶片上的光刻胶,在光刻胶中进行一系列化学处理以将图案雕刻到光刻胶下方的材料中。然后从图案化的衬底或晶片中获得多个测试测量22。
在步骤S210中,给出前述测试图案21和测试测量22作为模型校准的输入。在步骤S220中的模型校准工艺中,应用如图1中所示的模型模拟工艺来模拟具有输入测试图案的光学模型(步骤S221),并且针对测试测量22检查从模型模拟获得的模拟测量(即,模拟图像)以查看模型是否适合(步骤S222)。
如果测试测量22和模拟测量彼此相距较远,那么就意味着模型不稳健,并且工艺返回到步骤S221以校正并改变模型并再次执行模拟。重复执行模型校准,直到测试测量22和模拟测量彼此接近为止,并且因此,在步骤S230中输出校准模型23。
应注意,每当在步骤S221中改变光学模型时,光学器件的参数都会改变,并相应地再次执行模型模拟,这意味着需要计算新TCC内核,且因此计算是减慢模型校准工艺的关键因素。
在图2的下半部分中,示出掩模校正/优化工艺,并且在步骤S240中,给出在先前模型校准工艺中经过校准的模型23、掩模24和目标25作为步骤S250中的掩模校正/优化工艺中的多个模拟的输入。在掩模校正/优化工艺中,模拟中的一个或组合包含但不限于:在步骤S251中执行光学邻近校正(optical proximity correction,OPC)、逆光刻技术(inverselithography technology,ILT)、基于模型的辅助特征(model-based assist feature,MBAF)和源-掩模优化(source-mask optimization,SMO),并且在步骤S252中执行光刻工艺检查(lithography process check,LPC)以检查在校正/优化之后制造的掩模是否适合。
在一些实施例中,在步骤S252中,如果通过使用给定掩模24获得的结果图像经检查为与目标25(例如期望图案)不匹配或类似,那么再次执行步骤S251中的模拟以校正掩模。另一方面,在步骤S252中,如果通过使用给定掩模24获得的结果图像经检查为与目标25匹配或类似,那么在步骤S260中输出优化掩模26。
应注意,尽管步骤S251中的一些模拟(例如OPC/ILT)可以应用先前在步骤S220中产生的TCC内核,但是一些模拟可能需要新TCC内核。举例来说,在一些实施例中,针对源-掩模优化(SMO)的模拟,由于在模拟中优化光学器件,因此改变光学器件的参数并且期望计算新TCC内核。类似于步骤S221中的模拟,步骤S251中的新TCC内核的计算可减慢掩模校正/优化工艺并且导致掩模制造工艺的总时间增加。
鉴于上述情况,本公开开发一种用于计算光学模型模拟的特征内核的方法,以加快TCC内核的计算,从而缩短模型模拟工艺和掩模制造工艺的时间。
详细地说,图3示出根据本公开的实施例的用于计算光学模型模拟的特征内核的电子设备的框图。参看图3,本发明实施例的电子设备30是例如个人计算机、服务器、工作站或任何其它计算装置,并且可包含数据检索装置32、存储介质34以及处理器36,所述处理器耦合到数据检索装置32和存储介质34。
在一些实施例中,电子设备30在外部连接到光学成像系统(未绘示)并且配置以检索光学成像系统的多种特性并且扫描由对衬底或晶片执行热扫描的光学成像系统获得的图像,以便执行光学模型模拟。
在一些实施例中,数据检索装置32配置以连接光学成像系统并且检索光学成像系统的特性。数据检索装置32是例如任何有线或无线接口,例如USB、火线、雷电、通用异步接收器/发射器(universal asynchronous receiver/transmitter,UART)、串行外围接口总线(serial peripheral interface bus,SPI)、WiFi或蓝牙,但是本公开并不限于此。
存储介质34配置以存储由数据检索装置32检索到的特性。存储介质34是例如随机存取存储器(random access memory,RAM)、只读存储器(read-only memory,ROM)、闪存、硬盘、独立磁盘冗余阵列(redundant array of independent disks,RAID)、其它类似存储装置或其组合,但本公开并不限于此。
处理器36是例如中央处理单元(central processing unit,CPU)、其它可编程通用或专用微处理器、数字信号处理器(digital signal processor,DSP)、可编程控制器、专用集成电路(application specific integrated circuit,ASIC)、可编程逻辑装置(programmable logic device,PLD)、其它类似装置或其组合,但本公开并不限于此。
电子设备30配置以执行用于实施根据本公开的一些实施例的用于计算光学模型模拟的特征内核的方法的指令。具体来说,图4是示出根据本公开的实施例的用于计算光学模型模拟的特征内核的方法的流程图。本发明实施例的方法适用于图3的电子设备30,并且下文参看图3的电子设备30中的各种组件描述所述方法的详细步骤。
参看图4,在步骤S402中,处理器36识别在数学上描述光学成像系统的多种特性的特征矩阵。特征矩阵是例如在数学上描述光学系统的特性的TCC矩阵,例如光瞳函数、源强度、相干因子或偏振,但本公开并不限于此。
在步骤S404中,处理器36产生包括至少一个向量的取样矩阵,所述至少一个向量充当形成特征矩阵的低秩基底的输入。取样矩阵包括例如随机矩阵、伪随机矩阵或结构矩阵。在一些实施例中,当产生取样矩阵时,处理器36可执行过取样(例如将秩k放大到k+O(1))以便改善随机取样误差。在一些实施例中,同时对向量组(即,矩阵)起作用以利用计算机硬件高速缓存效率,但在一些实施例中,还可以导出每次单个向量变型(single-vector-at-a-time variant),这可以使计算结果更准确。
在步骤S406中,处理器36将取样矩阵迭代地乘以特征矩阵,并且根据数值稳定性自适应地重新按比例调整相乘结果直到满足收敛条件。重复应用相乘以改善频谱近似误差。
在一些实施例中,在不直接计算全特征矩阵的情况下,通过利用特性的子集来执行与特征矩阵的相乘。举例来说,通过执行与TCC矩阵T的相乘为Tx=P(P*x),计算可绕过T的构建(包含O(N3)个算子)并且至多包含O(N2)个算子,这给出算法的总O(N2)成本。还可以利用其它结构(例如低秩、稀疏或张量),并且本公开并不限于此。
在一些实施例中,当将多个源点铺设在共形均匀网格上时,通过使用快速傅里叶变换(fast Fourier transform,FFT)计算根据特性的光瞳矩阵与所述光瞳矩阵的共轭转置的卷积来执行与特征矩阵的相乘。光瞳矩阵的卷积结构允许使用FFT在O(N log N)时间内相乘,并给出算法的总O(N log N)。
在一些实施例中,当将任意分布的多个源点重新取样到均匀网格中时,通过使用FFT计算根据特性的光瞳矩阵与所述光瞳矩阵的共轭转置的卷积来执行与特征矩阵的相乘。可将源点重新取样到具有可控空中图像误差的均匀网格上,并且如上文所描述的一样给出算法的总O(N log N)。
在一些实施例中,当任意分布的多个源点不在均匀网格上时,通过使用非均匀FFT计算根据特性的光瞳矩阵与所述光瞳矩阵的共轭转置的卷积来执行与特征矩阵的相乘。即使对于任意分布的源点,也可通过使用非均匀FFT来实现算法的总O(N log N)。
在一些实施例中,除了明确考虑到具有提供快于O(N3)相乘的结构T=PP*的TCC矩阵之外,也可以应用具有更一般结构(例如对于对角不定(即,可为负的)矩阵D的T=PDP*)的TCC矩阵。
在步骤S408中,处理器36使用迭代结果来形成简化特征矩阵。
在一些实施例中,通过子空间迭代获得简化特征矩阵,使用所述矩阵来计算近似于主要特征空间的低秩基底。收敛可与执行固定数量的迭代相关或与监测子空间相似性的一些度量相关。也就是说,在一些实施例中,当迭代次数达到预定次数时确定满足收敛条件,而在一些实施例中,当监测自适应误差度量时确定满足收敛条件,但本公开并不限于此。因此,通过使用迭代的最终结果来形成简化特征矩阵。
在一些实施例中,通过块克雷洛夫方法获得简化特征矩阵,在所述简化特征矩阵中,中间结果以与传统兰索斯(Lanczos)方法类似的方式更好地用于构建可能比子空间迭代更大的基底。也就是说,通过初始化空基底并将每次迭代的中间结果累加到所述空基底来形成简化特征矩阵。与子空间迭代相比,近似误差严格来说并不太差,但代价是简化特征矩阵可能更大,且计算分解值的成本可能更高。
在步骤S410中,处理器36计算简化特征矩阵的分解值,并且在步骤S412中,处理器36从计算出的分解值中提取多个特征内核。在一些实施例中,通过使用块算法(例如特征值分解(EVD)算法或奇异值分解(SVD)算法)来计算简化特征矩阵的分解值,这些分解值分别在下文详细说明。
图5A是示出根据本公开的实施例的用于通过使用EVD算法计算光学模型模拟的特征内核的方法的流程图。图5B是示出根据本公开的实施例的通过使用EVD算法计算光学模型模拟的特征内核的实例的示意图。本发明实施例的方法适用于图3的电子设备30,并且下文参看图3的电子设备30中的各种组件描述所述方法的详细步骤。
参看图5A和图5B两者,在步骤S502中,处理器36获得光学系统的光瞳矩阵,以便计算在数学上描述光学系统的特性的TCC矩阵T。
在步骤S504中,借用来自随机数值线性代数的技术,处理器36发现大致保留相关内核的大小为O(k)的低秩基底Q。
在步骤S506中,处理器36将低秩基底Q乘以TCC矩阵T以形成简化的TCC矩阵在所述矩阵中,将低秩基底Q的共轭转置矩阵Q*乘以TCC矩阵T和低秩基底Q,以便将TCC矩阵T压缩成简化TCC矩阵/>
在步骤S508中,处理器36根据简化TCC矩阵计算简化EVD。
在步骤S510中,处理器36从计算出的简化EVD中提取内核。在一些实施例中,处理器36从计算出的简化EVD中提取简化TCC内核(即,),然后通过将低秩基底Q与简化TCC内核相乘来将所提取的简化TCC内核扩展到全空间,以便恢复TCC内核(即,/>)。
根据本发明实施例,发现低秩基底,并且使用所述低秩基底将TCC矩阵压缩成小简化矩阵(例如大小为k阶乘k阶的TCC矩阵,其中k是相关内核的数量)。在较小压缩空间中执行计算,并且将结果扩展回原始空间。因此,可显著降低计算复杂度,并且可加快TCC内核的产生。
在一些实施例中,可通过如上文所引入的子空间迭代或块克雷洛夫方法来执行用于发现低秩基底的方法。
详细地说,图6A是示出根据本公开的实施例的用于通过子空间迭代发现低秩基底的方法的流程图。图6B是示出根据本公开的实施例的用于通过块克雷洛夫方法发现低秩基底的方法的流程图。本发明实施例的方法适用于图3的电子设备30。
参看图6A,在步骤S602中,处理器36产生过取样随机矩阵。在一些实施例中,处理器36可通过对随机矩阵、伪随机矩阵或结构矩阵进行过取样来产生过取样随机矩阵。
在步骤S604中,处理器36将所产生的随机矩阵乘以TCC矩阵,并且在步骤S606中,根据数值稳定性自适应地重新按比例调整相乘结果。重新按比例调整被称为例如格莱姆-施密特正交化(Gram-Schmidt orthogonalization),其主要用于数值稳定性,但本公开并不限于此。
在步骤S608中,处理器36确定迭代是否完成。在一些实施例中,可基于设定迭代次数或一些自适应误差度量来使终止(即,收敛)固定。
如果不满足收敛条件,那么处理器36返回到步骤S604以将所产生的随机矩阵迭代地乘以TCC矩阵,并且再次重新按比例调整相乘结果直到满足收敛条件。一旦满足收敛条件,则处理器36继续进行到步骤S610以使重新按比例调整的相乘结果正交化并最终输出正交化的结果作为低秩基底。
参看图6B,在步骤S612中,处理器36初始化空基底,且随后在步骤S614中产生过取样随机矩阵。在一些实施例中,处理器36可通过对随机矩阵、伪随机矩阵或结构矩阵进行过取样来产生过取样随机矩阵。
在步骤S616中,处理器36将所产生的随机矩阵乘以TCC矩阵,并且在步骤S618中,根据数值稳定性自适应地重新按比例调整相乘结果,并将重新按比例调整的相乘结果与在步骤S612中进行初始化的空基底相加。重新按比例调整被称为例如格莱姆-施密特正交化,其主要用于数值稳定性,但本公开并不限于此。
在步骤S620中,处理器36确定迭代是否完成。在一些实施例中,可基于设定迭代次数或一些自适应误差度量来使终止(即,收敛)固定。
如果不满足收敛条件,那么处理器36返回到步骤S616以将所产生的随机矩阵迭代地乘以TCC矩阵,重新按比例调整相乘结果并且将重新按比例调整的相乘与累加基底相加直到满足收敛条件。一旦满足收敛条件,则处理器36继续进行到步骤S622以使累加基底正交化并最终输出正交化的基底作为低秩基底。
与子空间迭代相比,对于块克雷洛夫方法,累加中间结果以构建可能比子空间迭代更大的基底。因此,近似误差严格来说并不太差,但用于计算EVD的成本可能更昂贵,这是由于所产生的基底更大。
在图5A到图6B中所说明的前述方法中,使用EVD算法计算特征内核。然而,在一些实施例中,可使用SVD算法来计算特征内核。下文给出实施例以进行进一步说明。
图7A是示出根据本公开的实施例的用于通过使用SVD算法计算光学模型模拟的特征内核的方法的流程图。图7B是示出根据本公开的实施例的通过使用SVD算法计算光学模型模拟的特征内核的实例的示意图。本发明实施例的方法适用于图3的电子设备30,并且下文参看图3的电子设备30中的各种组件描述所述方法的详细步骤。
参看图7A和图7B两者,在步骤S702中,处理器36获得光学系统的光瞳矩阵P。
在步骤S704中,借用来自随机数值线性代数的技术,处理器36发现大致保留相关内核的大小为O(k)的低秩基底Q。
在步骤S706中,处理器36将低秩基底Q的共轭转置矩阵Q*乘以光瞳矩阵P以形成简化光瞳矩阵以便将光瞳矩阵P压缩成简化光瞳矩阵/>
在步骤S708中,处理器36根据简化光瞳矩阵计算简化SVD。对于具有特征分解/>的简化TCC矩阵,对于某个幺正矩阵/>对应的简化光瞳矩阵SVD是/>
在步骤S710中,处理器36从计算出的简化SVD中提取内核。在一些实施例中,处理器36从计算出的简化SVD中提取简化TCC内核(即,),然后通过将低秩基底Q与简化TCC内核相乘来将所提取的简化TCC内核扩展到全空间,以便恢复TCC内核(即,/>)。
在一些实施例中,可通过如上文所引入的子空间迭代或块克雷洛夫方法来执行用于发现低秩基底的方法。
具体地说,在一些实施例中,可通过子空间迭代获得简化特征矩阵,使用所述矩阵来计算近似于主要特征空间的低秩基底。收敛可以与执行固定数量的迭代相关或与监测子空间相似性的一些度量相关。因此,可通过将所产生的随机矩阵迭代地乘以TCC矩阵(即,通过光瞳矩阵和其共轭转置),重新按比例调整相乘结果并使重新按比例调整的相乘结果正交化来发现低秩基底。
在一些实施例中,可通过块克雷洛夫方法获得简化特征矩阵,在所述简化特征矩阵中,中间结果更好地用于构建可能比子空间迭代更大的基底。因此,可通过初始化空基底、将每次迭代的中间结果累加到所述空基底并使累加基底正交化来发现低秩基底。
在本发明实施例的光瞳矩阵P的快速低秩SVD算法中,由于光瞳矩阵P的频谱衰减仅是TCC矩阵T的频谱衰减的二分之一,取决于频谱衰减的误差可较大。
在上述实施例中,分别引入基于快速EVD的算法和基于快速SVD的算法以提取TCC内核,但在一些实施例中,相同构想还可用以产生优于传统分解的各种快速低秩因式分解,从而产生具有例如简化存储器复杂度的“非典型”TCC内核。
在一些实施例中,提供一种包括处理器可执行指令的非暂时性计算机可读记录介质,所述处理器可执行指令在经执行时执行如以上实施例中所说明的用于计算光学模型模拟的特征内核的方法。在一些实施例中,非暂时性计算机可读记录介质是其上编码有计算机可读数据的CD-R、DVD-R、快闪驱动器或硬盘驱动器的盘片等。计算机可读数据(例如包括多个零和一的二进制数据)又包括配置以根据本文中所阐述得原理中的一或多种操作的计算机指令集。在一些实施例中,例如,处理器可执行计算机指令配置以执行用于计算光学模型模拟的特征内核的方法,例如图4中说明的例示性方法中的至少一些。本领域普通技术人员可设计经配置以根据本文中所提出的技术操作的许多此类计算机可读记录介质。
根据一些实施例,提供一种适用于具有存储介质和处理器的电子设备的用于计算光学模型模拟的特征内核的方法。在所述方法中,通过处理器识别在数学上描述光学成像系统的多种特性的特征矩阵,所述多种特性从光学成像系统中检索并且存储在存储介质中。通过处理器产生包括至少一个向量的取样矩阵,所述至少一个向量充当形成特征矩阵的低秩基底的输入。由处理器将取样矩阵迭代地乘以特征矩阵,并且根据数值稳定性自适应地重新按比例调整相乘结果直到满足收敛条件。由处理器使用迭代结果来形成简化特征矩阵。通过处理器计算简化特征矩阵的分解值,并且从计算出的分解值中提取多个特征内核。
在相关实施例中,所述简化特征矩阵由所述迭代的最终结果形成。
在相关实施例中,通过初始化空基底并且将每次迭代的中间结果累加到所述空基底来形成所述简化特征矩阵。
在相关实施例中,所述取样矩阵包括随机矩阵、伪随机矩阵或结构矩阵,且产生所述取样矩阵的步骤还包括:执行过取样以改进随机取样误差。
在相关实施例中,所述收敛条件在迭代次数达到预定次数或自适应误差度量受监测时得以满足。
在相关实施例中,在不直接计算全特征矩阵的情况下,通过利用所述特性的子集来执行与所述特征矩阵的相乘。
在相关实施例中,当将多个源点铺设在共形均匀网格上时,通过使用快速傅里叶变换(FFT)计算根据所述特性的光瞳矩阵与所述光瞳矩阵的共轭转置的卷积来执行与所述特征矩阵的相乘。
在相关实施例中,当将任意分布的多个源点重新取样到均匀网格上时,通过使用快速傅里叶变换计算根据所述特性的光瞳矩阵与所述光瞳矩阵的共轭转置的卷积来执行与所述特征矩阵的相乘。
在相关实施例中,当任意分布的多个源点不在均匀网格上时,通过使用不均匀快速傅里叶变换计算根据所述特性的光瞳矩阵与所述光瞳矩阵的共轭转置的卷积来执行与所述特征矩阵的相乘。
在相关实施例中,所述计算所述简化特征矩阵的分解值的步骤包括:通过使用块算法来计算所述简化特征矩阵的所述分解值,所述块算法包括特征值分解(EVD)算法或奇异值分解(SVD)算法。
根据一些实施例,用于计算光学模型模拟的特征内核的设备包含数据检索装置、存储介质以及处理器。数据检索装置配置以检索光学成像系统的多种特性。存储介质配置以存储由数据检索装置检索到的特性。处理器耦合到数据检索装置和存储介质,并且配置以执行用以进行以下步骤的指令:识别在数学上描述光学成像系统的多种特性的特征矩阵;产生包括至少一个向量的取样矩阵,所述至少一个向量充当形成特征矩阵的低秩基底的输入;将取样矩阵迭代地乘以特征矩阵,并且根据数值稳定性自适应地重新按比例调整相乘结果直到满足收敛条件;使用迭代结果以形成简化特征矩阵;以及计算简化特征矩阵的分解值,并且从计算出的分解值中提取多个特征内核。
在相关实施例中,所述简化特征矩阵由所述迭代的最终结果形成。
在相关实施例中,通过初始化空基底并且将每次迭代的中间结果累加到所述空基底来形成所述简化特征矩阵。
在相关实施例中,所述取样矩阵包括随机矩阵、伪随机矩阵或结构矩阵,且所述处理器还执行过取样以改进随机取样误差。
在相关实施例中,所述收敛条件在迭代次数达到预定次数或自适应误差度量受监测时得以满足。
在相关实施例中,在不直接计算全特征矩阵的情况下,通过所述处理器利用所述特性的子集来执行与所述特征矩阵的相乘。
在相关实施例中,当将多个源点铺设在共形均匀网格上时,通过所述处理器使用快速傅里叶变换计算根据所述特性的光瞳矩阵与所述光瞳矩阵的共轭转置的卷积来执行与所述特征矩阵的相乘。
在相关实施例中,当将任意分布的多个源点重新取样到共形均匀网格上时,通过所述处理器使用快速傅里叶变换计算根据所述特性的光瞳矩阵与所述光瞳矩阵的共轭转置的卷积来执行与所述特征矩阵的相乘。
在相关实施例中,当任意分布的多个源点不在均匀网格上时,通过所述处理器使用不均匀快速傅里叶变换计算根据所述特性的光瞳矩阵与所述光瞳矩阵的共轭转置的卷积来执行与所述特征矩阵的相乘。
根据一些实施例,一种非暂时性计算机可读记录介质包含处理器可执行指令,所述处理器可执行指令在经执行时执行用于计算光学模型模拟的特征内核的方法。在方法中,识别在数学上描述光学成像系统的多种特性的特征矩阵。产生包括至少一个向量的取样矩阵,所述至少一个向量充当形成特征矩阵的低秩基底的输入。将取样矩阵迭代地乘以特征矩阵,并且根据数值稳定性自适应地重新按比例调整相乘结果直到满足收敛条件。使用迭代结果来形成简化特征矩阵。计算简化特征矩阵的分解值,并且从计算出的分解值中提取多个特征内核。
前文概述若干实施例的特征以使本领域的技术人员可更好地理解本公开的各方面。本领域的技术人员应了解,其可以易于使用本公开作为设计或修改用于进行本文中所介绍的实施例的相同目的和/或获得相同优势的其它方法和结构的基础。本领域的技术人员还应认识到,此类等效构造并不脱离本公开的精神和范围,且本领域的技术人员可在不脱离本公开的精神和范围的情况下在本文中进行各种改变、替代以及更改。
Claims (20)
1.一种用于计算光学模型模拟的特征内核的方法,适用于具有存储介质及处理器的电子设备,所述方法包括:
通过所述处理器识别在数学上描述光学成像系统的多种特性的特征矩阵,所述多种特性从所述光学成像系统中检索并且存储在所述存储介质中;
通过所述处理器产生包括至少一个向量的取样矩阵,所述至少一个向量充当形成所述特征矩阵的低秩基底的输入;
通过所述处理器将所述取样矩阵乘以所述特征矩阵以获得相乘结果,并且根据所述相乘结果的数值稳定性自适应地重新按比例调整所述相乘结果,以及将重新按比例调整的所述相乘结果迭代地乘以所述特征矩阵以获得当前迭代的相乘结果并且重新按比例调整所述相乘结果直到满足收敛条件;
通过所述处理器使用重新按比例调整的所述相乘结果以形成简化特征矩阵;以及
通过所述处理器计算所述简化特征矩阵的分解值,并且从计算出的所述分解值中提取多个特征内核。
2.根据权利要求1所述的方法,其中所述简化特征矩阵由迭代的最终结果形成。
3.根据权利要求1所述的方法,其中通过初始化空基底并且将每次迭代的中间结果累加到所述空基底来形成所述简化特征矩阵。
4.根据权利要求1所述的方法,其中所述取样矩阵包括随机矩阵、伪随机矩阵或结构矩阵,且产生所述取样矩阵的步骤还包括:
执行过取样以改进随机取样误差。
5.根据权利要求1所述的方法,其中所述收敛条件在迭代次数达到预定次数或自适应误差度量受监测时得以满足。
6.根据权利要求1所述的方法,其中在不直接计算全特征矩阵的情况下,通过利用所述特性的子集来执行与所述特征矩阵的相乘。
7.根据权利要求1所述的方法,其中当将多个源点铺设在共形均匀网格上时,通过使用快速傅里叶变换计算根据所述特性的光瞳矩阵与所述光瞳矩阵的共轭转置的卷积来执行与所述特征矩阵的相乘。
8.根据权利要求1所述的方法,其中当将任意分布的多个源点重新取样到均匀网格上时,通过使用快速傅里叶变换计算根据所述特性的光瞳矩阵与所述光瞳矩阵的共轭转置的卷积来执行与所述特征矩阵的相乘。
9.根据权利要求1所述的方法,其中当任意分布的多个源点不在均匀网格上时,通过使用不均匀快速傅里叶变换计算根据所述特性的光瞳矩阵与所述光瞳矩阵的共轭转置的卷积来执行与所述特征矩阵的相乘。
10.根据权利要求1所述的方法,其中计算所述简化特征矩阵的所述分解值的步骤包括:
通过使用块算法来计算所述简化特征矩阵的所述分解值,所述块算法包括特征值分解算法或奇异值分解算法。
11.一种用于计算光学模型模拟的特征内核的设备,包括:
数据检索装置,配置以检索光学成像系统的多种特性;
存储介质,配置以存储由所述数据检索装置检索到的所述特性;
处理器,耦合到所述数据检索装置及所述存储介质,并且配置以执行用以执行以下步骤的指令:
识别在数学上描述所述光学成像系统的所述多种特性的特征矩阵;
产生包括至少一个向量的取样矩阵,所述至少一个向量充当形成所述特征矩阵的低秩基底的输入;
将所述取样矩阵乘以所述特征矩阵以获得相乘结果,并且根据所述相乘结果的数值稳定性自适应地重新按比例调整所述相乘结果,以及将重新按比例调整的所述相乘结果迭代地乘以所述特征矩阵以获得当前迭代的相乘结果并且重新按比例调整所述相乘结果直到满足收敛条件;
使用重新按比例调整的所述相乘结果以形成简化特征矩阵;以及
计算所述简化特征矩阵的分解值,并且从计算出的所述分解值中提取多个特征内核。
12.根据权利要求11所述的用于计算光学模型模拟的特征内核的设备,其中所述简化特征矩阵由迭代的最终结果形成。
13.根据权利要求11所述的用于计算光学模型模拟的特征内核的设备,其中通过初始化空基底并且将每次迭代的中间结果累加到所述空基底来形成所述简化特征矩阵。
14.根据权利要求11所述的用于计算光学模型模拟的特征内核的设备,其中所述取样矩阵包括随机矩阵、伪随机矩阵或结构矩阵,且所述处理器还执行过取样以改进随机取样误差。
15.根据权利要求11所述的用于计算光学模型模拟的特征内核的设备,其中所述收敛条件在迭代次数达到预定次数或自适应误差度量受监测时得以满足。
16.根据权利要求11所述的用于计算光学模型模拟的特征内核的设备,其中在不直接计算全特征矩阵的情况下,通过所述处理器利用所述特性的子集来执行与所述特征矩阵的相乘。
17.根据权利要求11所述的用于计算光学模型模拟的特征内核的设备,其中当将多个源点铺设在共形均匀网格上时,通过所述处理器使用快速傅里叶变换计算根据所述特性的光瞳矩阵与所述光瞳矩阵的共轭转置的卷积来执行与所述特征矩阵的相乘。
18.根据权利要求11所述的用于计算光学模型模拟的特征内核的设备,其中当将任意分布的多个源点重新取样到共形均匀网格上时,通过所述处理器使用快速傅里叶变换计算根据所述特性的光瞳矩阵与所述光瞳矩阵的共轭转置的卷积来执行与所述特征矩阵的相乘。
19.根据权利要求11所述的用于计算光学模型模拟的特征内核的设备,其中当任意分布的多个源点不在均匀网格上时,通过所述处理器使用不均匀快速傅里叶变换计算根据所述特性的光瞳矩阵与所述光瞳矩阵的共轭转置的卷积来执行与所述特征矩阵的相乘。
20.一种非暂时性计算机可读记录介质,包括处理器可执行指令,所述处理器可执行指令在经执行时执行用于计算光学模型模拟的特征内核的方法,所述方法包括:
识别在数学上描述光学成像系统的多种特性的特征矩阵;
产生包括至少一个向量的取样矩阵,所述至少一个向量充当形成所述特征矩阵的低秩基底的输入;
将所述取样矩阵乘以所述特征矩阵以获得相乘结果,并且根据所述相乘结果的数值稳定性自适应地重新按比例调整所述相乘结果,以及将重新按比例调整的所述相乘结果迭代地乘以所述特征矩阵以获得当前迭代的相乘结果并且重新按比例调整所述相乘结果直到满足收敛条件;
使用重新按比例调整的所述相乘结果以形成简化特征矩阵;以及
计算所述简化特征矩阵的分解值,并且从计算出的所述分解值中提取多个特征内核。
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201862725271P | 2018-08-31 | 2018-08-31 | |
US62/725,271 | 2018-08-31 | ||
US16/423,190 | 2019-05-28 | ||
US16/423,190 US10809629B2 (en) | 2018-08-31 | 2019-05-28 | Method and apparatus for computing feature kernels for optical model simulation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110929366A CN110929366A (zh) | 2020-03-27 |
CN110929366B true CN110929366B (zh) | 2024-05-28 |
Family
ID=69640665
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910773551.4A Active CN110929366B (zh) | 2018-08-31 | 2019-08-21 | 用于计算光学模型模拟的特征内核的方法 |
Country Status (3)
Country | Link |
---|---|
US (2) | US10809629B2 (zh) |
CN (1) | CN110929366B (zh) |
TW (1) | TWI747036B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111856872B (zh) * | 2020-08-03 | 2021-11-30 | 中国科学院上海光学精密机械研究所 | 基于双采样率像素化掩模图形快速光学邻近效应修正方法 |
US11435670B1 (en) | 2021-02-26 | 2022-09-06 | Taiwan Semiconductor Manufacturing Co., Ltd. | Multi-component kernels for vector optical image simulation |
CN115390372B (zh) * | 2022-09-29 | 2023-04-11 | 武汉宇微光学软件有限公司 | 一种光刻系统掩模成像方法、装置及系统 |
CN116415562B (zh) * | 2023-06-06 | 2023-08-25 | 上海朝阳永续信息技术股份有限公司 | 用于解析金融数据的方法、设备和介质 |
CN117170175B (zh) * | 2023-11-03 | 2024-01-30 | 合肥晶合集成电路股份有限公司 | 光学临近效应修正建模方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1770164A (zh) * | 2005-09-14 | 2006-05-10 | 浙江大学 | 一种基于卷积核的集成电路光刻制造建模方法 |
CN101013195A (zh) * | 2007-01-10 | 2007-08-08 | 中国科学院光电技术研究所 | 基于线性相位反演复原技术的自适应光学系统 |
US7664713B1 (en) * | 2006-12-31 | 2010-02-16 | At&T Corp. | Method and apparatus for providing fast kernel learning on sparse data |
CN104977797A (zh) * | 2014-04-02 | 2015-10-14 | 中芯国际集成电路制造(上海)有限公司 | 光学邻近校正方法以及优化光学邻近校正模型的方法 |
CN107644131A (zh) * | 2017-09-13 | 2018-01-30 | 苏州珂晶达电子有限公司 | 基于极坐标采样的交叉传递函数快速分解方法 |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7030997B2 (en) * | 2001-09-11 | 2006-04-18 | The Regents Of The University Of California | Characterizing aberrations in an imaging lens and applications to visual testing and integrated circuit mask analysis |
US20050015233A1 (en) * | 2003-07-17 | 2005-01-20 | International Business Machines Corporation | Method for computing partially coherent aerial imagery |
US7079223B2 (en) * | 2004-02-20 | 2006-07-18 | International Business Machines Corporation | Fast model-based optical proximity correction |
JP5358780B2 (ja) * | 2004-09-11 | 2013-12-04 | ザ ボード オブ トラスティーズ オブ ザ レランド スタンフォード ジュニア ユニバーシティー | 光学導波管のモード特性をモデリングするための方法および装置 |
KR100634527B1 (ko) * | 2004-11-26 | 2006-10-16 | 삼성전자주식회사 | 층 기반 영상 처리 장치 및 방법 |
US7634137B2 (en) * | 2005-10-14 | 2009-12-15 | Microsoft Corporation | Unfolded convolution for fast feature extraction |
US7933471B2 (en) * | 2007-03-09 | 2011-04-26 | Synopsys, Inc. | Method and system for correlating physical model representation to pattern layout |
US10025198B2 (en) * | 2008-07-07 | 2018-07-17 | Asml Netherlands B.V. | Smart selection and/or weighting of parameters for lithographic process simulation |
TW201013339A (en) * | 2008-08-07 | 2010-04-01 | Nikon Corp | Substrate processing method, substrate processing equipment, exposure method and device fabrication method |
US8285030B2 (en) * | 2010-03-15 | 2012-10-09 | Synopsys, Inc. | Determining calibration parameters for a lithographic process |
US8875066B2 (en) * | 2013-03-15 | 2014-10-28 | Synopsys, Inc. | Performing image calculation based on spatial coherence |
US9658987B2 (en) * | 2014-05-15 | 2017-05-23 | International Business Machines Corporation | Regression using M-estimators and polynomial kernel support vector machines and principal component regression |
US10430543B2 (en) * | 2014-10-04 | 2019-10-01 | Synopsys, Inc. | Matrix reduction for lithography simulation |
KR102109059B1 (ko) * | 2014-11-26 | 2020-05-12 | 에이에스엠엘 네델란즈 비.브이. | 계측 방법, 컴퓨터 제품 및 시스템 |
WO2017060192A1 (en) * | 2015-10-08 | 2017-04-13 | Asml Netherlands B.V. | Method and apparatus for pattern correction and verification |
US20180113840A1 (en) * | 2016-10-25 | 2018-04-26 | Wisconsin Alumni Research Foundation | Matrix Processor with Localized Memory |
US20190122111A1 (en) * | 2017-10-24 | 2019-04-25 | Nec Laboratories America, Inc. | Adaptive Convolutional Neural Knowledge Graph Learning System Leveraging Entity Descriptions |
-
2019
- 2019-05-28 US US16/423,190 patent/US10809629B2/en active Active
- 2019-08-21 CN CN201910773551.4A patent/CN110929366B/zh active Active
- 2019-08-30 TW TW108131310A patent/TWI747036B/zh active
-
2020
- 2020-09-22 US US17/028,981 patent/US11003092B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1770164A (zh) * | 2005-09-14 | 2006-05-10 | 浙江大学 | 一种基于卷积核的集成电路光刻制造建模方法 |
US7664713B1 (en) * | 2006-12-31 | 2010-02-16 | At&T Corp. | Method and apparatus for providing fast kernel learning on sparse data |
CN101013195A (zh) * | 2007-01-10 | 2007-08-08 | 中国科学院光电技术研究所 | 基于线性相位反演复原技术的自适应光学系统 |
CN104977797A (zh) * | 2014-04-02 | 2015-10-14 | 中芯国际集成电路制造(上海)有限公司 | 光学邻近校正方法以及优化光学邻近校正模型的方法 |
CN107644131A (zh) * | 2017-09-13 | 2018-01-30 | 苏州珂晶达电子有限公司 | 基于极坐标采样的交叉传递函数快速分解方法 |
Non-Patent Citations (3)
Title |
---|
Digital equalization of optical nonlinearities in very high-speed optical communication systems;Armando N. Pinto 等;《2014 16th International Conference on Transparent Optical Networks(ICTON)》;20140814;1-5 * |
Fast Kernel Feature Extraction Method and Its Application;Xu Liang 等;《Computer Engineering》;20091231;第35卷(第24期);26-8, 32 * |
Selection of convolution kernel in no-uniform fast Fourier transform for Fourier domain optical coherence tomography;Kenny K.H. Chan 等;《Optics Express》;20111219;第19卷(第27期);26891-26904 * |
Also Published As
Publication number | Publication date |
---|---|
TW202011110A (zh) | 2020-03-16 |
US10809629B2 (en) | 2020-10-20 |
US20200073249A1 (en) | 2020-03-05 |
US20210003923A1 (en) | 2021-01-07 |
CN110929366A (zh) | 2020-03-27 |
US11003092B2 (en) | 2021-05-11 |
TWI747036B (zh) | 2021-11-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110929366B (zh) | 用于计算光学模型模拟的特征内核的方法 | |
KR102120523B1 (ko) | 프로세스 유도된 왜곡 예측 및 오버레이 에러의 피드포워드 및 피드백 보정 | |
US7921383B1 (en) | Photolithographic process simulation including efficient result computation for multiple process variation values | |
JP5960953B2 (ja) | 単一のプロセスウィンドウモデルを生成するためのシステム | |
US8285030B2 (en) | Determining calibration parameters for a lithographic process | |
TWI640934B (zh) | 基於預測模型化之聚焦誤差預測之方法及系統 | |
CN111627799B (zh) | 制造半导体元件的方法 | |
US20100119961A1 (en) | Methods and system for lithography calibration | |
CN110678961A (zh) | 在光学光刻中模拟近场图像 | |
US8136054B2 (en) | Compact abbe's kernel generation using principal component analysis | |
US10678971B2 (en) | Space exploration with Bayesian inference | |
WO2022105062A1 (zh) | 用于光学量测的方法、系统、计算设备和存储介质 | |
WO2020054443A1 (ja) | 学習装置、推論装置及び学習済みモデル | |
CN113238460B (zh) | 一种基于深度学习的用于超紫外的光学邻近校正方法 | |
KR20190039579A (ko) | 노광후 공정들의 모델링 | |
CN113779779A (zh) | 优化掩模的方法、设备和计算机可读存储介质 | |
CN114065693A (zh) | 超大规模集成电路结构布局的优化方法、系统和电子设备 | |
Kareem et al. | Layout pattern synthesis for lithography optimizations | |
US20240126179A1 (en) | Optical imaging method, device and system for photolithography system | |
KR102498694B1 (ko) | 전산 리소그래피 마스크 모델과 관련된 전자계를 결정하는 방법 | |
CN117452779B (zh) | 一种计算光刻非线性系统的建模、标定方法和装置 | |
Hordan et al. | Weisfeiler Leman for Euclidean Equivariant Machine Learning | |
WO2024100778A1 (ja) | 第一原理分子動力学法を計算するための量子情報処理方法 | |
Partlo et al. | Optimizing NA and sigma for subhalf-micrometer lithography | |
WO2024012800A1 (en) | Systems and methods for predicting post-etch stochastic variation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |