边坡稳定性计算方法 目前的边坡的侧压力理论，得出的计算结果，显然与实际情形不符。边坡稳定性计算，有直线法和圆弧法，当然也有抛物线计算方法，这些不同的计算方法，都做了不同的假设条件。 当然这些先辈拿出这些计算方法之前，也曾经困惑，不做假设简化，基本无法计算。而根据各种假设条件，是会得出理论上的结果，但与实际情况又不符。倒是有些后人不管这些假设条件，直接应用其计算结果，把这些和实际不符的公式应用到现有的规范和理论中。 瑞典条分法，其中的一个假设条件破裂面为圆弧，另一个条件为假设的条间土之间，没有相互作用力，这样的话，对每一个土条在滑裂面上进行力学分解，然后求和叠加，最后选取系数最小的滑裂面。从而得出判断结果。其实，那两个假设条件对吗？都不对！ 第一、土体的实际滑动破裂面，不是圆弧。第二、假设的条状土之间，会存在粘聚力与摩擦力。边坡的问题看似比较简单，只有少数的几个参数，但是，这几个参数之间，并不是线性相关。对于实际的边坡来讲，虽然用内摩擦角①和粘聚力C来表示，但对于不同的破裂面，破裂面上的作用力，摩擦力和粘聚力，都是破裂面的函数，并不能用线性的方法分别求解叠加，如果是那样，计算就简单多了。 边坡的破裂面不能用简单函数表达，但是，如果不对破裂面作假设，那又无从计算，直线和圆弧，是最简单的曲线，所以基于这两种曲线的假设，是计算的第一步，但由于这种假设与实际不符，结果肯定与实际相差甚远。
条分法的计算，是来源于微积分的数值计算方法，如果条间土之间，存在相互作用力，那对条状土的力学分解，又无法进行下去。 所以才有了圆弧破裂面的假设与忽略条间土的相互作用的假设。 其实先辈拿出这样与实际不符的理论，内心是充满着矛盾的。 实际看到的边坡的滑裂，大多是上部几乎是直线，下部是曲线形状，不能用简单函数表示，所以说，要放弃求解函数表达式的想法。计算还是可以用条分法，但要考虑到条间土的相互作用。 用微分迭代的方法求解，能够得出近似破裂面，如果每次迭代，都趋于收敛，那收敛的曲线，就是最终的破裂面。 参照图3，下面将介绍这种方法的求解步骤。
基于Matlab的MMSE的语音增强算法的研究 投递人发布于2013-06-27 17:05 评论(0)有 5 人阅读 本课题隶属于学校的创新性课题研究项目。2012年就已经做完了，今天一并拿来发表。 目录： --基于谱减法的语音信号增强算法 (1) 一：语音增强技术概述 (3) 二：语音增强的目的 (3) 三：语音信号的相关特性 (3) 1 语音特性 (3) 1. 1 语音信号具有短时平稳性 (3) 1.2.语音信号可以分为浊音和清音 (3) 1.3.语音信号可以利用统计分析特征描述 (4) 2 人耳感知特性 (4) 3 噪声特性 (4) 3.1周期性噪声 (4) 3.2脉冲噪声 (5) 3.3宽带噪声 (5) 3.4同声道语音干扰 (5) 3.5传输噪声 (5) 四：国内外有关抗噪声技术的解决方案 (5) 1 语音增强算法 (5) (1) 基于谱减法的语音增强 (6) (2) 自适应滤波法的语音增强 (6) (3) 短时对数谱的MMSE语音增强 (6) 2 寻找稳健的语音特征作为特征参数 (6)
3 基于模型参数自适应的噪声补偿算法 (6) 五：语音增强算法的三种具体算法分析与比较 (6) 1：谱减法 (6) (1) 谱减法算法的理论分析 (6) (2) 通过语音增强技术改善语音质量的过程 (7) (3) 谱减法的流程图 (7) (4)谱减法原理图 (8) 2：自适应噪声抵消法 (9) (1) 自适应滤波器原理 (9) (2) 结构框图 (10) 3：短时对数谱的MMSE语音增强算法 (11) 六：实验环境（matlab）简介 (13) 七：算法的实现及仿真结果 (14) 八：附件 (15) 九：参考文献 (17) 一：语音增强技术概述 在通信过程中语音受到来自周围环境、传输媒介引入的噪声，使接收到的语音信号并非纯净的原始语音信号，而是受噪声污染的带噪语音信号。这里的“噪音”定义为所需语音信号以外的所有干扰信号。 干扰信号可以是窄带的或宽带的、白噪声的或有色噪声的、声学的或电学的、加性的或乘性的，甚至可以是其它无关的语音。由噪声导致的语音质量的下降会使许多语音处理系统的性能急剧恶化。采用语音增强技术进行预处理，可有效地改善系统性能。 二：语音增强的目的 对收听人而言主要是改善语音质量，提高语音可懂度，减少疲劳感；对语音处理系统（识别器、声码器、手机）而言是提高系统的识别率和抗干扰能力。 三：语音信号的相关特性
二 抗剪强度折减系数法的理论 2．1抗剪强度折减系数法的概念 抗剪强度折减系数(SSRF ：Shear Strength Reduction Factor)定义为：在外荷载保持不变的情况下，边坡内土体所发挥的最大抗剪强度与外荷载在边坡内所产生的实际剪应力之比。这里定义的抗剪强度折减系数，与极限平衡分析中所定义的土坡稳定安全系数在本质上是一致的。 2．2抗剪强度折减系数法的具体内容 所谓抗剪强度折减技术就是将土体的抗剪强度指标C 和φ，用一个折减系数s F ，如式 (1)和(2) 所示的形式进行折减，然后用折减后的虚拟抗剪强度指标F C 和F φ，取代原来的抗剪强度指标C 和φ，如式（3）所示。 s F F C C /= （式1） )/)((tan tan 1s F F φφ-= （式2） F F fF C φστtan += （式3） 式中：F C 是折减后土体虚拟的粘聚力；F φ是折减后土体虚拟的内摩擦角；fF τ是折减后的抗剪强度。 折减系数s F 的初始值取得足够小，以保证开始时是一个近乎弹性的问题。然后不断增加s F 的值，折减后的抗剪强度指标逐步减小，直到某一个折减抗剪强度下整个土坡发生失稳，那么在发生整体失稳之前的那个折减系数值，即土体的实际抗剪强度指标与发生虚拟破坏时折减强度指标的比值，就是这个土坡的稳定安全系数。 2．3抗剪强度折减系数法的优点 结合有限差分法的抗剪强度折减系数法较传统的方法具有如下优点： （1）能够对具有复杂地貌、地质的边坡进行计算； （2）考虑了土体的本构关系，以及变形对应力的影响； （3）能够模拟土坡的边坡过程及其滑移面形状（通常由剪应变增量或者位移增量确定滑移面的形状和位置）； （4）能够模拟土体与支护结构（超前支护、土钉、面层等）的共同作用；
一、边坡稳定性计算方法 在边坡稳定计算方法中，通常采用整体的极限平衡方法来进行分析。根据边坡不同破裂面形状而有不同的分析模式。边坡失稳的破裂面形状按土质和成因不同而不同，粗粒土或砂性土的破裂面多呈直线形；细粒土或粘性土的破裂面多为圆弧形；滑坡的滑动面为不规则的折线或圆弧状。这里将主要介绍边坡稳定性分析的基本原理以及在某些边界条件下边坡稳定的计算理论和方法。 （一）直线破裂面法 所谓直线破裂面是指边坡破坏时其破裂面近似平面，在断面近似直线。为了简 化计算这类边坡稳定性分析采用直线破裂面法。能形成直线破裂面的土类包括：均质砂 性土坡；透水的砂、砾、碎石土；主要由内摩擦角控制强度的填土。 图9 －1 为一砂性边坡示意图，坡高H ，坡角β，土的容重为γ，抗剪 度指标为 c 、φ。如果倾角α的平面AC 面为土坡破坏时的滑动面，则可分析该滑 动体的稳定性。 沿边坡长度方向截取一个单位长度作为平面问题分析。 图9-1 砂性边坡受力示意图 已知滑体ABC重W ，滑面的倾角为α，显然，滑面AC 上由滑体的重量W= γ（ΔABC）产生的下滑力T 和由土的抗剪强度产生的 抗滑力Tˊ分别为： T=W ·sina 和 则此时边坡的稳定程度或安全系数可用抗滑力与下滑力来表示，即 为了保证土坡的稳定性，安全系数 F s 值一般不小于 1.25 ，特殊情况下可允许减小到 1.15 。对于C=0 的砂性土坡或是指边坡，其安全系 数表达式则变为 从上式可以看出，当α=β时，F s 值最小，说明边坡表面一层土最容易滑动，这时
当F s =1 时，β=φ，表明边坡处于极限平衡状态。此时β角称为休止角，也称安息角。 此外，山区顺层滑坡或坡积层沿着基岩面滑动现象一般也属于平面滑动类型。这类滑坡滑动面的深度与长度之比往往很小。当深长比小 于0.1 时，可以把它当作一个无限边坡进行分析。 图9-2 表示一无限边坡示意图，滑动面位置在坡面下H深度处。取一单位长度的滑动土条进 行分析，作用在滑动面上的剪应力为, 在极限平衡状态时，破坏面上的剪应 力等于土的抗剪强度，即 得 式中N s = c/ γH称为稳定系数。通过稳定因数可以确定α和φ关系。当c=0 时，即无 粘性土。α=φ，与前述分析相同。 二圆弧条法 根据大量的观测表明，粘性土自然山坡、人工填筑或开挖的边坡在破坏时，破裂面的形状多呈近似的圆弧状。粘性土的抗剪强度包括摩擦强 度和粘聚强度两个组成部分。由于粘聚力的存在，粘性土边坡不会像无粘性土坡一样沿坡面表面滑动。根据土体极限平衡理论，可以导出均质粘 这坡的滑动面为对数螺线曲面，形状近似于圆柱面。因此，在工程设计中常假定滑动面为圆弧面。建立在这一假定上稳定分析方法称为圆弧滑动 法和圆弧条分法。 1. 圆弧滑动法 1915 年瑞典彼得森（K.E.Petterson ）用圆弧滑动法分析边坡的稳定性，以后该法在各国得到广泛应用，称为瑞典圆弧法。 图9 － 3 表示一均质的粘性土坡。AC 为可能的滑动面，O 为圆心，R 为半径。 假定边坡破坏时，滑体ABC 在自重W 作用下，沿AC 绕O 点整体转动。滑动面AC 上的力系有：促使边坡滑动的滑动力矩M s =W ·d ；抵抗边坡滑动的抗滑力矩，它应该 包括由粘聚力产生的抗滑力矩M r =c ·AC ·R ，此外还应有由摩擦力所产生的抗滑力矩， 这里假定φ＝0 。边坡沿AC 的安全系数F s 用作用在AC 面上的抗滑力矩和下滑力 矩之比表示，因此有 这就是整体圆弧滑动计算边坡稳定的公式，它只适用于φ＝0 的情况。 图9-3 边坡整体滑动 2. 瑞典条分法
基于有限元强度折减法的边坡稳定性分析 报告 学院：土木工程与力学学院 专业：结构工程 姓名： 学号： 2016年7月
有限元强度折减法研究进展 摘要：在边坡稳定性分析中，相比于传统的极限平衡法、极限分析法等，有限元强度折减法具有明显的优势。这主要体现在其无须事先假定滑动面的形状和位置，只需通过不断降低边坡岩土体的强度参数，进而使边坡岩土体因抗剪强度不能抵抗剪切应力而发生破坏，并最终得到边坡的最危险滑动面及相应的安全系数。有限元强度折减法兼有数值计算方法和传统极限平衡方法的优点。本文介绍了有限元强度折减法的原理与主要研究现状，并对其中的一些重点问题进行了研究与总结。 关键词：强度折减法;有限元;边坡稳定 1 有限元强度折减法基本原理 所谓强度折减，就是在理想弹塑性有限元计算中将边坡岩土体抗剪切强度参数逐渐降低直到其达到破坏状态为止，程序可以自动根据弹塑性计算结果得到破坏滑动面(塑性应变和位移突变的地带)，同时得到边坡的强度储备安全系数ω, 于是有: ==。 '/,tan'tan/ c cω??ω 一般地，强度折减弹塑性有限元数值分析方法考察边坡稳定性的步骤是：首先对于某一给定的强度折减系数，通过逐级加载的弹塑性有限元数值计算确定边坡内的应力场、应变场或位移场，并且对应力、应变或位移的某些分布特征以及有限元计算过程中的某些数学特征进行分析，不断增大折减系数，直至根据对这些特征的分析结果表明边坡己经发生失稳破坏，将此时的折减系数定义为边坡的稳定安全系数。尽管强度折减有限元法在边坡稳定性分析中得到重视与发展，但其计算中需要采用一定的边坡失稳评判标准来确定边坡失稳的临界状态，但是，各种判据的选用至今并没有取得统一。 2 主要研究现状 强度折减概念由Zienkiewicz最早提出并用于边坡的稳定性分析，受限于当时数值计算和计算机水平而未能得到大的发展，直到近十几年来，随着数值计算和计算机技术的迅猛发展，强度折减法也得到了极大的发展，国内外许多学者在这方面做了大量的工作。 Ugai假定土体为理想的弹塑性材料，采用有限元强度折减法较为系统地分别对直立边坡、倾斜边坡、非均质边坡以及存在孔隙水压力的复杂边坡的稳定性进行了分析研究，并指出弹塑性强度折减有限元法具有较强的适应性和可行性。Matsui和San将强度折减技术与采用Duncan-Chang双曲线模型的非线性有限元法相结合，以剪应变作为边坡破坏评判指标，研究了人工填筑边坡和开挖边坡的稳定性，指出填筑边坡应采用总剪应变，而开挖边坡应采用局部剪应变增量作为失稳破坏标准，并将分析结果与极限平衡法进行了对比。Ugai和Leshchinsky 将强度折减技术引入弹塑性有限元法中进行边坡的三维稳定性分析，并与极限平衡法的计算结果进行了较全面的比较研究，指出尽管二者的理论基础、实现手段
边坡稳定性计算 一、编制依据 为保证挖方施工安全，施工现场做到“安全、文明”，满足施工进度要求，以下列法律、法规、标准、规范、规程、相关文件为强制性前提，进行边坡稳定性计算。 1、现有施工图设计； 2、《公路桥涵施工技术规范》（JTJ041-2000）； 3、《路桥施工计算手册》（人民交通出版社）； 4、《土力学与地基基础》； 二、工程概况及地质情况 岢岚至临县高速公路是《山西省高速公路网规划》“3纵11横11环”中西纵高速公路的重要组成部分，也是山西省西部把第四横（保德-五台长城岭）和第五横（平定杨树庄—佳县）高速公路窜连起来的重要路段。 项目区路线走廊带地形起伏极大，总体地势为东北高西南低，地貌主体为隆起的基岩中山与黄土梁峁，部分区域为海拔较低的河流沟谷及冲沟，。受构造活动和水流侵蚀作用的影响，本区地形切割剧烈，河谷发育，沟壑纵横，依据地貌成因类型及其显示特征，将本区划分为黄土丘陵区、侵蚀堆积河川宽谷区、山岭区、黄土覆盖中低山区四个地貌单元，岩性主要为第四系冲、坡积及风积粉土及粉质粘土等。 三、计算 本项目地形复杂，涵洞、桩基及路基施工作业面比较多。根据挖方路段在全线的分布情，选择有代表性路段进行分析计算。由于项目地质挖方为风积粉土及粉质粘土，是典型的黄土地貌。根据施工图纸给出的计算参数，对于黄土挖方路段，拟定边坡参数γ=19g/cm3，C=40 Kpa，φ=29°，采用瑞典条分法进行计算，稳定安全系数达到1.2以上。 3.1 瑞典条分法原理 如图所示边坡，瑞典条分法假定可能滑动面是一圆弧AD，不考虑条块两侧的作用力，即假设Ei和Xi的合力等于Ei+1和Xi+1的合力，同时它们的作用线
附件四：边坡稳定性计算书 1、汽机房区域边坡稳定性计算书（适用于基坑基底标高为-7.00m~-9.00m）H=8.5m 天然放坡支护 ---------------------------------------------------------------------- [ 基本信息 ] ---------------------------------------------------------------------- ---------------------------------------------------------------------- [ 放坡信息 ] ---------------------------------------------------------------------- ---------------------------------------------------------------------- [ 超载信息 ] ----------------------------------------------------------------------
---------------------------------------------------------------------- [ 土层信息 ] ---------------------------------------------------------------------- [ 土层参数 ] ---------------------------------------------------------------------- ---------------------------------------------------------------------- [ 整体稳定验算 ] ---------------------------------------------------------------------- 天然放坡计算条件: 计算方法：瑞典条分法 应力状态：总应力法 基坑底面以下的截止计算深度: 0.00m 基坑底面以下滑裂面搜索步长: 5.00m 条分法中的土条宽度: 1.00m 天然放坡计算结果:
基本原理： 强度折减法中边坡稳定的安全系数定义为：使边坡刚好达到临界破坏状态时，对岩、土体的抗剪强度进行折减的程度，即定义安全系数为岩土体的实际抗剪强度与临界破坏时的折减后剪切强度的比值。强度折减法的要点是公式1、2来调整岩土体的强度指标C 和φ(式中，F C 为折减后的粘结力，F φ为折减后的摩擦角，trial F 为折减系数），然后对边坡稳定性进行数值分析，不断地增加折减系数。反复计算，直至其达到临界破坏，此时得到的折减系数即为安全系数S F 。公式如下： trial F F C C /= （1） t a n =F φ-1)/)((tan trial F φ（2) 实现过程： 目前尚无统一的边坡失稳判据，现行的边坡失稳判据主要有以下几种： 1 以数值计算的收敛性作为失稳判据 2 以特征部位位移的突变性作为失稳判据 3 以塑性区的贯通性作为失稳判据 在FLAC3D 中求解安全系数时，单次安全系数的计算过程主要采用的是第一种失稳判据。假设数值计算模型所有非空区域都采用摩尔-库伦本构模型，便可使用命令Solve fos 来求解安全系数：首先，通过给粘结力设定一个大值来改变内部应力，以找到体系达到力平衡的典型时步r N ;接着，对于给定的安全系数s F ，执行r N 时步，如果体系不平衡力与典型内力比率R 小于10-3，则认为体系达到力平衡。如
果不平衡力比率R大于10-3，再执行r N时步，直至R小于10-3退出当前计算，开始新一轮折减计算过程。除上述以力不平衡比率小于10-3作为终止条件外，FLAC3D还采用： 1 前后典型时步运算结束时的不平衡力比率R差值小于10% 2 强度折减后的计算过程已运行了6个典型时步r N作为计算终止条件 计算过程中，只要满足上述三个标准中的任何一个，便退出当前计算。这样做的目的只要是为了控制整个强度折减法循环计算过程中的求解时间。 可以从这几个方面判断：边坡沿滑动面产生滑动、软弱面处产生的沿X方向的位移是否最大、剪切应变增量云图、安全系数、剪切应变增量云图、变形矢量图及速度矢量图、水平位移、竖直位移、垂直应力、最大不平衡力、在坡顶边缘和坡脚处设置监测点（水平应力竖直应力位移）。 FLAC是快速拉格朗日差分分析(Fast grangian Analysis of Continua)的简写。它是以岩石力学理论为基础，以介质物理力学参数和地质构造特性为计算依据，建立在客观反映原型和动态演化过程仿真力学效应基础上的一种新型数值方法。虽然其基本原理类同于离散元法，但却可与有限元一样适用于多种材料模式与边界条件非规则区域的连续问题求解。而且计算中利用的“混合离散化”技术可针对不同介质特
燕山大学 医学软件课程设计说明书 题目：基于MATLAB巴特沃斯滤波器的音频去噪的GUI设计 学院（系）：电气工程学院 年级专业： 13级生物医学工程 2 班 学号： 130103040041 学生姓名：魏鑫 指导教师：许全盛
目录 一、设计目的意义 (1) 1.1绪论 (1) 1.2设计目的 (1) 1.3意义 (1) 二、设计内容 (2) 2.1 设计原理 (2) 2.2 设计内容 (2) 三、设计过程及结果分析 (3) 3.1 设计步骤 (3) 3.2 MATLAB程序及结果 (3) 3.3 结果分析 (8) 四、总结 (9) 五、参考文献 (10)
一、设计目的意义 1.1 绪论 语音是语言的声学表现，是人类交流信息最自然、最有效、最方便的手段。随着社会文化的进步和科学技术的发展，人类开始进入了信息化时代，用现代手段研究语音处理技术，使人们能更加有效地产生、传输、存储、和获取语音信息，这对于促进社会的发展具有十分重要的意义，因此，语音信号处理正越来越受到人们的关注和广泛的研究。 1.2 设计目的 (1)掌握数字信号处理的基本概念，基本理论和基本方法。 (2)熟悉离散信号和系统的时域特性。 (3)掌握序列快速傅里叶变换方法。 (4)学会MATLAB的使用，掌握MATLAB的程序设计方法。 (5)掌握利用MATLAB对语音信号进行频谱分析。 (6)掌握滤波器的网络结构。 (7)掌握MATLAB设计IIR、FIR数字滤波器的方法和对信号进行滤波的方法。 1.3 意义 语音信号处理是一门比较实用的电子工程的专业课程，语音是人类获取信息的重要来源和利用信息的重要手段。通过语言相互传递信息是人类最重要的基本功能之一。语言是人类特有的功能，它是创造和记载几千年人类文明史的根本手段，没有语言就没有今天的人类文明。语音是语言的声学表现，是相互传递信息的最重要的手段，是人类最重要、最有效、最常用和最方便的交换信息的形式。 语音信号处理是研究用数字信号处理技术对语音信号进行处理的一门学科，它是一门新兴的学科，同时又是综合性的多学科领域和涉及面很广的交叉学科。
题目基于matlab谱减法音频降噪处理 班级 学号 姓名 指导 时间 景德镇陶瓷学院
数字信号处理课程设计任务书
目录 1、设计要求. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .1 2、设计原理. . . . . . . ……………………………………………….. . . . . . . . . . . .2 3、源程序清单. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .7 4、设计结果和仿真波形. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .11 5、参考文献. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15. 6、设计心得体会. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
1、设计要求 语言是人类最重要、直接、有效和便捷的交换信息的方式。随着近些年科学技术的飞速发展，人们也不满足于和计算机的信息交换方式，希望能够甩掉键盘和鼠标而实现用语言来对计算机进行控制。因此，语音信号处理技术便应运而生。语音信号处理是一门新兴的学科，同时也是综合多种学科和涉及面非常广泛的交叉学科。现在在一些职能系统中嵌入有语音处理系统，但它们只能在安静的环境中才能使用。然而，在语音信息的采集过程中难免会有各种噪声的干扰。噪声不仅降低了语音的可懂度和语音质量，还严重的影响语音处理的准确性，甚至使系统不能正常工作。本文将就对语音增强技术的原理和方法进行讨论，重点介绍语音增强的一种方法——谱减法及其改进算法。该方法能够有效消除平稳的加性噪声，其改进算法能够有效消除普通方法产生的“音乐噪声”，在很大程度上提高语音信号的信噪比。 目前，语言识别技术已经取得了重大进展，并开始进入实用阶段。但语音识别系统必须在相对比较安静的环境下运行，然而，在语言信息的采集中难免会有各种噪声的干扰，在较强的噪声背景下，语音识别系统的准确性会受到较大影响，甚至没法正常工作。所以在语音识别系统对语音信息处理前，应该对语音信息进行预处理，即背景噪声消除。语音背景噪声消除技术的出现使得语音识别技术更加稳定和精确，也使得语音信息的可懂度大大提高，使人们能够从较复杂的语音信息中提取到更多的有用信息。
... . . 土坡稳定性计算书 本计算书参照《建筑施工计算手册》江正荣编著中国建筑工业、《实用土木工程手册》第三版文渊编著人民教同、《地基与基础》第三版中国建筑工业、《土力学》等相关文献进行编制。 计算土坡稳定性采用圆弧条分法进行分析计算，由于该计算过程是大量的重复计算，故本计算书只列出相应的计算公式和计算结果，省略了重复计算过程。 本计算书采用瑞典条分法进行分析计算，假定滑动面为圆柱面及滑动土体为不变形刚体，还假定不考虑土条两侧上的作用力。 一、参数信息: 条分方法：瑞典条分法； 考虑地下水位影响； 基坑外侧水位到坑顶的距离(m)：1.56； 基坑侧水位到坑顶的距离(m)：14.000； 放坡参数： 序号放坡高度(m) 放坡宽度(m) 平台宽度(m) 条分块数 0 3.50 3.50 2.00 0.00 1 4.50 4.50 3.00 0.00 2 6.20 6.20 3.00 0.00 荷载参数： 土层参数：
序号土名称 土厚 度 （m） 坑壁土的重 度γ(kN/m3) 坑壁土的摩 擦角φ(°) 粘聚力 （kPa） 饱容重 (kN/m3) 1 粉质粘土15 20.5 10 10 20.5 二、计算原理: 根据土坡极限平衡稳定进行计算。自然界匀质土坡失去稳定，滑动面呈曲面，通常滑动面接近圆弧，可将滑裂面近似成圆弧计算。将土坡的土体沿竖直方向分成若干个土条，从土条中任意取出第i条，不考虑其侧面上的作用力时，该土条上存在着： 1、土条自重， 2、作用于土条弧面上的法向反力， 3、作用于土条圆弧面上的切向阻力。 将抗剪强度引起的极限抗滑力矩和滑动力矩的比值作为安全系数，考虑安全储备的大小，按照《规》要求，安全系数要满足>=1.3的要求。 将抗剪强度引起的极限抗滑力矩和滑动力矩的比值作为安全系数，考虑安全储备的大小，按照《规》要求，安全系数要满足>=1.3的要求。 三、计算公式:
处理宽带噪声的最通用技术是谱减法。它利用语音信号的短时平稳特性，从带噪语音的短时谱值中减去噪声的短时谱，从而得到纯净语音的频谱，达到语音增强得目的。谱减法包括幅度谱减法和功率谱减法：幅度谱减法就是在频域中从带噪语音的幅度谱上减去噪声的幅度谱作为语音信号的幅度谱；功率谱减法则是从带噪语音的功率谱中减去噪声的功率谱，得到纯净语音的功率谱估计，通过开方运算得到幅度谱。由于人耳对语音频谱分量的相位感知不敏感，因此这些算法都是在幅度上进行的修正，相位部分则保持不变，仍然使用带噪语音的相位。 谱减法通过从带噪语音的短时谱估值中减去噪声的短时谱来达到语音增强得目的，算法简单且容易实现。但在减去噪声谱后，还会有些较大功率谱分量的剩余部分，在频谱上呈现出随机出现的尖峰，在听觉上形成残留噪声。这种噪声具有一定的节奏性起伏感，被称之为“音乐噪声”。后来，Ephraim 等人又对谱减法进行了大量改进，部分解决了“音乐噪声”问题，但在带噪语言信噪比较低时其残余噪声还是很大，尤其是当信噪比小于5dB 的时候。因此，如何最大限度地消除谱减法中的“音乐噪声”，仍将是人们今后研究的重要课题。 以下为MATLAB源码的实现 %基本谱减法 clear; %[xx,fs]=wavread('E:\mywhisper\shu.wav'); %[xx,fs]=wavread('E:\speech\x\w1xun_01.wav'); %[xx,fs]=wavread('E:\speech\耳语音切割\b\w1ba_5'); [xx,fs]=wavread('MIC0.wav');% 读取音频文件yuan.wav，并返回采样数据给变量xx及采样率Fs [team,row]=size(xx);%将数组xx的行数赋给team,列数赋给row if row==2 x=(xx(:,1)+xx(:,2))/2; yy=x; %如果语音信号xx为2列，即信号为双声道，则将其转换成单声道信号，即取两列的平均值赋给x，并将x的值赋给yy else x=xx; yy=x; %若语音信号xx为单声道，则将xx的值赋给x，并将x的值赋给yy end x=x-mean(x)+0.1*rand(length(x),1); N=length(x);%将语音信号长度赋给变量N n=220;%对语音信号进行分帧，帧长为220 n1=160;%帧移为160 frame=floor((N-n)/(n-n1));%将分帧数赋给变量frame %frame=floor(N/n); for i=1:frame y1=x((i-1)*(n-n1)+1:(i-1)*(n-n1)+n).*hamming(n); %对每段分帧进行加窗处理
平面、折线滑动法边坡稳定性计算书计算依据： 1、《建筑基坑支护技术规程》JGJ120-2012 2、《建筑边坡工程技术规范》GB50330-2002 3、《建筑施工计算手册》江正荣编著 一、基本参数 边坡稳定计算方式折线滑动法边坡工程安全等级三级边坡边坡土体类型填土土的重度γ(KN/m3) 20 土的内摩擦角φ(°)15 土的粘聚力c(kPa) 12 边坡高度H(m) 11.862 边坡斜面倾角α(°)40 坡顶均布荷载q(kPa) 0.2 二、边坡稳定性计算 计算简图 滑动面参数 滑动面序号滑动面倾角θi(°)滑动面对应竖向土条宽度bi(m) 1 35 5.67 2 35 5.6 3 35 5.67 土条面积计算：
R1=(G1+qb1)cosθ1×tanφ+c×l1=(156.213+0.2×2.803)×cos(35°)×tg(15°)+12×6.922=117.474 kN/m T1=(G1+ qb1)sinθ1 =(156.213+0.2×2.803)×sin(35°)=89.922 kN/m R2=(G2+qb2)cosθ2×tanφ+c×l2=(131.759+0.2×0)×cos(35°)×tg(15°)+12×6.836=110.952 kN/m T2=(G2+ qb2)sinθ2 =(131.759+0.2×0)×sin(35°)=75.574 kN/m R3=(G3+qb3)cosθ3×tanφ+c×l3=(44.652+0.2×0)×cos(35°)×tg(15°)+12×6.922=92.865 kN/m T3=(G3+ qb3)sinθ3 =(44.652+0.2×0)×sin(35°)=25.611 kN/m K s=(∑R iψiψi+1...ψn-1+R n)/(∑T iψiψi+1...ψn-1+T n)，(i=1,2,3,...,n-1) 第i块计算条块剩余下滑推力向第i+1计算条块的传递系数为： ψi=cos(θi-θi+1)-sin(θi-θi+1)×tanφi K s=(∑R iψiψi+1...ψn-1+R n)/(∑T iψiψi+1...ψn-1+T n)=(117.474×1×1+110.952×1+92.865)/(89.922×1×1+75.574×1+25.611)=1.681≥1.25 满足要求！
匀质粘性土体边坡稳定性计算中 极限平衡法与强度折减法数值分析比较 1边坡的种类和滑面类型 边坡按形成的原因大致可以分为三类：自然边坡、人工开挖边坡和人工填筑边坡；边坡按土层的种类也大致可以分为三种：匀质黏土边坡（人工填筑边坡一般为匀质边坡）、非匀质粘性土边坡（自然的黏砂性土边坡，或者人工开挖黏土边坡）和存在潜在软弱面或层面强度差异较大的边坡（如存在既有滑面的滑坡和上层为黏土层下为岩层的边坡）。使边坡失稳的外因大致有：外荷载（地震、列车荷载、房屋和填土等）、重力和水的渗流、岩土的膨胀力等。边坡失稳时滑动及滑面类型主要取决于边坡的外部荷载和边坡土层的类别，匀质粘土边坡的失稳滑面主要为圆弧型，非匀质边坡的失稳滑面主要是曲面型（复合型、对数螺旋型等），存在潜在软弱面或层面强度差异较大的边坡一般沿既有软弱面或者沿强度差异较大的两层面层间滑动（一般多为折线型，也有直线型）。本将对匀质边坡进行有限差分强度折减法数值分析法和经典极限平衡法进行稳定性分析计算，并就安全系数的计算及粘性土边坡潜在滑动面的确定进行比较分析。 2匀质边坡的极限平衡法计算和数值模拟 2.1匀质粘土边坡的极限平衡法概述 粘性土边坡中,危险滑动面在土体的内部,常与圆弧面相似。经典的基于圆弧滑动面的边坡稳定分析方法称为圆弧滑动法,属极限平衡法[1]，本文对几种常用的经典分析方法进行简单概述如下。 （1)整体圆弧滑动法。又称瑞典圆弧法,用于分析均质黏性土边坡的稳定性,即只能分析内摩擦角υ=0时的边坡稳定问题。边坡稳定安全系数为抗滑力矩与滑动力矩之比, 一般公式为:11 (sin ) n i i r n s i i i cl M K M W α === = ∑∑ （2)简单圆弧条分法。又称瑞典条分法,仍假定滑动面为一个圆弧面。为分析摩擦角υ>0时粘性土边坡的稳定性,将土坡分成若干个条块,但只考虑作用在条块上的重力、滑弧面上的法向力和切向抗滑力,忽略条块侧面法向力和切向力的作用。一般公式为： 1 1 (cos tan ) (sin ) n i i i i r n s i i i W cl M K M W α φα ==+= = ∑∑ （3)简化毕肖普(Bishop)法。简化毕肖普法仍假设滑动面为圆弧,将滑动土体分为若干个条块。在进行第ｉ个条块的受力分析时,考虑条块侧面法向力的作用,但忽略切向力的作用。一般公式为：
有限元强度折减法和极限平衡法适用性研究 曾红丽 戚明军 陕西省建筑设计研究院有限责任公司 710003 西安 摘要：简述了有限元强度折减法原理和两种极限平衡理论进行边坡稳定性分析的方法，分别利用ANSYS 和SLOPE/W 程序对一均质土坡进行了稳定性分析。建立了同样尺寸的两种计算模型，得出了分别采用有限元强度折减法、Bishop 法和Janbu 法在未考虑地下水的情况下边坡安全系数。将计算结果进行比较分析，表明基于极限平衡理论的分析方法分析得到的结果偏小，而有限元强度折减法由于考虑了土体内部应力应变关系所得结果更切合工程实际。 关键词：有限元强度折减法；极限平衡法；适用性；稳定性 1引言 边坡工程是公路、铁路、水利水电以及矿山工程一个不可或缺的组成部分，因此边坡稳定性问题的研究就成为岩土工程界研究的热点问题。国内外学者，已经取得了大量研究成果。目前，研究边坡稳定性的传统方法主要有：极限平衡法，极限分析法，滑移线场法等。随着计算机技术的发展，有限元强度折减法 在边坡稳定性分析中已经得到了较好的应用[1～5] 。本文通过算例分析，比较分析有限元强度折减法和基于极限平衡理论的方法在边坡工程中适用性的差别。 2有限元强度折减法原理[6] 所谓强度折减，就是在理想弹塑性有限元计算中将边坡岩土体抗剪强度参数逐渐降低直到其达到破坏状态为止，程序可以自动根据弹塑性计算结果得到破坏滑动面，同时得到边坡的强度储备安全系数K ，也称为强度折减系数。 通常，边坡的稳定性安全系数定义为沿滑动面的抗剪强度与滑动面的实际剪力的比值，公式表示为 ??+=dA dA c K τ?σ)tan ( （1） 将式（1）两边同时除以K ，得 ????+=+=dA dA c dA dA K K c τ?στ?σ )'tan '()tan ( 1 （2） 式中，K c c ='，)arctan(tan 'K ??=，c 为粘聚力，?为内摩擦角，σ为滑动面上的法向应力，τ为滑动面上的抗剪强度。 通过逐步调整系数K ，得到不同的'c ，'?，将'c ，'?代入有限元程序，反复分析边坡，直到坡体达到临界状态，坡体达到临界状态时的K 值即作为边坡稳定性安全系数。 3极限平衡法 利用极限平衡法对边坡的稳定性进行了分析，该方法的基本特点是，只考虑静力平衡条件和土的Mohr-Coulomb 破坏准则，也就是说，通过分析土体的破坏那一刻的平衡来求得问题的解[7]。极限平衡理论的主要思想是将滑动土体进行条分，由极限状态下土条所受力和力矩的平衡来分析边坡稳定性。它是目前应用最多的一种分析方法。文中所用极限平衡理论的方法有以下几种[8]。 （1）Bishop 法：该方法考虑了土条间的作用力，这是对传统的瑞典条分法的重要改进。该方法忽略各土条之间的切向条间力，认为条间力的合力是水平的，同时假设破坏面是圆弧面，且定义边坡安全系数为沿整个滑动面上的抗剪强度与实际产生的剪应力的比值。该方法的计算结果比较接近实际，常用于土质、软岩质及碎岩边坡的稳定性分析。 （2）Janbu 法：Janbu 提出了同时满足力和力矩平衡的“通用条分法”。这一方法区别于其他方法的一个重要方面，就是通过假定土条侧向力的作用点而不是作用方向来求解安全系数的。采用非圆弧面，按条块滑动平衡确定条间力，按推力线确定法向力的作用点，简化计算条间切向力为零，然后再对稳定系数进行修正。 4算例分析
土坡稳定性计算书 本计算书参照《建筑施工计算手册》江正荣编著中国建筑工业出版社、《实用土木工程手册》第三版杨文渊编著人民教同出版社、《地基与基础》第三版中国建筑工业出版社、《土力学》等相关文献进行编制。 计算土坡稳定性采用圆弧条分法进行分析计算，由于该计算过程是大量的重复计算，故本计算书只列出相应的计算公式和计算结果，省略了重复计算过程。 本计算书采用瑞典条分法进行分析计算，假定滑动面为圆柱面及滑动土体为不变形刚体，还假定不考虑土条两侧上的作用力。 一、参数信息: 条分方法：瑞典条分法； 考虑地下水位影响； 基坑外侧水位到坑顶的距离(m)：1.56； 基坑内侧水位到坑顶的距离(m)：14.000； 放坡参数： 序号放坡高度(m) 放坡宽度(m) 平台宽度(m) 条分块数 0 3.50 3.50 2.00 0.00 1 4.50 4.50 3.00 0.00 2 6.20 6.20 3.00 0.00 荷载参数： 土层参数：
序号土名称 土厚 度（m） 坑壁土的重 度γ(kN/m3) 坑壁土的内 摩擦角φ(°) 粘聚力 （kPa） 饱容重 (kN/m3) 1 粉质粘土15 20.5 10 10 20.5 二、计算原理: 根据土坡极限平衡稳定进行计算。自然界匀质土坡失去稳定，滑动面呈曲面，通常滑动面接近圆弧，可将滑裂面近似成圆弧计算。将土坡的土体沿竖直方向分成若干个土条，从土条中任意取出第i条，不考虑其侧面上的作用力时，该土条上存在着： 1、土条自重， 2、作用于土条弧面上的法向反力， 3、作用于土条圆弧面上的切向阻力。 将抗剪强度引起的极限抗滑力矩和滑动力矩的比值作为安全系数，考虑安全储备的大小，按照《规范》要求，安全系数要满足>=1.3的要求。 将抗剪强度引起的极限抗滑力矩和滑动力矩的比值作为安全系数，考虑安全储备的大小，按照《规范》要求，安全系数要满足>=1.3的要求。 三、计算公式:
基于Matlab的谱减法语音增强的研究 【摘要】在实际应用中，待分析的语音信号一般是被噪声污染的语音，很大程度上影响了语音处理系统的性能。为此，需要我们对带噪语音进行语音增强处理，抑制噪声，恢复纯净语音。本文就谱减法对带噪语音进行处理，并利用matlab 行仿真，验证了谱减法的语音增强的效果。 【关键词】Matlab;语音增强;谱减法 1.引言 研究语音增强技术在实际中有重要价值。目前，语音增强己在很多方面得到广泛的应用，例如语音处理系统、通信、多媒体技术、数字化家电等领域。语音增强的一个主要目标，就是从带噪语音信号中提取尽可能纯净的原始语音。根据与输入语音信号的关系，噪声可分为加性噪声和非加性噪声两类。对某些非加性噪声而言，可以通过一定的变换转换成加性噪声。非加性噪声主要是残响和传送网络的电路噪声等。加性噪声通常分为宽带噪声、冲激噪声、语音干扰噪声、周期噪声等[1]。 2.谱减法基本原理 谱减法的基本思想是在假定加性噪声与短时平稳的语言信号相互独立的情况下，从带噪语音的功率谱中减去噪声功率谱，从而得到较为纯净的语音频谱[2]。如果设s（t）为纯净语音，n（t）为噪声信号，y（t）为带噪语音信号[3]，则有： y（t）=s（t）+n（t）（2-1） 用Y（）、S（）、N（）分别表示y（t）、s（t）、n（t）的傅里叶变换，则可得 Y（）=S（）+N（）（2-2） 由于假定语音信号与加性噪声是相互独立的，则可以得到： （2-3） 如果用、、分别表示y（t）、s（t）和n（t）的功率谱，则有： =+ （2-4） 而由于平稳噪声的功率谱在发声前和发声期间可以认为基本没有变化，这样可以通过发声前的所谓“寂静段”来估计噪声的功率谱，从而有：
基于改进型谱减法的语音增强 摘要 本文主要研究改进型谱减算法在语音增强中的应用，目的是增强语音质量，减少语音失真和提高其可度懂。我们首先介绍了语音增强的研究意义，然后介绍了语音信号的相关理论，进而阐述了语音增强的基本谱减法的原理，并在此基础上提出了一种改进型谱减算法。该算法通过语音激活检测(端点检测法)来确定“寂静段”（纯噪声段），从而对噪声功率谱进行重新估计。为了减小基音检测算法可能产生的检测误差，采用了组合递归平滑法来减小噪声谱估计的误差。 整个仿真实验中，我们对引入的加性噪声进行处理，其噪声谱估计的性能可在本文中的MATLAB仿真实验中体现。 仿真结果表明，该算法在去除背景噪声的同时，保证了较小的语音失真、提高了信噪比，达到了较好的测听效果。 最后，基于噪声与语音具有一定的相关性的实际情况，我们提出了算法的进一步改进设想，并对此思想做出了数学推导，得到了算法进一步改进的方向及可行性。关键词：语音增强；谱减法；噪声估计；端点检测；组合递归平滑；仿真实验；改进算法
第1章背景介绍 1.1研究背景 人们在语音通信过程中不可避免地会受到来自周围环境、传输媒介引入的噪声、通信设备内部电噪声乃至其他讲话者的干扰。这些干扰最终将使接收者接收到的语音已非纯净的原始语音信号，而是受噪声污染的带噪语音信号。例如，安装在汽车、飞机或舰船上的电话，街道、机场的公用电话，常受到很强背景噪声的干扰，严重影响通话质量。又如，室内会议电话的交混回响随同语音广播到每个会议地点，影响收听效果。再如深海潜水员在氦-氧面罩内讲话引起的失真，语言障碍残疾人的语音失真，有历史价值的旧唱片、旧录音带的噪声和失真等，都是带噪语音信号的例子。 环境噪声污染使许多语音处理系统的性能急剧恶化。例如，语音识别己取得重大进展，正在步入实用阶段。但目前的识别系统大都是在安静环境中工作的，在噪声环境中尤其是强噪声环境，语音识别系统的识别率将受到严重影响。低速率语音编码，特别是参数编码(如：声码器)，也遇到类似问题。由于语音生成模型是低速率参数编码的基础，当模型参数的提取受到混杂在语音中背景噪声严重干扰时，重建语音的质量将急剧恶化，甚至变得完全不可懂。在上述情况下，语音增强作为一种预处理手段，不失为解决噪声污染的一种有效途径。 在实际需求的推动下，早在上个世纪60年代语音增强这个研究课题就引起人们的注意，此后40多年人们一直锲而不舍地进行这方面的研究。随着数字信号处理理论的成熟，70年代曾形成一个理论研究高潮，取得了一些基础性成果，并使语音增强发展成为语音信号处理的一个重要分支。进入80年代后，DSP(数字信号处理)技术的发展和成熟为语音增强的实时实现提供了可能。 语音增强不但与语音信号数字处理理论有关，而且涉及到人的听觉感知和语音学范畴。再者，噪声的来源众多，随应用场合而异，它们的特性也各不相同。即使在实验室仿真条件下，也难以找到一种通用的语音增强算法能适用于各种噪声环境，所以必须针对不同噪声，采用不同的语音增强对策。目前，某些语音增强算法在实