[watcher]

函数tril(...)和triu(...)分别产生下三角和上三角矩阵，建议先看下帮助

可以这样写程序
b=triu(ones(4,1)*a,0);