如何用matlab绘制椭圆极化电场矢量随时间空间变化 的图像，请给出代码

以下是用MATLAB绘制椭圆极化电场矢量随时间空间变化的代码示例：

% 定义电场参数
E0 = 1; % 电场振幅
omega = 2*pi*1e9; % 角频率
k = 2*pi/800e-9; % 波矢
phi = 0; % 初相位
eps0 = 8.854e-12; % 真空介电常数
c = 1/sqrt(eps0); % 光速

% 定义空间和时间范围
x = linspace(-4e-6, 4e-6, 100); % x 轴范围
y = linspace(-4e-6, 4e-6, 100); % y 轴范围
t = linspace(0, 20e-12, 100); % 时间范围

% 生成网格
[X, Y, T] = meshgrid(x, y, t);

% 计算电场矢量
Ex = E0*cos(k*X - omega*T + phi); % x 方向分量
Ey = E0*cos(k*Y - omega*T + phi); % y 方向分量
Ez = 0; % z 方向分量

% 计算电场强度和极化方向
E = sqrt(Ex.^2 + Ey.^2 + Ez.^2); % 电场强度
theta = atan2(Ey, Ex); % 极化方向

% 绘制电场矢量图
quiver3(X, Y, T, Ex./E, Ey./E, Ez, 'color', 'b', 'linewidth', 1.5);

% 设置坐标轴和标题
xlabel('x (m)');
ylabel('y (m)');
zlabel('t (s)');
title('椭圆极化电场矢量随时间空间变化');

% 设置图形参数
axis tight;
grid on;
view(30, 30);

在上述代码中，我们首先定义了电场的参数，包括电场振幅、角频率、波矢等。然后定义了空间和时间范围，以及生成了网格。接着，我们根据电场公式计算了电场矢量的 x、y 方向分量，并设置了 z 方向分量为 0。然后，我们计算了电场强度和极化方向，并使用 `quiver3` 函数绘制了电场矢量图。最后，我们设置了坐标轴、标题和图形参数。

运行上述代码，即可得到椭圆极化电场矢量随时间空间变化的图像。