实验五 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帮助文件。
MATLAB程序设计实验指导书
MATLAB程序设计实验指导书MATLAB是集数值计算、符号运算及图形处理等强大功能于一体的科学计算语言。
作为强大的科学计算平台，它几乎能满足所有的计算需求。
在美国及其他发达国家的理工科院校里，MATLAB已经作为一门必修的课程；在科研院所、大型公司或企业的工程计算部门，MATLAB也是最普遍的计算工具之一。
有鉴于此，我院开设了《MATLAB程序设计》这门课程，它需要一定的理论基础，同时又具有很强的实践性。
如何加强理论课程的学习、加深学生对本课程中的基本理论知识及基本方法的理解，如何培养学生实践动手能力是教学的当务之急。
而MATLAB程序设计实验课程就是一种重要的教学手段和途径。
实验将MATLAB程序设计的基本方法灵活地运用在数学、电路等课程中，重点突出，内容丰富。
同时，注重理论分析与实际动手相结合，以理论指导实践，以实践验证基本原理，旨在提高学生分析问题、解决问题的能力及动手能力，使学生进一步巩固基本理论知识，建立比较全面的MATLAB程序设计的概念。
实验注意事项1、实验系统接通电源前请确保电源插座接地良好。
2、完成实验后请确保关闭电脑电源及插座电源。
实验一 矩阵及其运算一、实验目的1、熟悉MATLAB 工作环境2、掌握矩阵和数组的创建、寻访和运算 二、实验内容验证欧姆定律：iur ，其中i u r , ,分别是电阻（欧姆）、电压（伏特）、电流（安培）。
已知u=[0.89, 1.20, 3.09, 4.27, 3.62, 7.71, 8.99, 7.92, 9.70, 10.41]，i=[0.028, 0.040, 0.100, 0.145, 0.118, 0.258, 0.299, 0.257, 0.308, 0.345]。
三、实验器材PC 机 （装有 MA TLAB 软件 ） 1台 四、实验原理 4.1 Desktop 简介MATLAB R2006a 版的Desktop 操作桌面，是一个高度集成的MATLAB 工作界面。
实验五 用Matlab数据拟合
2. 已知观测数据点如表所示 x 1.6 2.7 1.3 4.1 3.6 2.3 y 17.7 49 13.1 189.4 110.8 34.5
0.6
4
4.9
409.1
3
65
2.4
36.9
求a, b, c的值, 使得曲线 f(x)=aex+bsin x+c lnx 与已知数据 点在最小二乘意义上充分接近.
例4 已知观测数据点如表所示
x 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 y 3.1 3.27 3.81 4.5 5.18 6 7.05 8.56 9.69 11.25 13.17 求三个参数 a, b, c的值, 使得曲线 f(x)=aex+bx2+cx3 与 已知数据点在最小二乘意义上充分接近.
1.55
500 106
2.47
2. 93
3. 03
2.89
1000 106 1500 106 2000 106 2375 106
3 1.953 103 1.517 103 1.219 10
/ / N / m 2 3.103 103 2.465 103
解: 描出散点图, 在命令窗口输入: t=[0:1:16] y=[30.0 29.1 28.4 28.1 28.0 27.7 27.5 27.2 27.0 26.8 26.5 26.3 26.1 25.7 25.3 24.8 24.0] plot(t,y,'*')
a=polyfit(t,y,1) a= -0.3012 29.3804 hold on
beta
= 3.0022 4.0304 0.9404
(完整版)Matlab实验5选择结构程序结构
实验五、选择与循环结构一、实验目的：1、 掌握建立和执行M 文件的方法。
2、 掌握利用if 语句实现选择结构的方法。
3、 掌握利用switch 语句实现多分支选择结构的方法。
4、 掌握try 语句的使用。
5、 掌握利用for 语句实现循环结构的方法。
6、 掌握利用while 语句实现循环结构的方法。
7、 熟悉利用向量运算来代替循环的操作方法。
二、实验内容：1、 列分段函数的值。
⎪⎩⎪⎨⎧--≠≠<≤+--≠<-+=其他且且,632,100,6530,6222x x x x x x x x x x x y要求：(1) 用if 语句实现，分别输出x =-0.5,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。
提示：x 的值从键盘输入，可以是向量。
%homework_5_1_1.mx=input('请输入x 的值：x=');if (x<0 & x~=-3)y= x.*x + x - 6elseif (x>=0 & x<10 & x~=2 & x~=3)y=x.*x-5.*x+6elsey=x.*x-x-6end>> homework_5_1请输入x 的值：x=[-0.5 -3.0 1.0 2.0 2.5 3.0 5.0]y =-5.2500 6.0000 -6.0000 -4.0000 -2.2500 0 14.0000(2) 用逻辑表达式实现上述函数。
%homework_5_1_2.mx=input('请输入x 的值：x=')y=(x<0 & x~=-3).*(x.*x+x-6)...+(x>=0 & x<10 &x~=2 &x~=3).*(x.*x-5.*x+6)...+(x>=10 | x==-3 | x==3 | x==2).*(x.*x-x-6)>> homework_5_1_2请输入x=[-0.5 -3.0 1.0 2.0 2.5 3.0 5.0]x =-0.5000 -3.0000 1.0000 2.0000 2.5000 3.0000 5.0000 y =-6.2500 6.0000 2.0000 -4.0000 -0.2500 0 6.00002、输入一个百分制成绩，要求输出成绩的等级为A，B，C，D，E。
MATLAB上机实验报告
MATLAB上机实验报告实验一、实验目的--------------------------------------------------------------------------------------------------------3二、实验内容--------------------------------------------------------------------------------------------------------3三、实验步骤--------------------------------------------------------------------------------------------------------31. M 文件的编辑----------------------------------------------------------------------------------------------32. 程序流程控制结构---------------------------------------------------------------------------------------31) for 循环结构-----------------------------------------------------------------------------------------42) while 循环结构----------------------------------------------------------------------------------------5练习-----------------------------------------------------------------------------------------------------------53）if-else-end分支结构---------------------------------------------------------------------------------54）switc-case结构----------------------------------------------------------------------------------------63. 子函数和参数传递------------------------------------------------------------------------------------------6练习-----------------------------------------------------------------------------------------------------------74. 局部变量和全局变量--------------------------------------------------------------------------------------8 MATLAB上机练习6.1 M 文件--------------------------------------------------------------------------------------------8（1）脚本文件-----------------------------------------------------------------------------------------------------8（2）函数文件-----------------------------------------------------------------------------------------------------9（3）函数调用-----------------------------------------------------------------------------------------------------9 6.2 MATLAB 的程序控制结构( 1) for 循环----------------------------------------------------------------------------------------------------9(2) while 循环----------------------------------------------------------------------------------------------------12(3)if语句---------------------------------------------------------------------------------------------------------12(4)switch 语句-------------------------------------------------------------------------------------------------13(5)try 语句-----------------------------------------------------------------------------------------------------14(6)程序流的控制--------------------------------------------------------------------------------------------14 6.3 数据的输入与输出-----------------------------------------------------------------------------------14（1）键盘输入语句------------------------------------------------------------------------------------------14（2）屏幕输出语句---------------------------------------------------------------------------------------------14 6.4 MATLAB文件操作(1)fopen 语句-----------------------------------------------------------------------------------------------------15(2) fclose 语句-----------------------------------------------------------------------------------------------------15 6.5 面向对象编程--------------------------------------------------------------------------------------------15(1)创建类目录----------------------------------------------------------------------------------------------------15(2) 建立类的数据结构-----------------------------------------------------------------------------------------15(3)创建类的基本方法---------------------------------------------------------------------------------------------16 6.6 MATLAB 程序优化-----------------------------------------------------------------------------------------17 6.7程序调试------------------------------------------------------------------------------------------------------17 6.8 习题----------------------------------------------------------------------------------------------------------17实验五MATLAB 程序设计一、实验目的掌握MATLAB 程序设计的主要方法，熟练编写MATLAB 函数。
MATLAB课程设计
一、课程目的初步学会使用MATLAB软件，掌握MATLAB基本的函数的使用，熟练MATLAB编程，以便为今后数字信号的处理打下基础。
二、课程内容实验一 MATLAB语言基础实验目的：基本掌握 MATLAB 向量、矩阵、数组的生成及其基本运算(区分数组运算和矩阵运算)、常用的数学函数。
了解字符串的操作。
实验内容:1、创建以下矩阵：A为3×4的全 1 矩阵、B为3×3的0矩阵、C 为3×3的单位阵、D为3×3的魔方阵、E由C和D纵向拼接而成、F抽取E的2～5行元素生成、G由F经变形为3×4的矩阵而得、以G为子矩阵用复制函数(repmat)生成6×8的大矩阵H。
源程序：A=ones(3,4)% A为3×4的全 1 矩阵B=zeros(3)% B为3×3的0矩阵C=eye(3)% C为3×3的单位阵D=magic(3)% D为3×3的魔方阵E=cat(1,C,D)% E由C和D纵向拼接而成F=E(2:5,:)% F抽取E的2～5行元素生G=reshape(F,3,4)% G由F经变形为3×4的矩阵而得H=repmat(G,2,2) %以G为子矩阵用复制函数(repmat)生成6×8的大矩阵H运行结果：A =1 1 1 11 1 1 11 1 1 1B =0 0 00 0 00 0 0C =1 0 00 1 0 0 0 1 D =8 1 6 3 5 7 4 9 2 E =1 0 0 0 1 0 0 0 1 8 1 6 3 5 7 4 92 F =0 1 0 0 0 1 8 1 6 3 5 7 G =0 3 1 1 0 1 5 6 8 0 0 7 H =0 3 1 1 0 3 1 1 0 1 5 6 0 1 5 6 8 0 0 7 8 0 0 7 0 3 1 1 0 3 1 1 0 1 5 6 0 1 5 6 8 0 0 7 8 0 0 72、（1）用矩阵除法求下列方程组的解 x=[x 1；x 2；x 3]；⎪⎩⎪⎨⎧-=---=++-=++73847523436321321321x x x x x x x x x (2) 求矩阵的秩；(3) 求矩阵的特征值与特征向量； (4) 矩阵的乘幂与开方； (5) 矩阵的指数与对数； (6) 矩阵的提取与翻转。
《MATLAB程序设计》实验指导书
三、
实验仪器和设备
1、 计算机一台。 2、 MATLAB7.0 以上集成环境。
四、
预习要求
2
《MATLAB 程序设计》实验指导书
1、 复习 MATLAB 的启动与Байду номын сангаас出，熟悉 MATLAB 运行环境。 2、 复习 MATLAB 中矩阵的生成以及矩阵运算的基本原理。
五、
实验内容及步骤
实验内容：
1、 求下列表达式的值 1） z1
5
《MATLAB 程序设计》实验指导书
2） 函数调用 函数文件编制好后，就可调用函数进行计算了。函数调用的一般格式为 [输出实参表]=函数名（输入实参表） 注意：函数调用时各实参出现的顺序、个数，应与函数定义时形参的顺序、个 数一致，否则会出错。函数调用时，实参先传递给形参，然后再执行函数功能。 5、 选择结构 1） if 语句 a） 单分支 if 语句 if 条件 语句块 end b） 双分支 if 语句 if 条件 语句块 1 else 语句块 2 end c） 多分支 if 语句 if 条件 1 语句块 1 elseif 条件 2 语句块 2 …… elseif 条件 n 语句块 n else 语句块 n+1 end 2） switch 语句 switch 表达式 case 结果表 1 语句块 1 case 结果表 2 语句块 2 ……
二、
实验原理
1、 M 文件 用 MATLAB 语言编写的程序，称为 M 文件，它们的扩展名均为.m。M 文件根据 调用方式的不同分为两类，命令文件（Script file）和函数文件（Function file） 。 2、 建立新的 M 文件 启动 MATLAB 文本编辑器有 3 种方法： 1） 单击工具栏上的“New M-File”命令按钮。 2） 从 MATLAB 主窗口的“File”菜单中选择“New”菜单项，再选择“M-file”命令。 3） 在 MATLAB 命令窗口输入命令“edit”。 3、 打开已有的 M 文件 1） 在当前目录窗口选中要打开的 M 文件，双击鼠标左键。 2） 单击 MATLAB 主窗口工具栏上的“Open File”命令按钮，再从弹出的对话框中 选择所需打开的 M 文件。 3） 从 MATLAB 主窗口的“File”菜单中选择“Open”命令，在“Open”对话框中选中 所需打开的文件。 4） 在 MATLAB 命令窗口输入命令“edit 文件名” 。 4、 函数文件 1） 函数文件的基本结构 函数文件由 function 语句引导，其基本结构为： function 输出形参表=函数名（输入形参表） 注释说明部分 函数体语句 注意：函数名的命名规则与变量名相同。当输出形参多于一个时，应用方括号 括起来。
MATLAB程序设计与应用课后实验答案
A-B+eye(3)
A*B
A.*B
A^3
A.^3
A/B
B\A
[A,B]
[A([1,3],:);B^2]
运算结果：
A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7];
A+6.*B
A-B+eye(3)
A*B
A.*B
A^3
A.^3
A/B
B\A
0 0 0 0
ans =
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
由ans,所以
2. 产生5阶希尔伯特矩阵H和5阶帕斯卡矩阵P，且求其行列式的值Hh和Hp以及它们的条件数Th和Tp，判断哪个矩阵性能更好。为什么？
解：M文件如下：
输出结果：
H =
P =
1 1 1 1 1
Columns 57 through 60
Column 61
z4 =
0 0
2. 已知：
求下列表达式的值：
(1) A+6*B和A-B+I（其中I为单位矩阵）
(2) A*B和A.*B
(3) A^3和A.^3
(4) A/B及B\A
(5) [A,B]和[A([1,3],:);B^2]
解：
M文件：
A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7];
实验四 循环结构程序设计
一、实验目的
1. 掌握利用for语句实现循环结构的方法。
2. 掌握利用while语句实现循环结构的方法。
实验五---M文件的编写与调试(答案)
MATLAB实验报告
专业班级
姓名
学号
成绩评定
电气与信息学院
和谐勤奋求是创新
实验教学考核和成绩评定方法
1．课内实验考核成绩，严格按照该课程教学大纲中明确规定的比重执行。
实验成绩不合格者，不能参加课程考试，待补做合格后方能参加考试。
2．单独设立的实验课考核按百分制评分，考核内容应包括基本理论、实验原理和实验。
3．实验考核内容包括：1〕实验预习；2〕实验过程〔包括实验操作、实验记录和实验态度、表现〕；3〕实验报告；权重分别为0.2 、0.4 、0.4；原则上根据上述三个方面进行综合评定。
学生未取得1〕和2〕项成绩时，第3〕项成绩无效。
4．实验指导教师应严格按照考核内容分项给出评定成绩，并及时批改实验报告，给出综合成绩，反馈实验中出现的问题。
实验成绩在教师手册中有记载。
实验报告主要内容
一．实验目的
二．实验仪器及设备
三．实验原理
四．实验步骤
五．实验记录及原始记录
六．数据处理及结论
七. 思考题
八．实验体会〔可选项〕
注：1. 为了节省纸张，保护环境，便于保管实验报告，统一采用A4纸，实验报告建议双面打印〔正文采用宋体五号字〕或手写,右侧装订。
2. 实验类别指验证、演示、综合、设计、创新〔研究〕、操作六种类型实验。
3. 验证性实验：是指为了使学生稳固课程基本理论知识而开设的强调演示和证明，注重实验结果〔事实、
概念或理论〕的实验。
4. 综合性实验：是指实验内容涉及本课程的综合知识或本课程相关的课程知识的实验。
5. 设计性实验：是指给定实验目的、要求和实验条件，由学生自行设计实验方案并加以实现的实验。
实验五 matlab基础知识(简单)
本次实验注意：《实验五MALTAB基础知识（简单）》《实验五基于Matlab的信号频谱分析（复杂）》选作一个即可实验五MALTAB基础知识（一）实验目的 (2)（二）实验设备 (2)（三）实验要求 (2)（四）实验内容 (2)1.1 MATLAB基础知识 (2)1.1.1 MATLAB程序设计语言简介 (2)1.1.2 MA TLAB界面及帮助 (2)1.2 MA TLAB基本运算 (4)1.2.1 MA TLAB内部特殊变量和常数 (4)1.2.2 变量类型 (4)1.2.3 内存变量管理 (5)1.2.4 MA TLAB常用数学函数 (5)1.2.5 MA TLAB矩阵生成 (5)1.2.6 MA TLAB矩阵运算 (8)1.2.7 MA TLAB中的矩阵分析 (10)1.3 MA TLAB程序设计 (10)1.3.1 M文件 (10)1.3.2 程序控制结构 (12)实验五MALTAB基础知识（一）实验目的●了解MA TLAB 程序设计语言的基本特点，熟悉MA TLAB软件运行环境●掌握创建、保存、打开m文件及函数的方法●掌握变量等有关概念，具备初步的将一般数学问题转化为对应的计算机模型并进行处理的能力（二）实验设备计算机，Matlab软件（三）实验要求本实验属于验证实验，请根据（四）实验内容的步骤，运行相应的指令或例子，并将仿真结果截图至文档（请自己新建一个word文档，注意，并不一定所有指令或例子的实验结果都要截图，截图数目大于等于5个即可，自己选择性截图，答案不唯一，自由发挥）请在页眉处填写班级、学号、姓名，并将实验报告命名为“实验五_学号_姓名”，并通过FTP上传至指定文件夹。
（四）实验内容1.1 MATLAB基础知识1.1.1 MATLAB程序设计语言简介MA TLAB，Matrix Laboratory的缩写，是由MathWorks公司开发的一套用于科学工程计算的可视化高性能语言，具有强大的矩阵运算能力。
matlab程序课程设计
matlab程序课程设计一、教学目标本课程的教学目标是使学生掌握MATLAB编程的基本知识和技能，能够运用MATLAB解决简单的数学和工程问题。
具体目标如下：1.理解MATLAB的基本概念，如变量、矩阵、数组等。
2.掌握MATLAB的基本语法，如运算符、函数、循环和条件语句等。
3.了解MATLAB的绘图功能，能够绘制基本的图形。
4.能够使用MATLAB进行简单的数学计算和数据分析。
5.能够编写MATLAB脚本程序，解决实际问题。
6.能够使用MATLAB的绘图功能，进行数据可视化。
情感态度价值观目标：1.培养学生的计算机编程思维，提高解决问题的能力。
2.培养学生团队合作的精神，提高沟通与协作能力。
3.培养学生对科学研究的兴趣，提高创新意识。
二、教学内容根据教学目标，本课程的教学内容主要包括以下几个方面：1.MATLAB概述：介绍MATLAB的发展历程、功能特点和应用领域。
2.MATLAB基本语法：变量、矩阵、数组、运算符、函数、循环和条件语句等。
3.MATLAB绘图功能：基本图形绘制、图形编辑和格式设置等。
4.MATLAB编程实践：数学计算、数据分析、实际问题解决等。
5.第1周：MATLAB概述和基本语法。
6.第2周：MATLAB绘图功能。
7.第3周：MATLAB编程实践（数学计算和数据分析）。
8.第4周：MATLAB编程实践（实际问题解决）。
三、教学方法为了达到教学目标，本课程将采用以下教学方法：1.讲授法：讲解MATLAB的基本概念、语法和绘图功能。
2.案例分析法：分析实际问题，引导学生运用MATLAB编程解决。
3.实验法：上机操作，让学生亲手编写MATLAB程序，巩固所学知识。
4.小组讨论法：分组完成项目任务，培养团队合作和沟通能力。
四、教学资源为了支持教学内容和教学方法的实施，本课程将采用以下教学资源：1.教材：《MATLAB程序设计》。
2.参考书：提供一些MATLAB编程的参考书籍，供学生自主学习。
实验五-MATLAB的程序设计
实验五MATLAB的程序设计一、实验目的1.熟悉MATLAB 的编程环境；2.熟悉和掌握脚本文件和函数文件建立与使用方法；3.掌握MATLAB程序设计的流程控制语句；4.掌握函数调用和参数传递；5.熟悉MATLAB 程序调试方法。
二、实验仪器1.计算机2.MATLAB 软件环境三、实验指导1. M 文本文件介绍一个M 文件通常包含五个部分，函数定义语句、H1 帮助行、帮助文本、函数体或者脚本文件语句和注释语句。
一个完整的M 文件的结构为：function f = fact(n) 函数定义语句% Compute a factorial value. H1 行% FACT(N) returns the factorial of N, 帮助文本% usually denoted by N!% Put simply, FACT(N) is PROD(1:N). 注释语句f = prod(1:n); 函数体函数定义语句只在函数文件中存在，定义函数名称、输入输出参数的数量和顺序。
函数定义语句的格式为：function y = fact(x)2. 函数流程控制顺序结构顺序结构是最简单的程序结构，系统在编译程序时，按照程序的物理位置顺序执行。
这种程序容易编制，但是结构单一，能够实现的功能有限。
判断语句（if…else…end）if…endif…else…endif…elseif…else…end1）if…end此时的程序结构如下：if 表达式执行代码块end这是最简单的判断语句，只有一个判断语句，其中的表达式为逻辑表达式，当表达式为真时，执行相应的语句，否则，直接跳到下一段语句。
2）if…else…end当程序有两个选择时，可以选择if…else…end 结构，此时程序结构为：if 表达式执行代码块1else执行代码块2end当判断表达式为真时，执行代码块1，否则执行代码块2。
3）if…elseif…else…end上面的两种形式中，分别包含一个选择和两个选择，当判断包含多个选择时，可以采用elseif 语句，结构为：if 表达式1执行代码块1elseif 表达式2执行代码块2elseif ..................else执行代码块end其中可以包含任意多个elseif 语句。
MATLAB程序设计与应用课后实验答案
实验一MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin851z e =+(2)21ln(2z x =，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+=++=-- (4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.52. 已知：求下列表达式的值：(1) A+6*B和A-B+I（其中I为单位矩阵）(2) A*B和A.*B(3) A^3和A.^3(4) A/B及B\A(5) [A,B]和[A([1,3],:);B^2]解：运算结果：3. 设有矩阵A和B(1) 求它们的乘积C。
(2) 将矩阵C的右下角3×2子矩阵赋给D。
(3) 查看MATLAB工作空间的使用情况。
解：. 运算结果：4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。
(2) 建立一个字符串向量，删除其中的大写字母。
解：(1) 结果：(2). 建立一个字符串向量 例如：ch='ABC123d4e56Fg9';则要求结果是：实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A O S +⎡⎤=⎢⎥⎣⎦。
解: M 文件如下；输出结果：由ans,所以22E R RS A OS +⎡⎤=⎢⎥⎣⎦ 2. 产生5阶希尔伯特矩阵H 和5阶帕斯卡矩阵P ，且求其行列式的值Hh 和Hp 以及它们的条件数Th 和Tp ，判断哪个矩阵性能更好。
为什么？解：M 文件如下：3. 建立一个5×5矩阵，求它的行列式值、迹、秩和范数。
- 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}。