超松弛迭代中“下限约束”怎样在matlab中实现

sunnyshore

在超松弛迭代程序中，有迭代初值。现在我想实现：给定一个下限值，每步迭代后若得到的数值小于下限值，则让其等于给定的下限！若得到的数值大于给定的下限值，则继续迭代计算。请问：怎样把约束下限加入到超松弛迭代程序中？

1. function y=sor(a,b,w,x0)
   
2. a=[2275 1224.4 474.53;1224.4 679.45 263.24;474.53 263.24 102.49];
   
3. b=[42422 23347 9063]';
   
4. x0=[0;0;0];
   
5. w=1.6496;
   
6. D=diag(diag(a));
   
7. U=-triu(a,1);
   
8. L=-tril(a,-1);
   
9. lw=(D-w*L)\((1-w)*D+w*U);
   
10. f=(D-w*L)\b*w;
    
11. y=lw*x0+f;n=1;
    
12. while norm(y-x0)>=1.0e-6
    
13. x0=y;
    
14. y=lw*x0+f;n=n+1;
    
15. end
    
16. y
    
17. n

复制代码

这是超松弛迭代程序，假设给定下限为2.0 ，5.0，11.0