您所在位置: 网站首页 / 北航数值分析大作业一.docx / 文档详情
北航数值分析大作业一.docx 立即下载
2025-08-26
约2.4万字
约52页
0
34KB
举报 版权申诉
预览加载中,请您耐心等待几秒...

北航数值分析大作业一.docx

北航数值分析大作业一.docx

预览

免费试读已结束,剩余 47 页请下载文档后查看

10 金币

下载文档

如果您无法下载资料,请参考说明:

1、部分资料下载需要金币,请确保您的账户上有足够的金币

2、已购买过的文档,再次下载不重复扣费

3、资料包下载后请先用软件解压,在使用对应软件打开

北航数值分析大作业一

第一篇:北航数值分析大作业一北京航空航天大学数值分析大作业一学院名称自动化专业方向控制工程学号ZY1403140学生姓名许阳教师孙玉泉日期2014年11月26日设有501501的实对称矩阵A,a1bcbAccbcba501其中,ai(1.640.024i)sin(0.2i)0.64e(i1,2,,501),b0.16,c0.064。矩阵A的特征值为i(i1,2,,501),并且有0.1i12501,|s|min|i|1i5011.求1,501和s的值。2.求A的与数k1k501140最接近的特征值ik(k1,2,,39)。3.求A的(谱范数)条件数cond(A)2和行列式detA。一方案设计求1,501和s的值。s为按模最小特征值,|s|min|i|。可使用反幂法求得。1i5011,501分别为最大特征值及最小特征值。可使用幂法求出按模最大特征值,如结果为正,即为501,结果为负,则为1。使用位移的方式求得另一特征值即可。求A的与数k1k501140最接近的特征值ik(k1,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)yk1uk1/max|uk1|uAyk1kksgn(max|uk1|)max|uk1|(2-1)终止迭代的控制理论使用|kk1|/|k|,实际使用||k||k1||/|k|(2-2)由于不保存A矩阵中的零元素,只保存主对角元素a[501]及b,c值。则上式中ukAyk1简化为: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(501u(i)cy(i2)by(i1)a(i)y(i)by(i1)cy(i2)(i3,,499)(2-3)反幂法使用如下迭代格式:(0)T任取非零向量u0(h1(0),,hn)yk1uk1/max|uk1|-1uAyk1kksgn(max|uk1|)max|uk1|(2-4)其中ukA1yk1Aukyk1,解方程求出uk。求解过程中使用LU分解,由于A为5对角矩阵,选择追赶法求取LU分解。求解过程如下:LUukyk1Lxkyk1Uukxkuk追赶法求LU分解的实现:a1bcbAccLUbcba501p1r21t1q1z3q499t500z501r501p5011由上式推出分解公式如下:p1a1,t1b/a1r2b,p2a2r2t1qic/pi,i1,...,499ti(briqi1)/pi,i2,...,500zic,i3,...,501ribcti2,i3,...,501piaicqi2riti1,i3,...,501推导出回代求解公式如下:x1y1/p1x2(y2r2x1)/p2xi(yizixi2rixi1)/pi,i3,...,501(2-5)(2-6)(2-7)u501x501u500x500t500u501uxtuqx,i499,...,1iii1ii2i(2-8)cond(A)2及A行列式求解cond(A)2|1|(2-9)|s|由式(2-5)可得:501detApii1三源程序#include#includedoubleep=1e-12,b=0.16,c=-0.064;intj=0;doublepower(doublea[501]);//幂法doubleinv
查看更多
单篇购买
VIP会员(1亿+VIP文档免费下)

扫码即表示接受《下载须知》

北航数值分析大作业一

文档大小:34KB

限时特价:扫码查看

• 请登录后再进行扫码购买
• 使用微信/支付宝扫码注册及付费下载,详阅 用户协议 隐私政策
• 如已在其他页面进行付款,请刷新当前页面重试
• 付费购买成功后,此文档可永久免费下载
全场最划算
12个月
199.0
¥360.0
限时特惠
3个月
69.9
¥90.0
新人专享
1个月
19.9
¥30.0
24个月
398.0
¥720.0
6个月会员
139.9
¥180.0

6亿VIP文档任选,共次下载特权。

已优惠

微信/支付宝扫码完成支付,可开具发票

VIP尽享专属权益

VIP文档免费下载

赠送VIP文档免费下载次数

阅读免打扰

去除文档详情页间广告

专属身份标识

尊贵的VIP专属身份标识

高级客服

一对一高级客服服务

多端互通

电脑端/手机端权益通用