雅克比法求矩阵特征值特征向量.docVIP

C语言课程设计报告 课程名称： 计算机综合课程设计 学 院： 土木工程学院 设计题目： 矩阵特征值分解 级 别： B 学生姓名： 学 号： 同组学生： 无 学 号： 无 指导教师： 2012年 9 月 5 日 C语言课程设计任务书 （以下要求需写入设计报告书） 学生选题说明: 以所发课程设计要求为准，请同学们仔细阅读； 本任务书提供的设计案例仅供选题参考；也可自选，但难易程度需难度相当； 鼓励结合本专业（土木工程、力学）知识进行选题，编制程序解决专业实际问题。 限2人选的题目可由1-2人完成（A级）；限1人选的题目只能由1人单独完成（B级）； 设计总体要求: 采用模块化程序设计； 鼓励可视化编程； 源程序中应有足够的注释； 学生可自行增加新功能模块（视情况可另外加分）； 必须上机调试通过； 注重算法运用，优化存储效率与运算效率； 需提交源程序(含有注释)及相关文件(数据或数据库文件)； (cpp文件、txt或dat文件等) 提交设计报告书，具体要求见以下说明。 设计报告格式： 目录 课程设计任务书（功能简介、课程设计要求）； 系统设计（包括总体结构、模块、功能等，辅以程序设计组成框图、流程图解释）； 模块设计（主要模块功能、源代码、注释（如函数功能、入口及出口参数说明，函数调用关系描述等）； 调试及测试：（调试方法，测试结果的分析与讨论，截屏、正确性分析）； 设计总结：（编程中遇到的问题及解决方法）； 心得体会及致谢； 参考文献  1.课程设计任务书 功能简介： 输入一个对称正方矩阵A，从文本文件读入； 对矩阵A进行特征值分解，将分解结果：即U矩阵、S矩阵输出至文本文件； 将最小特征值及对应的特征向量输出至文本文件； 验证其分解结果是否正确。 提示：A=USUT,具体算法可参考相关文献。 功能说明： 矩阵特征值分解被广泛运用于土木工程问题的数值计算中，如可用于计算结构自振频率与自振周期、结构特征屈曲问题等。 注：以三阶对称矩阵为例 2.系统设计 模块设计 #includestdio.h #includestdlib.h #includemath.h int main() { FILE *fp; int tezheng(double *a,int n,double *s,double *u,double eps,int itmax); //函数调用声明 int i,j,p,itmax=1000; //itmax为最大循环次数 double eps=1e-7,s[3][3],u[3][3]; //eps为元素精度，s为对角矩阵S，u为矩阵U double a[9];//a为待分解矩阵A i=tezheng(a,3,s,u,eps,1000); if(i0) //i对应函数中的返回值it { if((fp=fopen(juzhen.txt,w))==NULL) //打开待输入txt文件 { printf(无法打开文件.\n); return; } printf(U矩阵为:\n); //下几句分别向屏幕和txt文件输入矩阵U fprintf(fp,U矩阵为:\n); for(i=0;i3;i++) { for(j=0;j3;j++) { printf(%10.6f,u[i][j]); fprintf(fp,%10.6f,u[i][j]); } printf(\n); fprintf(fp,\n); } printf(S对角矩阵为:\n); //下几句分别向屏幕和txt文件输入矩阵S fprintf(fp,S对角矩阵为:\n); for(i=0;i3;i++) { for(j=0;j3;j++) { printf(