MATLAB实验四循环结构程序的设计说明
(完整版)MATLAB)课后实验[1]
![(完整版)MATLAB)课后实验[1]](https://img.taocdn.com/s3/m/88685371e53a580217fcfe63.png)
实验一 MATLAB 运算基础1. 先求以下表达式的值，尔后显示 MATLAB 工作空间的使用情况并保存全部变量。
(1)2sin 85 z1 21 e(2) 12z ln( x 1 x ) ，其中22 x2 1 2i5(3)ae e az sin( a 0.3) ln , a 3.0, 2.9, L , 2.9, 32 22t 0 t 1(4) 2z t 1 1 t 242t 2t 1 2 t 3，其中解：M 文件:z1=2*sin(85*pi/180)/(1+exp(2))x=[2 1+2*i;-.45 5];z2=1/2*log(x+sqrt(1+x^2))a=-3.0:0.1:3.0;z3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2)t=0:0.5:2.5;z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1)4. 完成以下操作：(1) 求[100,999] 之间能被 21 整除的数的个数。
(2) 建立一个字符串向量，删除其中的大写字母。
解：(1) 结果：m=100:999;n=find(mod(m,21)==0);length(n)ans =43(2). 建立一个字符串向量比方：ch='ABC123d4e56Fg9'; 那么要求结果是：ch='ABC123d4e56Fg9';k=find(ch>='A'&ch<='Z');ch(k)=[]ch =123d4e56g9实验二 MATLAB矩阵解析与办理1. 设有分块矩阵 A E R3 3 3 2O S2 3 2 2，其中 E、R、O、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试经过数值计算考据 2A E R RS2O S。
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编程【实验目的】1、了解Matlab运行方式2、熟悉Matlab程序流程与结构3、会编写简单的Matlab程序4、掌握M文件编辑器及程序的调试工具【实验内容】1、关系运算与逻辑运算2、m函数文件（function file）和m命令文件（script file）以及m文件的调试3、分支结构和循环结构4、相关语句或函数5、根据要求，完成实验报告中的内容【实验指导】1)关系运算与逻辑运算关系运算符>(大于) <(小于) >=(大于等于) <=(小于等于) ==(等于) ~=(不等于) 逻辑运算符& (与) |(或) ~(非) xor(异或)% 逻辑结果以“1”或者“非零数字”表示真, 以“0”表示假。
测试函数2)m函数文件（function file）基本数学函数自定义函数：以m文件形式出现✓格式：function [y1, y2, ..., yn]=fun_name(x1, x2, ..., xn)✓说明：◆函数名必须符合Matlab变量命名规则。
◆m文件名必须与函数名相同。
◆输出参数只有一个或者没有时，可以省略方括号。
◆输入参数可为0、1、2、...个，但不能省略圆括号。
◆M函数文件第一行必须以funtion开头。
◆函数有自已的工作空间，它与Matlab系统本身的工作空间是分开的。
◆当m文件到终点或者遇到return时结束调用。
◆函数可以调用其他函数，也可以调用自己。
◆一个m文件中可以定义多个函数，第一个出现的为主函数，其他的为子函数。
注意文件名与主函数名相同，子函数只能被主函数调用，而不能被该m文件外的函数或者m文件调用。
✓调用：[y1, y2, ..., yn]=fun_name(x1, x2, ..., xn)3)m命令文件（script file）该文件仅是一连串Matlab指令的集合，执行时按顺序执行。
它没有输入参数与输出参数，中间可以调用函数。
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. 已知：1234413134787,2033657327A B --⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求下列表达式的值：(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 和B1234530166789101769,111213141502341617181920970212223242541311A 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 OS +⎡⎤=⎢⎥⎣⎦。
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:: 解：2. 已知：1234413134787,2033657327A B --⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求下列表达式的值：(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 和B1234530166789101769,111213141502341617181920970212223242541311A 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 +⎡⎤=⎢⎥⎣⎦。
MATLAB)课后实验答案
实验一 MATLAB 运算基础1、 先求下列表达式得值，然后显示MATLAB 工作空间得使用情况并保存全部变量。
(1) 0122sin 851z 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 --+=++=--L (4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0、5:2、5 解：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 OS +⎡⎤=⎢⎥⎣⎦。
解: M 文件如下；5、 下面就是一个线性方程组：1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦(1) 求方程得解。
(2) 将方程右边向量元素b 3改为0、53再求解，并比较b 3得变化与解得相对变化。
(3) 计算系数矩阵A 得条件数并分析结论。
解： M 文件如下：实验三 选择结构程序设计1、 求分段函数得值。
2226035605231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且及其他用if 语句实现，分别输出x=-5、0,-3、0,1、0,2、0,2、5,3、0,5、0时得y 值。
实验报告4
实验四：循环结构程序设计班级：学生姓名：学号：一、实验目的1、理解循环的概念2、理解并掌握循环结构相关语句的含义、格式及使用3、学会循环的应用及控制，包括：①掌握使用循环输入多个数据的方法②掌握在多个数据中有选择地输出数据的方法③掌握在多个数据中对某种数据进行计数的方法④掌握求多个数据中最大值、最小值的方法⑤掌握使用break、continue语句终止循环4、掌握循环的嵌套二、知识要点1、循环变量、循环条件、循环体的概念2、三种循环语句的一般格式、执行过程3、理解选择结构与循环结构中“条件表达式”之不同含义4、二重循环的执行过程三、实验预习（要求做实验前完成）1、循环变量的主要用途是：2、用循环求多个数的和之前，先要把和的初始值赋为：3、用循环求多个数的乘积之前，先要把乘积的初始值赋为：4、字符变量能否作为循环变量？5、循环过程中，如果循环条件成立，但需要结束循环，可采用什么办法？6、什么叫循环的嵌套？四、实验内容(要求提供：①算法描述或流程图②源程序)1. 编程，利用循环计算以下表达式的值：（5＋52）*(4+42)*(3+32)*(2+22)*(1+12)*(1/2+1/3+1/4+1/5)（for循环）include <stdio.h>#include <stdlib.h>int main(){int a;double sum=1,sum1=0;for(a=1;a<=5;a++)sum=sum*(a+a*a);printf("结果为%lf\n",sum);for(a=2;a<=5;a++)sum1=sum1+(1.0/a);printf("%lf\n",sum1);printf("结果为%lf\n",sum*sum1);return 0;}2. 编程，从键盘输入若干个整数，当输入0时，输入停止。
matlab实验4_循环结构程序设计_参考解答
高教社刘卫国《MATLAB 程序设计与应用》（第二版）实验参考解答
实验四：循环结构程序设计
clear all; close all; clc; y=0; n=1; while y<3
y=y+1/(2*n-1); n=n+1; end n=n-1; y=y-1/(2*n-1) n=n-1
运行结果：
y= 2.9944
x1=a/(b+x0); y=abs(x1-x0); x0=x1;
MATLAB 语言
3
WangChengyou © 山东大学威海分校
高教社刘卫国《MATLAB 程序设计与应用》（第二版）实验参考解答
实验四：循环结构程序设计
n=n+1; end x1 r1=(-b+sqrt(b^2+4*a))/2 r2=(-b-sqrt(b^2+4*a))/2
clear all; close all; clc; prime=[]; for n=2:49
x=n*(n+1)-1; signx=0; for k=2:sqrt(x)
MATLAB 语言
5
WangChengyou © 山东大学威海分校
高教社刘卫国《MATLAB 程序设计与应用》（第二版）实验参考解答
根据
π2 6
=1 12
+
1 22
+1 32
+"+
1 n2
，求 π 的近似值。当
n
分别取
100、1000、10000
时，
结果是多少？
要求：分别用循环结构和向量运算（使用 sum 函数）来实现。 答：程序设计：
(1) 循环结构
matlab实验 循环结构
实验（四）项目名称：循环结构一、实验目的：1. 掌握利用for语句实现循环结构的方法。
2. 掌握利用while语句实现循环结构的方法。
3. 熟悉利用向量运算来代替循环操作的方法。
二、实验原理1.FOR 循环在for和end语句之间的{commands}按数组中的每一列执行一次。
在每一次迭代中，x被指定为数组的下一列，即在第n次循环中，x=array(:, n)。
2.WHILE循环只要在表达式里的所有元素为真，就执行while和end 语句之间的{commands}。
通常，表达式的求值给出一个标量值，但数组值也同样有效。
在数组情况下，所得到数组的所有元素必须都为真。
三、实验环境1.硬件：PC机2. 软件：Windows操作系统、matlab2015四、实验内容、步骤以及结果4.1.1实验要求：用while语句写一个程序，k=5，每循环一次，自动减1，并自动输出。
4.1.2实验步骤（1）启动matlab，新建一个M文件；（2）输入程序，如图1；（3）保存文件；（4）编译源程序，观察屏幕上显示的编译信息，修改出现的错误，直到编译成功；图1：程序输入4.1.3运行结果如下：图2：运行结果4.2.1实验要求：1、根据，求π 的近似值。
当n分别取100,1000,10000时，结果是多少？要求：分别用循环结构和向量运算（使用sum 函数）来实现4.2.2实验步骤：（1）启动matlab ，M 文件；（2）输入程序，如图3和4；（3）保存文件；4.2.3运行结果如下： 图3:for 循环结构图4:向量运算图5：运行结果五、实验总结通过本次实验，我慢慢理解了循环控制结构的执行过程，掌握了一些循环结构程序设计方法，但还不是很熟练。
我了解了一些MATLAB的编程方法和循环结构程序设计方法，但都不是特别的熟练，都有待加强。
MATLAB实验
MATLAB 实验报告班级:14通信1班 学号:201424124124 姓名:林启铭实验一 MATLAB 运算基础(一)一、实验目的1、掌握建立矩阵的方法。
2、掌握MATLAB 各种表达式的书写规则以及各种运算方法。
二、实验内容1、求下列表达式的值。
（1）20185sin 21ez += MATLAB 代码：z1=2*sin(85*pi/180)/(1+exp(2))%将角度化为弧度z1 =0.2375（2）()x x z ++=1ln 212，其中⎢⎣⎡-=45.02x ⎥⎦⎤+521i (可以分别对矩阵和元素群运算)MATLAB 代码：>> x=[2,1+2i;-0.45,5]x =2.0000 + 0.0000i 1.0000 + 2.0000i-0.4500 + 0.0000i 5.0000 + 0.0000i>> z2=1/2*log(x+sqrt(1+x))z2 =0.6585 + 0.0000i 0.6509 + 0.4013i-0.6162 + 0.0000i 1.0041 + 0.0000i（3）()3.0sin 232.03.0+⋅-=a e e z aa ， 0.3,9.2,8.2,...,8.2,9.2,0.3---=a （结果请用图形表示）（提示：利用冒号表达式生成a 向量；求各点的函数值时用点乘运算）MATLAB 代码：>> a=-3.0:0.1:3.0;%利用冒号表达式生成a 向量，加分号结尾避免大量数据刷屏 >> z3=(exp(0.3.*a)-exp(0.2.*a)).*sin(a+0.3)/2;>> plot(a,z3);%绘制出以a 为自变量，z3为因变量的曲线>>曲线图：2、已知⎢⎢⎢⎣⎡=33412A 65734⎥⎥⎥⎦⎤-7874 和 ⎢⎢⎢⎣⎡=321B 203-⎥⎥⎥⎦⎤-731 求下列表达式的值：（1）A+6*B 和A-B+I （其中I 为单位矩阵）。
matlab实验指导答案详解（非常详细正确）
matlab实验指导答案详解（⾮常详细正确）实验⼀ MATLAB ⼯作环境熟悉及简单命令的执⾏⼀、实验⽬的：熟悉MATLAB 的⼯作环境，学会使⽤MATLAB 进⾏⼀些简单的运算。
⼆、实验内容：MATLAB 的启动和退出，熟悉MATLAB 的桌⾯（Desktop ），包括菜单（Menu ）、⼯具条（Toolbar ）、命令窗⼝(Command Window)、历史命令窗⼝、⼯作空间(Workspace)等；完成⼀些基本的矩阵操作；学习使⽤在线帮助系统。
三、实验步骤：1、启动MATLAB ，熟悉MATLAB 的桌⾯。
2、在命令窗⼝执⾏命令完成以下运算，观察workspace 的变化，记录运算结果。
（1）（365-52?2-70）÷3 >>(365-52*2-70)/3 ans = 63.6667（2）>>area=pi*2.5^2 area = 19.6350（3）已知x=3，y=4，在MATLAB 中求z ：()232y x y x z -= >>x=3 >>y=4>>z = x ^2 * y ^3 / (x - y) ^2 z = 576（4）将下⾯的矩阵赋值给变量m1,在workspace 中察看m1在内存中占⽤的字节数。
m1=11514412679810115133216 执⾏以下命令>>m1 =[16 2 3 13 ; 5 11 10 8 ; 9 7 6 12 ; 4 14 15 1 ] >>m1( 2 , 3 ) ans = 10 >>m1( 11 ) ans = 6>>m1( : , 3 ) ans =3 10 6 15>>m1( 2 : 3 , 1 : 3 ) ans =5 11 10 9 7 6>>m1( 1 ,4 ) + m1( 2 ,3 ) + m1( 3 ,2 ) + m1( 4 ,1) ans = 34（5）执⾏命令>>help abs查看函数abs 的⽤法及⽤途，计算abs( 3 + 4i ) （6）执⾏命令>>x=0:0.1:6*pi; >>y=5*sin(x); >>plot(x,y)（6）运⾏MATLAB 的演⽰程序，>>demo ，以便对MATLAB 有⼀个总体了解。
MATLAB程序设计及应用(第二版)课后实验答案
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:: 解：2. 已知：1234413134787,2033657327A B --⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求下列表达式的值：(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 和B1234530166789101769,111213141502341617181920970212223242541311A 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 +⎡⎤=⎢⎥⎣⎦。
《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全部实验及答案
MATLAB全部实验及答案MATLAB全部实验及答案实验一、MATLAB基本操作实验内容及步骤4、有关向量、矩阵或数组的一些运算（1）设A=15；B=20；求C=A+B与c=a+b？（2）设A=[1 2 3;4 5 6;7 8 9]，B=[9 8 7;6 5 4;3 2 1]；求A*B 与A.*B？A*B就是线代里面的矩阵相乘A.*B是对应位置的元素相乘（3）设a=10，b=20；求i=a/b=0.5与j=a\b=2？（4）设a=[1 -2 3;4 5 -4;5 -6 7]请设计出程序，分别找出小于0的矩阵元素及其位置（单下标、全下标的形式），并将其单下标转换成全下标。
clear,clca=[1 -2 3;4 5 -4;5 -6 7];[x,y]=find(a<0);c=[];for i=1:length(x)c(i,1)=a(x(i),y(i));c(i,2)=x(i);c(i,3)=y(i);c(i,4)=(y(i)-1)*size(a,2)+x(i);endc（5）在MATLAB命令行窗口运行A=[1,2;3,4]+i*[5,6;7,8]；看结果如何？如果改成运行A=[1,2;3,4]+i[5,6;7,8]，结果又如何？前面那个是虚数矩阵，后面那个出错（6）请写出完成下列计算的指令：a=[1 2 3;3 4 2;5 2 3]，求a^2=？,a.^2=？a^2= 22 16 1625 26 2326 24 28a.^2=1 4 99 16 425 4 9（7）有一段指令如下，请思考并说明运行结果及其原因clearX=[1 2;8 9;3 6];X( : ) 转化为列向量（8）使用三元组方法，创建下列稀疏矩阵2 0 8 00 0 0 10 4 0 06 0 0 0方法一：clear,clcdata=[2 8 1 4 6];ir=[1 1 2 3 4 ];jc=[1 3 4 2 1];s=sparse(ir,jc,data,4,4);full(s)方法二：不用三元组法clear,clca=zeros(4,4);a(1,[1,3])=[2,8];a(2,4)=1;a(3,2)=4;a(4,1)=6;a（9）写出下列指令的运行结果>> A = [ 1 2 3 ]; B = [ 4 5 6 ];>> C = 3.^A>> D = A.^B5、已知+?=-334sin 234πt e y t 若需要计算t ∈[-1,1]，取间隔为0.01，试计算出相对应的y 值。
MATLAB实验四 选择结构和循环结构程序设计
实验四 选择结构和循环结构程序设计实验目的：1.掌握建立和执行M 文件的方法 2.掌握利用if 语句实现选择结构的方法 3.掌握利用switch 语句实现多分支选择结构的方法 4.掌握利用for 语句实现循环结构的方法 5.掌握利用while 语句实现循环结构的方法 6. 熟悉利用向量运算来代替循环操作的方法实验内容：1. 求分段函数的值。
3及2其它,1x 且50,653且0,6222≠⎪⎩⎪⎨⎧--≠<≤+--≠<-+=x x x x x x x x x x y用if 语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。
x=input('请输入x 的值:x=');if(x<0&x~=-3)y=x.*x+x-6elseif(x>=0&x<5&x~=2&x~=3)y=x.*x-5.*x+6elsey=x.*x-x-1end请输入x 的值:x=[-5.0 -3.0 1.0 2.0 2.5 3.0 5.0]y =29.0000 11.0000 -1.0000 1.0000 2.7500 5.0000 19.00002. 输入一个百分制成绩，要求输出成绩等级A 、B 、C 、D 、E 。
其中90分~100分为A ，80分~89分为B ，70分~79分为C ，60分~69分为D ，60分以下为E 。
要求：a) 分别用if 语句和switch 语句实现。
If 语句x=input('请输入分数：')if(x<=100&x>=90)disp('A')elseif(x>=80&x<=89)disp('B')elseif(x>=70&x<=79)disp('C')elseif(x>=60&x<=69)disp('D')elseif(x<60)disp('E')end请输入分数：95x =95Aswitch语句c=input('请输入分数:')switch ccase num2cell(90:100)disp('A')case num2cell(80:89)disp('B')case num2cell(70:79)disp('C')case num2cell(60:69)disp('D')case num2cell(0:59)disp('E')end请输入分数:83c =83Bb)输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。
实验4答案 Matlab程序设计1
实验4 Matlab 程序设计1实验目的：1、 掌握建立和执行M 文件的方法；2、 掌握实现选择结构的方法；3、 掌握实现循环结构的方法。
实验内容：1. 从键盘输入一个4位整数，按如下规则加密后输出。
加密规则：每位数字都加上7，然后用和除以10的余数取代该数字；再把第一位与第三位交换，第二位与第四位交换。
a = input('请输入一个四位整数:');a1 = fix(a/1000);a2 = rem(fix(a/100),10);a3 = rem(rem(fix(a/10),100),10);a4 = rem(a,10);a1 = rem(a1+7,10)a2 = rem(a2+7,10)a3 = rem(a3+7,10)a4 = rem(a4+7,10)b1 = a1;a1 = a3;a3 = b1;b2 = a2;a2 = a4;a4 = b2;b = a1*1000+a2*100+a3*10+a42. 求分段函数的值。
,x x x x y x x x x x x x ⎧+- <≠-⎪=-+ ≤<≠≠⎨⎪-- ⎩2226035605231且且及其他用if 语句实现，分别输出x=-5,-3,0,1,2,2.5,3,5时的y 值。
x = input(请输入X 的值:');if x <0 & x~=-3y = x.*x +x-6;elseif x>=0 & x<5 & x~=2 & x~=3y = x.*x-5*x+6;elsey = x.*x-x-1;endy3.输入一个百分制成绩，要求输出成绩等级A、B、C、D、E，其中90~100分为A，80~89分为B，70~79分为C，60~69分为D，60分以下为E。
要求：（1）分别用if语句和swich语句实现。
（2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。
实验4___MATLAB程序设计__4学时
程序2： 程序 ： 首先建立函数文件f2c.m。 首先建立函数文件 。 function c=f2c(f) c=5*(f-32)/9 然后在MATLAB的命令窗口调用该函数文件。 的命令窗口调用该函数文件。 然后在 的命令窗口调用该函数文件 clt temperature：'); ： x=f2c(y) 输出情况为： 输出情况为： Input Fahrenheit temperature：70 ： c= 21.1111 x= 21.1111
4.2.3 循环结构 1．for语句 ． 语句 for语句的格式为： 语句的格式为： 语句的格式为 for 循环变量 表达式 表达式 表达式 循环变量=表达式 表达式2:表达式 表达式1:表达式 表达式3 循环体语句 end 其中表达式1的值为循环变量的初值 表达式2的值 的值为循环变量的初值， 其中表达式 的值为循环变量的初值，表达式 的值 为步长，表达式3的值为循环变量的终值 步长为1时 的值为循环变量的终值。 为步长，表达式 的值为循环变量的终值。步长为 时， 表达式2可以省略 可以省略。 表达式2可以省略。
2．switch语句 ． 语句 switch语句根据表达式的取值不同，分别执行不同 语句根据表达式的取值不同， 语句根据表达式的取值不同 的语句，其语句格式为： 的语句，其语句格式为： switch 表达式 case 表达式 表达式1 语句组1 语句组 case 表达式 表达式2 语句组2 语句组 …… case 表达式 表达式m 语句组m 语句组 otherwise 语句组n 语句组 end
例4-8 一个三位整数各位数字的立方和等于该数本身则 称该数为水仙花数。输出全部水仙花数。 称该数为水仙花数。输出全部水仙花数。 程序如下： 程序如下： for m=100:999 m1=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 disp(m) end end 输出的水仙花数为： 输出的水仙花数为：153、370、371、407 、 、 、
(完整版)MATLAB)课后实验答案[1]
![(完整版)MATLAB)课后实验答案[1]](https://img.taocdn.com/s3/m/d27ebaa958fb770bf78a55b8.png)
实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin 851z 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 --+=++=--L (4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.5 解：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 文件如下；5. 下面是一个线性方程组：1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦(1) 求方程的解。
(2) 将方程右边向量元素b 3改为0.53再求解，并比较b 3的变化和解的相对变化。
(3) 计算系数矩阵A 的条件数并分析结论。
解： M 文件如下： 123d4e56g9实验三 选择结构程序设计1. 求分段函数的值。
2226035605231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且及其他用if 语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。
MATLAB程序设计与应用(第二版)课后实验答案
Matlab课后实验题答案实验一 MATLAB运算基础1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。
(1)0 122sin851ze =+(2)21ln( 2z x=+，其中2120.455i x+⎡⎤=⎢⎥-⎣⎦(3)0.30.330.3sin(0.3)ln, 3.0, 2.9,,2.9,3.0 22a ae e az a a--+=++=--(4)2242011122123t tz t tt t t⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t=0:0.5:2.52. 已知：1234413134787,2033657327A B --⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求下列表达式的值：(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 和B123453166789101769,111213141502341617181920970212223242541311A 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 OS +⎡⎤=⎢⎥⎣⎦。
程序设计基础实验报告4-循环结构
程序设计基础实验指导及实验报告专业班级：学号：姓名：时间：指导老师：魏银珍实验四循环结构程序设计一、实验目的1．熟练掌握三种循环语句的应用。
2．熟练掌握循环结构的嵌套。
3．掌握break和continue语句的使用。
4．练习调试与修改程序。
二、实验内容1、读下面程序并运行，结果是什么？2、读下面程序并运行，改正其中的错误。
main( ) main( ){{int m, k ,I , n =0; int s ;for(m=101;m<=200;m=m+2) float n ,t ,pi ;{ t=1 ;pi=0 ;n=1.0 ;s=-1 ;k=sqrt(m) ; while(fabs(t)>1e+6)for(i=2;i<=k; i++) {pi=pi+t ;if(m%i= =0)break; n=n+2 ;if(i>=k+1) s=-s ;{ t=s/n ;printf(“%d”,m); }n=n+1; pi=pi*4 ;} printf(“pi=%10.6d\n”, pi);if(n%10= =0)printf(“\n”); }}printf(“\n”);}3．用while语句、do-while语句编写程序，计算e≈1+1/1!+1/2!+...+1/n!。
使误差小于给定的ε=10-5。
输出e和总的项数n。
4、如果一个数恰好等于它的所有因子（包括1但不包括自身）之和，则称之为“完数”。
例如6的因子是1，2，3，且6=1+2+3，因此6是完数。
编写程序找出1000之内的所有完数及这些完数的和。
5．输入并运行下面的程序，观察程序的运行结果。
#include"stdio.h"main( ){int n;while (1){printf("Enter a number:");scanf("%d",&n);if (n%2==1){printf("I said");continue;}break;}printf("Thanks. I needed that!");6. 打印九九乘法表，用for语句实现1X1=11X2=2 2X2=4……1X9=9 …… 9X9=81教材 6.6，6.147. 5.提高题：打印出所有的“水仙花数”，所谓“水仙花数”是指一个3位数，其各位数字立方之和等于该数本身。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
程序设计：
x=1.043;x);
n=1;
while(abs(x-(a/x-b))>10e-5&n
v=500)
x=a/(b+x);
n=n+1;
end
x
r1=(-b+sqrt(bA2+4*a))/2
r2=(-b-sqrt(bA2+4*a))/2
运行结果：
x =
2
1.根据一丄A22,
6
求的近似值。当n分别取100、
1000、10000时，记过是多少？
n=100
程序设计：
n=100;
i=1: n;
f=1./i.A2;
x=sum⑴；
y=sqrt(6*x)
运行结果：
x =
3.1321
n=1000
程序设计：
n=1000;
i=1: n;
f=1./i.A2;
x=sum(f);
y=sqrt(6*x)
运行结果：
3.1406
n=10000
程序设计：
n=1000;
i=1: n;
f=1./i.A2;
x=sum(f);
y=sqrt(6*x)
运行结果：
x =
3.1415
2.根据y=1+--丄，求:
3 52n-1
(1)y<3时的最大n值。
(2)与(1)的n值对应的y值 程序设计：
y=1;n=1;
while(y<3)
n=n+2;
y=y+1/ n;
end
y=y-1/ n
n=( n+1)/2
运行结果:
y =
2.9944
57
3.
考虑以下迭代公式:
Xnl
a
bXn
其中a、b为正的常数。
（1）编写程序求迭代的结果，迭
代的终止条件为
Xn
5
10
迭代初值Xo=1.0，迭代次数不超
过500次
程序设计:
x=1.0;
运行结果：
x =
1.7016
r1 =
1.7016
r2 =
-4.7016
(a,b)取(10,0.1)
程序设计：
x=1.0;
a=10;
b=0.1;
x=a/(b+x);
a=2;
b=3; x=a/(b+x);
n=1;
while(abs(x-(a/x-b))>
10e-5&n v=500) x=a/(b+x); n=n+1;
end
x
n
运行结果:
x =
0.5616
n二
6
（2）如果迭代过程收敛于r，那
/~2
么r的准确值是b4a，当
2
（a,b）的值取（1,1）、（8,3）、
（10,0.1）时，分别对迭代结果 的准确值进行比较。
0.6181
r1 =
0.6180
r2 =
-1.6180
(a,b)取(8,3)
程序设计：
x=1.0;
a=8;
b=3;
x=a/(b+x);
n=1;
while (abs(x-(a/x-b))>10e-5&n
<=500)
x=a/(b+x);
n=n+1;
end
x
r1=(-b+sqrt(bA2+4*a))/2
r2=(-b-sqrt(bA2+4*a))/2