利用matlab根据非定常传热方程求沿x轴速度v移动的基模高斯激光辐照岩石温度场和应力场，

由于题目没有给出具体的非定常传热方程和岩石材料参数，以下是一种可能的求解方法，供参考。

1. 假设岩石材料是均匀、各向同性的，并且没有相变或化学反应等复杂的过程影响，可以采用热传导方程进行求解：

$$
\rho c_p \frac{\partial T}{\partial t} = \nabla \cdot (k \nabla T) + Q(x,t)
$$

其中，$T(x,y,z,t)$是温度场，$\rho$是岩石密度，$c_p$是比热容，$k$是热导率，$Q(x,t)$是激光辐照时产生的热源项，可以根据高斯激光功率密度公式计算：

$$
Q(x,t) = \frac{P_0}{\pi w^2} e^{-\frac{2(x-vt)^2}{w^2}} \delta(t)
$$

其中，$P_0$是激光功率，$w$是激光束腰半径，$\delta(t)$是热源项的脉冲函数。

2. 引入应力场的影响，可以将上述方程改为热弹性耦合方程：

$$
\rho c_p \frac{\partial T}{\partial t} = \nabla \cdot (k \nabla T) + Q(x,t) + \nabla \cdot \boldsymbol{q}
$$

$$
\boldsymbol{q} = -\lambda \nabla T + \boldsymbol{s}
$$

其中，$\boldsymbol{q}$是热流密度，$\lambda$是热导率，$\boldsymbol{s}$是应力场对热流的修正项，可以用胡克定律计算：

$$
\boldsymbol{s} = -\frac{\alpha E}{1-2\nu} \boldsymbol{\varepsilon} - 2\alpha \mu \boldsymbol{\varepsilon} + \alpha \beta (T-T_0) \boldsymbol{I}
$$

其中，$E$是弹性模量，$\nu$是泊松比，$\alpha$是热膨胀系数，$\mu$是剪切模量，$\beta$是岩石材料的体积膨胀系数，$T_0$是参考温度，$\boldsymbol{\varepsilon}$是应变张量，$\boldsymbol{I}$是单位张量。

3. 用matlab编写程序求解上述方程组。可以采用有限差分法或有限元法等数值方法进行离散化，然后用显式或隐式格式进行时间推进。需要设置合适的边界条件和初值条件，例如，可以将岩石表面设置为绝热边界，将岩石内部初始温度设为室温，将应力场设为零。根据求解结果，可以绘制出沿x轴的温度和应力分布图，以及时间演化图。

4. 进一步分析结果，可以计算岩石的热应力和热应变，评估岩石的热稳定性和破坏机制，以及激光加热对岩石物理性质的影响。还可以通过与实验数据对比，验证模型的准确性和可靠性，以及优化模型参数和求解方法，提高模拟效率和精度。