怎么用MATLAB实现下述计算-经管之家官网！

T0＝44.73，19.97，17.39，17.9
T1＝44.93，19.48，17.74，17.85
T2＝46.35，19.29，17.3，17.07
T3＝44.84，19.25，18.26，17.65
上述数据分别为02，03，04，05年四个地区所占比例的数据，现在要计算一步转移概率矩阵P，即T0*P=T1,T1*P=T2,T2*P=T3，其约束条件是：
依据残差平方和最小的思想，
第一步：即t年的状态j的概率＝t-1年状态i的概率*从状态i向状态j的转移概率，然后对所有的i求和；
第二步：用t年的j状态的比例减去上一步所求的和，对这个差求平方，然后对所有的j求和，即求平方和，让其最小。
同时，要求所求概率矩阵的各行元素和为1，各元素大于0小于1。
用一个笨方法来解的话是这样：
设转移概率矩阵为：P11 P12 P13 P14
P21 P22 P23 P24
P31 P32 P33 P34
P41 P42 P43 P44
t1年状态1的差为：44.93－(17.97*P21+17.39*P31+17.9*P41)
状态2的差为：19.48－（44.73*P12+17.39*P32+17.9*P42）
状态3的差为：17.74－（44.73*P13+19.97*P23+17.9*P43）
状态4的差为：17.85－（44.73*P14+19.97*P24+17.39*P34）
t2年状态1的差为：46.35－(19.48*P21+17.74*P31+17.85*P41)
状态2的差为：19.29－（44.93*P12+17.74*P32+17.85*P42）
状态3的差为：17.3－（44.93*P13+19.48*P23+17.85*P43）
状态4的差为：17.07－（44.93*P14+19.48*P24+17.74*P34）
t3年状态1的差为：44.84－(19.29*P21+17.3*P31+17.07*P41)
状态2的差为：19.25－（46.35*P12+17.3*P32+17.07*P42）
状态3的差为：18.26－（46.35*P13+19.29*P23+17.07*P43）
状态4的差为：17.65－（46.35*P14+19.29*P24+17.3*P34）
然后对上面所有的差分别求平方，再求和，让这个和最小，即残差平方和最小。
用这个约束条件进行求解，应该怎么运用MATLAB写语言呢？
正确结果应该是
P( 1, 1) 0.2569220
P( 1, 2) 0.1734723
P( 1, 3) 0.3918807
p( 1, 4) 0.1777250
P( 2, 1) 0.000000
p( 2, 2) 0.5165667
P( 2, 3) 0.000000
P( 2, 4) 0.4834333
P( 3, 1) 1.000000
P( 3, 2) 0.000000
P( 3, 3) 0.000000
P( 3, 4) 0.000000
P( 4, 1) 0.9227927
P( 4, 2) 0.7720731E-01
P( 4, 3) 0.000000
P( 4, 4) 0.000000

「经管之家」APP：经管人学习、答疑、交友，就上经管之家！
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP