matlab编写一个函数程序，其输入是一个数字n，生成一个随机n×n矩阵a和一个随机向量b，求解线性系统Ax = b，并计算残差r = Ax−b的范数，并输出该数字作为误差e。为n = 5、10、50...

首先，我们需要了解如何在 Matlab 中生成随机矩阵和向量。

我们可以使用 Matlab 内置的 rand 函数来生成随机矩阵。这个函数会生成一个 0 到 1 之间的随机数组。我们可以通过将其乘以一个常数来调整数组的大小，并使用 randn 函数来生成服从正态分布的随机数组。

例如，下面的代码会生成一个 3x3 的随机矩阵：

A = rand(3)

同样，我们可以使用 rand 函数生成一个随机向量，例如：

b = rand(3, 1)

这将生成一个 3x1 的随机向量。

接下来，我们需要求解线性系统 $Ax = b$。这可以使用 Matlab 中的 \ 运算符来实现，它会自动调用矩阵分解算法(如 LU 分解)来求解线性方程组。例如，我们可以使用下面的代码来求解：

x = A \ b

接下来，我们需要计算残差 $r = Ax - b$ 的范数，并输出这个数字作为误差 $e$。我们可以使用 Matlab 中的 norm 函数来计算向量的范数。例如，我们可以使用下面的代码来计算 $r$ 的范数：

r = A * x - b
e = norm(r)