myhilb.m

来自「控制系统计算机辅助设计——MATLAB语言与应用(源代码）」· M 代码 · 共 18 行

function A=myhilb(n, m)
%MYHILB  本函数用来演示 MATLAB 语言的函数编写方法。
%   A=MYHILB(N, M) 将产生一个 N 行 M 列的 Hilbert 矩阵 A；
%   A=MYHILB(N) 将产生一个 NxN 的方 Hilbert 阵 A;
%
%See also: HILB.

%  Designed by Professor Dingyu XUE, Northeastern University, PRC
%     5 April, 1995, Last modified by DYX at 30 July, 2001
if nargout>1, error('Too many output arguments.'); end
if nargin==1, m=n;   % 若给出一个输入，则生成方阵
elseif nargin==0 | nargin>2
   error('Wrong number of iutput arguments.');
end
for i=1:n, for j=1:m
   A(i,j)=1/(i+j-1);
end, end