=0即全体特征值为00,,00由需要P1AP=→PAP+7=1pAP+PP|=P(4+1)P=14+1→A+1=-1注(1)若AB(相似),则AB有相同特征值A,可引入记号:谱集(4)={2,2,…,λ}(全体特征值,含重复)A∽B→o()=o(B)(2)A∽B→1-A=1-B-(2-4元一2)…(-n),特征多项式PAP=B=A-A=p(1-A)P=A-B引理:若A0A2,则M-A|-|M1-4|-1-A1|2-A2→ar(4)=o(A)∪a(42k+1,Ak-2,…n1f(x2)设B,f(x)为多项式,则f(B)=of(,)引理:若n阶方阵A的谱集(4)=1,42,…},则)的全体特社值为)2,…,),x)为多项式Pf:由许尔定理,A∽B→f(4)∽f(B)f(x)的全体特征值为(A1)(42),,()},fx)为多项式例如:4为A的特征值→x为4的特征值。(x)=x)共113页矩阵理论A笔记第3页北京航空航天大学张京蕊工程系统工程系引理:令B,f(x)=x-B|=(x-41)(x-12)….(x-n)则fB)=(B-1D(B-21)…(B-A1D=0Pf:当n=2时,B=0x2f(x)=(x-1)(x-2)000→f(B)-(B-41)(B-21)(2-元)0(00∴得证★Cayley公式:设n阶方阵A的特征多项式为f(x)=|x-A|=a+a1x+…,+x则f4)=anl+a14+…,+4=0Pf:由许尔PAP=B=→P(4)P=fp3P)=f(B)=0(引理)定义若多项式x)使(4)=0,则称(x)为A的个零化式结论方阵A的特征多项式)=1x1-4为A的一个零化式g特征多项式fx)=x2可知:f(A)=A2+1=+I=00-1Hx)=|xI-A|=(x-)(x+i,(i=√-1,t2=-1)f(A)=(A-i)(4+i1=0也可取P=则PPAP=,对角形共113页矩阵理论A笔记第4页北京航空航天大学张京蕊工程系统工程系g:知A则A"=0Onxn由Cayley特征多项式:f(x)=x"→f(4)=4"=0Ex1.A=求P使得PP为对角阵,并验证Cayley定理2.A=cd/,求fx)=x1-4验证f4)-0补充知识(schur公式、Cayley公式)应用由A"=-(a0I+a1A+1A·AanA+a142+…+a.,A把①代入②→Am1=(-)+(+)4+…+(+)41可知:任何和(m≥n)都可写成,4,,A的线性组合任何多项式g(A),可写成lA,…,4的组合。Fg:若A|≠0,fx)=xI-A|=a0+a1x+…+x",ao=|-A|≠则A可用A的多项式表示∵a1A+a242+…+an21A-+A"--a072A(a1+a24+…+an-142+A)Aa1+…+an1A"2+A-1零化式定义:若g(x)=b+b1x+…+bnx,使得g(4)=bn+b14+…+bn4m=0,称g(x)为方阵A的零化式注:方阵A的零化式有无穷多个∴取特征多项式x)则4)=0任取式M(x),f(A(4)=0→f(x)(x)也是零化式极小式定义:在方阵A的零化式集合中,去次数最小的且首项系数为1的零化式m(x),称它为A的极小式共113页矩阵理论A笔记第5页北京航空航天大学张京蕊工程系统工程系注:极小式唯一性质:①极小式m(x)必为特征多项式fx)=|xI-A的因式。②特征多项式fx)=|x1-A的每个单因子(x-4)也是极小式的因子)f(x)=|x1-4=(x-x)(x-2)则极小式m(x)=(x-x)(x-2)y…(x-,),且1≤l1≤m1,1≤l2≤m2,…,1≤l≤n,41,A2…,n互不相同210EgA=020,B=020,求极小式mA(),m()解:(1)|xI-A|=(x-2)(x-1)极小式为:(x-2)(x-1)或(x-2)(x-1)计算:(4-2/)4-1)=000010k≠000000∴极小式为m4(x)=(x-2)2(x-1)(2)|-B|-=(x-2)2(x-1)00000计算:(B-2)B-1)=000010=000-1八000∴极小式为m(x)=(x-2)(x-1)Eg求下列极小式m(x)4604-60(1)A=-3-50,(2)B=2-303-6100210(3)C,(4)D=000010002000解:(1)特征多项式|x7-A|-(x-1)(x+2)极小式为:(x-1)(x+2)或(x-1)(x+2)共113页矩阵理论A笔记第6页北京航空航天大学张京蕊工程系统工程系验证:(4-D(A+2D=0∴极小式为m(x)-(x-1)(x+2)(3)解法如下引理:A1,A2的极小式为m1(x),m2(x)A10的极小式m(x)等丁m1(x),m2(x)的最小公倍式0A2(此引力可推广到A1,42,43)0100极小式为(x-1)2,0010极小式为(x-1)0取最小公倍式(x-1)2为C的极小式。460(5)F-/40,A1=020|,A00A0123-6101O引理;设D=,则D的极小式m(x)O验证:先证D的性质(右推公式)设A-(an)xn=(a1,2,…,n)则有AD=(0,01,a2,,.m1)AD2=(0,0,∞1,,x12)AD=(0,….0.,a1,,axn)单位向量技巧:∵AI=A(en,e2…,en)=(el,leAen)=A=(a1,a2,.a,)∴Ae1=01,Ae2=(2,.,A→AD=A(0,e1,e2,…,en-1)=(0,a1,a2…,an-)同理AD2=(AD)D=(0,.01,.12)可知:D-1-(D)Dy2-(0.,0,,e1)≠0D"=(D)D1=0,而特征多项式(x)=|x1-D|=x,极小式为某个x共113页矩阵理论A笔记第7页-IMDN开发者社群-imdn.cn">
=0即全体特征值为00,,00由需要P1AP=→PAP+7=1pAP+PP|=P(4+1)P=14+1→A+1=-1注(1)若AB(相似),则AB有相同特征值A,可引入记号:谱集(4)={2,2,…,λ}(全体特征值,含重复)A∽B→o()=o(B)(2)A∽B→1-A=1-B-(2-4元一2)…(-n),特征多项式PAP=B=A-A=p(1-A)P=A-B引理:若A0A2,则M-A|-|M1-4|-1-A1|2-A2→ar(4)=o(A)∪a(42k+1,Ak-2,…n1f(x2)设B,f(x)为多项式,则f(B)=of(,)引理:若n阶方阵A的谱集(4)=1,42,…},则)的全体特社值为)2,…,),x)为多项式Pf:由许尔定理,A∽B→f(4)∽f(B)f(x)的全体特征值为(A1)(42),,()},fx)为多项式例如:4为A的特征值→x为4的特征值。(x)=x)共113页矩阵理论A笔记第3页北京航空航天大学张京蕊工程系统工程系引理:令B,f(x)=x-B|=(x-41)(x-12)….(x-n)则fB)=(B-1D(B-21)…(B-A1D=0Pf:当n=2时,B=0x2f(x)=(x-1)(x-2)000→f(B)-(B-41)(B-21)(2-元)0(00∴得证★Cayley公式:设n阶方阵A的特征多项式为f(x)=|x-A|=a+a1x+…,+x则f4)=anl+a14+…,+4=0Pf:由许尔PAP=B=→P(4)P=fp3P)=f(B)=0(引理)定义若多项式x)使(4)=0,则称(x)为A的个零化式结论方阵A的特征多项式)=1x1-4为A的一个零化式g特征多项式fx)=x2可知:f(A)=A2+1=+I=00-1Hx)=|xI-A|=(x-)(x+i,(i=√-1,t2=-1)f(A)=(A-i)(4+i1=0也可取P=则PPAP=,对角形共113页矩阵理论A笔记第4页北京航空航天大学张京蕊工程系统工程系g:知A则A"=0Onxn由Cayley特征多项式:f(x)=x"→f(4)=4"=0Ex1.A=求P使得PP为对角阵,并验证Cayley定理2.A=cd/,求fx)=x1-4验证f4)-0补充知识(schur公式、Cayley公式)应用由A"=-(a0I+a1A+1A·AanA+a142+…+a.,A把①代入②→Am1=(-)+(+)4+…+(+)41可知:任何和(m≥n)都可写成,4,,A的线性组合任何多项式g(A),可写成lA,…,4的组合。Fg:若A|≠0,fx)=xI-A|=a0+a1x+…+x",ao=|-A|≠则A可用A的多项式表示∵a1A+a242+…+an21A-+A"--a072A(a1+a24+…+an-142+A)Aa1+…+an1A"2+A-1零化式定义:若g(x)=b+b1x+…+bnx,使得g(4)=bn+b14+…+bn4m=0,称g(x)为方阵A的零化式注:方阵A的零化式有无穷多个∴取特征多项式x)则4)=0任取式M(x),f(A(4)=0→f(x)(x)也是零化式极小式定义:在方阵A的零化式集合中,去次数最小的且首项系数为1的零化式m(x),称它为A的极小式共113页矩阵理论A笔记第5页北京航空航天大学张京蕊工程系统工程系注:极小式唯一性质:①极小式m(x)必为特征多项式fx)=|xI-A的因式。②特征多项式fx)=|x1-A的每个单因子(x-4)也是极小式的因子)f(x)=|x1-4=(x-x)(x-2)则极小式m(x)=(x-x)(x-2)y…(x-,),且1≤l1≤m1,1≤l2≤m2,…,1≤l≤n,41,A2…,n互不相同210EgA=020,B=020,求极小式mA(),m()解:(1)|xI-A|=(x-2)(x-1)极小式为:(x-2)(x-1)或(x-2)(x-1)计算:(4-2/)4-1)=000010k≠000000∴极小式为m4(x)=(x-2)2(x-1)(2)|-B|-=(x-2)2(x-1)00000计算:(B-2)B-1)=000010=000-1八000∴极小式为m(x)=(x-2)(x-1)Eg求下列极小式m(x)4604-60(1)A=-3-50,(2)B=2-303-6100210(3)C,(4)D=000010002000解:(1)特征多项式|x7-A|-(x-1)(x+2)极小式为:(x-1)(x+2)或(x-1)(x+2)共113页矩阵理论A笔记第6页北京航空航天大学张京蕊工程系统工程系验证:(4-D(A+2D=0∴极小式为m(x)-(x-1)(x+2)(3)解法如下引理:A1,A2的极小式为m1(x),m2(x)A10的极小式m(x)等丁m1(x),m2(x)的最小公倍式0A2(此引力可推广到A1,42,43)0100极小式为(x-1)2,0010极小式为(x-1)0取最小公倍式(x-1)2为C的极小式。460(5)F-/40,A1=020|,A00A0123-6101O引理;设D=,则D的极小式m(x)O验证:先证D的性质(右推公式)设A-(an)xn=(a1,2,…,n)则有AD=(0,01,a2,,.m1)AD2=(0,0,∞1,,x12)AD=(0,….0.,a1,,axn)单位向量技巧:∵AI=A(en,e2…,en)=(el,leAen)=A=(a1,a2,.a,)∴Ae1=01,Ae2=(2,.,A→AD=A(0,e1,e2,…,en-1)=(0,a1,a2…,an-)同理AD2=(AD)D=(0,.01,.12)可知:D-1-(D)Dy2-(0.,0,,e1)≠0D"=(D)D1=0,而特征多项式(x)=|x1-D|=x,极小式为某个x共113页矩阵理论A笔记第7页 - IMDN开发者社群-imdn.cn">
于 2020-12-09 发布
0 506
北京航空航天大学矩阵理论学习笔记,总结版,学霸总结,可以放心下载使用北京航空航天大学张京蕊工程系统工程系月录§0补充公式§1 Jordan(约当)标准形(简介)§2线性变换与矩阵.24§3欧式空间与QR分解.48§4常用矩阵分解●鲁D●●·,,,,,74§5范数与级数.81§6广义逆A..97§7直积拉直及应用105矩阵理论A笔记北京航空航天大学张京蕊工程系统工程系§0补充公式令A=(a)mxn∈C",风x)=4o+a1x0定义f(4)=a0+a1A+…+amAm,其中I=l若g(x)=bo+b1x+…+bkx,(x)g(x)=g(x)(x),则f(4)“g(A)=g(A)f(A)分块公式A10令A,A1,A2为方阵00 A(2)f(A),fx)为多项式令A=,A1,4为方阵AO(2)f(4)相似关系:A∽B,(PAP=B)则:(1)(P1AP)=P!AP,(k=0,1,2,(2)f(PAP)=PfA)P,f(x)为多项式许尔公式( schur):每个复方阼,A-(a)nxm都相似丁上三角形。共113页矩阵理论A笔记第1页北京航空航天大学张京蕊工程系统工程系即:P-1AP=其中41,,的次序可以任意指定Pf:用归纳法n=1时成立可以设为(n=1阶方阵成立对于n阶方阵A=(an)2×n设特征值为A,…,n取为对应的特征向量,记为a1≠0,A1=1ax1把a1扩展为可逆方阵Q=(a1,02,xn)22e又:g(a,a,…,.)=(Qa,Qba2,,Qan)其中Qe1,aQ0Q4=QA(a1a2,…an)2-I(Aa,,AAQ=(Qa,、+)…,(*)其中A1为(n-1所阶0人:0 A为由假设,对于A1必有(n-1)阶P,可推出PAPEg知n阶方阵A,适合A=0,则A+|=1共113页矩阵理论A笔记第2页北京航空航天大学张京蕊工程系统工程系Pf:A=0→任意特征值A=0→>=0即全体特征值为00,,00由需要P1AP=→PAP+7=1pAP+PP|=P(4+1)P=14+1→A+1=-1注(1)若AB(相似),则AB有相同特征值A,可引入记号:谱集(4)={2,2,…,λ}(全体特征值,含重复)A∽B→o()=o(B)(2)A∽B→1-A=1-B-(2-4元一2)…(-n),特征多项式PAP=B=A-A=p(1-A)P=A-B引理:若A0A2,则M-A|-|M1-4|-1-A1|2-A2→ar(4)=o(A)∪a(42k+1,Ak-2,…n1f(x2)设B,f(x)为多项式,则f(B)=o f(,)引理:若n阶方阵A的谱集(4)=1,42,…},则)的全体特社值为)2,…,),x)为多项式Pf:由许尔定理,A∽B→f(4)∽f(B)f(x)的全体特征值为(A1)(42),,()},fx)为多项式例如:4为A的特征值→x为4的特征值。(x)=x)共113页矩阵理论A笔记第3页北京航空航天大学张京蕊工程系统工程系引理:令B,f(x)=x-B|=(x-41)(x-12)….(x-n)则fB)=(B-1D(B-21)…(B-A1D=0Pf:当n=2时,B=0x2f(x)=(x-1)(x-2)000→f(B)-(B-41)(B-21)(2-元)0(00∴得证★ Cayley公式:设n阶方阵A的特征多项式为f(x)=|x-A|=a+a1x+…,+x则f4)=anl+a14+…,+4=0Pf:由许尔PAP=B=→P(4)P=fp3P)=f(B)=0(引理)定义若多项式x)使(4)=0,则称(x)为A的个零化式结论方阵A的特征多项式)=1x1-4为A的一个零化式g特征多项式fx)=x2可知:f(A)=A2+1=+I=00-1Hx)=|xI-A|=(x-)(x+i,(i=√-1,t2=-1)f(A)=(A-i)(4+i1=0也可取P=则PPAP=,对角形共113页矩阵理论A笔记第4页北京航空航天大学张京蕊工程系统工程系g:知A则A"=0Onxn由 Cayley特征多项式:f(x)=x"→f(4)=4"=0Ex 1. A=求P使得PP为对角阵,并验证 Cayley定理2.A=cd/,求fx)=x1-4验证f4)-0补充知识( schur公式、 Cayley公式)应用由A"=-(a0I+a1A+1A·AanA+a142+…+a.,A把①代入②→Am1=(-)+(+)4+…+(+)41可知:任何和(m≥n)都可写成,4,,A的线性组合任何多项式g(A),可写成lA,…,4的组合。Fg:若A|≠0,fx)=xI-A|=a0+a1x+…+x",ao=|-A|≠则A可用A的多项式表示∵a1A+a242+…+an21A-+A"--a072A(a1+a24+…+an-142+A)Aa1+…+an1A"2+A-1零化式定义:若g(x)=b+b1x+…+bnx,使得g(4)=bn+b14+…+bn4m=0,称g(x)为方阵A的零化式注:方阵A的零化式有无穷多个∴取特征多项式x)则4)=0任取式M(x),f(A(4)=0→f(x)(x)也是零化式极小式定义:在方阵A的零化式集合中,去次数最小的且首项系数为1的零化式m(x),称它为A的极小式共113页矩阵理论A笔记第5页北京航空航天大学张京蕊工程系统工程系注:极小式唯一性质:①极小式m(x)必为特征多项式fx)=|xI-A的因式。②特征多项式fx)=|x1-A的每个单因子(x-4)也是极小式的因子)f(x)=|x1-4=(x-x)(x-2)则极小式m(x)=(x-x)(x-2)y…(x-,),且1≤l1≤m1,1≤l2≤m2,…,1≤l≤n,41,A2…,n互不相同210EgA=020,B=020,求极小式mA(),m()解:(1)|xI-A|=(x-2)(x-1)极小式为:(x-2)(x-1)或(x-2)(x-1)计算:(4-2/)4-1)=000010k≠000000∴极小式为m4(x)=(x-2)2(x-1)(2)|-B|-=(x-2)2(x-1)00000计算:(B-2)B-1)=000010=000-1八000∴极小式为m(x)=(x-2)(x-1)Eg求下列极小式m(x)4604-60(1)A=-3-50,(2)B=2-303-6100210(3)C,(4)D=000010002000解:(1)特征多项式|x7-A|-(x-1)(x+2)极小式为:(x-1)(x+2)或(x-1)(x+2)共113页矩阵理论A笔记第6页北京航空航天大学张京蕊工程系统工程系验证:(4-D(A+2D=0∴极小式为m(x)-(x-1)(x+2)(3)解法如下引理:A1,A2的极小式为m1(x),m2(x)A10的极小式m(x)等丁m1(x),m2(x)的最小公倍式0A2(此引力可推广到A1,42,43)0100极小式为(x-1)2,0010极小式为(x-1)0取最小公倍式(x-1)2为C的极小式。460(5)F-/40,A1=020|,A00 A0123-6101O引理;设D=,则D的极小式m(x)O验证:先证D的性质(右推公式)设A-(an)xn=(a1,2,…,n)则有AD=(0,01,a2,,.m1)AD2=(0,0,∞1,,x12)AD=(0,….0.,a1,,axn)单位向量技巧:∵AI=A(en,e2…,en)=(el,leAen)=A=(a1, a2,. a,)∴Ae1=01,Ae2=(2,.,A→AD=A(0,e1,e2,…,en-1)=(0,a1,a2…,an-)同理AD2=(AD)D=(0,.01,.12)可知:D-1-(D)Dy2-(0.,0,,e1)≠0D"=(D)D1=0,而特征多项式(x)=|x1-D|=x,极小式为某个x共113页矩阵理论A笔记第7页
下载说明:请别用迅雷下载,失败请重下,重下不扣分!
-
小波包分解程序
最佳小波包分解程序
- 2020-11-29下载
- 积分:1
-
上位机(VB6.0)控制步进电机
【实例简介】实现功能:利用PC控制步进电机转动。可以控制转动方向、转动速度以及读取转动角度。上位机VB6.0编写,下位机ATMEGA16程序,均为源码。
- 2021-11-18 00:47:55下载
- 积分:1
-
SIFT tutorial(matlab版本的sift源码和详细介绍)
matlab版本的源码,有详细的sift特征提取和匹配的步骤介绍,具体sift是什么我就不多说了。这个matlab版本的源码就是供初学者学习之用
- 2020-12-05下载
- 积分:1
-
dsp tms320c 6713 的 例程源码
dsp 6713 的 例程源码 seed dec 6713
- 2020-12-04下载
- 积分:1
-
汽车噪声与振动-理论与应用
汽车NVH领域的经典入门读物,引领进入汽车振动噪声研究的最佳选择。第一章管道声学(1.12)时,声压幅值达到最大,反节点的位子是:(1.13)驻波是由频率相同的向右传播的入射波和向左传播的反射波迭加而成。驻波并不是运动的波,而是静止的,这是“驻”波名称的由来。波表示管道中的声音的模态。对於长度一定的管道来说,由于有许多频率的波,因此也就有很多驻波。这里所提到的驻波是假设管壁刚硬,所冇声波遇到管壁时全部被反射回来。可是实际上,管端壁不是完全刚性,因此反射波的声压不完全等於入射波声压,因此在节点处,入射波和反射波不可能完全抵消。但是这些点处的声压大部分被抵消,声压最低。第二节管道声阻抗阻抗是指当对媒质受到压力或者搾动力时,媒质会对传播产生阻碍。管道中的声学阻抗Z,定义声压与质点体积速度的比值,即(1.14)式中,u,U和S分别是管道中的速度,体积速度和截面积。体积速度与质点速度的关系为:L=SL。声吝在管道內传播,当管道的截面积发生变化的时候,声阻抗也发生变化。图1.3是截面积变化的管道,在变截面的地方,由于阻抗发生变化,一部分入射波就会被发射回原来的管道而另一部分入射波会在新的截面管道中继续传播。抗性消音器的工作原就是基于这种阻抗的变化。声波从发动机出来并在进气或者排气系统中传播,当遇到消音元件或者截面积变化时,入射声波被反射回发动机声源,从而抑制声音的传播。进排气系统中声阻抗不匹配的情况主要有截面积变化,主管道中插入了其他管道(如旁支消音器等),管道开口通往大气等等图1.3截面积变化的管道进排气系统中管道的长度都是有限的。图1.4表小一个长度为L的管道。假改管道两端的声阻抗分别已知,即在=处,声阻抗为,在=处,声阻抗为由公式(1.6)和(1.9),可以得到管道中仟一点的声阻抗为管道声学图1.4长度为L的管道将=代入公式(1.15)中,得到该处的声阻抗为:将三代入公式(1.15)中,得到该处的声阻抗:公式(17)可以重新写成下面的形式18将方程(1.16)代入到方程(1.18中,消除和,就得到输入声阻抗和输出声阻抗的关系,如下第一章管道声学(1.20第三节管口封闭与管口敞廾声波从管道入口端发射出来,传播到尾端。管道尾端通常有两种情况,一种是开口的,如进气管口,排气尾管口;另一种是封口的,如四分之一波长管。下面就来分析这两种尾端的声学特征。1.开∏-封闭管道图1.5表示管道尾端封闭状况。声音在管道里问石传播,当声波碰到刚性的封闭端时,声波被全部反弹冋来,再向左传播管口封闭图1.5开凵封闭管道对一个刚性的封闭口来说,其声阻抗为无穷大,即>0,根据公式(1.19),得到:1.21)声阻抗可以写成下面的形式:(1.22)式中R和粉别是阻抗的实部和虚部,R为声阻,称为声抗。声阻取决于结构的材料特性,而声抗则取决」结构的儿何特性。当声抗为零的时候,结构就发生共振。公式(1.21)中的声阻抗也可以写成公式(1.22)那样的形式,为(1.23)上式如果满足下亩的条件:(1.24)即,那么这个开口-封闭管道就发生共振,其固有频率为:(1.25)当n-=1,2,3,.,时,分别对应著管道第·阶、第二阶、第三阶,,.,等阶次频率图1.6是管道声波的第一阶和第阶模态。这个声波在封闭端时,声压达大最大值,然后发射第一章管道声学到入口处,使得入口端的声压为零,即在开口端形成驻波节点。四分之一波长管就是应用这个原理来工作的。图1.6管道声波的第一阶模态(A)和第二阶模态(B)公式(1.25)可以转变为管道长度与波长的关系,表达如下1.26)当n=1时,管道的长度是波长的四分之,即:。所以这种开∏封闭的管道通常叫著四分之一泼长管2.开口开口答道图1.7为一个尾端开口的管」。声波从入口端向右传播进入开口端时,声音与大气产生声耦合。大气的辐射声阻抗会将一部分声波返回管口敞开图1.7开口-开口管道声波在尾端的声阻抗为周围坏境的声阻抗,也就是说这个声阻抗不为零。为了使问题简化起见,我们先假设这个阻抗为零,然后再对所得到的结构进行修正。如果在x=处的声阻抗为零,那么由公式(1.19)可以得到下式(1.27)同样,当这个声阻抗中的声抗为零的吋候,管道就发生共振,这时必须满足:即:这时,开口-开口管道的共振频率为:当n=1,2,3,,时,分别对应著管道第一阶、第二阶、第三阶,.,等阶次频率第一章管道声学图1.8是开口-开口管道声波的第一阶和第二阶模态。图1.8开口-开口管道的第一阶模态(A)和第二阶模态(B)公式(1.30)可以转变为管道长度与波长的关系,衣达如下(1.31)3.开口管道的修正在推导尾端廾口公式时,我们假设了出口周围坯境的声阻抗为零,但是实际上这个阻抗不为零,因此必须对公式(1.27-1.31)的结论做修止。对图1.9这样的开∏终端,被称为自由自由开口。该开口处的声阻抗为:等效管图1.9自山开口-开口修正管道山于管道的直径非常小,因此和都远远小于1。山公式(1.27)和(1.32)得到:(1.33)这样,管道内的频率为34)管道长度与波长的关系为(1.35这样管道的长度比声阻抗为零的时候要短些,也就是说好像有一根等效的延长管与原来的管道相连接。管道的计算长度就是实际管子长度加上等效延长管长度△即第一章管道声学(1.36有时侯,在出口管处还会加类似与法兰的结构,如图1.10所示。这时,有效延长管的长度为△实际管子的长度为:△式中是管子的计算长度发等效管图1.10法兰开口-开口修正管道第四节四端网终分析进气系统或者排气系统都是有很多管道和消音元件组成。分析整个系统往往是非常复杂的,但是如果将系统分解到一些小的段落,那么分析起来就相对容易些。得到了每个段落或者是每个部件的分析结果,然后将之合成起来就得到了整个系统的结果。四端网络分析就是这种分析方法,在管道声学分析中得到了广泛的应用。对於管道中一小段质量(如图1.11)来说,动力方程可以写成如下:(1.39)式中,S是管道的截面积,是这个小质量段的长度,和分别是质量端两边的压力图1.11管道中一小段质量的受力分析公式(1.39)可以表达为(1.40)第一章管道声学对这一小段质量来说,假设两边的速度是相等的,即将这公式(1.40)和(1.41)写成矩阵形式,得到:(1.42)公式(1.42)建立起这段小质量块两边的压力和速度的关系。管道中小段质量块后端的压力和速度可以用它前端的压力和速度来表示。也就是说质量块后端与前端之间建立起来一种传递关系。同样对一个长度为L的管道(如图1.4所示)也可以得到管道两端的传递关系。在=处的压力和速度可以通过公式(1.6)和(1.9)分别求得(1.13)由以上两式可以得到和,如下:45(1.46根据公式(6)和(9),在处的压力和速度分别为将公式(1.45)和(1.46)中和的表达式代入公式(1.47)和(1.48)之中,就得到管道入冂与出∏之间声压和速度之间的关系,为:+49将公式(1.49)和(1.50)写成如下的矩阵形式第一章管道声学这样就得到了管道两边的压力和速度的传递关系。公式(1.51)可以简单地写成如下形式式中,被称为传递矩阵。如果管道的传递矩阵知道,那么只要知道管道端的压力和速度,就可以通过传递矩阵算出另一端的压力和速度。在传递矩阵两边分别是两个输入参数和两个输出参数。这四个参数的关系由传递矩阵来确定,因此这种表达方式称为四端网络法。上面介绍了小段质量和长度为L的管道的传递矩阵表达方法。这种方法可以推广到任何一个声学元件,其输入端和输出端的声压和速度都可以用四端网络米表示。图1.12代表某个声学元件i。图1.12一个管道元件的四端网终图这个元件两边的压力和速度关系为式中是传递矩阵,是传递矩阵系数。汽车的进气系统包括进气管道、空气过滤器、赫耳姆兹消音器、四分之波长管等。排气系统包括排气多支管、催化器、谐振器、消音器和管道等。一个系统如果由N个元件组成。而且每个元件的传递矩阵都知道,那么出声口的声压和速度就可以用声源的声压和速度来表示如下形式:(1.54)式中的L1是系统的传递矩阵,如下形式(1.55)
- 2020-12-05下载
- 积分:1
-
生日快乐 送女友源代码 小惊喜
在异性好友过生日的时候,特别制作了这个,全屏都在下蛋糕的小程序,蛋糕下降的过程中伴随生日快乐歌曲,文字闪烁标示出主题,绝对值得,自己原创,现将源码分享给大家,希望大家喜欢。
- 2020-11-29下载
- 积分:1
-
《MySQL实用教程》郑阿奇著实验报告完整版
的要求:(1)掌握MySQL服务器安装方法(2)掌握MySQL Administrator的基本使用方法(3)基本了解数据库及其对象实验准备:(1)了解MySQL安装的软硬件要求(2)了解MYSQL支持的身份验证模式(3)了解MySQL各组件的主要功能(4)基本了解数据库、表、数据库对象实验内容:1.安装MySQL服务器和MySQL界面工具安装MySQL界面工具。(插入安装好的界面工具截图即可)
- 2020-12-05下载
- 积分:1
-
DFA(Detrended fluctuation analysis) matlab 程序
DFA 是检测时间序列长程相关的范围的有效方法.进入包含着三个文件的目录,在Matlab命令窗口,直接输入runDFA,就可以得到结果。
- 2020-12-02下载
- 积分:1
-
Python实现循环神经网络RNN
基于Python的循环神经网络(RNN)实现
- 2020-11-28下载
- 积分:1
-
gpops工具箱+例子(详细注释)+官方手册+安装方法
gpops工具箱+例子(详细注释)+官方手册+安装方法
- 2020-12-06下载
- 积分:1