实验五 M文件和MATLAB程序设计
Matlab程序设计打印版
实验四 Matlab程序设计[实验目的]1．掌握字符串数组的创建和构造方法及常用字符串函数的使用。
2．熟练掌握MATLAB 控制流的使用方法。
3．熟悉M 脚本文件、函数文件的编写方法和技巧。
[实验原理]与数值数组相比，串数组在MATLAB 中的重要性较小，但不可缺少。
如果没有串数组及相应的操作，那么数据可视化、图形用户界面的制作将会遇到困难。
字符串与数值数组是两种不同的数据类，它们的创建方式也不同。
字符串的创建方式是：将待建的字符放在“单引号对”中。
注意，“单引号对”必须是在英文状态下输入，其作用是MATLAB 识别送来内容“身份”所必需的，如A=’This is an example!’；就创建了一个字符串A。
注意创建带单引号的字符串时，每个单引号符用“连续2 个单引号符”标识。
字符串的标识同数值数组同，而且也可以使用size 指令观察串数组的大小。
串数组的ASCII 码可以通过指令abs 和double 来获取，而用char 指令可以把ASCII 码变为串数组，另外，MATLAB 可以很好的支持中文字符串数组。
对于复杂串数组的创建，一是可以直接创建，但是要保证同一串数组的各行字符数相等，即保证各行等长，不推荐，太繁琐。
二是可以利用串操作函数创建多行数组，比如char, str2mat, strvcat 等，具体操作自己通过帮助体会。
另外还可以通过转化函数产生数码字符长，比如A_str=int2str(A) 就是把整数数组A 转换成串数组，如果是非整数将被四舍五入后再转换，类似的函数还有num2str（把非整数数组转换为串数组，常用于图形中数据点的标识）、mat2str（把数值数组转换成输入形态的串数组，常与eval 指令配用）。
假如想灵活运用MATLAB 去解决实际问题，想充分调动MATLAB——科学技术资源，想理解MATLAB 版本升级所依仗的基础，那么掌握M 脚本文件合函数的编写规则将十分有用。
MATLAB及其应用实验指导书
实验一 MATLAB 软件安装和使用初步一、 实验目的学习MATLAB 软件的安装过程，熟悉MA TLAB 软件界面的组成及基本使用方法。
二、 实验要求1、 掌握MATLAB 软件的启动与退出的方法。
2、 掌握MATLAB 软件的桌面环境。
3、 掌握MATLAB 软件的基本使用方法。
4、 掌握M 文件编辑器的使用。
5、 学会使用MATLAB 的帮助系统。
三、 实验内容1．通过桌面快捷图标或“开始—程序”菜单运行MATLAB ，熟悉MA TLAB 的桌面环境，如图1 所示。
图1 MATLAB 桌面环境2．浏览各个菜单项的内容，试着改变各个窗口的字体与字号。
3．已知1u =，3v =，用MA TLAB 分别执行下列语句。
a. 43u vb. ()222v u v -+ c. 333v v u - d.243v π 4. 分别运行以下命令并将运行结果与工作空间浏览器中的状态进行比较： >> who>> whos 5．画出函数0.22x y e -=的图形。
要求在M 文件编辑器中输入以下语句，然后保存，运行程序并查看显示的结果。
x = 0:0.1:10;y = 2 * exp(-0.2 * x); plot(x, y);查看工作空间浏览器workspce 中的当前状态，查看变量x 的内容，而后使用数字编辑器（array editor ）修改变量 x 的内容，而后在command window 中输入指令plot(x,y),观察figure window 中图形的变化。
6．通过工作目录浏览器（Current Directory Browser ）或者“File ”菜单中相应菜单项改变当前工作 目录。
改变当前路径如f:\mydir ，然后打开一个M 文件编辑器，输入如下指令： %Create an input array from -2*pi t0 2*pi t=-2*pi:pi/10:2*pi; x=abs(sin(t)); plot(t,x) 保存文件test1.m ，而后在command window 输入指令test1，观察结果。
实验五(线性方程组的数值解法和非线性方程求解)
1大学数学实验 实验报告 | 2014/4/5一、 实验目的1、学习用Matlab 软件数值求解线性代数方程组，对迭代法的收敛性和解的稳定性作初步分析；2、通过实例学习用线性代数方程组解决简化问题。
二、 实验内容项目一：种群的繁殖与稳定收获：种群的数量因繁殖而增加，因自然死亡而减少，对于人工饲养的种群（比如家畜）而言，为了保证稳定的收获，各个年龄的种群数量应维持不变。
种群因雌性个体的繁殖而改变，为方便起见以下种群数量均指其中的雌性。
种群年龄记作k=1,2,…,n ,当年年龄k 的种群数量记作x k ，繁殖率记作b k （每个雌性个体1年的繁殖的数量），自然存活率记作s k （s k =1−d k ，d k 为1年的死亡率），收获量记作ℎk ，则来年年龄k 的种群数量x ̌k 应该为x ̌k =∑b k n k=1x k , x ̌k+1=s k x k −ℎk , (k=1,2,…,n -1)。
要求各个年龄的种群数量每年维持不变就是要求使得x ̌k =x k , (k=1,2,…,n -1).(1) 如果b k , s k 已知，给定收获量ℎk ，建立求各个年龄的稳定种群数量x k 的模型（用矩阵、向量表示）.(2) 设n =5,b 1=b 2=b 5=0,b 3=5,b 4=3,s 1=s 4=0.4,s 2=s 3=0.6,如要求ℎ1~ℎ5为500,400,200,100,100，求x 1~x 5.(3) 要使ℎ1~ℎ5均为500，如何达到？问题分析：该问题属于简单的种群数量增长模型，在一定的条件（存活率，繁殖率等）下为使各年龄阶段的种群数量保持不变，各个年龄段的种群数量将会满足一定的要求，只要找到种群数量与各个参量之间的关系，建立起种群数量恒定的方程就可以求解出各年龄阶段的种群数量。
模型建立：根据题目中的信息，令x ̌k =x k ，得到方程组如下：{x ̌1=∑b k nk=1x k =x 1x ̌k+1=s k x k −ℎk =x k+1整理得到：{−x 1∑b k nk=1x k =0−x k+1+s k x k =ℎk2 大学数学实验 实验报告 | 2014/4/52写成系数矩阵的形式如下：A =[b 1−1b 2b 3s 1−100s 2−1…b n−1b n0000⋮⋱⋮000000000⋯00−10s n−1−1]令h =[0, ℎ1,ℎ2,ℎ3,…,ℎn−2,ℎn−1]Tx =[x n , x n−1,…,x 1]T则方程组化为矩阵形式：Ax =h ，即为所求模型。
实验5
实验五第3章MATLAB程序设计第二节程序控制结构程序控制结构有三种：顺序结构、旋择结构、循环结构.任何复杂的程序都可以由这三种基本结构构成.Matlab提供了实现控制结构的语句，利用这些语句可以编写解决实际问题的程序.一、顺序结构顺序结构是指按照程序中语句的排列顺序从上到下依次执行，直到程序的最后一个语句如例1的qiu文件和例2的dd1文件.这是最简单的一种程序结构.一般涉及数据的输入、计算或处理、数据的输出等内容.1. 数据的输入通过input命令来接收从终端输入的内容，它也可以显示文本和提示.命令集14 输入命令input(out ,in)在屏幕上显示出字符串out的文本并等待终端的输入.如果变量in是's'，则输入的内容以字符串的形式进行保存，通常MATL AB在保存前要尽可能地求出表达式的值.如果使用格式控制符号如'\n'，字符串out可以是若干行.2. 数据的输出可以通过简单地输入变量的名字来显示数字矩阵或者字符串向量的内容，结果将显示出变量的名字和内容.另一种显示变量的值就是使用命令disp.使用它只显示出变量的内容，这是有用的，特别是在字符串的应用中.命令集15显示命令disp(A)显示矩阵A的内容，如果A是字符串，则显示出它的文本.二、循环结构循环结构是按照给定的条件，重复执行指定的语句.Matlab用于实现循环结构的语句有for — end语句和while — end语句.1．for —end语句for-end语句的格式为：for 循环变量=表达式1 : 表达式2 : 表达式3循环体语句end其中表达式1的值为循环变量的初值，表达式2的值为步长，表达式3的值为循环变量的终值.步长为1时，表达式2可以省略.循环语句应用时应注意：● 循环语句可以嵌套使用；● 不能在 for — end 循环体内改变循环变量的值； ● 为了提高代码的运行效率，应尽可能提高代码的向量化程度，避免 for — end 循环的使用，如例8；● 如果预先就知道循环的次数，则可以采用 for — end 循环；否则，如果预先无法确定循环的次数，则可以使用 while — end 循环.2. while — end 循环while-end 语句的一般格式为：while (条件) 循环体语句 end 其执行过程为：若条件成立，则执行循环体语句，执行后再判断条件是否成立，如果不成立则跳出循环.■例1 使用while-end 结构计算1+2+3+…+100 . 解：■例2 计算 Matlab 中 eps 值 (exp1.m) 解：s=0;for i=1:100 s=s+i; end ss=0; i=1;while i<=100 s=s+i; i=i+1; end snum=0; EPS=1; while (1+EPS)>1 EPS=EPS/2; num=num+1; endEPS,num三、选择结构选择结构是根据给定的条件成立或不成立，分别执行不同的语句.Matlab 用于实现选择结构的语句有条件语句：if –end 和 开关语句：switch-end .1．单分支结构 if —end 语句 2．双分支结构 if —else —end 语句 if-end 语句的格式为： if-else-end 语句的格式为：3．多分支结构 if — elseif — else — end 语句 if — elseif — else — end 语句的格式为：if 条件表达式 语句组 endif 条件表达式 语句组1 else语句组2 endif 条件表达式1 语句组1 elseif 条件表达式 2 语句组2 ... ...elseif 条件表达式 m 语句组m else 语句组end■例3求函数1;00;1;x y x x -<⎧⎪= =0⎨⎪ >⎩ ，输入x 值，输出y 值.解：■例4一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数.输出全部水仙花数.解：程序如下运行结果：153 370 371 407function y=fuhao(x) if x<0 y=-1; elseif x==0 y=0; else y=1; endA=[];for m=100:999m1=fix(m/100); %求m 的百位数字 m2=rem(fix(m/10),10); %求m 的十位数字 m3=rem(m,10); %求m 的个位数字 if m==m1*m1*m1+m2*m2*m2+m3*m3*m3 A=[A,m];endend disp(A)4. switch-end 语句根据表达式的不同取值，分别执行不同的语句.其格式如下：■关于 switch-end 语句的几点注解：● Matlab 首先计算表达式的值，然后将它依次与各个case 指令后的检测值进行比较，当比较结果为真时，就执行相应的语句组，然后跳出 switch-end 结构● 如果所有的比较结果都为假，则执行 otherwise 后面的语句组，然后跳出 switch-end 结构 ● otherwise 指令可以不出现● switch 后面的表达式的值可以是一个标量或字符串 ● case 指令后的检测值超过一个时，应用｛｝括起来switch 表达式 case 值1 语句组1 case 值2 语句组2 ... ... case 值m 语句组m otherwise 语句组n end■例5 阅读程序method=input('请输入方法名：')switch methodcase {'linear','bilinear'}disp('Method is linear')case 'cubic'disp('Method is cubic')case 'nearest'disp('Method is nearest')otherwisedisp('Unknown method.')end■例6用switch-end语句编一程序，对于给定的一个百分制成绩，输出相应的五分制成绩，设：90分以上为…A‟, 80～89分为…B‟，70～79分为…C‟，60～69分为…D‟，60分以下为…E‟.function sctole(score)grade=fix(score/10);switch gradecase {10,9}grade='A'case 8grade='B'case 7grade='C'case 6grade='D'otherwisegrade='E'end四、其他流程控制语句●break语句●continue语句●return语句●try语句●echo命令●keyboard命令●pause命令1.break和continue语句与循环结构相关的语句还有break语句和continue语句.它们一般与if语句配合使用.break语句用于终止循环的执行.当在循环体内执行到该语句时，程序将跳出循环，继续执行循环语句后的下一语句.continue语句控制跳过循环体中的某些语句.当在循环体内执行到该语句时，程序将跳过循环体中所有剩下的语句，继续下一次循环. ■例7求[100，200]之间第一个能被21整除的整数.for n=100:200if rem(n,21)~=0continueendbreakendn for n=100:200if rem(n,21)==0breakendendn2.try语句语句格式为：try语句组1catch语句组2endtry语句先试探性执行语句组1，如果语句组1在执行过程中出现错误，则将错误信息赋给保留的lasterr变量，并转去执行语句组2.这种试探性执行语句是其他高级语言所没有的.■例8矩阵乘法运算要求两矩阵的维数相容，否则会出错.先求两矩阵的乘积，若出错，则自动转去求两矩阵的点乘.function C=rc(A,B)%A,B为矩阵tryC=A*B;catchC=A.*B;endlasterr %显示出错原因3.请求键盘输入命令keyboard暂停运行程序并处于等待键盘输入状态（便于调试及程序运行中修改变量），键入“return”程序继续运行.如运行如下：4. echo 命令（在命令窗口，程序调用前使用）一般情况下，M 文件执行时，文件中的命令不会显示在命令窗口.echo 命令可使文件命令在执行时可见.这对程序的调试及演示很有用.对命令文件和函数文件，echo 的作用稍微有些不同.●对命令文件，其格式如下：echo on 打开命令文件的回应命令 echo off 关闭回应命令●对函数文件，其格式如下：echo file on 使指定的file 文件的命令在执行中被显示出来 echo file off 关闭指定文件的命令在执行中的回应function y=f(x) keyboardy=(5*x-6)/((x-3)*(x-3)+2);>> f K>> x=4; K>> return ans = 4.6667>> y=f(4) K>> return y = 4.6667。
matlab实验五循环结构程序设计实验报告
实验五：循环结构程序设计实验报告一、实验目的1. 了解循环结构的基本概念和原理；2. 掌握使用Matlab进行循环结构程序设计的方法；3. 提高编程能力和问题解决能力。
二、实验内容1. 学习while循环和for循环的基本语法和用法；2. 设计并编写一些基于循环结构的Matlab程序；3. 分析程序运行结果并进行总结。
三、实验原理循环结构是程序设计中的重要组成部分，它可以让一段代码重复执行多次，从而简化程序的编写。
在Matlab中，循环结构主要有while 循环和for循环两种形式。
while循环在执行循环体前判断条件是否为真，只有条件为真时才会执行循环体；for循环则是指定循环的次数，每次迭代时执行一次循环体。
四、实验步骤与结果1. 编写一个使用while循环的程序，实现1~100的累加。
程序代码如下：```matlabsum = 0;i = 1;while i <= 100sum = sum + i;i = i + 1;enddisp(sum);```运行程序后，得到的结果为5050。
2. 编写一个使用for循环的程序，实现1~10的阶乘。
程序代码如下：```matlabresult = 1;for i = 1:10result = result * i;enddisp(result);```运行程序后，得到的结果为xxx。
五、实验总结通过本次实验，我深刻理解了循环结构的基本概念和原理，掌握了使用Matlab进行循环结构程序设计的方法。
在编写循环程序的过程中，我发现循环结构可以大大简化程序的编写，并且能够高效地处理重复性任务。
我也进一步提高了自己的编程能力和问题解决能力。
本次实验使我对Matlab中的循环结构有了更加深入的了解，我相信这对我的编程能力和日后的学习工作都将大有裨益。
六、实验感想和改进意见通过本次实验，我深刻认识到循环结构在程序设计中的重要性和灵活性。
循环结构能够帮助我们简化程序的编写，提高代码的重用性和可读性，因此在实际的程序设计中，合理地运用循环结构能够大大提高程序的效率和逻辑清晰度。
实验五 M文件及程序设计
一、实验目的1.练习函数的定义和书写及调用2.练习M文件下书写代码二、实验设备装有MATLAB的计算机。
三、实验内容1，分别编写命令（M）文件和函数（function）文件，在命令文件中可循环输入学生成绩，在函数文件判断成绩所属等级。
统计各分段的人数，并将各人的成绩按优、良、中、及格和不及格表示。
2.猜数游戏。
首先由计算机产生一个[1,100]之间的随机整数，然后由用户猜测所产生的随机数。
根据用户猜测的情况给出不同提示，如猜测的数大于产生的数，则显示“High”并请用户继续猜测，小于则显示“Low”并请用户继续猜测，等于则显示“You won”同时退出游戏。
用户最多可以猜7次。
如果第七次还未猜对，则显示“You fail!"并退出游戏。
四、实验步骤和结果1. 代码结果:输入正数开始5请输入学生成绩x=89请输入学生姓名l请输入学生成绩x=78请输入学生姓名h请输入学生成绩x=45请输入学生姓名k请输入学生成绩x=-1ans =优秀的人数0ans =良好的人数1ans =中等的人数1ans =及格的人数0ans =不及格的人数1姓名成绩等级l 89 良好h 78 中等k 45 不及格2.代码五、总结1.实验练习了函数的创建和使用，并且在M的文件中调用函数。
2.复习了单元数组和结构数组的知识。
3. 使用了字符串比较函数，在M文件中对各个等级的人数进行了统计。
4.猜数游戏中开始的代码，当第七次猜对后会输出”You Fail”,后对代码进行了改进。
MATLAB实验
MATLAB实验一：MATLAB语言基本概念实验实验目的：1. 熟悉MATLAB语言及使用环境；2.掌握MATLAB的常用命令；3．掌握MATLAB的工作空间的使用；4.掌握MATLAB的获得帮助的途径。
5.掌握科学计算的有关方法，熟悉MATLAB语言及其在科学计算中的运用；6.掌握MATLAB的命令运行方式和M文件运行方式；7.掌握矩阵在MATLAB中的运用。
实验方案分析及设计：本次实验主要目的是了解MATLAB的使用环境，以及常用的一些命令的使用；了解矩阵在MATLAB实验中的具体运用，以及相关的一些符号命令的使用。
实验器材：电脑一台，MATLAB软件实验步骤：打开MATLAB程序，将实验内容中的题目依次输入MATLAB中，运行得到并记录结果，最后再对所得结果进行验证。
实验内容及要求：1．熟悉MATLAB的菜单和快捷键的功能2．熟悉MATLAB的命令窗口的使用3．熟悉常用指令的使用format clc clear help lookfor who whos 4．熟悉命令历史窗口的使用5. 熟悉MATLAB工作空间的功能将工作空间中的变量保存为M文件，并提取该文件中的变量6.熟悉MATLAB获取帮助的途径将所有plot开头的函数列出来，并详细给出plotfis函数的使用方法1. 输入 A=[7 1 5;2 5 6;3 1 5]，B=[1 1 1; 2 2 2;3 3 3]，在命令窗口中执行下列表达式，掌握其含义：A(2, 3) A(:,2) A(3,:) A(:,1:2:3)A(:,3).*B(:,2) A(:,3)*B(2,:) A*BA.*BA^2 A.^2 B/A B./AA=[7 1 5;2 5 6;3 1 5]7 1 52 5 63 1 5>> B=[1 1 1; 2 2 2;3 3 3]1 1 12 2 23 3 3>> A(2, 3)6>> A(:,2)151>> A(3,:)3 1 5>> A(:,1:2:3)7 52 63 5>> A(:,3).*B(:,2)51215>> A(:,3)*B(2,:)10 10 1012 12 1210 10 10>> A*B24 24 2430 30 3020 20 20>> A.*B7 1 54 10 129 3 15>> A^266 17 6642 33 7038 13 46>> A.^249 1 254 25 369 1 25>> B/A0.1842 0.2105 -0.23680.3684 0.4211 -0.47370.5526 0.6316 -0.7105>> B./A0.1429 1.0000 0.20001.0000 0.4000 0.33331.0000 3.0000 0.60002．输入 C=1:2:20，则 C （i ）表示什么？其中 i=1,2,3, (10)1到19差为2，i 代表公差3. 试用 help 命令理解下面程序各指令的含义：cleart =0:0.001:2*pi;subplot(2,2,1);polar(t, 1+cos(t))subplot(2,2,2);plot(cos(t).^3,sin(t).^3)subplot(2,2,3);polar(t,abs(sin(t).*cos(t)))subplot(2,2,4);polar(t,(cos(2*t)).^0.5)4计算矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡897473535与⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡638976242之和。
实验五+MATLAB数值计算(含实验报告)
实验五 MATLAB 数值计算一、实验目的1．掌握求数值导数和数值积分的方法。
2．掌握代数方程数值求解的方法。
3．掌握常微分方程数值求解的方法。
二、实验的设备及条件计算机一台（带有MATLAB7.0以上的软件环境）。
设计提示1．参考本节主要内容，学习并理解相关函数的含义及调用方法。
三、实验内容1．线性系统方程：分别使用左除（\）和求逆（inv ）求解下面系统方程的解：⎪⎩⎪⎨⎧=+=+=++377251463c b b a c b a2. 数值积分：使用quad 和trapz 求解⎰-503/dx xe x 的数值积分，并与其解析解9243/5+--e 相比较；3. 请完成教材P154页中实验指导环节的实验内容第2题4. 请完成教材P155页中思考练习的第3题（1），并绘制解在该求解区间（即[0,5]）上的图像；。
5、请完成教材P164页实验指导环节的实验内容第5题。
（提示：该函数的符号导数，可以通过函数diff 求得。
首先定义符号变表达式，如求sin(x)的一阶符号导数，可以先定义f=’sin(x)’;df=diff(f);可求得df=cos(x)。
其中df 即为函数f 的一阶符号导数）。
四、实验报告要求（包含预习报告要求和最终报告要求）1.实验名称2.实验目的3.实验设备及条件4.实验内容及要求5.实验程序设计指程序代码。
6.实验结果及结果分析实验结果要求必须客观,现象。
结果分析是对实验结果的理论评判。
7.实验中出现的问题及解决方法8. 思考题的回答五、实验报告的提交方式Word文档，命名方式:实验号_你的学号_姓名例如本次实验：实验一_000000001_张三.doc（信息101提交报告邮箱）：E_mail: *******************(网络工程101提交作业邮箱)：E_mail: *******************（注意网络班的M是大写的)下一次课前提交,过期不收!六、参考文献参考教材和Matlab帮助文件。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验五 M文件和MATLAB程序设计一、实验目的matlab作为一种高级计算机语言，不仅可以命令行方式完成操作，也具有数据结构、控制流、输入输出等能力，本次实验通过熟悉和掌握m文件的建立与使用方法，以及函数与控制程序流程语句的使用，使学生具备一定的编程和程序调试能力。
1．掌握M文件的使用方法。
2．掌握if语句和switch语句的使用3. 掌握循环语句的使用4. 通过练习理解MATLAB编程方法。
二、实验原理1．m文件用matlab语言编写的程序，称为m文件。
M文件根据调用方式的不同分为两类，命令文件（Script file）和函数文件（Function file）。
区别?2．程序控制结构1）顺序结构2）选择结构（1）if语句a) 单分支if语句b) 双分支if语句c) 多分支if语句（2）switch 语句（3）try语句3）循环结构（1）for 语句（2）while语句（3）break语句、continue语句、return使用,区别?3．函数文件function 输出形参表＝函数名（输入形参表）注释说明部分函数体语句三、实验要求1．首先上机练习PPT中各种流程控制语句的有关实例。
2．然后上机练习下面的实验习题。
四、实验习题1．数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。
重复此过程，最终得到的结果为1。
如：2→13→10→5→16→8→4→2→16→3→10→5→16→8→4→2→1运行下面的程序，按程序提示输入n=1,2,3,5,7等数来验证这一结论。
%classic "3n+1" problem from number theory.while 1n=input('Enter n,negative quits:');if n<=0breakenda=n;while n>1if rem(n,2)==0n=n/2;elsen=3*n+1;enda=[a,n];endaendEnter n,negative quits:3a =3 10 5 16 84 2 12． 编程求满足∑=>m i i 1100002的最小m 值。
a=0;i=1; while (a<10000) a=a+pow2(i); i=i+1; end m=i-1; m 13 3． 编写一个函数，计算下面函数的值，给出x 的值，调用该函数后，返回y 的值。
function [y]=myfun1(x)⎪⎩⎪⎨⎧>+-≤<≤=3,630,0,sin )(x x x x x x x y选择一些数据测试你编写的函数。
function y=myfun1(x) if x<=0 y=sin(x);elseif x>3 y=-x+6; end 运行结果: >> y=myfun1(-pi/2) y = -1 >> y=myfun1(0)y = 0 >> y=myfun1(2) y = 2 >> y=myfun1(4) y = 24． 用如下迭代公式求a ，a 的值分别为：3，17，113。
迭代的终止条件为5n 1n 10x x -+≤-，迭代初值0.1x 0=,迭代次数不超过100次。
分别对迭代结果和准确值进行比较，并统计迭代次数。
迭代公式：11122--++=n n n x a x xfunction [x, n]=sqrt_a(a) x=1.0; for k=1:100 m=x; x=x/2+a/(2*x); if abs(x-m)<=10^(-5) break end end x; n=k; s=(x-sqrt(a)); if s<=10^(-5) disp('正确'); else disp('错误'); end 下面调用举例：运行结果:>> [x ,n]=sqrt_a(3) 正确x = 1.7321 n = 5 >> [x, n]=sqrt_a(17) 正确 x = 4.1231 n =正确x =10.6301n =85．定义一个函数文件，求给定复数的指数、对数、正弦和余弦，并在命令文件中调用该函数文件。
程序设计：function [e ln s c]=num(x)e=exp(x)ln=log(x)s=sin(x)c=cos(x)end运行结果：>> num(5i)e =0.2837 - 0.9589iln =1.6094 + 1.5708is =0 +74.2032ic =74.2099ans =0.2837 - 0.9589i6． 设 f(x)=01.011.01)3()2(42+++--x x ，编写一个MATLAB 函数文件fx.m ，使得调用f(x)时，x 可用矩阵代入，得出的f(x)为同阶矩阵。
程序设计：函数文件fx.m:function A=fx(x)A=1./((x-2).^2+0.1)+1./(((x-3).^4)+0.01)end调用fx.m 的命令文件：>> A=fx([1 2;2 3;4 3])运行结果:A =0.9716 10.990110.9901 100.90911.2340 100.90917． 已知y=)20()30()40(f f f 当f(n)=n+10ln(n 2+5)时，求y 的值。
程序设计：函数文件fn.m:function x=fn(n)x=n+10*log(n^2+5)end调用fn.m 的命令：y=fn(40)/(fn(30)+fn(20))运行结果:y =0.63908．将下列脚本式m 文件改写成函数式m 文件：（1）输入一个字符，若为大写字母，则输出其对应的小写字母；若为小写字母，则输出其对应的大写字母；若为数字字符则输出其对应的数值，若为其他字符则原样输出。
c=input('请输入一个字符','s'); %a 是字符，'a'是字符串if c>='A' & c<='Z'disp(setstr(abs(c)+abs('a')-abs('A'))); %加32elseif c>='a'& c<='z'disp(setstr(abs(c)- (abs('a')-abs('A'))));elseif c>='0'& c<='9'disp(abs(c)-abs('0'));elsedisp(c);end %必须与if 配对使用编写一个函数式m 文件，实现程序的功能c=input('请输入一个字符','s');Function M=disp(c);if c>='A' & c<='Z'M=setstr(abs(c)+abs('a')-abs('A'))elseif c>='a'& c<='z'M=setstr(abs(c)- (abs('a')-abs('A')))elseif c>='0'& c<='9'M=abs(c)-abs('0')end（2）某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price 来表示)：price<200 没有折扣200≤price<500 2%折扣500≤price<1000 5%折扣1000≤price<2500 10%折扣2500≤price<5000 25%折扣5000≤price 50%折扣输入所售商品的价格，求其实际销售价格。
程序如下：price=input('请输入商品价格');switch fix(price/100) %向零的方向取整case {0,1} %价格小于200rate=0;case {2,3,4} %价格大于等于200但小于500rate=2/100;case num2cell(5:9) %价格大于等于500但小于1000%为简便采用元胞数组rate=5/100;case num2cell(10:24) %价格大于等于1000但小于2500rate=10/100;case num2cell(25:49) %价格大于等于2500但小于5000rate=25/100;otherwise %价格大于等于5000rate=50/100;endprice=price*(1-rate) %输出商品实际销售价格编写一个函数式m文件，计算顾客实际付款数目。
Function price=f(rate)switch fix(price/100)case {0,1}。