机械振动大作业——简支梁的各情况分析2

合集下载

机械振动大作业——简支梁的各情况分析

机械振动大作业——简支梁的各情况分析

机械振动大作业姓名:徐强学号:SX1302106专业:航空宇航推进理论与工程能源与动力学院2013年12月简支梁的振动特性分析题目:针对简支梁、分别用单、双、三、十个自由度以及连续体模型,计算其固有频率、固有振型。

单、双、三自由度模型要求理论解;十自由度模型要求使用李兹法、霍尔茨法、矩阵迭代法、雅可比法、子空间迭代法求解基频;连续体要求推导理论解,并通过有限元软件进行数值计算。

解答:一、 单自由度简支梁的振动特性如图1,正方形截面(取5mm ×5mm )的简支梁,跨长为l =1m ,质量m 沿杆长均匀分布,将其简化为单自由度模型,忽略阻尼,则运动微分方程为0=+••kx x m ,固有频率ωn =eqeq m k ,其中k 为等效刚度,eq m 为等效质量。

因此,求出上述两项即可知单自由度简支梁的固有频率。

根据材料力学的结果,由于横向载荷F 作用在简支梁中间位置而引起的变形为)(224348EI F -)(x l x x y -=(20l x ≤≤), 48EI F -3max l y =为最大挠度,则: eq k =δF=348EIl梁本身的最大动能为:)(224348EI F -)(x l xx y -==)(223max43x l l x y -T max =2×dx x y l m l 220)(21⎭⎬⎫⎩⎨⎧•⎰=2max 351721•y m )(如果用eq m 表示简支梁的质量等效到中间位置时的大小,它的最大动能可表示为:T max =2max21•y m eq所以质量为m 的简支梁,等效到中间位置的全部质量为: m m eq 3517=故单自由度简支梁横向振动的固有频率为:ωn =eqeq m k =3171680mlEImk图1 简支梁的单自由度模型二、 双自由度简支梁的振动特性如图2,将简支梁简化为双自由度模型,仍假设在简支梁中间位置作用载荷,根据对称性,等效质量相等,因此只要求出在3/l 处的等效质量即可。

振动力学(梁的横向振动)分析

振动力学(梁的横向振动)分析

EI 2 2 求出后得到固有频率 i i a i , (i 1, 2 ) A
Φ( x) C1 sin x C2 cos x C3 sh x C4 ch x 振型为
sin l sh l C sin x sh x (cos x ch x) ch l cos l
u ( x, t ) Φ( x)q(t )
代入方程得到
a x2
2 2
d Φ( x) d q(t ) q(t ) dx2 Φ( x) dt 2
2 2
写为
2 d 2Φ( x) d 2 q(t ) 2 2 2 x dx 2 dt a 2 Φ ( x) q(t )
l
l d d 2 i d d i j EI EI dx j 2 2 0 dx dx dx dx 0
2
l
d d i d i d j j EI EI 2 2 dx 0 dx dx 0 dx
2 2
d d i EI 2 0 i dx 2 dx
l 2 2
l d i 2 dx EI dx Mi i 2 0 dx
弹性体的振动
梁在激励力作用下的响应
和一维波动方程一样,用振型叠加法求响应
u ( x, t ) qi (t )Φ ( x)
u ( x, t ) u ( x, t ) EI k0 2 x x x l
2
x l
3u ( x, t ) EI x3
ku (l , t )
x l
弹性体的振动

(完整版)机械振动课后习题和答案第二章习题和答案

(完整版)机械振动课后习题和答案第二章习题和答案

2.1 弹簧下悬挂一物体,弹簧静伸长为δ。

设将物体向下拉,使弹簧有静伸长3δ,然后无初速度地释放,求此后的运动方程。

解:设物体质量为m ,弹簧刚度为k ,则:mg k δ=,即:n ω==取系统静平衡位置为原点0x =,系统运动方程为: δ⎧+=⎪=⎨⎪=⎩&&&00020mx kx x x (参考教材P14)解得:δω=()2cos n x t t2.2 弹簧不受力时长度为65cm ,下端挂上1kg 物体后弹簧长85cm 。

设用手托住物体使弹簧回到原长后无初速度地释放,试求物体的运动方程、振幅、周期及弹簧力的最大值。

解:由题可知:弹簧的静伸长0.850.650.2()m =-=V所以:7(/)n rad s ω=== 取系统的平衡位置为原点,得到:系统的运动微分方程为:20n x x ω+=&& 其中,初始条件:(0)0.2(0)0x x=-⎧⎨=⎩& (参考教材P14) 所以系统的响应为:()0.2cos ()n x t t m ω=- 弹簧力为:()()cos ()k n mg F kx t x t t N ω===-V因此:振幅为0.2m 、周期为2()7s π、弹簧力最大值为1N 。

2.3 重物1m 悬挂在刚度为k 的弹簧上并处于静平衡位置,另一重物2m 从高度为h 处自由落到1m 上而无弹跳,如图所示,求其后的运动。

解:取系统的上下运动x 为坐标,向上为正,静平衡位置为原点0x =,则当m 有x 位移时,系统有: 2121()2T E m m x =+& 212U kx =由()0T d E U +=可知:12()0m m x kx ++=&& 即:12/()n k m m ω=+系统的初始条件为:⎧=⎪⎨=-⎪+⎩&2020122m gx k m x gh m m (能量守恒得:221201()2m gh m m x =+&) 因此系统的响应为:01()cos sin n n x t A t A t ωω=+其中:ω⎧==⎪⎨==-⎪+⎩&200021122n m g A x k x m g ghk A k m m即:ωω=-2()(cos )n n m g x t t t k2.4 一质量为m 、转动惯量为I 的圆柱体作自由纯滚动,圆心受到一弹簧k 约束,如图所示,求系统的固有频率。

单自由度系统自由振动——简支梁

单自由度系统自由振动——简支梁

单自由度系统自由振动(简支梁)一、 实验目的 1、测定简支梁的等效弹簧常数k ; 2、记录简支梁的自由振动曲线,用分析仪测定系统的有阻尼时的固有频率d ω及相对阻尼系数ζ; 3、用附加质量法测定简支梁的等效质量m ; 4、初步了解振动测试的一些仪器设备及测试方法。

  二、 实验装置及原理 1、 实验装置 一根均匀的、截面为矩形的简支梁,其简图如图1所示。

这个系统可看作如图2所示的,有阻尼的单自由度弹簧质量系统,有阻尼时的振动微分方程为: 0=++kx x c xm &&& (1) 令m c n =2,mk n =2ω (2) 则(1)式为:022=++x x n x n ω&&& (3) 再令nn ωζ= (4) 则式(3)为:022=++x x x n n ωςω&&& (5) 其中: m :为简支梁系统的等效质量; k :为简支梁系统对于跨度中点的等效弹簧常数; c :为简支梁下的阻尼常数,n 称为衰减系数,ζ称为相对阻尼系数; n ω:为简支梁系统固有频率,n n f πω2=,d ω为系统的有阻尼固有频率,d d f πω2=。

 2、 实验原理 (1) 等效弹簧常数的测定 由于梁在弹性范围内的挠度与梁所受载荷成正比,因此只要在简支梁的跨中点加载,同时图2用百分表读出该点的挠度值,即可测出等效弹簧常数。

 (2)记录简支梁系统的自由振动曲线 在简支梁跨度中点贴应变片作用是使梁在振动时的应变量变化转化成电阻量的变化,再将应变片按半桥接法接到动态应变仪上,把电阻量的变化信号放大,并转化成电压量的变化信号,输出到示波器或分析仪,这样即可观察和记录波形。

测试系统框图如图3所示。

(3)附加质量法测等效质量 根据式(2),因为()222n n f m k πω==,21ζωω−=n d ,d d f πω2=要测出简支梁的等效质量m ,只要在原来的简支梁上附加一个已知质量∆,再次求得带有附加质量∆时的固有频率2∆n ω,然后通过下式计算得到m : ()()()()22222222∆∆∆==∆+=n n n n n n f f f f m m ππωω (6) ()()1111222222−∆=−−−∆=∆∆∆d d d d f f f f m &ζζ (7)  三、 实验步骤 1、 测定简支梁系统的等效弹簧常数 在简支梁跨中点处用砝码加载(i=1,2, …., 5),同时用百分表读出该点相对应的挠度值,并记录表1中,按公式算出。

简支梁的变形与振动分析

简支梁的变形与振动分析

简支梁的变形与振动分析简支梁是一种常见的结构形式,广泛应用于桥梁、楼板等工程中。

在实际工程项目中,我们需要对简支梁的变形和振动进行分析,以确保结构的安全性和稳定性。

本文将从数学模型到应用实例,全面深入地探讨简支梁的变形与振动分析。

一、简支梁的基本理论简支梁是在两端支座的约束下,承受集中力或均布力作用下的一种结构形式。

为了研究其变形和振动特性,我们需要建立数学模型。

1. 简支梁的受力分析在进行简支梁的变形和振动分析前,首先需要了解其受力情况。

在两端支座的约束下,简支梁主要受到弯矩和剪力的作用。

通过弯矩和剪力的分析,可以得出简支梁的受力公式,进而计算结构在承受力作用下的变形。

2. 简支梁的变形分析简支梁在受力作用下会发生一定的变形。

根据梁的假设和力学原理,可以建立简支梁的弹性变形方程。

通过求解弹性变形方程,可以得到简支梁在各个位置的变形情况。

3. 简支梁的振动分析在实际工程中,简支梁还可能受到外力的激励,导致振动现象的发生。

为了分析简支梁的振动特性,我们可以建立简支梁的振动微分方程,并求解得到简支梁的振动模态。

二、简支梁的应用实例1. 桥梁工程简支梁在桥梁工程中得到广泛应用。

为了确保桥梁在运行过程中的安全性和稳定性,需要进行简支梁的变形与振动分析。

通过分析得到的变形和振动数据,可以对桥梁的结构参数进行优化,提高桥梁的工作性能。

2. 建筑结构在楼板、屋顶等建筑结构中,简支梁也扮演着重要的角色。

在设计建筑结构时,需要对简支梁进行变形与振动分析,以确保结构的稳定性和安全性。

通过合理调整支座位置或增加梁的截面尺寸,可以改善简支梁的变形和振动特性。

三、总结简支梁的变形与振动分析对于工程项目的设计和施工至关重要。

通过建立数学模型,进行受力分析和变形分析,可以预测结构在实际工况下的变形情况。

同时,通过振动分析,可以了解简支梁的振动特性,为结构的稳定性提供参考。

在实际工程中,我们还可以利用现代软件进行简支梁的有限元分析,获得更加准确的变形和振动数据。

机械振动大作业——简支梁的各情况分析2

机械振动大作业——简支梁的各情况分析2

图2 简支梁的双自由度模型
三、三自由度简支梁的振动特性
如图3,将简支梁简化为三自由度模型,按照双自由度类似的等效思想,可得等效质量:
因此,质量矩阵为:
由机械振动中文教材例6.6可知,系统的柔度矩阵为:
其中, 。
动力矩阵:
令特征行列式为零,得到频率方程为:
其中, ,将上式整理得:
其中, 。
利用Matlab软件,求解上述方程的根为:
2.2447
0.0629
0.1282
-0.0388
1.8642
0.0834
0.0932
-0.0347
1.3333
0.0768
0.0568
-0.0237
0.6941
0.046
0.0258
-0.011
将各列分别归一化得:
求得 和 分别为:
再由李兹法得特征值问题为:
解出:
, , ,
其中, ,相应的主振型为:
机械振动大作业
*****
学号:SX*******
专业:航空宇航推进理论与工程
能源与动力学院
2013年12月
简支梁的振动特性分析
题目:针对简支梁、分别用单、双、三、十个自由度以及连续体模型,计算其固有频率、固有振型。单、双、三自由度模型要求理论解;十自由度模型要求使用李兹法、霍尔茨法、矩阵迭代法、雅可比法、子空间迭代法求解基频;连续体要求推导理论解,并通过有限元软件进行数值计算。
0.0306
0.0158
0.0309
0.0443
0.0552
0.0628
0.0664
0.0650
0.0579
0.0443
0.0240

简支梁实验

简支梁实验

简支梁模态参数测定之一—测定固有频率与振型一、实验目的1、加深对系统固有频率和主振型的理解;2、掌握振动系统固有频率及主振型的一种测量方法(共振法);3.了解压电式传感器及与它相配的测量系统的工作原理,掌握正确使用的方法;4、了解激振系统的工作原理。

二、实验装置框图图1 表示实验装置系统框图图1 实验装置系统框图三、实验原理试验模态分析法是确定结构固有频率的有效方法,在结构分析中应用广泛,而简支梁也是桥梁结构中一种常见的模型,现代桥梁中依然存在不少采用简支梁模型的桥梁结构。

所以本事通过试验模态法得到简支梁的固有频率和振型,也是桥梁结构分析中一种常用的方法很有实际意义,实验所用的均质等截面简支梁模型,属于小阻尼和连续的无限自由度的振动系统。

本实验模型是一矩形截面简支梁,它是一无限自由度系统。

理论上说,它应有无限个固有频率和主振型,在一般情况下,梁的振动是无穷多个主振型的迭加。

如果给梁施加一个合适大小的激振力,且该力的频率正好等于梁的某阶固有频率,就会产生共振,对应于这一阶固有频率而确定的振动形态叫做这一阶主振型,这时其它各阶振型的影响小得可以忽略不计。

用共振法确定梁的各阶固有频率及振型,具体步骤是首先得找到梁的各阶固有频率,并让激扰力频率等于某阶固有频率,使梁产生共振,然后,测定共振状态下梁上各测点的振动加速度值,从而确定前三阶主振型。

振型:即振动形态,即梁上各个测量点和振幅的关系图。

如图所示为一阶,二阶和三阶的振型图。

在正弦激励下振幅的比值等于加速度的比值。

所以本次试验测量加速度与位置之间的关系就能正确画出振型,大致如图2所示。

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 10 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.80.9 10 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.80.9 1图2 前三阶振型图根据材料力学理论下简支梁固有频率的计算:2012f l ππ⎛⎫= ⎪⎝⎭E 为材料的弹性模量,查表取E=210Gpa 测量得简支梁b=0.05m h=0.15m l=1m312bh I =s 为梁的横截面积37850kgm ρ=2201135.1622f Hzl l ππππ⎛⎫⎛⎫=== ⎪ ⎪⎝⎭⎝⎭10f f =214140.6f f Hz == 319316.4f f Hz==四、实验方法1、 激振器安装把激振器安装在支架上,将激振器和支架固定在实验台基座上,并保证激振器顶杆对简支梁有一定的预压力(不要超过激振杆上的红线标识),用专用连接线连接激振器和DH1301输出接口。

振动实验报告讲解

振动实验报告讲解

振动实验报告讲解振动与控制系列实验姓名:李⽅⽴学号:201520000111电⼦科技⼤学机械电⼦⼯程学院实验1 简⽀梁强迫振动幅频特性和阻尼的测量⼀、实验⽬的1、学会测量单⾃由度系统强迫振动的幅频特性曲线。

2、学会根据幅频特性曲线确定系统的固有频率f 0和阻尼⽐。

⼆、实验装置框图图3.1表⽰实验装置的框图图3-1 实验装置框图KCX图3-2 单⾃由度系统⼒学模型三、实验原理单⾃由度系统的⼒学模型如图3-2所⽰。

在正弦激振⼒的作⽤下系统作简谐强迫振动,设激振⼒F 的幅值B 、圆频率ωo(频率f=ω/2π),系统的运动微分⽅程式为:扫频信号源动态分析仪计算机系统及分析软件打印机或绘图仪简⽀梁振动传感器激振器⼒传感器质量块M或 M F x dt dxdt x d M F x dt dx n dtx d FKx dt dx C dtx d M /2/222222222=++=++=++ωξωω(3-1)式中:ω—系统固有圆频率ω =K/Mn ---衰减系数 2n=C/M ξ---相对阻尼系数ξ=n/ωF ——激振⼒ )2sin(sin 0ft B t B F πω== ⽅程①的特解,即强迫振动为:)2sin()sin(0?π?ω-=-=f A A x (3-2)式中:A ——强迫振动振幅--初相位20222024)(/ωωωn M B A +-=(3-3)式(3-3)叫做系统的幅频特性。

将式(3-3)所表⽰的振动幅值与激振频率的关系⽤图形表⽰,称为幅频特性曲线(如图3-3所⽰):3-2 单⾃由度系统⼒学模型 3-3 单⾃由度系统振动的幅频特性曲线图3-3中,Amax 为系统共振时的振幅;f 0为系统固有频率,1f 、2f 为半功率点频率。

振幅为Amax 时的频率叫共振频率f 0。

在有阻尼的情况下,共振频率为:221ξ-=f f a (3-4) 当阻尼较⼩时,0f f a =故以固有频率0f 作为共振频率a f 。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

机械振动大作业姓名:徐强学号:SX1302106专业:航空宇航推进理论与工程能源与动力学院2013年12月简支梁的振动特性分析题目:针对简支梁、分别用单、双、三、十个自由度以及连续体模型,计算其固有频率、固有振型。

单、双、三自由度模型要求理论解;十自由度模型要求使用兹法、霍尔茨法、矩阵迭代法、雅可比法、子空间迭代法求解基频;连续体要求推导理论解,并通过有限元软件进行数值计算。

解答:一、单自由度简支梁的振动特性如图1,正方形截面(取5mm×5mm,质量m沿杆长均匀分布,将其简化为单自由度模型,忽略阻尼,则运动固有频率ωn其中k等效质量。

因此,求出上述两项即可知单自由度简支梁的固有频率。

根据材料力学的结果,由于横向载荷F作用在简支梁中间位置而度,则:梁本身的最大动能为:T max=2大动能可表示为:T max所以质量为m的简支梁,等效到中间位置的全部质量为: 故单自由度简支梁横向振动的固有频率为:ωn图1 简支梁的单自由度模型二、 双自由度简支梁的振动特性如图2,将简支梁简化为双自由度模型,仍假设在简支梁中间位置作用载荷,根据对称性,等效质量相等,去小量得:所以,质量矩阵为:双自由度简支梁的柔度矩阵:在所以,柔度矩阵为:动力矩阵:令特征行列式为零,得到频率方程为:解上述方程的根为:第一、二阶主振型分别为:图2 简支梁的双自由度模型三、 三自由度简支梁的振动特性如图3,将简支梁简化为三自由度模型,按照双自由度类似的等效思想,可得等效质量:因此,质量矩阵为:由机械振动中文教材例6.6可知,系统的柔度矩阵为:动力矩阵:令特征行列式为零,得到频率方程为:利用Matlab 软件,求解上述方程的根为:第一、二、三阶主振型分别为:图3 简支梁的三自由度模型四、十自由度简支梁的数值方法将简支梁简化为十自由度模型(如图4)。

图4 简支梁的十自由度模型通过在一点施加单位力,计算其余点的挠度,可得柔度矩阵:0.0137 0.0240 0.0306 0.0339 0.0344 0.0324 0.0284 0.0227 0.0158 0.0081 0.0240 0.0443 0.0579 0.0650 0.0664 0.0628 0.0552 0.0443 0.0309 0.0158 0.0306 0.0579 0.0787 0.0904 0.0934 0.0891 0.0787 0.0633 0.0443 0.0227 0.0339 0.0650 0.0904 0.1071 0.1131 0.1093 0.0973 0.0787 0.0552 0.0284 0.0344 0.0664 0.0934 0.1131 0.1229 0.1212 0.1093 0.0891 0.0628 0.0324 0.0324 0.0628 0.0891 0.1093 0.1212 0.1229 0.1131 0.0934 0.0664 0.0344 0.0284 0.0552 0.0787 0.0973 0.1093 0.1131 0.1071 0.0904 0.0650 0.0339 0.0227 0.0443 0.0633 0.0787 0.0891 0.0934 0.0904 0.0787 0.0579 0.0306 0.0158 0.0309 0.0443 0.0552 0.0628 0.0664 0.0650 0.0579 0.0443 0.0240 0.0081 0.0158 0.0227 0.0284 0.0324 0.0344 0.0339 0.0306 0.0240 0.0137表1十自由度简支梁为十个集中质量的振动模型,每个质量都近似等动力矩阵为:下面,用如下几种方法计算十自由度简支梁的固有频率与振型。

1、邓克莱法利用邓克莱法求基频(比准确值小):2、瑞利法(1)瑞利第一商柔度矩阵求逆得刚度矩阵:2。

2.0433 -1.9003 0.7778 -0.2649 0.1647 0.0356 -0.2817 0.2834 -0.0446 -0.0926-1.9003 2.9228 -2.4675 1.4375 -0.6862 0.0299 0.5193 -0.6226 0.3193 -0.0446 0.7778 -2.4675 3.8387 -3.3468 1.6331 -0.1417 -0.6684 0.8588 -0.6226 0.2834-0.2649 1.4375 -3.3468 4.2339 -2.9123 0.778 0.4309 -0.6684 0.5193 -0.2817 0.1647 -0.6862 1.6331 -2.9123 3.4193 -2.2932 0.778 -0.1417 0.0299 0.0356 0.0356 0.0299 -0.1417 0.778 -2.2932 3.4193 -2.9123 1.6331 -0.6862 0.1647 -0.2817 0.5193 -0.6684 0.4309 0.778 -2.9123 4.2339 -3.3468 1.4375 -0.2649 0.2834 -0.6226 0.8588 -0.6684 -0.1417 1.6331 -3.3468 3.8387 -2.4675 0.7778 -0.0446 0.3193 -0.6226 0.5193 0.0299 -0.6862 1.4375 -2.4675 2.9228 -1.9003 -0.0926 -0.0446 0.2834 -0.2817 0.0356 0.1647-0.2649 0.7778 -1.9003 2.0433表2假设力作用在简支梁中间位置而得到各点的静变形,可以表示为:因此,可以假设振型:(2)瑞利第二商同样假设力作用在简支梁中间位置,由瑞利第二商公式:可得3、兹法将十自由度简支梁缩减为三自由度,假设振型为:则可求出:由式,得:以及:所以系统的前三阶主振型的近似为:4、矩阵迭代法单位力作用在简支梁中间位置得到各点的挠度变形,将首项化72.96 0 -0.6 0 23.06 1.2 -0.6 1.2 180.1279 0 0.13540 4.3927 -1.2322 0.1354 -1.2322 4.2842一,得:因此可假设振型:利用矩阵迭代求第一阶固有频率和主振型:由上,仅3的第一阶主振型为:第一阶固有频率为:5、雅可比法根据雅可比法原理,依次找出上三角非对角线上(考虑对称性)的最大元素,代入旋转矩阵,可得:0.349 -0.056 -0.417 0.215 0.396 -0.394 -0.217 0.415 0.058 -0.3470.396 0.323 -0.134 -0.420 -0.212 0.212 0.420 0.134 -0.324 -0.395 0.398 -0.369 -0.032 0.334 -0.185 -0.215 0.340 0.002 -0.436 0.455 0.321 0.424 0.230 -0.122 -0.387 -0.387 -0.122 0.231 0.423 0.3220.161 -0.324 0.419 -0.346 0.087 0.244 -0.460 0.454 -0.282 0.102 0.232 0.388 0.422 0.322 0.120 -0.120 -0.322 -0.422 -0.388 -0.2320.104 -0.231 0.354 -0.442 0.466 -0.424 0.351 -0.263 0.155 -0.053 0.120 0.231 0.322 0.388 0.422 0.422 0.388 0.322 0.231 0.120 0.439 0.101 -0.383 -0.219 0.304 0.337 -0.236 -0.391 0.127 0.415 -0.404 0.453 -0.158 -0.184 0.324 -0.242 0.044 0.217 -0.448 0.399则:0.0005 0 0 0 0 0 0 0 0 00 0.0027 0 0 0 0 0 0 0 0 0 0 0.0003 0 0 0 0 0 0 0 0 0 0 0.0084 0 0 0 0 0 00 0 0 0 0.0001 0 0 0 0 00 0 0 0 0 0.0424 0 0 0 0 0 0 0 0 0 0 0.0001 0 0 0 0 0 0 0 0 0 0 0.6775 0 0 0 0 0 0 0 0 0 0 0.0011 0 0 0 0 0 0 0 0 0 0 0.0002由此得到系统的固有频率:6、子空间迭代法取假设振型:由动力矩阵迭代得到:0.6941 -0.046 0.0341 -0.01761.3333 -0.0768 0.0706 -0.02871.8642 -0.0834 0.108 -0.03122.2447 -0.0629 0.138 -0.03042.4421 -0.0234 0.1537 -0.03182.4421 0.0234 0.1506 -0.03612.2447 0.0629 0.1282 -0.03881.8642 0.0834 0.0932 -0.03471.3333 0.0768 0.0568 -0.02370.6941 0.046 0.0258 -0.011 将各列分别归一化得:再由兹法得特征值问题为:解出:所以:0.2850.0527 -0.0249 -0.5397 0.5477 0.0679 -0.0244 -0.9033 0.7661 0.0332 0.0048 -0.985 0.9229 -0.0198 0.0311 -0.7463 1.0042 -0.0589 0.0222 -0.2785 1.0042 -0.0607 -0.0147 0.279 0.9229 -0.0196 -0.0367 0.7477 0.7661 0.037 -0.0146 0.9855 0.5477 0.0665 0.0277 0.901 0.2850.04950.0368 0.5365各列分别归一化后,得:0.283808 0.776141-0.67663 -0.54764 0.545409 1-0.66304 -0.91659 0.762896 0.488954 0.130435 -0.99949 0.91904-0.29161 0.845109 -0.75728 1 -0.86745 0.603261 -0.2826 1 -0.89396 -0.39946 0.283105 0.91904 -0.28866 -0.99728 0.7587010.7628960.544919-0.3967410.545409 0.979381 0.752717 0.9142570.283808 0.729013 1 0.544394 重复上述过程进行第二次迭代。

相关文档
最新文档