z=21，m=1.5，a=20°，s= linspace(13.125, 13.96,30)，用matlab计算g=（0.5dsin a）的平方-(0.5*d)的平方+s的平方，并绘制出曲线

可以使用MATLAB的向量化操作，快速计算并绘制出曲线。具体实现代码如下：

z = 21;
m = 1.5;
a = 20;
s = linspace(13.125, 13.96, 30);

d = m * z ./ s;  % 计算 d
g = (0.5 * d .* sin(deg2rad(a))).^2 - (0.5 * d).^2 + s.^2;  % 计算 g

plot(s, g);  % 绘制曲线
xlabel('s');
ylabel('g');

解释一下代码的实现过程：

首先，根据公式 $d = \frac{mz}{s}$，计算出向量 `d`。其中，`m`、`z`、`s` 都是向量，所以 `d` 也是向量。

然后，根据公式 $g = (\frac{1}{2}d\sin a)^2 - (\frac{1}{2}d)^2 + s^2$，计算出向量 `g`。其中，`d`、`a`、`s` 都是向量，所以 `g` 也是向量。

最后，使用 `plot` 函数绘制出曲线，横轴为 `s`，纵轴为 `g`。