重庆大学--数学实验-温罗生

重重 庆庆 大大 学学 学学 生生 实实 验验 报报 告告 实验课程名称实验课程名称 数学实验数学实验 开课实验室开课实验室 DS1404 学学 院院 2014 年级年级 电子科学与技术电子科学与技术 专业班专业班 01 学学 生生 姓姓 名名 学学 号号 开开 课课 时时 间间 2016 学年第学年第 2 学期学期 总总 成成 绩绩 教师签名教师签名 数数 学学 与与 统统 计计 学学 院院 制制 开课学院 实验室 开课学院 实验室 DS1404 实验时间实验时间 2016 年年 3 月月 10 日日 实验项目类型实验项目类型 课程课程 名称名称 数学实验数学实验 实验项目实验项目 名名 称称 验证验证演示演示综合综合设计设计其他其他 指导指导 教师教师 温老师温老师成成 绩绩 1 代码 x 3 14 y1 sign x y2 x 3 y3 sin x y4 cos x y5 tan x 2 y6 2 tan x y7 3 x y8 log abs x x x y9 log2 x log2 3 y10 fix x y11 asin x 2 3 pi y12 x 1 6 y13 sqrt 2 x 2 4 x log x exp x sin x 2 x 2 y y1 y2 y3 y4 y5 y6 y7 y8 y9 y10 y11 y12 y13 y 结果 y 1 0000 30 9591 0 0016 1 0000 0 4645 0 0032 31 4891 3 5025 1 0415 3 0000 0 4498 1 2101 4 6822 2 代码 y rand 5 10 max y max y 2 min y max y 2 sum y sort y length y size y 说明 rand x 随机产生矩阵 max y 求取矩阵 y 每一行最大值 max y 2 求取矩阵 y 每一 列最大值 其中 min y 类似 sum y 相应的与 max y 相似的用法 表示求矩阵每一行或 每一列的和 这几个函数在用法上大体一致 只是实现的功能有所不同而已 结果 y 0 4229 0 6999 0 5309 0 9686 0 7788 0 2810 0 5181 0 6761 0 2548 0 7805 0 0942 0 6385 0 6544 0 5313 0 4235 0 4401 0 9436 0 2891 0 2240 0 6753 0 5985 0 0336 0 4076 0 3251 0 0908 0 5271 0 6377 0 6718 0 6678 0 0067 0 4709 0 0688 0 8200 0 1056 0 2665 0 4574 0 9577 0 6951 0 8444 0 6022 0 6959 0 3196 0 7184 0 6110 0 1537 0 8754 0 2407 0 0680 0 3445 0 3868 ans 0 6959 0 6999 0 8200 0 9686 0 7788 0 8754 0 9577 0 6951 0 8444 0 7805 ans 0 9686 0 9436 0 6718 0 9577 0 8754 ans 0 0942 0 0336 0 4076 0 1056 0 0908 0 2810 0 2407 0 0680 0 2240 0 0067 ans 0 9686 0 9436 0 6718 0 9577 0 8754 ans 2 2825 1 7604 3 1313 2 5417 1 7132 2 5810 3 2978 2 4001 2 3355 2 4515 ans 0 0942 0 0336 0 4076 0 1056 0 0908 0 2810 0 2407 0 0680 0 2240 0 0067 0 4229 0 0688 0 5309 0 3251 0 1537 0 4401 0 5181 0 2891 0 2548 0 3868 0 4709 0 3196 0 6544 0 5313 0 2665 0 4574 0 6377 0 6718 0 3445 0 6022 0 5985 0 6385 0 7184 0 6110 0 4235 0 5271 0 9436 0 6761 0 6678 0 6753 0 6959 0 6999 0 8200 0 9686 0 7788 0 8754 0 9577 0 6951 0 8444 0 7805 ans 10 3 代码 n 阶行列式 function u v DET x y A rand x u det A v 0 for i 1 x a A a y a i v A y i det a 1 y i v end end 下面任意选取 n 阶行列式进行验证 例如四阶行列式 结果 u v DET 4 2 u 0 0261 v 0 0261 4 代码 x 10 0 1 10 y exp x 2 2 4 subplot 2 2 1 plot x y theta pi 0 1 pi rou sin 3 theta subplot 2 2 2 plot theta rou t 10 0 1 10 y1 3 t 1 t 3 y2 3 t 2 1 t 3 subplot 2 2 3 plot y1 y2 subplot 2 2 4 ezplot 9 x 2 2 4 y 6 2 36 6 2 4 10 图形 5 代码 x y meshgrid 5 0 1 5 z sin pi sqrt x 2 4 y 2 mesh x y z 图形 6 代码 H n 1 for a 1 9 for b 0 9 for c 0 9 h 100 a 10 b c j a 3 b 3 c 3 if h j H n j n n 1 end end end end H 结果 H 153 370 371 407 7 代码 1 function f judge x f 1 if x 4 else for i 2 x 1 k mod x i if k 0 f 0 end end end end 2 n 0 H for i 2 20000 f judge i if f 1 n n 1 H n i end end H L size H hist H 10 等距区间内质数的个数的直方图 以 1000 为间距 8 代码 function f1 f2 f3 PI x h H 3 1 3 for i 10 x 1 10 x 1 j fix i 3 k fix i 3 3 for l k j h i l i l if abs h 3 pi abs H 3 pi H h end end end f1 H 1 f2 H 2 f3 H 3 end 结果 a b x PI 2 a 22 b 7 x 3 1429 a b x PI 3 a 355 b 113 x 3 1416 a b x PI 4 a 1065 b 339 x 3 1416 说明 a b x PI y 其中 a 代表分子 b 代表分母 y 为位数 附加 验证克莱姆法则 n 阶 程序 function X cramer a b A det a k size a if A 0 disp 输入有错 解不唯一 else for i 1 k 1 D a D i b Di det D X i Di A end end X X End 结果验证 如验证 4 阶 a rand 4 b 1 2 3 4 X cramer a b c a X Ycramer X 73 6009 6 6966 69 4126 3 4100 c 1 0000 2 0000 3 0000 4 0000 总结与体会 第一次写实验报告是初次将学到的东西进行运用 也将平时老师讲的一些知识点进行连 接 实验有助于巩固自己学到的知识 一二题是老师上课讲的一些比较基础的知识 三题需 要用到矩阵的知识 四五题则是相应的作图知识 尤其是在四