科学和社会的不断发展带来很多新的现象和问题，经典的数学模型及其分析方法表现出滞后和不足。一方面，这促使一些原本不被熟知的数学理论及描述方法被引入到各个应用领域中，另一方面，自然界及科学领域中的随机现象无处不在，在建模和分析中有必要考虑随机因素。本文的研究对象是含分数阶导数项的随机系统和含时滞项的随机系统，重点分析其中导数的非整数阶阶数及时滞对系统响应的影响，主要的工作和结论如下：针对分数阶导数的定义对含分数阶导数的随机系统提出数值模拟的实用方法。在Caputo定义下，基于边界单元法的思想，引入一组线性的、非耦合的随机分数阶模拟方程，借助Laplace变换，得到含多个分数阶导数随机系统数值求解的半解析方法；对Grünwald-Letnikov定义，给出GL系数序列的概念，利用其随项数迅速递减的特点对某一时刻分数阶导数的数值计算进行合理截断，减低其对历史数据的长期依赖性，结合经典算法得到含分数阶导数随机系统响应分析的数值方法。基于频域分析对随机系统中的分数阶导数建立全面认识。细致地讨论了Gauss白噪声激励下含分数阶导数线性系统响应的功率谱密度，将分数阶导数分解为阻尼与刚度的线性组合，给出组合系数随导数阶数的演化过程，阐明分数阶导数所描述的粘弹性力在系统中同时表现为阻尼和刚度。应用随机平均法，得到响应的Markov近似，建立概率密度函数的Fokker-Planck-Kolmogorov（FPK）方程并求出稳态解。利用改进的统计线性化法对Gauss白噪声激励下含分数阶导数的Duffing系统进行稳态响应分析。借助广义谐波平衡技术得到等效线性系统，并将分数阶导数分解为阻尼和刚度的线性组合。条件功率谱密度与随机平均法给出的振幅稳态概率密度完成响应的功率谱密度估计，结合分数阶导数分解的组合系数进行频域分析，说明分数阶导数所描述的粘弹性力同时蕴含阻尼和刚度，且等效线性化过程保留了Duffing系统非线性刚度的影响。对关联Gauss白噪声激励下的多时滞系统做出稳态响应分析。推广单时滞情形的结果，利用泛函计算，借助Novikov定理，给出支配响应概率密度演化的时滞FPK方程的推导过程，应用小时滞近似技术给出时滞FPK方程的稳态解，详细分析参数变化对稳态概率密度的影响，说明关联Gauss白噪声激励下多时滞系统中噪声关联强度与时滞对系统响应产生的不同影响。