对称三对角矩阵特征值的二分法.doc_第1页
对称三对角矩阵特征值的二分法.doc_第2页
对称三对角矩阵特征值的二分法.doc_第3页
对称三对角矩阵特征值的二分法.doc_第4页
对称三对角矩阵特征值的二分法.doc_第5页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

数学软件实验任务书课程名称数学软件实验班级数0901实验课题对称三对角矩阵特征值的二分法实验目的熟悉对称三对角矩阵特征值的二分法实验要求运用Matlab/C/C+/Java/Maple/Mathematica等其中一种语言完成实验内容对称三对角矩阵特征值的二分法成绩教师实验1对称三对角矩阵特征值的二分法1 实验原理对于对称三对角矩阵：设，记特征矩阵的左上角的阶子式为，设，利用行列式的展开式，可得的递推公式：为的特征多项式，个零点为矩阵的个特征值2 实验程序程序：int ebstq(int n,double b,double c,double q,double eps,int l)int i,j,k,m,it,u,v; double d,f,h,g,p,r,e,s; cn-1=0.0;d=0.0;f=0.0; for (j=0; jd) d=h; m=j; while (md) m=m+1; if (m!=j) do if (it=l) printf(failn); return(-1); it=it+1; g=bj; p=(bj+1-g)/(2.0*cj); r=sqrt(p*p+1.0); if (p=0.0) bj=cj/(p+r); else bj=cj/(p-r); h=g-bj; for (i=j+1; i=j; i-) g=e*c; h=e*p; if (fabs(p)=fabs(c) e=c/p; r=sqrt(e*e+1.0); ci+1=s*p*r; s=e/r; e=1.0/r; else e=p/c; r=sqrt(e*e+1.0); ci+1=s*c*r; s=1.0/r; e=e/r; p=e*b-s*g; bi+1=h+s*(e*g+s*b); for (k=0; kd); bj=bj+f; for (i=0; i=n-1; i+) k=i; p=b; if (i+1=n-1) j=i+1; while (j=n-1)&(bj=p) k=j; p=bj; j=j+1; if (k!=i) bk=b; b=p; for (j=0; j=n-1; j+
人人文库> 全部分类> 教育资料 > 课件下载

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论