第10章节数学软件使用初涉3MATLAB.ppt

* 中山大学计算机科学系 10.3.4 编写 MATLAB 程序 三、循环结构程序 包括 for 语句和 while 语句 例10.31 求 1～100 内所有素数之和。 for 实现 （S1031.m） while 实现 （S1031while.m） Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. * 中山大学计算机科学系 10.3.4 编写 MATLAB 程序 四、定义 M 函数 M 函数文件具有如下基本结构： function 输出形参表=函数名(输入形参表) %注释说明部分 函数体 其中，以 function 开头的一行为引导行，表示该 M 文件是一个函数文件；当输出形参有多个时，则用方括号[]括起来；第一行注释行是可以由命令 lookfor 搜索的文本行。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. * 中山大学计算机科学系 10.3.4 编写 MATLAB 程序 四、定义 M 函数 例10.32 定义一个求圆面积的 M 函数。程序如下： function area=S1032(r) %求圆面积 area=pi*r^2; %确定函数返回值 注意以下两点： 在函数体中，为输出形参赋值可以控制函数的返回值。 M 函数文件名应当与函数名一致，从而保持函数调用的一致性。例如，若要计算半径为3的圆面积，则输入： S1032(3) ans = 28.2743 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. * 中山大学计算机科学系 10.3.4 编写 MATLAB 程序 四、定义 M 函数 例10.33 定义求圆柱体体积的 M 函数，程序如下： function v=S1033(r,h) %（主函数）求圆柱体体积（圆半径为r，高为h） area=circleArea(r); %调用子函数circleArea v=area*h; function area=circleArea(r) %（子函数）求圆面积 area=pi*r^2; Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. * 中山大学计算机科学系 10.3.4 编写 MATLAB 程序 四、定义 M 函数 在一个 M 函数文件中，可以定义多个函数，但只有第1个函数（称为主函数）可被外部调用，而：其他函数（称为子函数）只能在该文件中被调用。如调用形式 S1033(3,4) 全局变量和局部变量 局部变量：按赋值方式定义的变量，存贮于特定工作空间，即 MATLAB 工作空间和函数工作空间 全局变量：使用 global 关键字声明的变量，存贮于 MATLAB 全局工作空间 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. * 中山大学计算机科学系 10.3.5 绘图 一、绘制二维图形 例10.34 在同一个直角坐标图中，绘制正弦函数 sin(x) 和余弦函数 cos(x) 在 [0,4π] 区间的曲线图。 x=0:pi/80:4*pi; plot(x,sin(x)); %绘制 sin 曲线 hold on %保持原图形 plot(x,cos(x),‘:’); %参数‘:’指定数据间连线为点线 legend(sin(x) ,cos(x),0); %自动放置图例注释 title(sin 与 cos 曲线); %设置图形标题 xlabel(X 轴); %置 X 轴标签 ylabel(Y 轴); %置 Y 轴标签 grid on; %显示网格 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011