实验三(3-1)

数学与信息科学学院实验报告

学生姓名： 李双喆 实验地点：11-326、11-327 实验项目序号：3-1 学号：10040342013 班级： 信本二班 实验日期：2011年10 月20 指导教师：曹慧荣 日 实验类型：验证型 学时：2 一、实验项目名称：matlab基本编程—MATLAB编程介绍及循环结构与应用 二、实验目的与要求：

利用Matlab语言编写程序，掌握Matlab文件操作；循环结构与应用 三、实验原理：

M文件的建立、保存与运行；disp语句与input语句的使用；循环结构与应用。 四、实验方案设计：

五、测试数据与实验结果

测试数据1: 采用以下方式绘制球面x2?y2?z2?r2的图形

（1） 采用input语句输入球的半径，编写m文件绘制给定半径时球面图

形；

r=input('r='); t=0:pi/10:2*pi; s=0:pi/10:2*pi; [T,S]=meshgrid(t,s); x=r.*sin(T).*cos(S); y=r.*sin(T).*sin(S); z=r.*cos(T); surf(x,y,z)

（2） 编写函数文件qiu.m,其中球的半径为输入变量，在命令窗口输入

qiu(4),qiu(3),qiu(9)时保存对应的球面图形；

420-2-4420-2-4-4-2204

3210-1-2-3420-2-4-4-2204

1050-5-101050-5-10-10-55010

测试数据2: 参考教材52页例题3-2，解决例题后提出的【进一步问题】 n=input('n=') p=1; for i=1:n; p=p*i

fprintf('i=%.0f,p=%.0f\\n',i,p) end

测试数据3: 参考教材52页例题3-3，解决例题后提出的【进一步问题】 e=input('e=：') n=1;i=1;p=1;s=1; while e<3/(p*(i+1))

p=p*i s=s+1/p i=i+1

fprintf('i=%.0f,s=%.8f\\n',i,s) end

测试数据4: 参考教材53页例题3-4，解决例题后提出的【进一步问题】 n=0; s=0;

while s<=1000 n=n+1; s=s+sqrt(n); if s>1000 s=s-sqrt(n) n=n-1 break end end 六、实验总结

七、部分参考代码（可附页或提交电子版）

测试数据1的Matlab代码：qui r=2 >> qui r=4 >> qui

r=3 >> qui r=9 >>

测试数据2的Matlab代码： n=input('n=') p=1; for i=1:n; p=p*i

fprintf('i=%.0f,p=%.0f\\n',i,p) end liti32 n=2 n = 2 p = 1 i=1,p=1 p = 2 i=2,p=2

测试数据3的Matlab代码：

e=input('e=：') n=1;i=1;p=1;s=1; while e<3/(p*(i+1)) p=p*i s=s+1/p i=i+1

fprintf('i=%.0f,s=%.8f\\n',i,s) end li e=：0.1 e = 0.1000 p = 1 s = 2 i = 2

i=2,s=2.00000000 p = 2 s =

2.5000 i = 3

i=3,s=2.50000000 p = 6 s = 2.6667 i = 4

i=4,s=2.66666667 >>

测试数据4的Matlab代码： n=0; s=0;

while s<=1000 n=n+1; s=s+sqrt(n); if s>1000 s=s-sqrt(n) n=n-1 break end end

genhao s = 993.6487 n = 130 >>

本文来源：https://www.bwwdw.com/article/0r6x.html