数值分析计算实习作业一
11页1、数 值 分 析计算实习题一学号:姓名: 院系: 2015年11月5日数值分析计算实习作业一一、分析1.1算法分析题目要求求出:1)特征值从小到大排列的最小特征值和最大特征值。2)特征值中模最小的特征值。3)靠近一组数的一组特征值。4)矩阵A的条件数cond(A)2。5)行列式detA。解决方法:1)若将所有行列式按模的大小排列则模最大的特征值一定是和中的一个,因此利用幂法求出模最大的特征值。然后利用带原点平移的幂法,将系数矩阵变为即将所有特征值都减去,则特征值按大小顺序排列的次序不变,模最大的特征值依然在整个排列的两端,再用一次幂法得到模最大的特征值,其中为带原点平移的幂法求出的特征值,最后两个特征值、比较大小,大的为,小的为。2)因为为按模最小的特征值,因此用反幂法可求的其特征值。3)因为靠近数,因此一定是所有的中模最小的,因此可利用带原点平移的反幂法求出特征值,此时的系数矩阵变为。4)条件数cond(A)2为模最小的特征值与模最大的特征值的比的绝对值,因此利用1和2中求出的和可解出条件数。5)可对矩阵A进行LU分解,即则,又因为矩阵L对角线元素为1,则=1,所以,U为上三角阵,行列
2、式为对角线元素的乘积,因此可得A的行列式。1.2程序分析1.2.1 因为A为拟三角阵,储存时零元素不储存,因此将矩阵A压缩为5*501的矩阵CA的带内元素=C中的元素程序中A5501即为压缩后的矩阵。1.1.2 程序中的B5501为过渡矩阵,在幂法迭代、反幂法迭代以及LU分解中均用矩阵B来计算,计算之间对B进行适当的赋值。先令B=A,调用幂法函数求;再令B=A-I,调用幂法函数求出;再令B=A,调用反幂法函数求;在令B=,调用反幂法函数求,最后令B=A,将B进行LU分解,计算A的行列式。1.2.3幂法求解过程:1)取初始向量u0=1,1,1;2)进行k次迭代,k=1,2,33)对每一次迭代,计算上一次迭代中的的模(2范数),将单位化后赋值给,即,计算矩阵B与向量的乘积,计算与的内积。4)如果两次相邻两次迭代的满足:(允许误差),则结束迭代,的值可认为是矩阵B对应的模最大的特征值。如果不满足误差条件,则重复3)的迭代直到达到误差允许值。1.2.4 反幂法求解过程:1)取初始向量u0=1,1,1;对矩阵B进行LU分解。2)进行k次迭代,k=1,2,33)对每一次迭代,计算上一次迭代中的的模
3、(2范数),将单位化后赋值给,即,利用LU分解后的矩阵B,求线性方程组,计算与的内积。4)如果两次相邻两次迭代的满足:(允许误差),则结束迭代,的值可认为是矩阵B对应的模最大的特征值。如果不满足误差条件,则重复3)的迭代直到达到误差允许值。1.2.5 拟三角矩阵的LU分解过程可参照数值分析书26页的算法得到。二、程序整个数值分析在c语言中的程序如下:#include #include double mifa();double fanmifa();void fenjie();/声明三个函数,mifa是用幂法求模最大的特征值,fanmifa是用反幂法求模最小的特征值,fenjie是对系数矩阵做LU分解double B5501=0;double A5501=0;/定义两个系数矩阵为全局变量,B作为求解过程中的系数矩阵使用double sigma=1e-12;/定义误差允许值/*主程序*/void main() int i,j,t; double lamk,lamm,min,max,miu,lam,x,lams; double cond2; double det=1; for(i=0;i=50
《数值分析计算实习作业一》由会员suns****4568分享,可在线阅读,更多相关《数值分析计算实习作业一》请在金锄头文库上搜索。
土地管理与地籍测量---第八章界址点测量
人机工程学案例分析(2)
工程安全培训_201303
第9章房地产投资决策分析
第2章房地产经纪制度
ACM程序设计-东北林业大学acm05
《亲爱的汉修先生》读书交流会
中原_深圳新世界尖岗山项目市场汇报_40P_2012年_别墅_项目分析_量价走势
五年级数学质量分析演示文稿
人工智能小镇-智慧小镇建设20180525
景观基本知识及发展历程
建设工程信息管理(2)
机电驱动技术第二章步进驱动技术
工程力学-第9章圆轴扭转时的应力变形分析与强度刚度设计
第一章第二节幼儿园文化环境建设的原则
第一章检测技术的基础知识
第一章__现代表面工程技术
第六章钢结构工程
第9节项目试运行管理
班主任工作经验交流课件(4)
2023-02-24 13页
2023-01-08 279页
2023-12-29 4页
2023-02-05 14页
2023-11-19 3页
2023-06-29 3页
2022-09-22 9页
2024-02-17 14页
2023-01-19 8页
2022-10-22 7页