第4章MATLAB程序设计.pptVIP

全局变量a的定义形式为：global a 需要注意的是，每一个全局变量在函数第一次使用之前必须声明，如果在本地工作区中已经被创建，那么声明为再次声明全局变量将会产生错误。为了避免这种错误，在函数中的初始注释行之后和第一个可执行性语句之前声明全局变量。 第4章MATLAB程序设计全文共52页，当前为第31页。 4.1 M文件 4.2 程序控制结构 4.3 函数文件 4.4 程序举例 4.5 程序调试 4.6 本章小结 第4章MATLAB程序设计全文共52页，当前为第32页。 例1 编写M函数实现：求一个数是否为素数，再编写一脚本文件，要求通过键盘输入一个整数，然后判断其是否为素数。 MATLAB代码如下： function [t] = sushu(n) % 判断yield数是否是素数 % 对根号n向下取整 k=floor(sqrt(n)); for i=2:k % 判断余数 if mod(n,i)==0 t=0; %当被i整除时，不是素数，终止循环 break; else t=1; end end MATLAB的脚本文如下： a = input(输入数据); if (sushu(a)==1) disp(a 是素数); else disp(a不是素数); end MATLAB程序的运行结果如下： 输入数据4 a不是素数 输入数据7 a 是素数 第4章MATLAB程序设计全文共52页，当前为第33页。 例2 编写求解方程ax2+bx+c=0的根的函数，这里应根据b2-4ac的不同取值分别处理，并输入几组典型值加以检验。 clc;clear; %输入参数值 a=input(请输入a的值：); b=input(请输入b的值：); c=input(请输入c的值：); disp([方程：,num2str(a),*x^2+,num2str(b),*x+,num2str(c),’=0’]) % 条件1 if a==0 % 语句1 x = -b/c; disp([根为：x=,num2str(x)]) 第4章MATLAB程序设计全文共52页，当前为第34页。 % 条件2 elseif b^2-4*a*c0 % 语句2 x1 = (-b+sqrt(b^2-4*a*c))/(2*a); x2 = (-b-sqrt(b^2-4*a*c))/(2*a); disp([根为：x1=,num2str(x1),;x2=,num2str(x2)]) % 条件3 elseif b^2-4*a*c==0 % 语句3 x=-b/(2*a); disp([根为：x1=x2=,num2str(x)]) % 条件4 else %语句4 disp(无解) end 输入a为4，b为8，c为4时，MATLAB程序的运行结果如下： 请输入a的值：4 请输入b的值：8 请输入c的值：4 方程：4*x^2+8*x+4=0 根为：x1=x2=-1 输入a为0，b为7，c为7时，MATLAB程序的运行结果如下： 请输入a的值：0 请输入b的值：7 请输入c的值：7 方程：0*x^2+7*x+7=0 x =-1 根为：x=-1 第4章MATLAB程序设计全文共52页，当前为第35页。 例3 利用menu函数输入选择参数ch。当ch=1时，产生[-10,10]之间均匀分布的随机数；当ch=2时，产生[-5,5]之间均匀分布的随机数；当ch=3时，产生[-1,1]之间均匀分布的随机数；当ch=4时，产生均值为0，方差为1的正态分布随机数。（要求使用switch函数） 第4章MATLAB程序设计全文共52页，当前为第36页。 MATLAB程序如下： %建立一个菜单 f = menu(ch,1,2,3,4); a = []; %情况选择 switch f case 1 a = 10-rand(2)*20 case 2 a = 5-rand(2)*10 case 3 a = 1-rand(2)*2 case 4 a = randn(2); otherwise disp(error!) end 单击按钮1的结果如下： a = -6.2945 7.4603 -8.1158 -8.2675 单击按钮2的结果如下： a = -1.3236 2.2150 4.0246 -0.4688 单击按钮3的结果如下： a = -0.9150 0.6848 -0.9298 -0.9412 单击按钮4的结果如下： a = 0.7254 0.7147 -0.0631 -0.2050 第4章M