采用回溯法,利用動態搜索樹的數據結構來構造一棵解樹,實現3著色問題。 本程序使用迭代回溯算法
標簽: 回溯法
上傳時間: 2017-09-26
上傳用戶:bruce5996
共軛梯度法為求解線性方程組而提出。后來,人們把這種方法用于求解無約束最優化問題, 使之成為一種重要的最優化方法。 共軛梯度法的基本思想是把共軛性與最速下降方法相結合, 利用已知點處的梯度構造一組共 軛方向, 并沿這組方向進行搜索, 求出目標函數的極小點。 根據共軛方向的基本性質, 這種 方法具有二次終止性。 在各種優化算法中, 共軛梯度法是非常重要的一種。 其優點是所需存 儲量小,具有步收斂性,穩定性高,而且不需要任何外來參數。 共軛方向 無約束最優化方法的核心問題是選擇搜索方向 . 在本次實驗中 , 我們運用基于共軛方向的一種 算法 — 共軛梯度法 三.算法流程圖: 四.實驗結果: (1). 實驗函數 f=(3*x1-cos(x2*x3)-1/2)^2+(x1^2-81*(x2+0.1)+sin(x3)+1.06)^2+(exp(-x1*x2)+20*x3+ 1/3*(10*3.14159-3))^2; 給定初始點 (0,0,0) , k=1 ,最 大迭代次數 n ? ? d 確定搜索方向 進 退 法 確 定 搜 索 區 間 分割法確定最 優步長
上傳時間: 2016-05-08
上傳用戶:saren11
function [alpha,N,U]=youxianchafen2(r1,r2,up,under,num,deta) %[alpha,N,U]=youxianchafen2(a,r1,r2,up,under,num,deta) %該函數用有限差分法求解有兩種介質的正方形區域的二維拉普拉斯方程的數值解 %函數返回迭代因子、迭代次數以及迭代完成后所求區域內網格節點處的值 %a為正方形求解區域的邊長 %r1,r2分別表示兩種介質的電導率 %up,under分別為上下邊界值 %num表示將區域每邊的網格剖分個數 %deta為迭代過程中所允許的相對誤差限 n=num+1; %每邊節點數 U(n,n)=0; %節點處數值矩陣 N=0; %迭代次數初值 alpha=2/(1+sin(pi/num));%超松弛迭代因子 k=r1/r2; %兩介質電導率之比 U(1,1:n)=up; %求解區域上邊界第一類邊界條件 U(n,1:n)=under; %求解區域下邊界第一類邊界條件 U(2:num,1)=0;U(2:num,n)=0; for i=2:num U(i,2:num)=up-(up-under)/num*(i-1);%采用線性賦值對上下邊界之間的節點賦迭代初值 end G=1; while G>0 %迭代條件:不滿足相對誤差限要求的節點數目G不為零 Un=U; %完成第n次迭代后所有節點處的值 G=0; %每完成一次迭代將不滿足相對誤差限要求的節點數目歸零 for j=1:n for i=2:num U1=U(i,j); %第n次迭代時網格節點處的值 if j==1 %第n+1次迭代左邊界第二類邊界條件 U(i,j)=1/4*(2*U(i,j+1)+U(i-1,j)+U(i+1,j)); end if (j>1)&&(j U2=1/4*(U(i,j+1)+ U(i-1,j)+ U(i,j-1)+ U(i+1,j)); U(i,j)=U1+alpha*(U2-U1); %引入超松弛迭代因子后的網格節點處的值 end if i==n+1-j %第n+1次迭代兩介質分界面(與網格對角線重合)第二類邊界條件 U(i,j)=1/4*(2/(1+k)*(U(i,j+1)+U(i+1,j))+2*k/(1+k)*(U(i-1,j)+U(i,j-1))); end if j==n %第n+1次迭代右邊界第二類邊界條件 U(i,n)=1/4*(2*U(i,j-1)+U(i-1,j)+U(i+1,j)); end end end N=N+1 %顯示迭代次數 Un1=U; %完成第n+1次迭代后所有節點處的值 err=abs((Un1-Un)./Un1);%第n+1次迭代與第n次迭代所有節點值的相對誤差 err(1,1:n)=0; %上邊界節點相對誤差置零 err(n,1:n)=0; %下邊界節點相對誤差置零 G=sum(sum(err>deta))%顯示每次迭代后不滿足相對誤差限要求的節點數目G end
標簽: 有限差分
上傳時間: 2018-07-13
上傳用戶:Kemin
機翼極限環振蕩(LCO)是典型的非線性氣動彈性問題,嚴重的會造成機翼的結構破壞。為了精確捕捉極限環振蕩初始臨界點,準確預測極限環的幅值,為機翼的設計提供準確的數據參考,本文綜合考慮了氣動與結構非線性的影響,提出了一種松耦合氣動彈性仿真方法。在子迭代過程中分別采用LUSGS雙時間推進和多步推進法交替求解氣動和結構動力學方程;一種高效的插值技術應用于耦合界面數據的映射與傳遞;采用精確動網格技術模擬氣體的非定常流動。對標準模型切尖三角翼的跨音速極限環振蕩的計算與分析,表明相比同類仿真方法,通過此方法得到結果與實驗值吻合更好;證明了結構幾何非線性與氣動非線性是誘發LCO的重要原因。耦合仿真方法保真度高,能為強非線性結構的強度設計提供重要依據。
上傳時間: 2013-10-22
上傳用戶:m62383408
文中詳細地介紹了正交投影子空間跟蹤算法(OPAST),它是一種基于最優化問題的方法,保證了每次迭代時權向量的正交性,并具有和PAST算法一樣的線性復雜度,以及與自然冪法(NP)一樣的全局收斂性。然而將其應用于盲多用戶檢測時,在迭代一定次數后,會出現誤碼率突然增大現象,這就導致了算法性能的下降,為了解決這一問題,文中提出一種方法,并通過仿真結果,證明它是行之有效的。
上傳時間: 2014-11-11
上傳用戶:xaijhqx
根據ID 號找出不同照片中的同名編碼點, 順次對照片進行相對定向和絕對定向。然后,根據多幅圖像的多極線幾何約束, 實現非編碼點的匹配, 消除誤匹配。采用前方交會法重建標志點的三維坐標, 利用光束平差對計算出的結果和內外部參數做迭代修正.
上傳時間: 2013-10-19
上傳用戶:love_stanford
CT重建算法,迭代算法ARTII,乘型ART算法,該算法法比初始值選擇不能為0
標簽: 重建算法
上傳時間: 2014-01-11
上傳用戶:zhouli
線性方程組的數值解法,包括高斯消去,三角分解,追趕法和Jacobi等迭代方法
上傳時間: 2013-12-26
上傳用戶:問題問題
優點:1.對于信噪比高的信號濾波效果好; % 2.對于邊沿的保護強過閾值濾波,不會產生閾值濾波情況下的過于平滑與Gibbs現象。 %缺點:1.由于對邊沿信號沒做任何處理,所以邊沿可能會有脈沖噪聲保留下來; % 2.計算相關系數中,如果計算出來的小波系數點位置偏差大,則相關系數計算受影響; % 3.需要迭代運算,迭代的噪聲能量閾值選取很重要,這里以開始段無信號處估計噪聲; % 4.需要迭代運算,所以運算量比閾值法大; % 5.受分解層次影響,在大尺度上小波系數點位置偏差更大,相關系數計算不準確。 %需要具體調整的地方:1.分解的尺度;
上傳時間: 2013-12-26
上傳用戶:wlcaption
數值線性代數的Matlab應用程序包 共13個程序函數,每個程序函數有相應的例子函數一一對應,以*Example.m命名 程序名稱 用途 Method 方法 GrmSch.m QR因子分解 classical Gram-Schmidt orthogonalization 格拉母-斯密特 MGrmSch.m QR因子分解 modified Gram-Schmidt iteration 修正格拉母-斯密特 householder.m QR因子分解 Householder 豪斯霍爾德QR因子分解 ZXEC.m 最小二乘擬合 polynomial interpolant 最小二乘插值多項式 NCLU.m LU因子分解 Gaussian elimination 不選主元素的高斯消元 PALU.m LU因子分解 partial pivoting Gaussian elimination 部分選主元的高斯消元 cholesky.m 楚因子分解 Cholesky Factorization 楚列斯基因子分解 PwItrt.m 求最大特征值 Power Iteration 冪迭代 Jacobi.m 求特征值 Jacobi iteration 按標準行方式次序的雅可比算法 Anld.m 求上Hessenberg Arnoldi Iteration 阿諾爾迪迭代 zuisu.m 解線性方程組 Steepest descent 最速下降法 CG.m 解線性方程組 Gradients 共軛梯度 BCG.m 解線性方程組 Biconjugate Gradients 雙共軛梯度
上傳時間: 2016-05-17
上傳用戶:小鵬