




如果您无法下载资料,请参考说明:
1、部分资料下载需要金币,请确保您的账户上有足够的金币
2、已购买过的文档,再次下载不重复扣费
3、资料包下载后请先用软件解压,在使用对应软件打开
北航数值分析大作业一 第一篇:北航数值分析大作业一北京航空航天大学数值分析大作业一学院名称自动化专业方向控制工程学号ZY1403140学生姓名许阳教师孙玉泉日期2014年11月26日设有501501的实对称矩阵A,a1bcbAccbcba501其中,ai(1.640.024i)sin(0.2i)0.64e(i1,2,,501),b0.16,c0.064。矩阵A的特征值为i(i1,2,,501),并且有0.1i12501,|s|min|i|1i5011.求1,501和s的值。2.求A的与数k1k501140最接近的特征值ik(k1,2,,39)。3.求A的(谱范数)条件数cond(A)2和行列式detA。一方案设计求1,501和s的值。s为按模最小特征值,|s|min|i|。可使用反幂法求得。1i5011,501分别为最大特征值及最小特征值。可使用幂法求出按模最大特征值,如结果为正,即为501,结果为负,则为1。使用位移的方式求得另一特征值即可。求A的与数k1k501140最接近的特征值ik(k1,2,...,39)。题目可看成求以k为偏移量后,按模最小的特征值。即以k为偏移量做位移,使用反幂法求出按模最小特征值后,加上k,即为所求。求A的(谱范数)条件数cond(A)2和行列式detA。矩阵A为非奇异对称矩阵,可知,cond(A)2|max|min(1-1)其中max为按模最大特征值,min为按模最小特征值。detA可由LU分解得到。因LU均为三角阵,则其主对角线乘积即为A的行列式。二算法实现幂法使用如下迭代格式:(0)T任取非零向量u0(u1(0),,un)yk1uk1/max|uk1|uAyk1kksgn(max|uk1|)max|uk1|(2-1)终止迭代的控制理论使用|kk1|/|k|,实际使用||k||k1||/|k|(2-2)由于不保存A矩阵中的零元素,只保存主对角元素a[501]及b,c值。则上式中ukAyk1简化为:u(1)a(1)y(1)by(2)cy(3)u(2)by(1)a(2)y(2)by(3)cy(4))u(500)cy(498)by(499)a(500)y(500)by(501)cy(499)by(500)a(501)y(501)u(501u(i)cy(i2)by(i1)a(i)y(i)by(i1)cy(i2)(i3,,499)(2-3)反幂法使用如下迭代格式:(0)T任取非零向量u0(h1(0),,hn)yk1uk1/max|uk1|-1uAyk1kksgn(max|uk1|)max|uk1|(2-4)其中ukA1yk1Aukyk1,解方程求出uk。求解过程中使用LU分解,由于A为5对角矩阵,选择追赶法求取LU分解。求解过程如下:LUukyk1Lxkyk1Uukxkuk追赶法求LU分解的实现:a1bcbAccLUbcba501p1r21t1q1z3q499t500z501r501p5011由上式推出分解公式如下:p1a1,t1b/a1r2b,p2a2r2t1qic/pi,i1,...,499ti(briqi1)/pi,i2,...,500zic,i3,...,501ribcti2,i3,...,501piaicqi2riti1,i3,...,501推导出回代求解公式如下:x1y1/p1x2(y2r2x1)/p2xi(yizixi2rixi1)/pi,i3,...,501(2-5)(2-6)(2-7)u501x501u500x500t500u501uxtuqx,i499,...,1iii1ii2i(2-8)cond(A)2及A行列式求解cond(A)2|1|(2-9)|s|由式(2-5)可得:501detApii1三源程序#include#includedoubleep=1e-12,b=0.16,c=-0.064;intj=0;doublepower(doublea[501]);//幂法doubleinv

猫巷****傲柏
实名认证
内容提供者


最近下载
一种基于双轨缆道的牵引式雷达波在线测流系统.pdf
一种基于双轨缆道的牵引式雷达波在线测流系统.pdf
一种胃肠道超声检查助显剂及其制备方法.pdf
201651206021+莫武林+浅析在互联网时代下酒店的营销策略——以湛江民大喜来登酒店为例.doc
201651206021+莫武林+浅析在互联网时代下酒店的营销策略——以湛江民大喜来登酒店为例.doc
用于空间热电转换的耐高温涡轮发电机转子及其装配方法.pdf
用于空间热电转换的耐高温涡轮发电机转子及其装配方法.pdf
用于空间热电转换的耐高温涡轮发电机转子及其装配方法.pdf
用于空间热电转换的耐高温涡轮发电机转子及其装配方法.pdf
用于空间热电转换的耐高温涡轮发电机转子及其装配方法.pdf