代码:
%得到如下形式的矩阵
% [ 1 2 3 4
% 1/2 1 2 3
% 1/3 1/2 1 2
% 1/4 1/3 1/2 1 ]
%EDIT BY PSL@CSU
%QQ:547423688
%Email:[email protected]
function M=get_M(n)
%M=get_M(4)可得到如上矩阵
M=zeros(n);
for i=1:n
for j=1:i-1
M(i,j)=1/(i+1-j);
end
M(i,i:n)=[1:n-i+1];
end
%EDIT BY PSL@CSU
%QQ:547423688
%Email:[email protected]