數(shù)值計(jì)算模擬導(dǎo)論

講座模擬導(dǎo)論問題講座模擬導(dǎo)論問題實(shí)例及基本JacobDeepakRamaswamy,MichalLucaDaniel,ShihhsienKuoandKarenSMA-HPC?2003模擬的應(yīng)——模擬的應(yīng)——————課程原問題實(shí)——————SMA-HPC?2003方近來發(fā)方近來發(fā)SMA-HPC?2003方——近來發(fā)方——近來發(fā)——SMA-HPC?2003方——機(jī)械單元（板、梁、殼）——方——機(jī)械單元（板、梁、殼）——近來發(fā)——SMA-HPC?2003方近來方近來發(fā)——SMA-HPC?2003方——近來發(fā)方——近來發(fā)——SMA-HPC?2003方——近方——近來發(fā)SMA-HPC?2003股票價(jià) 價(jià)格選 方股票價(jià) 價(jià)格選 方近來發(fā)——SMA-HPC?2003套期保值基金的股票定價(jià)方——近來發(fā)方——近來發(fā)——SMA-HPC?2003方——近來發(fā)方——近來發(fā)——SMA-HPC?2003方——近來發(fā)——方——近來發(fā)——SMA-HPC?2003白質(zhì)分子（子SMA-HPC?2003SMA-HPC?2003對一些現(xiàn)代技術(shù)的回對一些現(xiàn)代技術(shù)的回為什么證明定理SMA-HPC?2003 主電源線 至少有3V 主電源線 至少有3V通過ALU嗎SMA-HPC?2003SMA-HPC?2003車空間框架在載荷作用下車空間框架在載荷作用下下垂太多SMA-HPC?2003地面固貨梁接下圖示為將貨物（紅色）圖示為將貨物（紅色）SMA-HPC?2003發(fā)動機(jī)變得太熱嗎發(fā)動機(jī)變得太熱嗎SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003選擇拓?fù)浣Y(jié)構(gòu)和選擇拓?fù)浣Y(jié)構(gòu)和支柱的長度與寬度，滿足SMA-HPC?2003選擇形狀，滿足選擇形狀，滿足SMA-HPC?2003給定拓?fù)浣Y(jié)構(gòu)和金給定拓?fù)浣Y(jié)構(gòu)和金屬的長寬，確定SMA-HPC?2003IBM,Motorola,TI,Intel,Compaq,IBM,Motorola,TI,Intel,Compaq,Sony,SMA-HPC?2003一些大公司SMA-HPC?2003上千個(gè)小型公小型公司應(yīng)用上千個(gè)小型公小型公司應(yīng)用于磁盤驅(qū)動器電路、圖形加速卡、CD放器，蜂窩電話非功能原型的成本是多少SMA-HPC?2003SMA-HPC?200315638156386914MentorGraphics這些公司通過改善分析效率進(jìn)行競爭SMA-HPC?2003誰在設(shè)計(jì)VLSI工具？SMA-HPC?2003SMA-HPC?2003們可以構(gòu)建電路來確定各功能模塊的電壓，并確定VLSI電路是否SMA-HPC?2003電一電壓電路電 電路電 SMA-HPC?2003電一電壓電路電 電路電 SMA-HPC?2003SMA-HPC?2003功能模塊看作是電壓源電路單 SMA-HPC?2003功能模塊看作是電壓源電路單 SMA-HPC?2003SMA-HPC?2003金屬線看作電 物理符 電路模 長 電阻 金屬線看作電 物理符 電路模 長 電阻 SMA-HPC?2003線的長度（電流需進(jìn)一步通過），（電流有更多的路徑選擇線的長度（電流需進(jìn)一步通過），（電流有更多的路徑選擇 SMA-HPC?2003綜合在一起結(jié)果如示意圖所綜合在一起結(jié)果如示意圖所SMA-HPC?2003SMA-HPC?2003實(shí)例的簡化如圖所示SMA-HPC實(shí)例的簡化如圖所示SMA-HPC?2003負(fù)接壓螺SMA-HPC?2003載荷看成力SMA-HPC?2003載荷看成力SMA-HPC?2003力的力的模型可看作沿著Y負(fù)向的拉力（Y是垂直方向，X是水平方向）SMA-HPC?2003梁看作支柱 非受拉長 截面 梁看作支柱 非受拉長 截面 材料特 SMA-HPC?20031.2.1.2.不彎 SMA-HPC?20033.3.為了確定K，考慮梁伸長量為ΔL的力是）將（I）（II）（III）SMA-HPC?2003綜合在一SMA-HPC?2003綜合在一SMA-HPC?2003兩種類型的未知量兩種類型的未知量電路－每一節(jié)點(diǎn)的電流和為支柱－每一連接點(diǎn)的合力為基本方程SMA-HPC?2003熱問題：沿著桿的溫度分布是什么 SMA-HPC?2003輸入熱問題：沿著桿的溫度分布是什么 SMA-HPC?2003輸入熱離散表示SMA-HPC?2003熱離散表示SMA-HPC?2003熱離散表示 通過每一截面的熱 當(dāng)截面足夠小時(shí)取極 SMA-HPC?2003熱離散表示 通過每一截面的熱 當(dāng)截面足夠小時(shí)取極 SMA-HPC?2003熱離散表示兩個(gè)相鄰截 流入控制體積的純熱 SMA-HPC?2003控制體熱離散表示兩個(gè)相鄰截 流入控制體積的純熱 SMA-HPC?2003控制體輸入熱離散表示 流入控制體積的純熱 當(dāng)截面足夠小時(shí)取極 SMA-HPC?2003熱離散表示 流入控制體積的純熱 當(dāng)截面足夠小時(shí)取極 SMA-HPC?2003輸入熱電路模擬SMA-HPC?2003熱電路模擬SMA-HPC?2003支柱－每一連接點(diǎn)的合力為桿－進(jìn)入控制體積的熱流和SMA-HPC?2003電路實(shí)例識別未知量SMA-HPC?2003電路實(shí)例識別未知量SMA-HPC?2003給定一電給定一電路示意圖，問題是要確定節(jié)點(diǎn)電壓和單元電流。開始，節(jié)點(diǎn)編碼為0具有特殊意義，它是參考點(diǎn)。雖然電壓不是絕對值，但必須相對某一參考點(diǎn)度為了更好理解這一點(diǎn)，考慮一個(gè)簡單的電流源和電阻實(shí)例。為了使流過電阻的電流為1安培，V1－V0必須是1伏特。但究竟是＝0SMA-HPC?2003電路實(shí)例識別未知量SMA-HPC?2003電路實(shí)例識別未知量SMA-HPC?2003SMA-HPC?2003電路實(shí)例識別未知量SMA-HPC?2003電路實(shí)例識別未知量SMA-HPC?2003致）SMA-HPC?2003電路實(shí)例識別未知量SMA-HPC?2003電路實(shí)例識別未知量SMA-HPC?2003且此時(shí)還應(yīng)該注意電流的方向，如果電流通過左節(jié)點(diǎn)到右節(jié)，則左節(jié)點(diǎn)電壓應(yīng)比有節(jié)點(diǎn)高RI電路實(shí)例小電路實(shí)例的未知單元電流（除了電流源電路實(shí)例小電路實(shí)例的未知單元電流（除了電流源實(shí)例的方注意：方程數(shù)＝未知量SMA-HPC?2003模擬的許多應(yīng)模擬的許多應(yīng)兩種類型的方電路－每一節(jié)點(diǎn)的電流和為SMA-HPC?2003講座模擬導(dǎo)論講座模擬導(dǎo)論Jacob感謝DeepakRamaswamy,MichalRewienski,KarenSMA-HPC?2003從示意圖形成方從從示意圖形成方從示意圖構(gòu)建矩兩種表達(dá)方SMA-HPC?2003支柱實(shí)例識別未知量建立每一接頭的位置支柱實(shí)例識別未知量建立每一接頭的位置坐標(biāo)X，Y，將某一接頭設(shè)為零點(diǎn)。SMA-HPC?2003該符號通常用來表示一個(gè)固定結(jié)構(gòu)（如象水泥墻）SMA-HPC?2003支柱實(shí)例識別未知量給每一支柱設(shè)置X支柱實(shí)例識別未知量給每一支柱設(shè)置X，Y力分量SMA-HPC?2003支柱實(shí)例識別未知量SMA-HPC?2003支柱實(shí)例識別未知量SMA-HPC?2003如圖所示，支柱的拉力如圖所示，支柱的拉力f如圖示將局部坐標(biāo)映射到總統(tǒng)坐標(biāo)系計(jì)算fx和fy。從圖示幾何關(guān)系得到從ffxfy支柱實(shí)例守恒定律在每一接頭的X方向合力＝0支柱實(shí)例守恒定律在每一接頭的X方向合力＝0SMA-HPC?2003。如下圖，如果支柱伸長?力f力fafb大小相等方向相反。這是由于fa在X的正方向，而fb在X的反方向。其它力+fa=0。盡管右邊接頭的平衡方程是：其它力fb=其它力-fa0fafb作為獨(dú)立變量。而代之以所選力和建立支柱另一邊力間恒律中f3符號為正，但接頭2守恒律中為具有負(fù)號的反向SMA-HPC?2003支柱實(shí)例守恒定律SMA-HPC?2003支柱實(shí)例守恒定律SMA-HPC?2003于翻轉(zhuǎn)支柱的x1-x2>0或者對應(yīng)于壓縮支SMA-HPC?2003支柱實(shí)例總SMA-HPC支柱實(shí)例總SMA-HPC?2003沿著X軸的兩個(gè)支結(jié)點(diǎn)沿著X軸的兩個(gè)支結(jié)點(diǎn)1：f1xf2x節(jié)點(diǎn)2：-f2xf=SMA-HPC?2003 沿著X軸的兩個(gè) 沿著X軸的兩個(gè)支 基本方 SMA-HPC?2003沿著X軸的兩個(gè)支沿著X軸的兩個(gè)支簡化（節(jié)點(diǎn)）SMA-HPC?2003沿著X軸的兩個(gè)支節(jié)沿著X軸的兩個(gè)支節(jié)點(diǎn)方程的解fL=10(沿著X正方向的力SMA-HPC?2003沿著X軸的兩個(gè)支沿著X軸的兩個(gè)支注意力的符f2x=10(x正向的力f1x=?10(x負(fù)向的力SMA-HPC?2003上次的實(shí)例SMA-HPC?2003上次的實(shí)例SMA-HPC?2003兩類未知兩類未知兩類方電路－每一節(jié)點(diǎn)電流和SMA-HPC?2003電路實(shí)電路實(shí)SMA-HPC?2003二維支柱實(shí)二維支柱實(shí)JS對矩陣行表示單元基本方程（線性SMA-HPC?2003守恒方程電路實(shí)例SMA-HPC?2003守恒方程電路實(shí)例SMA-HPC?2003SMA-HPC?2003守恒方程電路實(shí)例SMA-HPC?2003守恒方程電路實(shí)例SMA-HPC?2003守恒方程電路實(shí)例 右邊對應(yīng)源電 矩陣A非方SMA-HPC守恒方程電路實(shí)例 右邊對應(yīng)源電 矩陣A非方SMA-HPC?2003一行表示每一守恒方程電路實(shí)例電阻對每一矩陣如何產(chǎn)生作守恒方程電路實(shí)例電阻對每一矩陣如何產(chǎn)生作A的每列不超過兩個(gè)非零元素SMA-HPC?2003當(dāng)電阻一端當(dāng)電阻一端連在參考點(diǎn)（0節(jié)點(diǎn)）時(shí)，矩陣會是怎SMA-HPC?2003守恒方程電路實(shí)例每一電流源對右邊的貢獻(xiàn)怎守恒方程電路實(shí)例每一電流源對右邊的貢獻(xiàn)怎SMA-HPC?2003守恒方程電路實(shí)例守恒矩陣方程產(chǎn)生算對每一電阻守恒方程電路實(shí)例守恒矩陣方程產(chǎn)生算對每一電阻令I(lǐng)s為零向?qū)τ诿恳浑娏鱏MA-HPC?2003守恒方程電路實(shí)例SMA-HPC?2003守恒方程電路實(shí)例SMA-HPC?2003基本方程電路實(shí)例首先確定電阻的電壓（支路電壓SMA-HPC?2003基本方程電路實(shí)例首先確定電阻的電壓（支路電壓SMA-HPC?2003電阻的電電阻的電壓是V1-V2，通過電阻的電流是：注意如果V1大于V2，則iSMA-HPC?2003基本方程電路實(shí)例SMA-HPC?2003基本方程電路實(shí)例SMA-HPC?2003由于有B個(gè)支路電壓和N個(gè)節(jié)點(diǎn)，所以矩陣有BSMA-HPC?2003基本方程電路實(shí)例SMA-HPC?2003基本方程電路實(shí)例SMA-HPC?2003對守恒律，如下矩陣所示支路k對對矩陣第k列兩個(gè)非零因素有貢注意到支路k的電壓是VlVm因此，正如下列矩陣，第k不能看出每一支路單元對影響矩陣的一列有貢SMA-HPC?2003基本方程電路實(shí)例基本方程電路實(shí)例SMA-HPC?2003基本方程電路實(shí)例SMA-HPC?2003基本方程電路實(shí)例SMA-HPC?2003節(jié)點(diǎn)支路關(guān)系電路實(shí)例SMA-HPC?2003節(jié)點(diǎn)支路關(guān)系電路實(shí)例SMA-HPC?2003例一對列表示每一未知——J——例一對列表示每一未知——J——S一對列表示每一方——J——SSMA-HPC?2003下列方法類似于電應(yīng)用A下列方法類似于電應(yīng)用ATSMA-HPC?2003守恒方程支柱實(shí)例SMA-HPC?2003守恒方程支柱實(shí)例SMA-HPC?2003SMA-HPC?2003守恒方程支柱實(shí)例 沖壓 SMA-HPC?2003守恒方程支柱實(shí)例 沖壓 SMA-HPC?2003力右邊表示力在接頭查看矩陣將產(chǎn)對于j1守恒方程的貢獻(xiàn)需要用專用下標(biāo)，我們用j1xj1y表示兩行，SxSySMA-HPC?2003守恒方程支柱實(shí)例如果（j1不固定如果（j2不固定如果（j1不固定守恒方程支柱實(shí)例如果（j1不固定如果（j2不固定如果（j1不固定A最多有兩個(gè)非零列SMA-HPC?2003基本方程支柱實(shí)例 如果x1,y1近似于x0,基本方程支柱實(shí)例 如果x1,y1近似于x0, SMA-HPC?2003如上所示，通過支柱的力如上所示，通過支柱的力，且xy其如果x和y偏離x0和y0則由于Fx(x0,y0)=同理類似表達(dá)式對y成立注意到就是不拉伸支柱，其旋轉(zhuǎn)將破壞小擾動條件，由于指向不正確的方向A-PC00很好逼近力基本方程支柱實(shí)例SMA-HPC?2003基本方程支柱實(shí)例SMA-HPC?2003基本方程支柱實(shí)例SMA-HPC?2003基本方程支柱實(shí)例SMA-HPC?2003基本方程節(jié)點(diǎn)支柱關(guān)系S＝J＝SMA-HPC?2003基本方程節(jié)點(diǎn)支柱關(guān)系S＝J＝SMA-HPC?2003比支柱實(shí)例SMA-HPC?2003比支柱實(shí)例SMA-HPC?2003電路實(shí)例產(chǎn)生矩陣計(jì)算節(jié)點(diǎn)數(shù)量，其中一個(gè)節(jié)點(diǎn)作為對每一節(jié)點(diǎn)寫出守恒律。除了0SMA-HPC?2003電路實(shí)例產(chǎn)生矩陣計(jì)算節(jié)點(diǎn)數(shù)量，其中一個(gè)節(jié)點(diǎn)作為對每一節(jié)點(diǎn)寫出守恒律。除了0SMA-HPC?2003基本矩陣電路實(shí)例SMA-HPC?2003基本矩陣電路實(shí)例SMA-HPC?2003節(jié)點(diǎn)n1的節(jié)點(diǎn)n1的節(jié)點(diǎn)n2的因此和Rk相關(guān)的矩陣是SMA-HPC?2003電路實(shí)例產(chǎn)生矩陣SMA-HPC?2003電路實(shí)例產(chǎn)生矩陣SMA-HPC?2003產(chǎn)生矩陣SMA-HPC?2003產(chǎn)生矩陣SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003G矩陣的特征SMA-HPC?2003G矩陣的特征SMA-HPC?2003節(jié)電支路形式節(jié)點(diǎn)支路公不對稱和對角節(jié)電支路形式節(jié)點(diǎn)支路公不對稱和對角占矩陣是(n+b)x(SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003問題單元電壓SMA-HPC?2003問題單元電壓SMA-HPC?2003問題單元電壓應(yīng)用電壓源能形成節(jié)點(diǎn)支路問題單元電壓應(yīng)用電壓源能形成節(jié)點(diǎn)支路基本方程SMA-HPC?2003問題單元電壓基本方 乘以 守恒 問題單元電壓基本方 乘以 守恒 形式：守恒量＝F（節(jié)點(diǎn)電壓SMA-HPC?2003問題單元?jiǎng)傂詣傂許MA-HPC?2003問題單元?jiǎng)傂詣傂許MA-HPC?2003與稀疏矩陣比較實(shí)例問題SMA-HPC?2003與稀疏矩陣比較實(shí)例問題SMA-HPC?2003與稀疏矩陣比較實(shí)例問題100x10SMA-HPC?2003與稀疏矩陣比較實(shí)例問題100x10SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003節(jié)點(diǎn)支路關(guān)節(jié)點(diǎn)支路關(guān)節(jié)點(diǎn)公SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003數(shù)值模擬數(shù)值模擬導(dǎo)論-感謝DeepakRamaswamyMichalRewienski,KarenVeroyandJacobWhite解的存在和唯解的存在和唯一高斯消元LU適應(yīng)條·無電源或剛性支承·無電源或剛性支承·對角和嚴(yán)格對角占優(yōu)矩陣·n×n方陣SMA-HPC?2003要求加權(quán)變量x要求加權(quán)變量x，使得矩陣M各列的加權(quán)和等于右邊的bSMA-HPC?2003GGGx1M1x2M2...xNGGGx1M1x2M2...xNMN疑問解答·－這個(gè)方程是否有解？－解是否唯一？·看是否有解存在一組變量x1,…..xn，使得：GGGx1M1x2M2...xNMN由此看出：只有當(dāng)b在由M各列組成的向量空間內(nèi)，解才存在。SMA-HPC?2003GGGy1M1y2M2...yNMNGGGy1M1y2M2...yNMN疑問解答續(xù)·解是否唯一？y1,…,y，的變量Gy1M1y2M2...yNMNM(x+y)=b又如果，則由此可見：當(dāng)且僅當(dāng)矩陣M的各列向量線性無關(guān)，方程解唯一。SMA-HPC?2003疑問解答續(xù)方·給定，其中M為N*N方陣疑問解答續(xù)方·給定，其中M為N*N方陣各列向量相互線性不相關(guān)的方陣稱之為非奇SMA-HPC?2003重要工具M(jìn)x用高斯消元法解線性方程·重要工具M(jìn)x用高斯消元法解線性方程·一種“直接”的方有限步求準(zhǔn)確解（不計(jì)舍入誤差）·求解增廣矩陣的精確解·計(jì)算所消耗的機(jī)時(shí)。SMA-HPC?2003舉例說明3*3SMA-HPC?2003舉例說明3*3SMA-HPC?2003舉例說明解題思路用矩陣的第一行消去第二和第三行的舉例說明解題思路用矩陣的第一行消去第二和第三行的SMA-HPC?2003舉例說明矩陣形式SMA-HPC?2003乘舉例說明矩陣形式SMA-HPC?2003乘舉例說明SMA-HPC?2003乘舉例說明SMA-HPC?2003乘舉例說明符號簡化表示SMA-HPC?2003舉例說明符號簡化表示SMA-HPC?2003舉例說明消去第三行的SMA-HPC?2003乘舉例說明消去第三行的SMA-HPC?2003乘舉例說明形成的三角陣SMA-HPC?2003舉例說明形成的三角陣SMA-HPC?2003舉例說明右邊向量的變化SMA-HPC?2003舉例說明右邊向量的變化SMA-HPC?2003舉例說明組合各部分SMA-HPC?2003舉例說明組合各部分SMA-HPC?2003Mx解方第一第二Mx解方第一第二消元過程解方程Ly回代過程第三Ux解方SMA-HPC?2003Mx解方第一第二Mx解方第一第二消元過程解方程Ly回代過程第三Ux解方SMA-HPC?2003解三角矩陣矩陣特點(diǎn)第一個(gè)方程只有y1一個(gè)未知量解三角矩陣矩陣特點(diǎn)第一個(gè)方程只有y1一個(gè)未知量第個(gè)方程只有y1和y2兩個(gè)未知SMA-HPC?2003解三角矩陣算SMA-HPC?2003解三角矩陣算SMA-HPC?2003解三角矩陣方程是直接的但又是費(fèi)時(shí)的。y1可以用一個(gè)除法運(yùn)算直接求得，y2可以用一乘法，一個(gè)減法，一個(gè)除法求得。一旦求得yk－1解三角矩陣方程是直接的但又是費(fèi)時(shí)的。y1可以用一個(gè)除法運(yùn)算直接求得，y2可以用一乘法，一個(gè)減法，一個(gè)除法求得。一旦求得yk－1，yk可以由k－1步乘法，k－2步加法，一N除法＋0加/減＋1加/減＋、、、、N-1加/減＋0乘法＋1乘法＋………N-1乘法求求＝（N-1）（N-2）加/減＋（N-1）（N-2）乘＋N＝N2分圖SMA-HPC?2003分圖SMA-HPC?2003用同樣的方法處理第二行。計(jì)算消去第三行和第四行中x2的乘子，并且用這些乘分算{i=1j=i+1每一{每一要消去的目標(biāo)MMMii為對角{對角元后面分算{i=1j=i+1每一{每一要消去的目標(biāo)MMMii為對角{對角元后面的元nk=i+1MMM乘}}}SMA-HPC?2003分對角元為零第i第i第jMM做一下分對角元為零第i第i第jMM做一下簡單變形（部分對角元MiijMii第j行和第i行交換SMA-HPC?2003分零對角元兩重要定理1）只有當(dāng)M為非奇異矩陣時(shí)分零對角元兩重要定理1）只有當(dāng)M為非奇異矩陣時(shí)列主元法（交換行）能有效2）對嚴(yán)格對角占優(yōu)矩陣LU分解將不會產(chǎn)生零對角元SMA-HPC?2003不會產(chǎn)生零對角元證明：1）求出第一步消元后的矩陣2）考察（n－1）×（n－1）的次矩陣仍然是完全對角占優(yōu)矩陣第一第一步消元后的第二由此得數(shù)值問題小對角元特我們能夠解釋這種現(xiàn)象嗎SMA-HPC數(shù)值問題小對角元特我們能夠解釋這種現(xiàn)象嗎SMA-HPC?2003出出數(shù)值問題小對角元浮點(diǎn)運(yùn)算的一個(gè)性雙精度主要問題結(jié)論數(shù)值問題小對角元浮點(diǎn)運(yùn)算的一個(gè)性雙精度主要問題結(jié)論避免大小數(shù)之間相加SMA-HPC?2003有效數(shù)位符號數(shù)值問題小對角元回過頭來看這個(gè)例子SMA-HPC?2003數(shù)值問題小對角元回過頭來看這個(gè)例子SMA-HPC?2003數(shù)值問題小對角元MMj那么第iMSMA-HPC?2003數(shù)值問題小對角元MMj那么第iMSMA-HPC?2003得是得是數(shù)值問題小對角元如果在LU分解過程中，矩陣是對角占數(shù)值問題小對角元如果在LU分解過程中，矩陣是對角占優(yōu)或用了列主元法來減小舍入誤差那么存在以下的點(diǎn)2）LU因子各位上的最大值將小于等于原始矩陣SMA-HPC?2003得是得是實(shí)多項(xiàng)式插值數(shù)據(jù)用一個(gè)N次的多項(xiàng)式來近似表示數(shù)據(jù)實(shí)多項(xiàng)式插值數(shù)據(jù)用一個(gè)N次的多項(xiàng)式來近似表示數(shù)據(jù)ftt ...2n012nSMA-HPC?2003實(shí)矩陣形SMA-HPC?2003實(shí)矩陣形SMA-HPC?2003實(shí)SMA-HPC?2003實(shí)SMA-HPC?2003在在我們用一個(gè)高次的多項(xiàng)式函數(shù)來逼近曲線t時(shí)。當(dāng)數(shù)據(jù)需要一個(gè)100次的多項(xiàng)式擾動分析 誤差范數(shù)向量的范n2xi2nx1擾動分析 誤差范數(shù)向量的范n2xi2nx1maxxiSMA-HPC?2003正方單位矩陣的范xxAx由矩陣A引起誤差的范數(shù)的最大值等于x被A放大的最大倍數(shù)矩陣的范xxAx由矩陣A引起誤差的范數(shù)的最大值等于x被A放大的最大倍數(shù)nmaxA1jA2nmaxAijSMA-HPC?2003擾動分析誤差范數(shù)擾動分析擾動計(jì)MMxxMxMxM擾動分析擾動計(jì)MMxxMxMxMxMx實(shí)際MxMxxxM1Mxx1xMMMMMxMSMA-HPC?2003由線性代數(shù)由線性代數(shù)學(xué)我們知道解x的變化范圍與原始的變化范圍之間有一個(gè)與矩陣M有關(guān)的倍數(shù)，這個(gè)倍M奇異矩陣超出本課程的考慮范圍，如果想了解請咨詢Trefethen或擾動分析更清晰的幾何逼近法GGM1GMx1M1x2擾動分析更清晰的幾何逼近法GGM1GMx1M1x2MM當(dāng)用向量來逼近時(shí)很難確定M1和與之對應(yīng)的M2的值SMA-HPC?2003幾何分析多項(xiàng)是插值多項(xiàng)式值得級數(shù)接線幾何分析多項(xiàng)是插值多項(xiàng)式值得級數(shù)接線SMA-HPC?2003問行的縮放比問行的縮放比例是否減小等比級行的縮放比例是否減小了條件條件定理：如果使用浮點(diǎn)運(yùn)算，當(dāng)進(jìn)行行縮放時(shí)在某種意義上不會減少等比級數(shù)沒有舍入誤解的存在和唯解的存在和唯一高斯消元LU適應(yīng)條數(shù)值模擬數(shù)值模擬導(dǎo)論-Luca感謝DeepakRamaswamyMichalRewienski,KarenVeroyandJacobWhite??回顧LU分解法稀疏矩陣—??回顧LU分解法稀疏矩陣—珩架和節(jié)點(diǎn)，電阻網(wǎng)，3d熱三角矩陣分解―一般的稀疏矩陣分―填充和重排—圖表逼稀疏矩陣數(shù)據(jù)結(jié)構(gòu)—散??SMA-HPC?2003分圖SMA-HPC?2003分圖SMA-HPC?2003U的x11，1，1變?yōu)榱恪Ｔ儆帽壤蜃樱ㄓ址Q之為乘子）代替這些零11/1，1位2341和11和1也被他們的乘子所代替。在這一過程中第三行其余的位置的值也會隨之改變，因此也將他們變?yōu)樗{(lán)色。用同樣的方法處理第二行。計(jì)算消去第三行和第四行中x2的乘子，并且用這矩陣分解算n-i=1{每一n{每一要消去的目標(biāo)到MMM到對角n矩陣分解算n-i=1{每一n{每一要消去的目標(biāo)到MMM到對角n{對角元后的元jiMikMjkM乘}}}SMA-HPC?2003矩陣分解對角占優(yōu)矩陣的性對一個(gè)對角占優(yōu)的矩陣進(jìn)行LU分解時(shí)不會產(chǎn)生零對角元。矩陣分解對角占優(yōu)矩陣的性對一個(gè)對角占優(yōu)的矩陣進(jìn)行LU分解時(shí)不會產(chǎn)生零對角元。。置上的值增加不會超過2n1)SMA-HPC?2003不會產(chǎn)生零對角元證明：1）求出第一步消元后的矩陣2）考察（n－1）×（n－1）的次矩陣仍然是完全對角占優(yōu)矩陣第一第一步消元后的第二由此得應(yīng)空間珩架空間珩節(jié)點(diǎn)矩未知方程：節(jié)點(diǎn)位合應(yīng)空間珩架空間珩節(jié)點(diǎn)矩未知方程：節(jié)點(diǎn)位合力SMA-HPC?2003應(yīng)電阻未知方程SMA-HPC?2003節(jié)點(diǎn)電電流應(yīng)電阻未知方程SMA-HPC?2003節(jié)點(diǎn)電電流和電阻網(wǎng)電阻網(wǎng)是一種特殊情況，它的數(shù)學(xué)模型是偏微分方程。（我們將在以后學(xué)習(xí)到）我們現(xiàn)在考察一下這個(gè)節(jié)點(diǎn)矩陣并注意矩陣中的非零數(shù)。從一個(gè)如下4的非零數(shù)應(yīng)電阻100×10電阻網(wǎng)矩陣的非零SMA-HPC?2003應(yīng)電阻100×10電阻網(wǎng)矩陣的非零SMA-HPC?2003應(yīng)四方體中的溫度場SMA-HPC?2003電模應(yīng)四方體中的溫度場SMA-HPC?2003電模三角帶矩陣矩陣形式SMA-HPC?2003三角帶矩陣矩陣形式SMA-HPC?2003三角帶矩陣高斯算法 i=1n-{每一到n{每一要消去的目標(biāo)MM對角到{三角帶矩陣高斯算法 i=1n-{每一到n{每一要消去的目標(biāo)MM對角到{對角元后的元MjiMik乘}}M需要N步運(yùn)}SMA-HPC?2003矩陣的填充例電阻例節(jié)點(diǎn)矩此矩陣為稱且對角矩陣的填充例電阻例節(jié)點(diǎn)矩此矩陣為稱且對角優(yōu)矩SMA-HPC?2003回想第二章回想第二章的內(nèi)容，節(jié)點(diǎn)矩陣的維數(shù)可以通過計(jì)算電阻器來確定，請看下面：這個(gè)電阻占據(jù)節(jié)點(diǎn)矩陣的如下也可以理解為：Gii就等于解點(diǎn)i處的電阻率（一個(gè)以上電阻）矩陣的填充例矩陣的非零數(shù)據(jù)結(jié)LU分解后的矩陣形矩陣的填充例矩陣的非零數(shù)據(jù)結(jié)LU分解后的矩陣形X都不為SMA-HPC?2003在LU分解的在LU分解的過程中，基準(zhǔn)行以外的各行要減去基準(zhǔn)行乘以一個(gè)因子。請通過下面的簡單例子，思考LU分解矩陣填充例填充的傳第一步的填充影響第二步的矩陣填充例填充的傳第一步的填充影響第二步的填SMA-HPC?2003行乘以一個(gè)系數(shù)，填充了第三列和第四列。當(dāng)用第三行和第四行減去第二乘以一個(gè)系數(shù)后，那么第三行和第四就產(chǎn)生了二級填充矩陣的填充重排SMA-HPC?2003未填填矩陣的填充重排SMA-HPC?2003未填填矩陣的填充重排在什么地方可以產(chǎn)生填充SMA-HPC?2003矩陣的填充重排在什么地方可以產(chǎn)生填充SMA-HPC?2003已經(jīng)分矩陣的填充重排MarkowitzFori=1到找到最小的Markowitz乘積的對角矩陣的填充重排MarkowitzFori=1到找到最小的Markowitz乘積的對角元j。交換第j行和第i行，第j列和第i列交換分解新的第i行，并確定填充結(jié)理想算SMA-HPC?2003第三步就是分解重新排列的矩陣并且填充矩陣。如果矩陣本來非常那么第三步將會很省接下來又必須計(jì)算分解后的矩陣的Markowitz乘積，這有又需要（N-1）步運(yùn)算。繼續(xù)，依此類推，只計(jì)算Markowitz乘積就需要約因此，我們要改善這種情況，方法便是計(jì)算矩陣每經(jīng)過一次分Makowitz矩陣的填充重排為什么只討論對角元與在節(jié)點(diǎn)公式中節(jié)矩陣的填充重排為什么只討論對角元與在節(jié)點(diǎn)公式中節(jié)點(diǎn)重排列相協(xié)調(diào)???減少搜索耗保留矩陣原有性－對角占－對SMA-HPC?2003矩陣的填充填充之后矩陣的樣SMA-HPC?2003稠矩陣的填充填充之后矩陣的樣SMA-HPC?2003稠非常稀非常稀矩陣的填充未分解的任意矩陣SMA-HPC?2003矩陣的填充未分解的任意矩陣SMA-HPC?2003矩陣的填充未分解的任意矩陣SMA-HPC?2003矩陣的填充未分解的任意矩陣SMA-HPC?2003矩陣意對稱矩陣的結(jié)構(gòu)和據(jù)陣矩陣每一行一個(gè)矩陣意對稱矩陣的結(jié)構(gòu)和據(jù)陣矩陣每一行一個(gè)節(jié)每一個(gè)非對角元對應(yīng)一條邊線SMA-HPC?2003當(dāng)且僅），矩陣當(dāng)且僅），矩陣如果矩陣在結(jié)構(gòu)上是對稱的與一個(gè)簡單的圖表相聯(lián)在圖中矩陣的每一行有一個(gè)如果0，在節(jié)點(diǎn)i和j之間有一個(gè)邊矩陣圖有兩個(gè)重要的性1)節(jié)點(diǎn)次數(shù)的平方就是Markowitz乘變化矩陣MarkowitZ乘Markowitz乘積＝（節(jié)點(diǎn)次數(shù)的平方SMA-HPC?2003矩陣MarkowitZ乘Markowitz乘積＝（節(jié)點(diǎn)次數(shù)的平方SMA-HPC?2003第第i個(gè)節(jié)點(diǎn)次數(shù)的平方等于Markowitz乘積，它可以由圖中很容易看出。節(jié)點(diǎn)矩陣矩陣分解LU分解的第一·刪掉與對角行相對矩陣矩陣分解LU分解的第一·刪掉與對角行相對應(yīng)的節(jié)·連接矩陣圖的邊SMA-HPC?2003每一步每一步LU分解需要一步浮點(diǎn)運(yùn)算并產(chǎn)生一個(gè)簡化矩陣，如完成了第i步分解之后，未分解的矩陣就被縮小為1必須加入與填充對應(yīng)的矩陣圖邊那么在矩陣圖中節(jié)點(diǎn)i也必須除去。另外，所用與節(jié)點(diǎn)i相鄰的節(jié)（相鄰節(jié)點(diǎn)通過邊界相聯(lián)）將會通過增加必要的邊界來使他們彼相矩陣?yán)齅arkowitz＝（節(jié)點(diǎn)次數(shù)SMA-HPC?2003據(jù)陣矩陣?yán)齅arkowitz＝（節(jié)點(diǎn)次數(shù)SMA-HPC?2003據(jù)陣矩陣?yán)粨Q第2和第1SMA-HPC?2003矩陣?yán)粨Q第2和第1SMA-HPC?2003更加平行的對角矩陣的另一順矩陣稀疏矩陣電阻網(wǎng)的例子未知：節(jié)電電方程：電矩陣稀疏矩陣電阻網(wǎng)的例子未知：節(jié)電電方程：電流和等于SMA-HPC?2003））的例子中我們可以很容易看出矩陣的特點(diǎn)。一個(gè)4×4系統(tǒng)它的節(jié)點(diǎn)矩陣如下零數(shù)是由于rowcoupling與三角帶的耦合產(chǎn)生的。矩陣稀疏矩陣網(wǎng)格的例子分解一個(gè)M×M的矩陣需要花矩陣稀疏矩陣網(wǎng)格的例子分解一個(gè)M×M的矩陣需要花費(fèi)多少機(jī)時(shí)推想中間列最后是否被消去SMA-HPC?2003要分解一個(gè)與要分解一個(gè)與類似由電阻器序列產(chǎn)生的M×M個(gè)刪格相對應(yīng)的雖然要得到這個(gè)問題的精確解已經(jīng)超出了本課程的內(nèi)容但我們可以步運(yùn)算，所以這表明分解這個(gè)由M×M步運(yùn)算，既然分解這個(gè)稠密矩陣刪格組成的矩陣也用非數(shù)值的元素用非數(shù)值的元素代替原來的矩對角占優(yōu)或?qū)ΨQ正用矩陣圖確定矩陣的次要用到許多矩陣圖操作的技必須認(rèn)真的組織計(jì)SMA-HPC?2003稀疏矩陣的數(shù)據(jù)結(jié)SMA-HPC?2003行指針向每稀疏矩陣的數(shù)據(jù)結(jié)SMA-HPC?2003行指針向每一行的數(shù)據(jù)排為了有效的為了有效的儲存一個(gè)稀疏矩陣，需要一種只記錄矩陣非零元素的這個(gè)矩陣的數(shù)據(jù)結(jié)我們可以看到這里沒有儲存任何的零稀疏矩陣的數(shù)據(jù)結(jié)冗余數(shù)從目標(biāo)行j中消去原始行i稀疏矩陣的數(shù)據(jù)結(jié)冗余數(shù)從目標(biāo)行j中消去原始行i考察j行所有位，我們發(fā)現(xiàn)其中3個(gè)位置與i行相對應(yīng)SMA-HPC?2003第j第i為了有效的為了有效的儲存一個(gè)稀疏矩陣，需要一種只記錄矩陣非零元素對應(yīng)。請看下面的例子這個(gè)矩陣的數(shù)據(jù)結(jié)我們可以看到這里沒有儲存任何的零稀疏矩陣的數(shù)據(jù)結(jié)構(gòu)冗余數(shù)據(jù)每一個(gè)多余的數(shù)據(jù)稀疏矩陣的數(shù)據(jù)結(jié)構(gòu)冗余數(shù)據(jù)每一個(gè)多余的數(shù)據(jù)都會占用不必要的內(nèi)存空SMA-HPC?2003為了有效的為了有效的儲存一個(gè)稀疏矩陣，需要一種只記錄矩陣非零元素對應(yīng)。請看下面的例子這個(gè)矩陣的數(shù)據(jù)結(jié)我們可以看到這里沒有儲存任何的零稀疏矩陣的數(shù)據(jù)結(jié)構(gòu)避免冗余數(shù)據(jù)的擴(kuò)散考察j稀疏矩陣的數(shù)據(jù)結(jié)構(gòu)避免冗余數(shù)據(jù)的擴(kuò)散考察j行所有的元素，將他們擴(kuò)展成為長度為n的向用標(biāo)定數(shù)組儲存需要的元SMA-HPC?2003第jLU分解和對角占LU分解和對角占－無數(shù)據(jù)元的因稀疏矩－壓桿，電阻網(wǎng)，3D熱對角矩陣的分－O(N)稀疏矩陣分解概－Markowitz重排步以達(dá)到最小填基于表格的方－分解和填－對降低用高斯消去解稀疏矩陣的復(fù)雜性是有幫助SMA-HPC?2003數(shù)值模擬數(shù)值模擬導(dǎo)論-感謝DeepakRamaswamy,MichalRewienski,KarenVeroyandKarenVeroy奇異矩陣?yán)覮U分解的不足之處雖然上面圖形的節(jié)點(diǎn)矩陣奇異矩陣?yán)覮U分解的不足之處雖然上面圖形的節(jié)點(diǎn)矩陣是一個(gè)奇異矩陣，但是仍舊存在一種解SMA-HPC?2003負(fù)載節(jié)拉奇異矩陣?yán)覮U分解的不足之處雖然上面圖形的節(jié)奇異矩陣?yán)覮U分解的不足之處雖然上面圖形的節(jié)點(diǎn)矩陣是一個(gè)奇異矩陣，但是存在一種解SMA-HPC?2003奇異矩陣?yán)踊仡櫨仃嚫髁械募訖?quán)和，并觀察下面奇異矩陣?yán)踊仡櫨仃嚫髁械募訖?quán)和，并觀察下面的等式GGx1M1x2M2...xNMN雖然矩陣M是奇異矩陣但是向量B在矩陣M的列向量組成的向量空間內(nèi)。SMA-HPC?2003正交如果M有正交列向量GiMiMj用第i列乘以加權(quán)列向量得正交如果M有正交列向量GiMiMj用第i列乘以加權(quán)列向量得GGGGMi(x1M1x2M2...xNMN利用正交向量將方程簡化為)MiGGG b Miiii(MMiiSMA-HPC?2003正交矩陣M正交的幾何意義G如果GGjMjM正交矩陣M正交的幾何意義G如果GGjMjMj則矩陣M正Mij二維向量的幾何意非正正SMA-HPC?2003正交QR法的基本思想原始矩帶有正交列向量正交QR法的基本思想原始矩帶有正交列向量的QybyQT怎么來完成這一步變換SMA-HPC?2003正交推導(dǎo)公式GGGM2r12M1M2，求滿GGG正交推導(dǎo)公式GGGM2r12M1M2，求滿GGGGGM1MM1GG1G即必有M11SMA-HPC?2003正交標(biāo)準(zhǔn)如果我們將向量標(biāo)準(zhǔn)化，公式將會變得簡單，因此我們先來將向量Q1標(biāo)準(zhǔn)化：GGGGG1正交標(biāo)準(zhǔn)如果我們將向量標(biāo)準(zhǔn)化，公式將會變得簡單，因此我們先來將向量Q1標(biāo)準(zhǔn)化：GGGGG1M1Q1Q1 M1rGGGGGQ2M22Q2現(xiàn)在要以便滿Gr12Q1MG G1G最后求GGrQ2SMA-HPC?2003正交2*2矩陣的變化過程既然Mx等于Qy，那么正交2*2矩陣的變化過程既然Mx等于Qy，那么我們可以找到x與y之間得關(guān)系SMA-HPC?2003正交2*2矩陣的分解過程標(biāo)準(zhǔn)正交如果給出QR分解后，需要正交2*2矩陣的分解過程標(biāo)準(zhǔn)正交如果給出QR分解后，需要兩步進(jìn)行求QRxbRxQTb第一Rx回代第二SMA-HPC?2003上三正交普通矩陣的QR分解3×3矩陣的情為了正交普通矩陣的QR分解3×3矩陣的情為了確保第三列正交使SMA-HPC?2003正交分解3*3矩陣必須解程求系SMA-HPC正交分解3*3矩陣必須解程求系SMA-HPC?2003正交解方程求系數(shù)將第N個(gè)向量正交N2項(xiàng)內(nèi)正交解方程求系數(shù)將第N個(gè)向量正交N2項(xiàng)內(nèi)積需要N3步運(yùn)SMA-HPC?2003正交使用正交向量3×3矩陣的情為了保證正交使用正交向量3×3矩陣的情為了保證第三列正交SMA-HPC?2003基本運(yùn)算法則改進(jìn)的Gram-SchmidnGGrM標(biāo)準(zhǔn)GiiN需要2N2N2步操基本運(yùn)算法則改進(jìn)的Gram-SchmidnGGrM標(biāo)準(zhǔn)GiiN需要2N2N2步操G j＝i＋1GnGNi2NNjGG MjSMA-HPC?2003基本運(yùn)算法則用圖形表示SMA-HPC?2003基本運(yùn)算法則用圖形表示SMA-HPC?2003基本運(yùn)算法則用圖形表示SMA-HPC?2003基本運(yùn)算法則用圖形表示SMA-HPC?2003基本運(yùn)算法則零如果有一列元素全為零該怎么基本運(yùn)算法則零如果有一列元素全為零該怎么辦此矩陣肯定是奇異矩陣SMA-HPC?2003基本運(yùn)算法則零列（序QR分解結(jié)果SMA-HPC?2003基本運(yùn)算法則零列（序QR分解結(jié)果SMA-HPC?2003奇異矩陣舉例回顧矩陣各列的加權(quán)和，并觀察下面的奇異矩陣舉例回顧矩陣各列的加權(quán)和，并觀察下面的等式當(dāng)M為奇異矩陣時(shí)會出現(xiàn)以下兩種情況：GG{M1MNb屬于向量空間{Q1,...,QNG{M1,...,MNSMA-HPC?2003最小值法求解公式求x使下式獲最小值NR最小值法求解公式求x使下式獲最小值NRx2Rx i等價(jià)于，如果b屬于向量空間{cols(M)}推出Mx＝bminxRRxT最小值法擴(kuò)展到非奇異或非平方的情SMA-HPC?2003最小值法一維空間最小值法GMxxMGxxx最小值法一維空間最小值法GMxxMGxxx 1 一維空間最小值正交標(biāo)SMA-HPC?2003最小值法一位空間最小值法（圖示一維最小值法產(chǎn)生的結(jié)果和b最小值法一位空間最小值法（圖示一維最小值法產(chǎn)生的結(jié)果和bSMA-HPC?2003最小值法二維空間最小值法xxGxMxxM最小值法二維空間最小值法xxGxMxxMGxMG 12用殘向量最小值耦合SMA-HPC?2003最小值法二維向量最小值法從更一般的角度考GGMp2T最小值法二維向量最小值法從更一般的角度考GGMp2TMSMA-HPC?2003最小值法構(gòu)建MTM正交最小化方向第i次搜索方向等于MTM正交化單位向GGGGGjpi最小值法構(gòu)建MTM正交最小化方向第i次搜索方向等于MTM正交化單位向GGGGGjpirpT MMpjjMGGpjjSMA-HPC?2003最小值法搜索方向上的最小值單獨(dú)做去藕最小值求的最小GGGTT2vbii將其對v求GGGT 最小值法搜索方向上的最小值單獨(dú)做去藕最小值求的最小GGGTT2vbii將其對v求GGGT 2bTMiiG求得TbMpivGiTMMSMA-HPC?2003最小最小化算法對每一目標(biāo)i=1GGpieforj到iGGrpjTTMpiGr最小最小化算法對每一目標(biāo)i=1GGpieforj到iGGrpjTTMpiGrGpi MMiGGrxxviiSMA-HPC?2003最小化法和QR分解SMA-HPC?2003正交最小化法和QR分解SMA-HPC?2003正交搜索方向正交單位向量——,G,...,G{G,搜索方向正交單位向量——,G,...,G{G,G p2NMTM12搜索單位可以使用其它的初始向,...,GG,b,Mb,Mkrylov子空MTM正交12搜索方為什么SMA-HPC?2003QR算——————修改QR算——————修改Gram-SchmidtQR和奇異矩零QR分解的最小化簡單提一下關(guān)于調(diào)整搜索方SMA-HPC?2003數(shù)值模擬數(shù)值模擬導(dǎo)論-感謝DeepakRamaswamy,MichalRewienski,KarenVeroyandJacobWhite常規(guī)的子空間常規(guī)的子空間極小化算——回顧學(xué)過的正交化和投射定GCR算－krylov-子空－對稱矩陣的簡－收斂條回顧特征值和特征向－范數(shù)和譜半－譜映射定G,w0k選擇一個(gè)kGG,w0k選擇一個(gè)kGG可以近似為向量GkiikSMA-HPC?2003GGG 最小殘向量b·殘向量定義kkGGGGG 最小殘向量b·殘向量定義kkGGb如果i ii殘向量最小化的思路為：取s，最小化式iGGkk2TbMwbkriiii2i0i0SMA-HPC?2003最小殘向量 算MGGG，那 或i22們將很容易計(jì)算的最小值bii22GG最小殘向量 算MGGG，那 或i22們將很容易計(jì)算的最小值bii22GGMGMT或非正MpjGG那么要建立一組正交的向量p,,使向pk0G并GG0，i,,wk＝向量空TMp0SMA-HPC?2003最小殘向量 運(yùn)算步驟G給定M，b并且一組搜索方,,wk01）通過將Mws正交化生GsjjMpkforj＝0to最小殘向量 運(yùn)算步驟G給定M，b并且一組搜索方,,wk01）通過將Mws正交化生GsjjMpkforj＝0tojipjpMpMpTiii2）解xk計(jì)算r的最小r0riMpMpkk iikxppi0iTTMp MpMpiiiiSMA-HPC?2003最小殘向量 圖示計(jì)算步驟1）正最小殘向量 圖示計(jì)算步驟1）正交2）解計(jì)算r的最小SMA-HPC?2003最小化算法bforj＝0topjfori＝0toj- p最小化算法bforj＝0topjfori＝0toj- pMpMp正交搜索方j(luò)jjii p標(biāo)準(zhǔn)化向pjjjTjxjrjjjMpxj更新結(jié)jjMprj更新殘向jjSMA-HPC?2003子空間的選擇 標(biāo)GG選的標(biāo),,GG對所有在空間{w0wk子空間的選擇 標(biāo)GG選的標(biāo),,GG對所有在空間{w0wk1}中Mk的值都很bMxkbiGG在向量空間{w0,wkA1b當(dāng)k??Nx Ge,G向量空1如果k＝N進(jìn)行QR分如果k<N情況會很糟SMA-HPC?2003子空間的選擇傳統(tǒng)方法fx1xTMxxTb的最小值，其中假求子空間的選擇傳統(tǒng)方法fx1xTMxxTb的最小值，其中假求MMT（對2矩陣）并xTMxxfxMxbx推導(dǎo)出x最小化bGw0,fxk1fx0,,取向量空xx這便是f的最快下降方向，但f并不是殘余值這種方法不能用于非對成矩陣，和不滿足xTMx的情SMA-HPC?2003子空間的選krylov子空向量空間＝向量空間fx,,xfxk子空間的選krylov子空向量空間＝向量空間fx,,xfxk注意,, 0k0G＝向量空間,,rk1}w0,,k那r0iik1r0{r0,,rk1}＝向量并且向量空,Mr0,,krylov子空SMA-HPC?2003 GCR算 krylovGCP的第krkMp求解第k步搜索方向的步kkkkrk 更新結(jié)k GCR算 krylovGCP的第krkMp求解第k步搜索方向的步kkkkrk 更新結(jié)kkkk更新殘向rkrkkkMrk1Tk計(jì)算新的正交搜索方rk1pkjp jTjjjSMA-HPC?2003 GCR算 krylovkrk 向量內(nèi)積，O(n)k kkkxk GCR算 krylovkrk 向量內(nèi)積，O(n)k kkkxkxkrkpO(n)如果是稀疏矩陣需要O(n) 向量加rkkkrkrkO(K)內(nèi)積，總共需要O(nk)kk如果M是稀疏矩陣，當(dāng)用k步來逼近n時(shí)，總共SMA-HPC?2003 GCR算 krylov我們會發(fā)現(xiàn)下面的情況那MpjrkMrk1Mrk1kj GCR算 krylov我們會發(fā)現(xiàn)下面的情況那MpjrkMrk1Mrk1kjjk k一步完成正pppkjkkTMpMpkkn，那么矩陣對稱，稀疏，GCR需要運(yùn)如O(n2kSMA-HPC?2003“krylov吸近端溫遠(yuǎn)端溫將棒離散“krylov吸近端溫遠(yuǎn)端溫將棒離散節(jié)點(diǎn)平衡方SMA-HPC?2003krylovSMA-HPC?2003krylovSMA-HPC?2003krylov近端溫遠(yuǎn)端溫離散節(jié)點(diǎn)平krylov近端溫遠(yuǎn)端溫離散節(jié)點(diǎn)平衡方SMA-HPC?2003krylov節(jié)點(diǎn)平衡方SMA-HPC?2003krylov節(jié)點(diǎn)平衡方SMA-HPC?2003殘余誤迭代次反殘余誤迭代次反復(fù)迭代后的log（殘余誤差）對SMA-HPC?2003(Rhs＝－1，＋(Rhs＝－1，＋反復(fù)迭代后的log（殘余誤差）對SMA-HPC?2003收斂性分析Krylovw,....wr0,Mr0,....Mk0kk xk收斂性分析Krylovw,....wr0,Mr0,....Mk0kk xk M0kik次多項(xiàng)kMr0Irkr0 i kw,....wr0,Mr0 kr00kSMA-HPC?2003收斂性分析Krylovj0，在GCRjr0,Mr0 kr0p,p k收斂性分析Krylovj0，在GCRjr0,Mr0 kr0p,p k2(M是k次多項(xiàng)式，這個(gè)多項(xiàng)式可以最rkkxkk2Ir0MMbMxk1r03rkkk22(M)r0是（k＋1）次多項(xiàng)式，要想最rk這里kk1(0)SMA-HPC?2003收斂性分析KrylovGCR多項(xiàng)式的22GCR最優(yōu)化的rk(Mk收斂性分析KrylovGCR多項(xiàng)式的22GCR最優(yōu)化的rk(Mk22kk1(0)2rk2SMA-HPC?2003基本定義Mi為iiM基本定義Mi為iiMiIi是矩陣M如果(MiI)GSMA-HPC?2003基本定義SMA-HPC?2003基本定義SMA-HPC?2003注意問題幾乎所有的N*N矩陣都有注意問題幾乎所有的N*N矩陣都有N個(gè)線性獨(dú)立的特征向量。矩陣M的一組特征值被稱之為MSMA-HPC?2003注意問題幾乎所有的N*N矩陣都有N個(gè)線性獨(dú)立注意問題幾乎所有的N*N矩陣都有N個(gè)線性獨(dú)立的特征向量?？梢缘扔诓⒉灰馕吨鳰SMA-HPC?2003譜半矩陣M所有矩陣MSMA-HPC?2003譜半矩陣M所有矩陣MSMA-HPC?2003熱流例子SMA-HPC?2003熱流例子SMA-HPC?2003熱流例子（續(xù)SMA-HPC?2003熱流例子（續(xù)SMA-HPC?2003熱流例子（續(xù)SMA-HPC?2003熱流例子（續(xù)SMA-HPC?2003譜半徑定理fx01x譜半徑定理fx01x...ppfM01M...pp譜半徑fMf譜半徑MSMA-HPC?2003譜半徑定理證明MMUU1UU1U2U UpU譜半徑定理證明MMUU1UU1U2U UpUp UU...UpU101ppU I...01pI... U p01pSMA-HPC?2003譜分xGGGN譜分xGGGN用矩陣M來代替xMxM(GGG... N)GGu...u SMA-HPC?2003收斂性分析Krylov1）GCR運(yùn)算法則在最多n步內(nèi)收斂于真實(shí)解。x收斂性分析Krylov1）GCR運(yùn)算法則在最多n步內(nèi)收斂于真實(shí)解。xx1x2....xn其中iM證明：令0，因此r 0n2）如果M只有q個(gè)特征值，那么GCR運(yùn)算法則最多以需SMA-HPC?2003····GCR－QSMA-HPC?2003數(shù)值模擬數(shù)值模擬導(dǎo)論-感謝DeepakRamaswamy,MichalRewienski,KarenVeroyandJacobWhite·回顧－Krylov·回顧－Krylov·回顧特征值和·收斂速度的－Chebychev多項(xiàng)－近似LU預(yù)處G,w0kbforj＝0topjrfori＝0toG,w0kbforj＝0topjrfori＝0toj- pMpMpjjjii ppjjMpMpxjrjjjMpxjjjMprjjjSMA-HPC?2003標(biāo)準(zhǔn)圖示運(yùn)算步1）正交2）xk標(biāo)準(zhǔn)圖示運(yùn)算步1）正交2）xk計(jì)算rSMA-HPC?2003GGG 開始的幾步rGGG 開始的幾步rb· ，0r0·p00 1r·p1.0 Mr11p SMA-HPC?2003GGG GGG x1r11· ·bMx2r0Mr0rM2r0 p Mr12pp SMA-HPC?2003GGG GCR的第kkGGG GCR的第kkTpkrkkMr jjj kkkrMpTk第k步運(yùn)算的最佳kxkxk rkrkkkSMA-HPC?2003GGG 多項(xiàng)式梗概的有j0，那么jGGG 多項(xiàng)式梗概的有j0，那么j1）向量空間p,p,...,p向量空間r0Mr0Mrk kMr02rk是kkk2bMxk1r0M IMr0 Mr003）rkkkk01這里是（k＋1）次多項(xiàng)式kk2rk2SMA-HPC?2003殘向量最小值多項(xiàng)式梗概,Mr0,...,Mrkxk2如最小：rk2Mr021xkrk是殘向量最小值多項(xiàng)式梗概,Mr0,...,Mrkxk2如最?。簉k2Mr021xkrk是k次多項(xiàng)式，最小kk2r IMMb002）rkkkk1）次多項(xiàng)式，如0 是kk那么他可以最小2rk2這里多項(xiàng)式作為解題工具，只有一個(gè)作用。那就最小化殘向量SMA-HPC?2003“與外界物熱交Krylov的例子吸近端溫遠(yuǎn)端“與外界物熱交Krylov的例子吸近端溫遠(yuǎn)端溫離散節(jié)點(diǎn)平衡方SMA-HPC?2003Krylov近端溫離散遠(yuǎn)端溫節(jié)點(diǎn)Krylov近端溫離散遠(yuǎn)端溫節(jié)點(diǎn)平衡方SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003(Rhs＝－1，＋(Rhs＝－1，＋SMA-HPC?2003殘向量最小值多項(xiàng)式最優(yōu)化MMrkk殘向量最小值多項(xiàng)式最優(yōu)化MMrkkkk1是任意的k次順序多項(xiàng)式，所以有krkSMA-HPC?2003矩陣的放大倍數(shù)問題假設(shè)y，那么y矩陣的放大倍數(shù)問題假設(shè)y，那么y比xSMA-HPC?2003回顧向量的范L2L1L回顧向量的范L2L1LSMA-HPC?2003矩陣L階范1MlxxlxllnM1j矩陣L階范1MlxxlxllnM1jnMjSMA-HPC?2003矩陣L階范nx1,1,...矩陣L階范nx1,1,...maxM1nmaxjMM2通過前面的介紹我們知道，解X的值的變動同一個(gè)因子有關(guān)，這個(gè)因子是關(guān)于A的。這個(gè)因子是，A (Cond(A)min(SMA-HPC?2003譜半徑理論fx1x...譜半徑理論fx1x...ppfMM...M01pspectrumfMfspectrumMSMA-HPC?2003收斂性分析矩陣多項(xiàng)式的標(biāo)M特征空間的條圖中英文為：矩陣收斂性分析矩陣多項(xiàng)式的標(biāo)M特征空間的條圖中英文為：矩陣MSMA-HPC?2003收斂性分析矩陣多項(xiàng)式的標(biāo)SMA-HPC?2003收斂性分析矩陣多項(xiàng)式的標(biāo)SMA-HPC?2003收斂性分析重要發(fā)1）一個(gè)Krylov子空間殘向量最小化的運(yùn)算，最多N步收斂于x收斂性分析重要發(fā)1）一個(gè)Krylov子空間殘向量最小化的運(yùn)算，最多N步收斂于xx1x2....xn，這里iM令M0ininxx1x2..xqqSMA-HPC?2003對稱矩陣的收斂性多項(xiàng)式的殘余如果M1）M有對稱矩陣的收斂性多項(xiàng)式的殘余如果M1）M有標(biāo)準(zhǔn)正交的特征如果M正定MSMA-HPC?2003導(dǎo)熱棒矩陣的多項(xiàng)式殘余量圖導(dǎo)熱棒矩陣的多項(xiàng)式殘余量圖SMA-HPC?2003kkSMA-HPC?2003導(dǎo)熱棒矩陣的多項(xiàng)式殘余量圖對稱矩陣的收斂Mmin,max,對稱矩陣的收斂Mmin,max,0pk（小的多項(xiàng)式性能比較好）多項(xiàng)式的最大最小值問題已經(jīng)被ChebyshevSMA-HPC?2003對稱矩陣的收斂Chebyshev方法：Chebyshev多項(xiàng)對稱矩陣的收斂Chebyshev方法：Chebyshev多項(xiàng)SMA-HPC?2003Chebyshev多項(xiàng)式Chebyshev多項(xiàng)式的最小化超出了SMA-HPC?2003對稱矩陣的收斂KrylovChebyshev的范SMA-HPC?2003對稱矩陣的收斂KrylovChebyshev的范SMA-HPC?2003對稱矩陣的收斂KrylovChebyshev的結(jié)SMA-HPC?2003對稱矩陣的收斂KrylovChebyshev的結(jié)SMA-HPC?2003前處KrylovSMA-HPC?2003前處KrylovSMA-HPC?2003前處Krylov讓M=D+M其中DD1Mx前處Krylov讓M=D+M其中DD1MxID1MndxSMA-HPC?2003系統(tǒng)離散SMA-HPC?2003系統(tǒng)離散SMA-HPC?2003SMA-HPCSMA-HPC?2003前處理矩陣特征SMA-HPC?2003前處理矩陣特征SMA-HPC?2003熱流動的例子維密稀疏熱流動的例子維密稀疏SMA-HPC?2003前處LU的前處L為下三角矩陣，ULU前處LU的前處L為下三角矩陣，ULUMxLU1b中Mx形式就等LUy。SMA-HPC?2003前處LU前處理（續(xù)在真實(shí)的LU填充這個(gè)LUSMA-HPC?2003前處LU前處理（續(xù)在真實(shí)的LU填充這個(gè)LUSMA-HPC?2003SMA-HPC?2003分解SMA-HPC?2003分解2前處KrylovLU近似分解前處理（續(xù)）放棄LUSMA-HPC?2003前處KrylovLU近似分解前處理（續(xù)）放棄LUSMA-HPC?2003－Krylov－Krylov－Chebychev－近似LUSMA-HPC?2003數(shù)值模擬分析導(dǎo)論數(shù)值模擬分析導(dǎo)論——第八講Jacob感謝ThankstoDeepakRamaswamyJaimePeraireRewienski,andKarenSMA-HPC?2003非線性問——非線性問——理查森和線性收——牛頓——————————SMA-HPC?2003壓桿實(shí)例SMA-HPC?2003壓桿實(shí)例SMA-HPC?2003壓桿實(shí)例提示：求壓桿力SMA-HPC?2003壓桿實(shí)例提示：求壓桿力SMA-HPC?2003壓桿實(shí)例SMA-HPC?2003壓桿實(shí)例SMA-HPC?2003壓桿實(shí)例為什么是非線性的增大對壓桿的拉桿受力壓桿實(shí)例為什么是非線性的增大對壓桿的拉桿受力大小和方向都發(fā)生改SMA-HPC?2003電路實(shí)例SMA-HPC?2003電路實(shí)例SMA-HPC?2003迭代求解fx0迭代求解fx0Wxkxkfxk1迭代收斂于正確的解迭代收斂的速SMA-HPC?2003定理查森迭代的定fxkx定理查森迭代的定fxkxkxk迭代收斂于方程的一個(gè)xkfxkxkx解SMA-HPC?2003實(shí)例SMA-HPC?2003實(shí)例SMA-HPC?2003實(shí)例SMA-HPC?2003實(shí)例SMA-HPC?2003實(shí)例SMA-HPC?2003實(shí)例SMA-HPC?2003收設(shè)xkfxxkx收設(shè)xkfxxkxfxxkxfxkfx需要xk1SMA-HPC?2003收均值定理SMA-HPC?2003收均值定理SMA-HPC?2003收均值定理的應(yīng)用xxkxfxfxxk收均值定理的應(yīng)用xxkxfxfxxk1xxkxfxkfxx f1kxSMA-HPC?2003收理查森法則SMA-HPC?2003線性收或那并如收理查森法則SMA-HPC?2003線性收或那并如實(shí)例SMA-HPC?2003實(shí)例SMA-HPC?2003問只是線性x,ff-——在電問只是線性x,ff-——在電路中表示位移f表示施加-——在桿件問但是算法fx——SMA-HPC?2003牛頓SMA-HPC?2003通過泰勒公式牛頓SMA-HPC?2003通過泰勒公式展開求圖形表示SMA-HPC?2003圖形表示SMA-HPC?2003例SMA-HPC?2003二漸漸例例SMA-HPC?2003二漸漸例例SMA-HPC?2003例SMA-HPC?2003收斂SMA-HPC?2003收斂SMA-HPC?2003收斂續(xù)SMA-HPC?2003兩邊同除以dfdx收斂續(xù)SMA-HPC?2003兩邊同除以dfdx收斂實(shí)例SMA-HPC?2003收斂實(shí)例SMA-HPC?2003收斂實(shí)例SMA-HPC?2003注收斂實(shí)例SMA-HPC?2003注收斂例子SMA-HPC?2003收斂例子SMA-HPC?2003收斂SMA-HPC?2003證收斂那收斂SMA-HPC?2003證收斂那收斂定df若L有界收斂定df若L有界總是收斂么SMA-HPC?2003收斂實(shí)收斂性依賴于給定恰當(dāng)?shù)某跏糞MA-HPC收斂實(shí)收斂性依賴于給定恰當(dāng)?shù)某跏糞MA-HPC?2003收斂收斂性檢測需要一個(gè)“”來檢驗(yàn)是否產(chǎn)生收斂收斂性檢測需要一個(gè)“”來檢驗(yàn)是否產(chǎn)生錯(cuò)誤的收SMA-HPC?2003收斂收斂性檢測fx”來檢驗(yàn)是否產(chǎn)生收斂收斂性檢測fx”來檢驗(yàn)是否產(chǎn)生錯(cuò)誤的收同樣需要一個(gè)SMA-HPC?2003非線性問——非線性問——理查森和線性收——牛頓——————————SMA-HPC?2003數(shù)值模擬數(shù)值模擬導(dǎo)論-雅克比·感謝DeepakRamaswamy,JaimeMichalRewienski,andKaren簡單回顧一維牛頓——收斂性簡單回顧一維牛頓——收斂性檢多維牛頓——基本算——雅可比矩陣的描——方程多維收斂——證明局部收斂——收斂性的改fx*x*fx2xffx*x*fx2xfx*xxx*fx2fxxfxSMA-HPC?2003初始給定值，kfxk初始給定值，kfxkxkfxkkkfxk1xk1SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003SMA-HPC?2003收斂性檢驗(yàn)SMA-HPC?2003收斂性檢驗(yàn)x”來檢測是否產(chǎn)生錯(cuò)誤的收斂同樣需要一個(gè)x”來檢測是否產(chǎn)生錯(cuò)誤的收斂同樣需要一個(gè)fSMA-HPC?2003收斂性檢驗(yàn)SMA-HPC?2003局部收斂SMA-HPC?2003局部收斂Fxff0GFxx0FyL或xllF00Ll 22xyyFyll0 lllFxff0GFxx0FyL或xllF00Ll 22xyyFyll0 llll E0c0l0xxlllfx0lylylfy0llSMA-HPC?2003桿件和節(jié)點(diǎn)問題節(jié)點(diǎn)分析1i2節(jié)點(diǎn)分析1i2gv3gv1v2SMA-HPC?2003Fx*問題：求解出使得xFx*問題：求解出使得xRN且FRN利用泰勒級數(shù)展開式RNFx*xxxH.OFxJ雅可比矩陣若近似于精確解則J xSMA-HPC?2003xFx*節(jié)點(diǎn)分析xR2且FR2xllF00Llxy節(jié)點(diǎn)分析xR2且FR2xllF00LlxyF00LlyG?Jx?FSMA-HPC?2003桿件和節(jié)點(diǎn)問題xR2且FR2i12GxR2且FR2i12G在節(jié)點(diǎn)2處：i3i2FGgvgvvv1312GJxFSMA-HPC?2003非線性電阻器問JFxxFxxFFxFJFxxFxxFFxF"1N1JFxx#%"##FxNN11SMA-HPC?2003 x為奇數(shù)FxFx"1N1JF x為奇數(shù)FxFx"1N1JFxx#x%"##FFxNN11SMA-HPC?2003奇數(shù)實(shí)例x0＝初始計(jì)算Fxk,xkF解方程JFxx0＝初始計(jì)算Fxk,xkF解方程JFxkk1xxFx求得kkkkkFxk1}xk1,SMA-HPC?2003計(jì)算雅可比矩陣和函數(shù)考慮節(jié)點(diǎn)之間的一個(gè)非線性電阻的作用gvb 2Fnvg計(jì)算雅可比矩陣和函數(shù)考慮節(jié)點(diǎn)之間的一個(gè)非線性電阻的作用gvb 2Fnvg n2上： n212g vvgn對節(jié)點(diǎn)n求偏微分121121SMA-HPC?2003計(jì)算雅可比矩陣和函數(shù)一個(gè)電阻器上受到的電#####gvgvnnv# 計(jì)算雅可比矩陣和函數(shù)一個(gè)電阻器上受到的電#####gvgvnnv# g####11##ngvgv # #2##JF##2#FSMA-HPC?2003進(jìn)一步計(jì)算牛頓算法x0k＝初始值重復(fù)Fx進(jìn)一步計(jì)算牛頓算法x0k＝初始值重復(fù)Fxk,F將JF和FxkxkxkFxk求得kkFFxk1xk1xk直,足夠小為}實(shí)例：直通棒中的熱流什么是雅可比行實(shí)例：直通棒中的熱流什么是雅可比行矩陣？SMA-HPC?2003多維收斂定理主要定1kJFJF多維收斂定理主要定1kJFJFxJFylxyz定理陳多維收斂定理重要輔助定JFxJFylxy多維收斂定理重要輔助定JFxJFylxyzl2FxFy yxF2這里沒有多維的均值定理SMA-HPC?2003通過牛頓迭代的定義和假定雅可比矩陣逆的極限值得FxkFxkxk通過牛頓迭代的定義和假定雅可比矩陣逆的極限值得FxkFxkxk1FkJx再次利用牛頓迭代的定義FxkFxk1kxk1xkxF0最后利用輔助定理22xkxkxkSMA-HPC?2003等式重xk1xkxkxkxk等式重xk1xkxkxkxk2x1若2xkx0xkkxk1xk則一維必須設(shè)法限定X的變一維必須設(shè)法限定X的變化范SMA-HPC?2003用牛頓算法求x0重復(fù)＝初始值，k計(jì)算Fx用牛頓算法求x0重復(fù)＝初始值，k計(jì)算FxkxkF解方程JFxFx求得kkkkxklimxk1xkkk1Fxk1xk直,足夠小}迭代方向正xkxkiiik表示迭代方向正xkxkiiik表示1xi迭代方向錯(cuò)li kkx mi,kx由此推斷，牛頓迭代不能保證全局SMA-HPC?2003一般阻尼定解方程JFxxFx 求解kkkk一般阻尼定解方程JFxxFx 求解kkkkkxkxk主要思想：線性搜Fxkkxk12找出k取極使2FxkFxkFxkkxk1kxk122該法在牛頓收斂方向執(zhí)行一維搜SMA-HPC?2003FkJxJFxyFkJxJFxyx(導(dǎo)數(shù)為Lipschitz存在一些列k0,1FFFxk1其中每進(jìn)行一步迭代則減小了F——全局收斂性SMA-HPC?2003x0k＝重復(fù)計(jì)算FxkkF解方程JFxx0k＝重復(fù)計(jì)算FxkkF解方程JFxxFx求得kkkkFxk找出k使kxk取極小xkxkkxkkkFxk1xk}直,足夠小為如何求阻尼系數(shù)嵌套迭阻尼牛頓奇異雅可比矩陣阻尼牛頓阻尼牛頓奇異雅可比矩陣阻尼牛頓法“推動”迭代趨向局部極小找出雅可比矩陣的奇異SMA-HPC?2003-——收斂性-——收斂性-第十?dāng)?shù)值-第十?dāng)?shù)值模擬導(dǎo)改進(jìn)的牛頓雅克比·懷感謝DeepakRamaswamy,JaimeMichalRewienski,andKaren阻尼牛頓定—若雅可阻尼牛頓定—若雅可比矩陣是非奇異矩陣，則全局收—奇異雅可比矩陣收斂非常困介紹連續(xù)定—源/載荷步問牛頓算法牛頓算法求初始給定值，k重復(fù)計(jì)算Fxk牛頓算法牛頓算法求初始給定值，k重復(fù)計(jì)算FxkxkFxkxkxkFx求得kFxklimxk1xkkk1fxk1}xk足夠直SMA-HPC?2003阻尼牛頓定律一般阻尼定 xxFx求解kk阻尼牛頓定律一般阻尼定 xxFx求解kkkkkxkxk主要思想：線性搜Fxkkxk12找出k使取極2Fxkkxk1Fxkkxk1TFxkkxk122該法在牛頓收斂方向執(zhí)行一維搜SMA-HPC?2003阻尼牛頓收斂定律如FkJx(反之不成立xy阻尼牛頓收斂定律如FkJx(反之不成立xyx(導(dǎo)數(shù)為Lipschitz那存在一些列 0,1使kFxk1Fxkkxk1Fxk其中每進(jìn)行一步迭代則減小了F——全局收斂SMA-HPC?2003阻尼牛頓1IrrVd Is1Vt阻尼牛頓1IrrVd Is1Vt節(jié)點(diǎn)方程v211016ev201f210SMA-HPC?2003實(shí)阻尼牛問解答續(xù)SMA-HPC?2003阻尼牛問解答續(xù)SMA-HPC?2003嵌套迭x0＝初始值，k重復(fù)計(jì)算Fxk,xkF解方程嵌套迭x0＝初始值，k重復(fù)計(jì)算Fxk,xkF解方程JFxFx求得kkkkFxk取極小找出 使kxkkkxkxkkkFxk1}xk,如何求阻尼系數(shù)SMA-HPC?2003阻尼牛頓定定理證通過定義牛頓迭F1kxkF牛頓方向多維均值迭l2Fx阻尼牛頓定定理證通過定義牛頓迭F1kxkF牛頓方向多維均值迭l2FxFyyxy2xF綜合以2