点安全系数与强度折减法对比

基于matlab的MMSE

基于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) 一：语音增强技术概述 在通信过程中语音受到来自周围环境、传输媒介引入的噪声，使接收到的语音信号并非纯净的原始语音信号，而是受噪声污染的带噪语音信号。这里的“噪音”定义为所需语音信号以外的所有干扰信号。 干扰信号可以是窄带的或宽带的、白噪声的或有色噪声的、声学的或电学的、加性的或乘性的，甚至可以是其它无关的语音。由噪声导致的语音质量的下降会使许多语音处理系统的性能急剧恶化。采用语音增强技术进行预处理，可有效地改善系统性能。 二：语音增强的目的 对收听人而言主要是改善语音质量，提高语音可懂度，减少疲劳感；对语音处理系统（识别器、声码器、手机）而言是提高系统的识别率和抗干扰能力。 三：语音信号的相关特性

机械设计中的安全系数选择问题

工程中的材料强度、刚度、稳定性。 强度-构件在确定的外力作用下，不发生破坏或过量塑性变形的能力。 杆-拉杆与压杆。 工程中承受拉伸的杆件统称为拉杆，受压的杆件成为杆或柱，承受扭转的杆件称为轴，承受弯曲的杆件统称为梁。 在工程力学中，把一些杆轴交汇于一点的工程结构称为桁架结构，这种结构受力特征是内力只有轴力，没有弯矩和剪力。如：井架的主体桁架、建筑脚手架、三角形屋架梁等。 许用应力与安全系数 最近听到对于建井结构安全的一些言论，有的说安全凭经验即可，我原来怎样用的，现在怎样用是没有问题的；有的说，计算是什么结果，应该遵守。 用伟人毛泽东的哲学思想是“实践—理论—实践”， 我们正常工作中选用的钢丝绳安全系数、钢材安全系数许用应力和安全系数都是比较成熟的，是规范推荐值或强制值。 在非标准或特殊情况下，安全应由自己评估。许用应力与安全系数常常应由自己选取决策。强度—在确定的外力作用下，不发生破坏的能力。 刚度—在确定的外力作用下，变形或位移在工程允许的范围内。 稳定性—在可能的外力作用下不会发生突然转变的能力。例如：建筑施工脚手架，强度、刚度能满足，但由于局部结构不稳定，使整个脚手架倾覆或塌陷。 材料名称屈服点σs抗拉强度σb抗剪强度τ单位材料使用地点 Q235 235 375 MPa或N/mm^2 普通结构 45 355 600 轴类件 30CrMnTi 1470 60Si2CrVA 1678 1865 钢丝 安全系数S应该综合荷载确定的准程度、材料性能数据的可靠性、所有计算方法的合理性、加工装配精度以及所设计的零件的重要性来确定。各行各业都有一些经验的安全系数，目前均偏于保守。目前，流行的安全系数法是部分系数法，他将各个对安全系数有影响的因素分别用一个分系数如：S1、S2、……标示，这些系数的乘积即即为安全系数：S=S1?S2?S3。。。。在实际应用中，取大取小带有一定主观性，即一般取大值或中间值，考虑的因素越多，系数值越大。 名称 S 抗疲劳计算系数 1.5～3 抗变形计算系数 1.2～2 抗断裂计算系数 2～4 抗不稳定计算系数 3～5 工作重要性系数 1.0～1.3 计算误差系数 1.2～1.3 轧制工艺可靠性系数 1.05～1.1 锻造工艺可靠性系数 1.05～1.1 铸造工艺可靠性系数 1.15～1.2 使用磨损系数 1.15～1.25 锈蚀系数 1.15～1.2 钢丝绳结构系数 1.217 案例：凿井提升钩头的安全系数S怎样确定？

强度折减法的原理

二 抗剪强度折减系数法的理论 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）能够模拟土体与支护结构（超前支护、土钉、面层等）的共同作用；

基于有限元强度折减法的边坡稳定性分析报告

基于有限元强度折减法的边坡稳定性分析 报告 学院：土木工程与力学学院 专业：结构工程 姓名： 学号： 2016年7月

有限元强度折减法研究进展 摘要：在边坡稳定性分析中，相比于传统的极限平衡法、极限分析法等，有限元强度折减法具有明显的优势。这主要体现在其无须事先假定滑动面的形状和位置，只需通过不断降低边坡岩土体的强度参数，进而使边坡岩土体因抗剪强度不能抵抗剪切应力而发生破坏，并最终得到边坡的最危险滑动面及相应的安全系数。有限元强度折减法兼有数值计算方法和传统极限平衡方法的优点。本文介绍了有限元强度折减法的原理与主要研究现状，并对其中的一些重点问题进行了研究与总结。 关键词：强度折减法;有限元;边坡稳定 1 有限元强度折减法基本原理 所谓强度折减，就是在理想弹塑性有限元计算中将边坡岩土体抗剪切强度参数逐渐降低直到其达到破坏状态为止，程序可以自动根据弹塑性计算结果得到破坏滑动面(塑性应变和位移突变的地带)，同时得到边坡的强度储备安全系数ω, 于是有: ==。 '/,tan'tan/ c cω??ω 一般地，强度折减弹塑性有限元数值分析方法考察边坡稳定性的步骤是：首先对于某一给定的强度折减系数，通过逐级加载的弹塑性有限元数值计算确定边坡内的应力场、应变场或位移场，并且对应力、应变或位移的某些分布特征以及有限元计算过程中的某些数学特征进行分析，不断增大折减系数，直至根据对这些特征的分析结果表明边坡己经发生失稳破坏，将此时的折减系数定义为边坡的稳定安全系数。尽管强度折减有限元法在边坡稳定性分析中得到重视与发展，但其计算中需要采用一定的边坡失稳评判标准来确定边坡失稳的临界状态，但是，各种判据的选用至今并没有取得统一。 2 主要研究现状 强度折减概念由Zienkiewicz最早提出并用于边坡的稳定性分析，受限于当时数值计算和计算机水平而未能得到大的发展，直到近十几年来，随着数值计算和计算机技术的迅猛发展，强度折减法也得到了极大的发展，国内外许多学者在这方面做了大量的工作。 Ugai假定土体为理想的弹塑性材料，采用有限元强度折减法较为系统地分别对直立边坡、倾斜边坡、非均质边坡以及存在孔隙水压力的复杂边坡的稳定性进行了分析研究，并指出弹塑性强度折减有限元法具有较强的适应性和可行性。Matsui和San将强度折减技术与采用Duncan-Chang双曲线模型的非线性有限元法相结合，以剪应变作为边坡破坏评判指标，研究了人工填筑边坡和开挖边坡的稳定性，指出填筑边坡应采用总剪应变，而开挖边坡应采用局部剪应变增量作为失稳破坏标准，并将分析结果与极限平衡法进行了对比。Ugai和Leshchinsky 将强度折减技术引入弹塑性有限元法中进行边坡的三维稳定性分析，并与极限平衡法的计算结果进行了较全面的比较研究，指出尽管二者的理论基础、实现手段

第三章疲劳强度计算练习题

第三章机械零件的疲劳强度设计 一、选择题 3－1 45钢的持久疲劳极限σ-1＝270MPa,，设疲劳曲线方程的幂指数m=9，应力循环基数N0=5×106次，当实际应力循环次数N=104次时，有限寿命疲劳极限为____________MPa。 （1）539 （2）135 （3）175 （4）417 3－2 有一根阶梯轴，用45钢制造，截面变化处过渡圆角的疲劳缺口系数Kσ=1.58，表面状态系数β＝0.28，尺寸系数εσ＝0.68，则其疲劳强度综合影响系数KσD=____________。 （1）0.35 （2）0.88 （3）1.14 （4）2.83 3－3 形状、尺寸、结构和工作条件相同的零件，采用下列不同材料制造：a）HT200；b）35钢；c）40CrNi钢。其中设计零件的疲劳缺口系数最大和最小的分别是____________。 （1）a）和b）（2）c）和a）（3）b）和c） （4）b）和a）（5）a）和c）（6）c）和b） 3－4 零件的截面形状一定，如绝对尺寸（横截面尺寸）增大，疲劳强度将随之____________。 （1）增高（2）不变（3）降低 3－5 零件的形状、尺寸、结果相同时，磨削加工的零件与精车加工相比，其疲劳强度____________。 （1）较高（2）较低（3）相同 3－6 零件表面经淬火、渗氮、喷丸、滚子碾压等处理后，其疲劳强度____________。 （1）增高（2）降低（3）不变（4）增高或降低视处理方法而定 3－7 影响零件疲劳强度的综合影响系数KσD或KτD与____________等因素有关。 （1）零件的应力集中、加工方法、过载 （2）零件的应力循环特性、应力集中、加载状态 （3）零件的表面状态、绝对尺寸、应力集中 （4）零件的材料、热处理方法、绝对尺寸。 3－8 已知设计零件的疲劳缺口系数Kσ=1.3、尺寸系数εσ＝0.9、表面状态系数βσ＝0.8。则疲劳强度综合影响系数KσD为____________。 （1）0.87 （2）0.68 （3）1.16 （4）1.8 3－9 已知零件的极限应力σr=200MPa，许用安全系数[S]=2，影响零件疲劳强度的系数为Kσ=1.2，εσ＝0.83，βσ＝0.90。则许用应力为[σr]___________MPa。 （1）160.6 （2）106.7 （3）62.25 （4）110.7 3－10 绘制设计零件的σm－σa极限应力简图时，所必须的已知数据是___________。 （1）σ－1，σ0，σs,Kσ（2）σ－1，σ0，σs, KσD （3）σ－1，σs, ψσ，Kσ（4）σ－1，σ0，ψσ, KσD 3－11 在图示设计零件的σm－σa极限应力简图中，如工作应力点M所在的ON线与横轴间夹角θ＝45o，则该零件受的是___________。

疲劳强度的计算

摘要：零件的疲劳强度是一个值得深刻探讨的问题，在众多领域有着至关重要 的地位，零件的疲劳强度决定了其疲劳寿命，也就决定了对零件的选择和对这个器件的设计。本论文在参考多方资料，以及在平日学习中积累总结的经验之后，对零件疲劳强度的计算有了一些结论，得出影响导致零件疲劳的原因有破坏应力与循环次数之间量的变化影响，静应力的影响，应力集中的影响，零件绝对尺寸的影响，表面状态与强化的影响等方面。在分析零件疲劳产生原因之后，得出许多关系变化图与计算方法。运用这些计算方法，对零件疲劳极限进行了计算上的确定。并总结出疲劳强度在一些条件下的相关计算方法，如在简单应力状态，复杂应力状态下的不同。对疲劳强度安全系数的确定也进行了一系列分析，最后，尝试建立了疲劳强度的统计模型。 Abstract：The fatigue strength of parts is a worthy of deep discussion, have a vital role in many fields, the fatigue strength of parts determines its fatigue life, also decided on the part of the selection and the device design.This paper in reference to various data, and after the usual study accumulation experience, calculation of the fatigue strength of parts have some conclusion, that caused damage should change between force and the number of cycles of the causes of fatigue parts, the influence of static stress, effect of stress concentration, affects the absolute size, surface state and strengthening effect etc.. After the analysis of fatigue causes, draw many relationship graph and calculation method. Using the calculation method of fatigue limit, determined the calculation. And summarizes the related calculation under some conditions the method of fatigue strength, as in the simple stress state, the complex stress state under the different. Determination of the fatigue strength safety factor is also carried out a series of analysis, finally, try to establish a statistical model of fatigue strength. 关键词：零件疲劳寿命疲劳强度 Key word:Spare parts Fatigue life Fatigue strength

边坡强度折减法

基本原理： 强度折减法中边坡稳定的安全系数定义为：使边坡刚好达到临界破坏状态时，对岩、土体的抗剪强度进行折减的程度，即定义安全系数为岩土体的实际抗剪强度与临界破坏时的折减后剪切强度的比值。强度折减法的要点是公式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音频降噪课程设计报告.doc

燕山大学 医学软件课程设计说明书 题目：基于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谱减法音频降噪处理

题目基于matlab谱减法音频降噪处理 班级 学号 姓名 指导 时间 景德镇陶瓷学院

数字信号处理课程设计任务书

目录 1、设计要求. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .1 2、设计原理. . . . . . . ……………………………………………….. . . . . . . . . . . .2 3、源程序清单. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .7 4、设计结果和仿真波形. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .11 5、参考文献. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15. 6、设计心得体会. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16

1、设计要求 语言是人类最重要、直接、有效和便捷的交换信息的方式。随着近些年科学技术的飞速发展，人们也不满足于和计算机的信息交换方式，希望能够甩掉键盘和鼠标而实现用语言来对计算机进行控制。因此，语音信号处理技术便应运而生。语音信号处理是一门新兴的学科，同时也是综合多种学科和涉及面非常广泛的交叉学科。现在在一些职能系统中嵌入有语音处理系统，但它们只能在安静的环境中才能使用。然而，在语音信息的采集过程中难免会有各种噪声的干扰。噪声不仅降低了语音的可懂度和语音质量，还严重的影响语音处理的准确性，甚至使系统不能正常工作。本文将就对语音增强技术的原理和方法进行讨论，重点介绍语音增强的一种方法——谱减法及其改进算法。该方法能够有效消除平稳的加性噪声，其改进算法能够有效消除普通方法产生的“音乐噪声”，在很大程度上提高语音信号的信噪比。 目前，语言识别技术已经取得了重大进展，并开始进入实用阶段。但语音识别系统必须在相对比较安静的环境下运行，然而，在语言信息的采集中难免会有各种噪声的干扰，在较强的噪声背景下，语音识别系统的准确性会受到较大影响，甚至没法正常工作。所以在语音识别系统对语音信息处理前，应该对语音信息进行预处理，即背景噪声消除。语音背景噪声消除技术的出现使得语音识别技术更加稳定和精确，也使得语音信息的可懂度大大提高，使人们能够从较复杂的语音信息中提取到更多的有用信息。

高强度螺栓疲劳校核

16、轮盘连接高强度螺栓疲劳强度校核 说明： 轮盘在设备的设计使用寿命期限内，始终处于受压状态，其三根弦杆承受压力作用，轮盘的整体弯矩由内、外弦杆的压力调幅来平衡，弦杆法兰连接的高强度螺栓承受的、由单独弦杆的弯矩引起的交变力很小。 由于法兰结合面的载荷全部为压力载荷，故螺栓的工作应力都小于其预紧力，故螺栓的拉力载荷总在预紧力一下某一范围波动。对螺栓而言，保证法兰结合面不松开，其压力载荷越大，螺栓残余预紧力就越小，螺栓的拉力就越小。本文的计算模型转变为较小圆角过度的阶梯轴拉伸（如图一），校核过渡截面的疲劳应力。 观览车的运行速度很慢，每周循环的时间为20分钟，考虑50年的使用寿命期，每年300天，每天工作8小时，共运行300000次循环，选小于结构钢S-N曲线的转折点的循环次数，且本文的计算载荷为正常满载+15m/s风载的载荷情况，故计算结果有一定的保守性。 疲劳设计方法是一门以试验为基础的设计方法，本计算选取的疲劳性能数据选自国内公开的《机械设计手册》数据。 图一：计算模型

附：螺栓无限寿命校核说明书 一、螺栓参数和预紧力： 螺栓直径：M30x160 性能等级：10.9级 过渡圆角：r=0.5mm 螺栓材料的破断强度：1000MPa 螺栓副连接的相对刚度：m b b C C C +=0.25 选用的单个螺栓预紧力矩：Nm T 1600= 则预紧力：kN N d T Q p 2671067.2030 .02.016002.05=?=?== 二、螺栓组载荷： 主管法兰圆周应力分布及载荷谱： 530*30螺栓组主管件轴力， 六点方位N=-4729kN ，七点半N=-4487kN ，九点N=-3785kN ，十点半N=-3181kN ，十一点N=-2961kN ，十二点N=-2300kN ，一点N=-2960kN ，一点半N=-3253kN ，三点N=-3891kN ，四点半N=-4552kN 。 最大压力：kN F a 4729-= 换算到单个螺栓的最大压力载荷：kN F F a 39412/472912/-=-== 螺栓最小拉力：kN F F C C C Q Q m b b p 1680.25267min =+=++ = 最小压力：kN F a 2300-=

语音增强--谱减法原理及源代码实现

处理宽带噪声的最通用技术是谱减法。它利用语音信号的短时平稳特性，从带噪语音的短时谱值中减去噪声的短时谱，从而得到纯净语音的频谱，达到语音增强得目的。谱减法包括幅度谱减法和功率谱减法：幅度谱减法就是在频域中从带噪语音的幅度谱上减去噪声的幅度谱作为语音信号的幅度谱；功率谱减法则是从带噪语音的功率谱中减去噪声的功率谱，得到纯净语音的功率谱估计，通过开方运算得到幅度谱。由于人耳对语音频谱分量的相位感知不敏感，因此这些算法都是在幅度上进行的修正，相位部分则保持不变，仍然使用带噪语音的相位。 谱减法通过从带噪语音的短时谱估值中减去噪声的短时谱来达到语音增强得目的，算法简单且容易实现。但在减去噪声谱后，还会有些较大功率谱分量的剩余部分，在频谱上呈现出随机出现的尖峰，在听觉上形成残留噪声。这种噪声具有一定的节奏性起伏感，被称之为“音乐噪声”。后来，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); %对每段分帧进行加窗处理

有限元强度折减法和极限平衡法适用性研究

有限元强度折减法和极限平衡法适用性研究 曾红丽 戚明军 陕西省建筑设计研究院有限责任公司 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算例分析

基于Matlab的谱减法语音增强的研究

基于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） 而由于平稳噪声的功率谱在发声前和发声期间可以认为基本没有变化，这样可以通过发声前的所谓“寂静段”来估计噪声的功率谱，从而有：

安全系数算法

3 安全度分析 根据标准图的设计说明，隧道按照喷锚构筑法原理，衬砌结构由初支和二次衬砌组成，支护参数主要以工程类比为主，并辅以结构数值分析检算。计算时，初期支护为主要承载结构。Ⅱ~Ⅲ级围岩二次衬砌作为安全储备，按承受围岩荷载的30% 检算；Ⅳ~Ⅴ级围岩二次衬砌作为承载结构，分别按承受围岩荷载的50%~70% 检算，得出荷载与结构安全系数。 3.1 围岩压力计算 衬砌荷载根据隧道的地形和地质条件、埋置深度、结构特征和施工方法等因素，按有关公式计算或按工程类比确定，主要考虑围岩压力、结构自重、围岩约束衬砌变形的弹性反力等，不考虑列车活载、冻胀力、地下水压等附加荷载。当施工发现其与设计不符时，应及时修正。对复杂地质条件的隧道，必要时应通过实地量测确定荷载的计算值及其分布规律，本图考虑在浅埋地段的隧道视具体情况采用加强衬砌。 3.1.1 深埋隧道围岩压力计算 计算深埋隧道衬砌时，围岩压力按松散压力考虑，其垂直及水平匀布压力可按下列规定确定。 （1）竖直压力 10.452S q h γγω-=?=??? （3-1） 式中： q ——围岩垂直匀布压力（kPa ）； γ——围岩重度（kN/m3）； h ——围岩压力计算高度（m ）； S ——围岩级别； ω——宽度影响系数，1(5)i B ω=+-； B ——坑道宽度（m ）； i ——坑道宽度每增减1m 时的围岩压力增减率。当B<5m 时，取i ＝0.2， B>5m 时，可取i =0.1。 （2）侧压力 水平匀布压力可按下式计算确定。

e q λ=? （3-2） 式中：λ——侧压力系数，其取值参照围岩级别分别取值。 3.1.2 浅埋隧道围岩压力计算 地面基本水平的浅埋隧道，所受的荷载具有对称性。其计算为： （1）竖直压力 tan 1h q h B γθγ?? =- ?? ? （3-3） [] θ?θ?ββ?βλtan tan )tan (tan tan 1tan tan tan c c c +-+-= （3-4） θ ????βtan tan ) tan()1(tan tan tan 2-++=c c c c （3-5） a h h 5.2= （3-6） 10.452S a h ω-=?? （3-7） ()10.10.5B ω=+?- （3-8） （2）侧压力 λγi i h e = （3-9） 式中: q ——垂直压力（N/m 2）； γ——围岩重度（N/m3）； h ——洞顶地面高度(m)； θ——洞顶土柱两侧摩擦角（°）； λ——侧压力系数，按照围岩级别分别取值； h i ——内外侧任意点至地面的距离(m)； c ?——围岩计算摩擦角（°）； β——产生最大推力时的破裂角（°）； a h ——深埋隧道垂直荷载计算高度（m ）； S ——围岩级别； ω——深埋隧道的宽度影响系数； B ——隧道开挖跨度（m ）。

第三章疲劳强度计算练习题dayin

第三章机械零件的疲劳强度设计 三、设计计算题 3－47 如图所示某旋转轴受径向载荷F=12kN作用，已知跨距L=1.6m，直径d=55mm，轴的角速度为ω，求中间截面上A点的应力随时间变化的表达式，并求A点的σmax、σmin、σa和σm。 3－48 一内燃机中的活塞连杆，当气缸发火时，此连杆受压应力σmax=－150MPa，当气缸进气开始时，连杆承受拉应力σmin=50MPa，试求：（1）该连杆的平均应力σm、应力幅σa 和应力比r；（2）绘出连杆的应力随时间而变化的简图。 3－49 一转动轴如图所示，轴上受有轴向力F a=1800N，径向力F r=5400N，支点间的距离L=320mm，轴是等截面的，直径d=45mm。试求该轴危险截面上的循环变应力的σmax、σmin、σm、σa和r。 题3－49图题3－50图 3－50 某一转轴的局部结构如图所示，轴的材料为Q235普通碳钢，精车制成。若已知直径D=120mm，d=110mm，圆角半径r=5mm，材料的力学性能为：σb=450MPa，σs=220MPa，试求截面变化处的疲劳强度综合影响系数KσD和KτD。 3－51 由脆性材料制成的受弯平板的平面尺寸如图所示，板厚30mm。A、B两处各有一个直径5mm的穿透小孔，弯矩M=20kN·m。试分别计算Ⅰ、Ⅱ两截面上的最大应力值。疲劳缺口系数查题3－53附图。 3－52 一转轴的各段尺寸及其受载情况如图所示。所有圆角半径均为r=3mm。试分别计算Ⅰ—Ⅰ至Ⅶ—Ⅶ各截面上的最大弯曲应力的名义值和实际值。疲劳缺口系数查题3－53附图。

题3－51图题3－52图 3－53 用高强度碳钢制成的构件 的平面尺寸如图所示，厚8mm，受拉力 F=50kN。该构件的Ⅰ、Ⅱ、Ⅲ截面上分别 有φ15mm的圆孔、R7.5mm的半圆缺口 和R7.5mm的圆角。试分别计算这三个截 面上的最大应力。 题3－53附图 附注：这三种结构的疲劳缺口系数值可从上图曲线中查得。 3－54 题3－53中如载荷F在25～85kN之间做周期性的变化，材料改为20CrMnTi，其力学性能为σs=835MPa，σ-1=345MPa，σ0=615MPa。危险截面上的疲劳缺口系数Kσ=1.45，尺寸系数εσ＝0.75，表面状态系数βσ＝0.9，按无限寿命考虑。试画出σm－σa极限应力图，并用图解法和解析法确定安全系数Sσ。 3－55 用题3－54的条件画出σm－σmax和σmin极限应力图，并用图解法和解析法确定安全系数。可参阅[5]。 3－56 按题3－54的条件，除载荷F变为在－32～64kN之间作周期性变化外，其余条件不变。试画出σmin－σmax极限应力简图，并用图解法和解析法确定安全系数。可参阅[5]。 3－57 一阶梯轴轴肩处的尺寸为D=60mm，d=50mm，r=4mm，如材料的力学性能为：σb=650MPa，σs=360MPa，σ-1=200MPa，σ0=320MPa。试绘制此零件的简化极限应力线图。 3－58 如上题中危险截面处的平均应力σm=30MPa，应力幅σa=45MPa，试分别按（1）r=c；（2）σm=c求出该截面上的计算安全系数Sσ。 3－59 一转轴的危险截面上作用有周期性波动的载荷：弯矩M=100~200N·m，转矩T=0~100N·m。轴的材料为45钢，力学性能：σs=400MPa，σ-1=270MPa，σ0=480MPa，τs=216MPa，τ-1=156MPa，τ0=300MPa。若截面直径d=25mm，疲劳缺口系数Kσ=1.78，Kτ=1.45，尺寸系数εσ=0.9，ετ=0.93，表面状态系数βσ=0.91，βτ=0.95。试确定安全系数S。计算时可

基于改进型谱减法的语音增强技术研究

基于改进型谱减法的语音增强 摘要 本文主要研究改进型谱减算法在语音增强中的应用，目的是增强语音质量，减少语音失真和提高其可度懂。我们首先介绍了语音增强的研究意义，然后介绍了语音信号的相关理论，进而阐述了语音增强的基本谱减法的原理，并在此基础上提出了一种改进型谱减算法。该算法通过语音激活检测(端点检测法)来确定“寂静段”（纯噪声段），从而对噪声功率谱进行重新估计。为了减小基音检测算法可能产生的检测误差，采用了组合递归平滑法来减小噪声谱估计的误差。 整个仿真实验中，我们对引入的加性噪声进行处理，其噪声谱估计的性能可在本文中的MATLAB仿真实验中体现。 仿真结果表明，该算法在去除背景噪声的同时，保证了较小的语音失真、提高了信噪比，达到了较好的测听效果。 最后，基于噪声与语音具有一定的相关性的实际情况，我们提出了算法的进一步改进设想，并对此思想做出了数学推导，得到了算法进一步改进的方向及可行性。关键词：语音增强；谱减法；噪声估计；端点检测；组合递归平滑；仿真实验；改进算法

第1章背景介绍 1.1研究背景 人们在语音通信过程中不可避免地会受到来自周围环境、传输媒介引入的噪声、通信设备内部电噪声乃至其他讲话者的干扰。这些干扰最终将使接收者接收到的语音已非纯净的原始语音信号，而是受噪声污染的带噪语音信号。例如，安装在汽车、飞机或舰船上的电话，街道、机场的公用电话，常受到很强背景噪声的干扰，严重影响通话质量。又如，室内会议电话的交混回响随同语音广播到每个会议地点，影响收听效果。再如深海潜水员在氦-氧面罩内讲话引起的失真，语言障碍残疾人的语音失真，有历史价值的旧唱片、旧录音带的噪声和失真等，都是带噪语音信号的例子。 环境噪声污染使许多语音处理系统的性能急剧恶化。例如，语音识别己取得重大进展，正在步入实用阶段。但目前的识别系统大都是在安静环境中工作的，在噪声环境中尤其是强噪声环境，语音识别系统的识别率将受到严重影响。低速率语音编码，特别是参数编码(如：声码器)，也遇到类似问题。由于语音生成模型是低速率参数编码的基础，当模型参数的提取受到混杂在语音中背景噪声严重干扰时，重建语音的质量将急剧恶化，甚至变得完全不可懂。在上述情况下，语音增强作为一种预处理手段，不失为解决噪声污染的一种有效途径。 在实际需求的推动下，早在上个世纪60年代语音增强这个研究课题就引起人们的注意，此后40多年人们一直锲而不舍地进行这方面的研究。随着数字信号处理理论的成熟，70年代曾形成一个理论研究高潮，取得了一些基础性成果，并使语音增强发展成为语音信号处理的一个重要分支。进入80年代后，DSP(数字信号处理)技术的发展和成熟为语音增强的实时实现提供了可能。 语音增强不但与语音信号数字处理理论有关，而且涉及到人的听觉感知和语音学范畴。再者，噪声的来源众多，随应用场合而异，它们的特性也各不相同。即使在实验室仿真条件下，也难以找到一种通用的语音增强算法能适用于各种噪声环境，所以必须针对不同噪声，采用不同的语音增强对策。目前，某些语音增强算法在实