*- 使用stata计算过程
   
    reg z x y
       
        reg z x y, nocons
       
        constraint def 1 x + y = 1
        constraint def 2 x > 0
        constraint def 3 x < 1       
        constraint def 4 y > 0
        constraint def 5 y < 1          
   
        cnsreg z x y, nocons constraint(1-5)   
   
   
   *- 结果说明：目前stata的constraint命令目前不支持定义系数的范围，因此只能用优化算法。
   
*======================================================================
   
   *- 优化算法的基本步骤
   
   
     clear
     set obs 1000

     gen beta = .    // x前的系数
         gen rho = .     // M系列与Z的相关系数

     local range = 1
         dis `range'
     local j = 1

     forvalues beta = 0(0.001)`range'{
       qui replace beta = `beta' in `j'
          
           preserve
            use dd.dta, clear
                gen m = `beta' * x + (1 - `beta') * y
                corr m z
                dis `r(rho)'
                drop m
           restore  
          
       qui replace rho = `r(rho)'  in `j++'
  
     }
         *
     * 测算结果：发现x前面系数为1时，M和Z的相关系数最大。
         
          gsort -rho  //请看第一行！
          
         * 原始数据x与z高度相关!! 请看下面命令的结果。
         
          use dd.dta, clear
          corr x y z
          pwcorr_a x y z

求受限的最大相关系数的excel数据

Grid search

通过哪个软件，如何做呢?非常感谢！

可以试一下这个方法：
用Eviews软件，做以下OLS回归： ls z=c(1)*x+(1-c(1))*y