《数学实验方法》课程实验报告.docVIP

本科生实验报告 实验课程 数学实验方法 学院名称 管理科学学院 专业名称 信息与计算科学 学生姓名 学生学号 指导教师 王茂芝 实验地点 6C402 实验成绩 二〇 一五 年 九 月 24 二〇 一五 年 十 月 填写说明 适用于本科生所有的实验报告（印制实验报告册除外）； 专业填写为专业全称，有专业方向的用小括号标明； 格式要求： 用A4纸双面打印（封面双面打印）或在A4大小纸上用蓝黑色水笔书写。 打印排版：正文用宋体小四号，1.5倍行距，页边距采取默认形式（上下2.54cm，左右2.54cm，页眉1.5cm，页脚1.75cm）。字符间距为默认值（缩放100%，间距：标准）；页码用小五号字底端居中。 具体要求： 题目（二号黑体居中）； 摘要（“摘要”二字用小二号黑体居中，隔行书写摘要的文字部分，小4号宋体）； 关键词（隔行顶格书写“关键词”三字，提炼3-5个关键词，用分号隔开，小4号黑体)； 正文部分采用三级标题； 第1章 ××(小二号黑体居中，段前0.5行) 1.1 ×××××小三号黑体×××××（段前、段后0.5行） 1.1.1小四号黑体（段前、段后0.5行） 参考文献（黑体小二号居中，段前0.5行），参考文献用五号宋体，参照《参考文献著录规则（GB/T 7714－2005）》。  实验一 矩阵(方阵)特征值和特征向量计算 1 实验内容 （1）熟悉幂法计算方阵特征值和特征向量； （2）利用C语言实现，Matlab对计算结果进行验证； （3）以AHP中的正互反矩阵作为测试数据 2 数据结构与算法描述 1）变量及函数的定义 变量/函数名 类 型 说 明 A[N][N] Double[][] 用于计算的矩阵 v0[N] Double[] 初始向量 max_element double 最终的特征值 2）程序流程图 3 实验数据与实验结果（可用文字描述或贴图的方式进行说明） 1）测试数据 矩阵A[N][N]= 1 1 0.5 1 1 0.25 0.5 0.25 2.0 初始特征向量v0[N]={1,1,1}; （2）实验结果 4 程序代码清单（可直接将可运行源代码粘贴在下面的方框中） #includestdio.h #includemath.h #define N 3 void show(double v1[]) { printf((); for(int i=0;iN;i++) printf(%f ,v1[i]); printf()\n); } void main() { double A[N][N]={{1,1,0.5},{1,1,0.25},{0.5,0.25,2.0}}; double v0[N]={1,1,1}; double v1[N]={0,0,0}; double dx=0; double max_element=0; int i,j,k; printf(迭代过程:\n); printf(==============================\n\n); for(k=0;k=100;k++) { for(i=0;iN;i++) { v1[i]=0; for(j=0;jN;j++) v1[i]+=A[i][j]*v0[j]; } max_element=0; for(i=0;iN;i++) if(fabs(v1[i])=max_element) max_element=v1[i]; for(i=0;iN;i++) v1[i]=v1[i]/max_element; dx=0; for(i=0;iN;i++) dx+=fabs(v1[i]-v0[i]); for(i=0;iN;i++) v0[i]=v1[i]; show(v1); if(dx0.001) break; } printf(\n==============================); if(k==300) { printf(没求出); } printf(\n特征向量是：); show(v1); printf(\n特征值=%f\n,max_element); } 实验二 Matlab外部编程接口 1 实验内容 利用计算引擎技术通过C调用Matlab实现下述画图效果。 提示：在matlab中用到的命令包括： [X,Y] = meshgrid(-2:.2: