用Matlab求方程组x'=A*x+f(t)满足初始条件的特解g(t)，其中，矩阵A=[1 0 3;8 1 -1;5 1 -1],f(t)=[sint;cost;t],g(0)=[-1;1;-1],写出源代码

以下是使用Matlab求解方程组x'=A*x+f(t)满足初始条件的特解g(t)的代码：

% 定义矩阵A和向量f(t)
A = [1 0 3; 8 1 -1; 5 1 -1];
f = @(t) [sin(t); cos(t); t];

% 定义初始条件g(0)
g0 = [-1; 1; -1];

% 定义ODE函数
odefun = @(t, x) A * x + f(t);

% 求解ODE方程
[t, x] = ode45(odefun, [0, 10], g0);

% 绘制结果
plot(t, x(:,1), t, x(:,2), t, x(:,3));
legend('x1', 'x2', 'x3');
xlabel('t');

解释：首先定义矩阵A和向量f(t)，然后定义初始条件g(0)。接着定义ODE函数，其中x是状态向量，t是时间，使用MATLAB内置的ode45函数求解ODE方程，并绘制结果图像。