matlab如何跳转程序,如何用matlab编写如下程序？急求

谢谢版主的肯定，我再补充一下，K=6,也就是求6阶一步转移概率矩阵的36个元素，N也即我选取研究的基金的种类，目前我暂定为10，Si,j(n)我也通过数据数出来了，取值是0到36，且Si,1(n)+Si,2(n)+........Si,6(n)=36.  而当Si,j(n)=0时， 说明(πi,j) ^ Si,j(n)=1,即为常数项。但要求36个πi,j,我想Si,j(1)+Si,j(2+...........Si,j(6)不能等于0，否则该πi,j在目标max函数中就为常数项，或者由约束条件，对于每个i=k(k属于1,2,3,4,5,6)时，其中最多允许有一个πk,j在目标max函数的系数为0，不然无法由约束条件每行之和=1解出。大家认为捏？

我数出了Si,j(n)构成的矩阵，可能是由于我选取的基金太少，只有10组，而且特点鲜明，很多都是待在一个状态好多时期都不转移的。还有时间跨度可能不够，我只选取了3年，以每个月有一个标准时间段，状态也可以告诉大家，即每只基金的评级星级，5星=6,4星=5,3星=4，。。。。。0星=1.导致我有好几个Si,j(1)+Si,j(2)+.......Sij(6)=0，我想这很失败。目标函数的表达式没有这么多的未知元πi,j,光由2个约束条件是无法求出这些系数为0的项的值的。我想再多加机组，扩大一下时间段，争取36个未知元的系数都不为0.

由于还要时间继续收集数据，我先举个例子让大家想想该怎么实现

比如MAX函数的表达式为πi,j=argmax[(π1,1^1+π1,1^3)+(π1,2^4+5+π1,2^7)+(π1,3^x+.........)+(π1,4^y+.......)+.........+(π6,6^z+.......+)]。

再由我的2楼上附件图片中的约束条件约束，该如何用SQP方法求解捏？我知道matlab中好像某个函数中有SQP选项。

另外我会尽快附上我的真实的Si,j(n)构成的矩阵，代入就得到了一个确定的目标函数表达式，可能n和各个未知量的系数次数不同，但是我给大家举的例子就是想请教一该如何用matlab实现这一过程的大致构想。

非常感谢各位关注的朋友和思考解答我的问题的朋友。