一種基于廣義傅里葉級數構建改進微分求積法的方法技術

本發明專利技術公開了一種基于廣義傅里葉級數構建改進微分求積法的方法，該方法包括：將需要求解的函數表示成傅里葉級數，并在傅里葉級數的基礎上加上附加項；根據邊界條件，確定附加項系數和傅里葉級數系數之間的關系；將求解域離散成獨立點，把帶有附加項的傅里葉級數代入各個獨立點，形成各傅里葉級數系數之間的關系矩陣；根據所述關系矩陣求解出傅里葉級數系數，確定需要求解的函數。該方法提高求解結果準確性。

【技術實現步驟摘要】
一種基于廣義傅里葉級數構建改進微分求積法的方法
本專利技術涉及數值求解
，特別是涉及一種基于廣義傅里葉級數構建改進微分求積法的方法。
技術介紹
目前，在工程和科學領域里，求解線性、非線性方程或方程組是一項很重要的任務。由于只有少數情況能得到解析解，所以大家把重點均放在求其近似解或數值解上。微分求積法在1971年提出，由于其較高的計算精度和效率，越來越受到人們的重視。但是微分求積法存在一些缺陷，即在處理結構振動方程的彈性邊界時存在處理不精確，導致求解誤差，求解結果準確性較低。在處理結構彈性邊界的方法中，2000年美國韋恩州立大學Li.W.L.教授在雜志《JournalofSoundandVibration》中第237卷中發表文章《Freevibrationsofbeamswithgeneralboundaryconditions》提出了廣義傅里葉級數，即在傳統傅里葉級數的基礎上加上附加項，用于滿足結構邊界處的不連續，此方法成功解決了彈性邊界條件下梁的橫向自由振動振動，2007年，杜敬濤在雜志《JournalofSoundandVibration》第306卷中發表文章《Ananalyticalmethodforthein-planevibrationanalysisofrectangularplateswithelasticallyrestrainededges》，此文章利用廣義傅里葉級數方法解決了彈性邊界矩形板的振動問題。目前還沒有廣義傅里葉級數和微分求積法相結合的技術出現以及相關文獻報道。
技術實現思路
本專利技術的目的是提供一種基于廣義傅里葉級數構建改進微分求積法的方法，以實現提高求解結果準確性。為解決上述技術問題，本專利技術提供一種基于廣義傅里葉級數構建改進微分求積法的方法，該方法包括：將需要求解的函數表示成傅里葉級數，并在傅里葉級數的基礎上加上附加項；根據邊界條件，確定附加項系數和傅里葉級數系數之間的關系；將求解域離散成獨立點，把帶有附加項的傅里葉級數代入各個獨立點，形成各傅里葉級數系數之間的關系矩陣；根據所述關系矩陣求解出傅里葉級數系數，確定需要求解的函數。優選的，附加項系數由邊界條件決定。優選的，所述需要求解的函數包括結構彈性邊界約束條件的微分方程。優選的，所述邊界條件為彈性邊界條件。優選的，所述將求解域離散成獨立點，包括：對所求解的函數區間進行離散，選取獨立點。優選的，所述把帶有附加項的傅里葉級數代入各個獨立點，形成各傅里葉級數系數之間的關系矩陣，包括：把帶有附加項的傅里葉級數代入各個獨立點，得到各個獨立點對應的方程，形成方程組；將方程組整理成矩陣形式，得到矩陣方程。優選的，所述矩陣方程中的未知系數為傅里葉級數系數。本專利技術所提供的一種基于廣義傅里葉級數構建改進微分求積法的方法，將需要求解的函數表示成傅里葉級數，并在傅里葉級數的基礎上加上附加項；根據邊界條件，確定附加項系數和傅里葉級數系數之間的關系；將求解域離散成獨立點，把帶有附加項的傅里葉級數代入各個獨立點，形成各傅里葉級數系數之間的關系矩陣；根據所述關系矩陣求解出傅里葉級數系數，確定需要求解的函數。可見，將微分求積法在求解微分方程時的精確性和高效性同廣義傅里葉級數在處理結構彈性邊界條件的靈活性相結合，所構造的改進微分求積法將在保證計算精度的同時很好的解決結構彈性邊界微分方程的求解，實現提高求解結果準確性。附圖說明為了更清楚地說明本專利技術實施例或現有技術中的技術方案，下面將對實施例或現有技術描述中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本專利技術的實施例，對于本領域普通技術人員來講，在不付出創造性勞動的前提下，還可以根據提供的附圖獲得其他的附圖。圖1為本專利技術所提供的一種基于廣義傅里葉級數構建改進微分求積法的方法的流程圖；圖2為彈性邊界梁示意圖。具體實施方式本專利技術的核心是提供一種基于廣義傅里葉級數構建改進微分求積法的方法，以實現提高求解結果準確性。為了使本
的人員更好地理解本專利技術方案，下面將結合本專利技術實施例中的附圖，對本專利技術實施例中的技術方案進行清楚、完整地描述，顯然，所描述的實施例僅僅是本專利技術一部分實施例，而不是全部的實施例。基于本專利技術中的實施例，本領域普通技術人員在沒有做出創造性勞動前提下所獲得的所有其他實施例，都屬于本專利技術保護的范圍。請參考圖1，圖1為本專利技術所提供的一種基于廣義傅里葉級數構建改進微分求積法的方法的流程圖，該方法包括：S11：將需要求解的函數表示成傅里葉級數，并在傅里葉級數的基礎上加上附加項；S12：根據邊界條件，確定附加項系數和傅里葉級數系數之間的關系；S13：將求解域離散成獨立點，把帶有附加項的傅里葉級數代入各個獨立點，形成各傅里葉級數系數之間的關系矩陣；S14：根據所述關系矩陣求解出傅里葉級數系數，確定需要求解的函數。可見，將微分求積法在求解微分方程時的精確性和高效性同廣義傅里葉級數在處理結構彈性邊界條件的靈活性相結合，所構造的改進微分求積法將在保證計算精度的同時很好的解決結構彈性邊界微分方程的求解，實現提高求解結果準確性。基于上述方法，具體的，附加項系數由邊界條件決定。其中，在傅里葉級數的基礎上加上附加項就構成了廣義傅里葉級數。進一步的，所述需要求解的函數包括結構彈性邊界約束條件的微分方程。其中，所述邊界條件為彈性邊界條件。進一步的，步驟S13中，所述將求解域離散成獨立點的過程具體為：對所求解的函數區間進行離散，選取獨立點。進一步的，步驟S13中，把帶有附加項的傅里葉級數代入各個獨立點，形成各傅里葉級數系數之間的關系矩陣的過程具體為：把帶有附加項的傅里葉級數代入各個獨立點，得到各個獨立點對應的方程，形成方程組；將方程組整理成矩陣形式，得到矩陣方程。其中，所述矩陣方程中的未知系數為傅里葉級數系數。矩陣方程為微分方程對應的矩陣方程。具體的，將微分方程中所求函數表述成傅里葉級數加上附加項的形式，其中附加項的系數由邊界條件決定。對所求區間進行離散，選取獨立點。將廣義傅里葉級數帶入各離散獨立點，得到各個獨立點對應的方程。將方程組整理成矩陣形式，得到微分方程對應的矩陣方程，方程中未知系數即傅里葉級數系數。其中，根據彈性邊界條件控制方程確定附加項系數同傅里葉級數系數間的關系。將附加項級數用傅里葉級數的系數表示，并將所求函數表示成只含傅里葉級數系數的形式。將所求函數區間離散成各個獨立點。將廣義傅里葉級數帶入各離散獨立點，帶進各離散獨立點，得到關于傅里葉級數系數的方程組。求解方程組，得到傅里葉系數，從而確定所求函數。詳細的，提供本方法應用于兩端彈性約束梁控制微分方程求解的例子。基于本方法，需要求解的函數為梁的振動控制方程。例如，圖2所示梁的振動控制方程是：邊界條件是：Dy”＝K0y'Dy”'＝-k0yx＝0(2a)Dy”＝K1y'Dy”'＝-k1yx＝L(2b)其中D＝EI，設梁的位移函數為：其中ξi(x)為容許函數，吸收傅里葉級數方法邊界不連續，加快傅里葉級數的收斂性。將梁的位移函數帶進邊界條件，得到如下式子：其中參數C1…C8為：將式(5)寫成如下的矩陣形式：其中，Qm和A的表達式如下：A＝{a0a1a2…aM}T(9)令R如下：則得到：從式(11)可以看出，系數b和a不是獨立的變量，系數b可以由系本文檔來自技高網...

【技術保護點】
一種基于廣義傅里葉級數構建改進微分求積法的方法，其特征在于，包括：將需要求解的函數表示成傅里葉級數，并在傅里葉級數的基礎上加上附加項；根據邊界條件，確定附加項系數和傅里葉級數系數之間的關系；將求解域離散成獨立點，把帶有附加項的傅里葉級數代入各個獨立點，形成各傅里葉級數系數之間的關系矩陣；根據所述關系矩陣求解出傅里葉級數系數，確定需要求解的函數。

【技術特征摘要】
1.一種基于廣義傅里葉級數構建改進微分求積法的方法，其特征在于，包括：將需要求解的函數表示成傅里葉級數，并在傅里葉級數的基礎上加上附加項；根據邊界條件，確定附加項系數和傅里葉級數系數之間的關系；將求解域離散成獨立點，把帶有附加項的傅里葉級數代入各個獨立點，形成各傅里葉級數系數之間的關系矩陣；根據所述關系矩陣求解出傅里葉級數系數，確定需要求解的函數。2.如權利要求1所述的方法，其特征在于，附加項系數由邊界條件決定。3.如權利要求1所述的方法，其特征在于，所述需要求解的函數包括結構彈性邊界約束條件的微分方...

【專利技術屬性】
技術研發人員：王剛，
申請(專利權)人：蘇州大學，
類型：發明
國別省市：江蘇,32