采用回溯法,利用动态搜索树的数据结构来构造一棵解树,实现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
上传用户:小鹏