已知
R是一个k×k阶的矩阵，其中每个元素Pij可以表示为：
Pij=Ф[((j*d)-(1-r)*(i-0.5)*d)/(r)]- Ф[((j-1)*d-(1-r)*(i-0.5)*d)/(r)] ；Ф为标准正态概率分布函数，其中，i,j=1,2,…,k;r=0.1
Y=P(I-R)-11
其中，P=[0,…,1,…,0]1×k,P(1,k/2)=1;I为k阶单位阵，1为所有元素均为1的k×1的列向量.
现在已知y的值为500，
对Y进行矩阵计算得到Rk×k=I-1/ARL*Ek×1*P1×k
这样我得到了R矩阵，已知里面任意一个元素，再根据Pij的计算表达式返回去求解d的值，是不是就转化成 这样一个表达式（m1,m2和a分别为常数），我想用通过二分法等求解非线性方程的方法求出d的值，但是首先得找到关于d的目标函数，麻烦各位帮忙看看，我对于matlab了解不多，不知道这个思路对不对？恳请各位指点！！！