2022年用MATLAB编程求π值定义 .pdf

上海电力学院高级程序设计（ C）课程设计报告题目：用 MATLAB 编程求 院系：计算机科学与技术学院专业年级：信息安全 2012 级学生姓名：涂桂花学号：20123333 指导教师：魏为民2015年 3 月 14 日名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 6 页 - - - - - - - - - 目录一、实验目的 . 1 二、实验内容和步骤 . 1 1. 操作环境 . 1 2. 系统配置 . 1 3. 操作步骤 . 1 4. 程序源代码 . 3 三、实验结果 . 3 四实验小结 . 4 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 6 页 - - - - - - - - - 1 上 海 电 力 学 院实 验 报 告课程名称实验项目姓名学号班级专业同组人姓名指导教师魏为民实验日期注意：报告文件名格式：学号_姓名 _实验项目 .doc ，电子版发邮件到，打印版按要求黑白打印。一、实验目的1. 了解 MATLAB 编译器的功能并学会使用；2. 学会使用 MATLAB 语言编写简单的程序；3. 用下列公式求值，并了解 的其他公式。二、实验内容和步骤1. 操作环境电脑型号华硕 X550EP 笔记本电脑操作系统 Windows 7 旗舰版 64 位 SP1 ( DirectX 11 ) 2. 系统配置处理器AMD E1-2100 APU with Radeon HD Graphics 双核主板华硕 X550EP ( AMD 未知 (1536h) ) 内存4 GB ( DDR3 1600MHz ) 主硬盘日立 HGST HTS545050A7E680 ( 500 GB / 5400 转/ 分 ) 显卡ATI Radeon HD 8210 ( 512 MB / 华硕 ) 显示器奇美 CMO15A7 ( 15.7 英寸 ) 光驱松下 DVD-RAM UJ8C2 S DVD刻录机声卡瑞昱 High Definition Audio 网卡瑞昱 RTL8168/8111/8112 Gigabit Ethernet Controller / 华硕3. 操作步骤3.1 打开 MATLAB 编辑器，如图1-1 所示。信息安全用 MA TLAB 编程求 涂桂花20123333 2012252 信息安全无2015.3.14 7566534431222名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 6 页 - - - - - - - - - 2 图 1-1 3.2 打开 editor编辑框编写程序，在“Command Window ”窗口里输入edit ，然后回车。结果如图1-2所示。图 1-2 3.3 在上图的editor里编写代码，编写完成之后，按 F5，然后在Command Window 里输入你要求的值p, 然后回车，显示出值结果。结果如图1-3 所示。图 1-3 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 6 页 - - - - - - - - - 3 3.4 改变 N的值，分析得到的不同值的大小与N 的关系。PI=3. 141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825 （1）当 N=10*104 时， =3.14151411868195662435709891724400222301483154296875 （2）当 N=10*108 时， =3.141592643066262180440162410377524793148040771484375 4. 程序源代码digits(10000); N=100000000;p=1; for m=1:N;m1=m*m;m2=2*m; n=(4*m1)/(m2-1)*(m2+1); p=p*n; end p=vpa(p*2); 三、实验结果1.测试数据（1）N=10*104 ， =3.14151411868195662435709891724400222301483154296875 （2）N=10*108 ， =3.141592643066262180440162410377524793148040771484375 （3）N=10*109 ， = 3.141592643066262180440162410377524793148040771484375 （4）N=10*1010, =3.141592643066262180440162410377524793148040771484375 2. 测试结果（1）当 N=10*104 时， =3.14151411868195662435709891724400222301483154296875 （2）当 N=10*108 时， =3.141592643066262180440162410377524793148040771484375 （5）当 N=10*109 ，= 3.141592643066262180440162410377524793148040771484375 （6）当 N=10*1010, =3.141592643066262180440162410377524793148040771484375 3. 截图名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 6 页 - - - - - - - - - 4 四实验小结如结果分析、心得体会、参考文献，以及与老师交流或需要老师特别注意的地方等。在这次实验中，我对MATLAB 有一些了解，发现它用作画图工具是很方便的。另外，还可以做一些计算方面的工作，语句简单易于理解。可是在使用MATLAB 编译器的时候，有时候总是操作导致Command Window里 的 不见了，导致后面无法将更改后的程序执行出结果，这时候就要重新启动MATLAB ，重新执行，比较麻烦。另外我发现当N值越大的时候，PI 的值精确度越高，然而N 大到一定范围后得出结果的过程也就变得非常慢，并且在10*108 后每增加10 倍后 PI 没有很大的变化，并且运算速度也变得非常慢了。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 6 页 - - - - - - - - -