蒙特卡洛法解最优化

moutqvw

1. %MonteCarlo法解最优化的例
   
2. clear;
   
3. %设定一个初值
   
4. vmax=-inf;
   
5. %specrnd生成给定分布律的随机数
   
6. %r=specrnd(x,p)返回一个来自分布律P(x)=p的随机数
   
7. %r=specrnd(x,p,m,n)返回m*n随机数矩阵
   
8. %p的默认值为等概率
   
9. %例如 分布律
   
10. %      ξ         7      8      9     10
    
11. %      p         0.1   0.2   0.4    0.3
    
12. %命令    specrnd([7 8 9 10],[0.1 0.2 0.4 0.3],10,10)
    
13. %产生100个模拟数据  
    
14. x2=specrnd(10:20,[],1,5);
    
15. x3=specrnd(-5:16,[],1,5);
    
16. for i=1:5
    
17.    for j=1:5
    
18.       if x2(i)+2*x3(j)>=10&3*x2(i)+2*x3(j)<=62,
    
19.          v=(x2(i)+10)*x2(i)*x3(j);
    
20.          if v>vmax,
    
21.             vmax=v;x20=x2(i);x30=x3(j);
    
22.          end
    
23.       end
    
24.    end
    
25. end
    
26. x1=x20+10,x2=x20,x3=x30
    
27. vmax
    

复制代码