Matlab求最大特征根

程序一

function dx=Lorenz(t,x); dx(1,1)=10*(x(2)-x(1)); dx(2,1)=x(1)*(30-x(3))-x(2); dx(3,1)=x(1)*x(2)-8/3*x(3); dx(4,1)=0; dx(5,1)=0; dx(6,1)=0;

function lambda_1=lyapunov_wolf1(data,N,m,tau,P)

% 该函数用来计算时间序列的最大Lyapunov 指数--Wolf 方法 % m: 嵌入维数 % tau:时间延迟 % data:时间序列 % N:时间序列长度

% P:时间序列的平均周期,选择演化相点距当前点的位置差，即若当前相点为I，则演化相点只能在|I－J|>P的相点中搜寻

% lambda_1:返回最大lyapunov指数值

%************************************************************************** % ode计算整数阶系统的时间序列

%*****************************************************************

用MATLAB编程求π值

上海电力学院

高级程序设计（C） 题 目： 院 系：专业年级：学生姓名：指导教师：课程设计报告

用MATLAB编程求π

计算机科学与技术学院 信息安全2012级 涂桂花 学号： 20123333 魏为民

2015年3月14日

用MATLAB编程求π值

目录

一、实验目的 .......................................................................... 1 二、实验内容和步骤 .............................................................. 1 1. 操作环境 ........................................................................ 1 2. 系统配置 ........................................................................ 1 3. 操作步骤 ...............

摘要

摘 要

随着计算机技术的发展，语音交互已经成为人机交互的必要手段,语音特征参数的精确度直接影响着语音合成的音质和语音识别的准确率。因此语音信号参数提取是语音信号处理研究中一项非常有意义的工作。

要提取语音信号特征参数，首先必须对语音信号进行数字化处理，本文采用Matlab录音工具录音，由计算机自动完成模拟信号到数字信号的转换。其次，进行预处理，包括预加重、加窗和分帧，本文讨论了预处理中各种参数的选择，以使信号特征提取更加准确。第三，讨论了各种时域特征参数的算法，包括短时能量、短时过零率、短时自相关函数和短时平均幅度差函数，比较了浊音、清音特征参数的区别。最后，利用算得的特征参数，完成了清音、浊音、静音的判别和浊音基音周期的提取。

本文采用Matlab软件提取参数，仿真判别清、浊音的不同算法，并从准确性和高效性两方面，结合软件可视化功能，对算法进行了评价。运用Matlab实现语音信号特征参数提取具有方法简单、速度快、误差小等优点，能够大大提高研究精度和效率。

关键词：语音信号；特征参数；提取；Matlab

I

Abstract

ABSTRACT

With the development of computer technology, spe

求最大公因数

说出下列各组数的最大公因数。(口答)第一组： 12和4 4 第二组： 1和10 1

13和39 1399和98 1

8和64 8 7和13 1

倍数关系：最大公因数是小的数。

互质关系：最大公因数是1。

第三组： 12和18 6

30和20 10

24和36 12

一般关系：??????

1、口答。

24的因数有： 1、2、3、4、6、8、12、24 。 36的因数有： 1、2、3、4、6、9、12、18、36。 24和36的最大公因数是： 12 。

2、用短除法分解质因数。 24 36 24=2×2×2×3 相同的质因数相 乘的积，就是他 36=2×2×3×3

们的最大公因数。

填 一 填

甲数= 2×3×5 乙数= 2×5×7 甲数和乙数的最大公因数是：。

短除法:

2

12 18

……用公因数2去除 ……用公因数3去除 ……除到公因数只有1为止

3 6 2

9 3

12和18的最大公因数是：2×3=6

(六)方法： 一、列举法：

找 最 大 公 因 数

1.先找各个数的因数。 2.找出两个数公有的因数。 3.确定最大公因数。 二、用倍数关系找： 如果两个数是倍数关

求最大公因数

说出下列各组数的最大公因数。(口答)第一组： 12和4 4 第二组： 1和10 1

13和39 1399和98 1

8和64 8 7和13 1

倍数关系：最大公因数是小的数。

互质关系：最大公因数是1。

第三组： 12和18 6

30和20 10

24和36 12

一般关系：??????

1、口答。

24的因数有： 1、2、3、4、6、8、12、24 。 36的因数有： 1、2、3、4、6、9、12、18、36。 24和36的最大公因数是： 12 。

2、用短除法分解质因数。 24 36 24=2×2×2×3 相同的质因数相 乘的积，就是他 36=2×2×3×3

们的最大公因数。

填 一 填

甲数= 2×3×5 乙数= 2×5×7 甲数和乙数的最大公因数是：。

短除法:

2

12 18

……用公因数2去除 ……用公因数3去除 ……除到公因数只有1为止

3 6 2

9 3

12和18的最大公因数是：2×3=6

(六)方法： 一、列举法：

找 最 大 公 因 数

1.先找各个数的因数。 2.找出两个数公有的因数。 3.确定最大公因数。 二、用倍数关系找： 如果两个数是倍数关

摘要

摘 要

随着计算机技术的发展，语音交互已经成为人机交互的必要手段,语音特征参数的精确度直接影响着语音合成的音质和语音识别的准确率。因此语音信号参数提取是语音信号处理研究中一项非常有意义的工作。

要提取语音信号特征参数，首先必须对语音信号进行数字化处理，本文采用Matlab录音工具录音，由计算机自动完成模拟信号到数字信号的转换。其次，进行预处理，包括预加重、加窗和分帧，本文讨论了预处理中各种参数的选择，以使信号特征提取更加准确。第三，讨论了各种时域特征参数的算法，包括短时能量、短时过零率、短时自相关函数和短时平均幅度差函数，比较了浊音、清音特征参数的区别。最后，利用算得的特征参数，完成了清音、浊音、静音的判别和浊音基音周期的提取。

本文采用Matlab软件提取参数，仿真判别清、浊音的不同算法，并从准确性和高效性两方面，结合软件可视化功能，对算法进行了评价。运用Matlab实现语音信号特征参数提取具有方法简单、速度快、误差小等优点，能够大大提高研究精度和效率。

关键词：语音信号；特征参数；提取；Matlab

I

Abstract

ABSTRACT

With the development of computer technology, spe

第15章 求矩阵特征值和特征向量

幂 法

幂法规范化算法

1. 输入矩阵A、初始向量u，误差eps 2. k?1

3. 计算V(k) ?Au(k-1)

4. mk ?max(V), mk-1 ?max(V) 5. uk ? V(k)/mk

(1)

6. 如果| mk - mk-1|

注：如上算法中的符号max(V)表示取向量V中绝对值最大的分量。本算法使用了数据规范化处理技术以防止计算过程中出现益出错误。

(k)

(k-1)

(0)

规范化幂法程序

Clear[a,u,x];

a=Input[\系数矩阵A=\；

u=Input[\初始迭代向量u(0)=\； n= Length[u];

eps= Input[\误差精度eps =\；

nmax=Input[“迭代允许最大次数nmax=”]; fmax[x_]:=Module[{m=0,m1,m2},

Do[m1=Abs[x[[k]]];

If[m1>m,m2=x[[k]];m=m1], {k,1,Length[x]}]; m2] v=a.u; m0=fmax[

2016届毕业生 毕业论文

题 目： 非线性方程求重根方法研究 院系名称： 理学院 专业班级： 学生姓名： 学 号： 指导教师： 教师职称：

2016年05月20日

摘 要

随着科学技术的发展，在现代科学和工程技术中，经常会遇到大量而复杂的数学计算问题。这些问题常常归结为非线性方程求根的问题。求解非线性方程的单根已经具有了比较成熟和丰富的构造技术手段。例如，其中在工程和其他领域的科学计算中的广泛应用迭代算法，它从某个初始点出发，由迭代格式生成一种收敛于方程根的序列。这些方法在面对非线性方程单根的时候可以很好的解决问题，然而这些方法在求解非线性方程的重根时，构造的算法显得相当的复杂甚至是无效的。举一个简单的例子就是平时我们经常研究的经典的牛顿迭代法。它对方程的单根二阶收敛，但是对于于方程的重根只能线性收敛，并且收敛速度变慢。因此非线性方程重根的高阶，尤其是最优解的迭代格式如何构造是一项具有挑战性的工作。直到

利用Matlab绘制系统的根轨迹

本章前面的内容介绍了控制系统根轨迹的绘制以及利用系统大致的根轨迹图分析系统性能的方法，若要由根轨迹获得系统在某一特定参数下准确的性能指标或者准确的闭环极点，需要依据幅值条件精确地作图。如果利用MATLAB工具箱中函数，则可方便、准确地作出根轨迹图，并利用图对系统进行分析。

MATLAB工具箱中，求系统根轨迹的几个常用函数有rlocus, rlocfind, sgrid，下面通过具体的例子来说明这些函数的应用。

例4-13 控制系统的开环传递函数为 G(s)H(s)=

绘制系统的根轨迹图。

解 利用函数rlocus函数可直接作出系统的根轨迹图，程序如下： % example4-13 % num=[1,5];

dun=[1,6,11,6,0]; rlocus(num,dun) 执行该程序后，可得到如图4-20所示的根轨迹。

图4-20 例4-13题根轨迹图

利用函数rolcus可画出系统的根轨迹图后，可用rlocfind函数在根轨迹上选择任意极点，得到相应的开环增益 和其它闭环极点。

例4-14 控制系统

clear clc

s=[0 0 0 0 0 0 35 50 2 40 50 160 75 150 4 80 150 200 115 250 14 180 475 300 150 350 24 280 800 400 250 420 36 565 1600 800 ];

x=input('请输入浓度值：') ];

for i=1:6 for j=1

if x(i)

elseif s(j,i)<=x(i)&&x(i)<=s(j+1,i)

r(j,i)=(s(j+1,i)-x(i))/(s(j+1,i)-s(j,i)) else

r(j,i)=0 end end

for j=2:4

if s(j-1,i)<=x(i)&&x(i)<=s(j,i)

r(j,i)=(x(i)-s(j-1,i))/(s(j,i)-s(j-1,i)) elseif s(j,i)<=x(i)&&x(i)<=s(j+1,i) r(j,i)=(s(j+1,i)-x(i))/(s(j+1,i)-s(j,i)) else

r(j,i)=0 end end for j=5

if x(i)>s(j,i) r(j,