实验四 MATLAB程序设计
河北工业大学MATLAB实验四
2014秋2012级《MATLAB 程序设计》实验报告一、实验目的：1、掌握MATLAB 二维图形绘制命令及其图形控制；2、熟悉MATLAB 三维图形绘制命令及其图形控制；3、熟悉特殊二维图形、三维图形的绘制方法。
二、实验内容：1、在同一坐标系下绘制下面三个函数在[0,40]t ∈的图象，要求三种曲线采用不同颜色以及不同的线型，同时对每条曲线利用函数进行标注，并添加图例。
120.134sin()t y t y y e t π-===2、编写程序，选择合适的步距，绘制下面函数在区间[-6,6]中的图象，并对分段的曲线进行标注，同时添加x 轴和y 轴的说明。
sin ,0(),036,3x x y x x x x x ≤⎧⎪=<≤⎨⎪-+>⎩3、利用subplot 函数在同一绘图窗口中用不同颜色和线型绘制以下两个函数在t ∈[-2π,2π]范围内的图象。
0.50.21222t t y y e -==4、某学校有30位教师,其中教授5人,副教授8人,讲师12人,助教5人.试画出职称比例结构饼图，并强调图形的教授部分同时添加标注。
5、生成1×10维的随机数向量a ，分别用红、黄、蓝、绿色绘出其连线图、脉冲图、阶梯图和条形图，并分别标出标题“连线图”、“脉冲图”、“阶梯图”、“条形图”。
6、使用subplot 函数，把图形窗口分成两个部分，分别绘制sin2t 和3cos3t 曲线,t 范围:0-5.要求前者用红色实线，数据点形状为五角星,有网格线，x 轴加标注“x ”，y 轴加标注“y=sin2t ”,加题目“y=sin2t 的曲线”；后者用蓝色虚线，数据点形状为圆形，无网格线，x 轴加标注“x ”，y 轴加标注“y=3cos3t ”,加题目“y=3cos3t 的曲线”；7、绘制下列三维曲线:(1)/20/20cos sin ,02t t x e ty e t t z t π--⎧=⎪=≤≤⎨⎪=⎩ (2) 23,01x t y t t z t =⎧⎪=≤≤⎨⎪=⎩8、绘制下列曲面图，并调整三维图的视角、背景色、着色以及透视效果。
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.熟悉各种绘图函数的使用；2.掌握图形的修饰方法和标注方法；3.了解MATLAB 中图形窗口的操作。
二、实验内容：1．用图形表示连续调制波形Y=sin(t)sin(9t)及其包络线。
程序代码如下：包络线：2．x=［－2π，2π］，y1=sinx、y2=cosx、y3=sin2x、y4=cos 2x①用MATLAB语言分四个区域分别绘制的曲线，并且对图形标题及横纵坐标轴进行标注。
程序：结果：②另建一个窗口，不分区，用不同颜色、线型绘出四条曲线，并标注图例注解。
程序：结果：③绘制三维曲线：⎪⎩⎪⎨⎧=≤≤==)cos()sin()200()cos()sin(t t t z t t y t x π程序：结果：3．绘制极坐标曲线ρ=asin(b+nθ),并分析参数a、b、n对曲线形状的影响。
（1）a=1;b=1;n=1（2）a=10;b=1;n=1（3）a=10;b=10;n=1 （4）a=10;b=10;n=10参数a、b、n对曲线形状的影响：由上面绘制的图形可知：a决定图形的大小，当a为整数时，图形半径大小就是a；b决定图形的旋转角度，图形的形状及大小不变；n决定图形的扇叶数，当n 为奇数时，扇叶数为n，当n为偶数时，扇叶数为2n。
三、结论本次实验用到了曲线绘图、三位曲线绘图的知识，与老师上课的内容一致，让我学的matlab绘图的知识得到了巩固，我还学会了如何使用title、subplot、plot、axis等函数。
在做实验的过程复习了hold on指令是覆盖函数继续绘图的意思。
MATLAB程序设计与应用实验报告4
《MATLAB程序设计与应用》实验报告专业：机械设计制造及其自动化班级：机械三班学号：17428320姓名：刘羽新2019年上期实验四：MATLAB 绘图基础一、实验目的通过本实验掌握MATLAB 二维、三维作图的基本操作，了解绘制图形的辅助操作。
二、实验要求1. 根据实验内容，编写相应的MATLAB 程序，并将程序及结果放置于相应位置。
2. 实验内容采用小四号宋体、1.5倍行距进行排版。
3. 实验报告以word 文档的.doc 格式保存；文件文件名为：学号_姓名，如“15312345_张三”。
三、实验内容1. 绘制下列曲线 (1) 21001y x =+, 55x -≤≤ （用绿色的虚线绘制，添加图形说明“二维曲线绘图”） x=-5:0.001:5;y=100./(1+x.*x);plot(x,y,'b:');title('二维曲线绘图');(2)235x ty t⎧=⎪⎨=⎪⎩, 1010x-≤≤（用红色的双划线绘制，添加图例“参数曲线绘图”）t=-10:0.001:10;x=t.*t;y=5*t.*t.*t;plot(x,y,'r--');legend('参数曲线绘图');(3)cossinx ty tz t=⎧⎪=⎨⎪=⎩, 02xπ≤≤（用黑色实线绘制，添加坐标轴名称）t = 0:0.001:2*pi; x=cos(t);y=sin(t);plot3(x,y,z,'k-'); xlabel('X');ylabel('Y');zlabel('Z');2. 绘制下列分段函数的曲线，并在各段曲线附近分别添加文本说明“21y x x =++”与“21sin y x π=+”：221,013sin ,12x x x y x x π⎧++≤<⎪=⎨+≤≤⎪⎩ t_1=0:0.001:1;t_2=1:0.001:2;y_1=x_1.*x_1+x_1+1;y_2=3+(sin(pi*x_2)).^2;plot(x_1,y_1,x_2,y_2);gtext('y=x^2+x+1');gtext('y=1+(sin(pi*x))^2');3. 已知241(11)1xy xx+=-≤≤+，sin3cos(0)sin3sinx t tty t tπ=⎧≤≤⎨=⎩，以子图形式绘制两条曲线。
实验四 MATLAB的绘图
西安邮电大学《Matlab程序设计基础》课内实验报告题目：实验四MATLAB的绘图院系名称：计算机学院专业名称：计算机科学与技术班级：计科1502班内序号：05141107(29)学生姓名：赵阳指导教师：张老师实验四MATLAB的绘图一、实验目的：掌握基本的绘图函数plot。
二、实验内容：1、在同一坐标系下绘制下面三个函数在[0，4pi]的图象。
代码如下：t=0:pi/100:4*pi;y1=t;y2=sqrt(t);y3=4*pi*exp(-0.1*t).*sin(t);plot(t,y1,t,y2,t,y3)2、编写程序，选择合适的步距，绘制下面函数在区间[-6，6]中的图象。
代码如下：x=linspace(-6,6,100);y=[];for x0=xif x0<=0y=[y,sin(x0)];elseif x0<=3y=[y,x0];elsey=[y,-x0+6];endendplot(x,y)3、用compass函数画下面相量图ua=1;ub=cos(-2*pi/3)+sin(-2*pi/3)*i;uc=cos(2*pi/3)+sin(2*pi/3)*i;compass([ua,ub,uc,ua-ub,ub-uc,uc-ua])代码如下：ua=1;ub=cos(-2*pi/3)+sin(-2*pi/3)*i;uc=cos(2*pi/3)+sin(2*pi/3)*i;compass([ua,ub,uc,ua-ub,ub-uc,uc-ua])4、三维空间曲线绘制z=0:0.1:4*pi;x=cos(z);y=sin(z);plot3(x,y,z)代码如下：z=0:0.1:4*pi;x=cos(z);y=sin(z);plot3(x,y,z)5、用mesh或surf函数，绘制下面方程所表示的三维空间曲面，x和y的取值范围设为[-3，3]。
代码如下：[x,y]=meshgrid(-3:0.5:3);z=-x^2/10+y^2/10;mesh(x,y,z)三、思考题在同一坐标系下，用不同颜色和线型绘制以下两个函数在tÎ[-2p，2p]范围内的图象。
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实现拉格朗日插值、分段线性插值
实验四用MATLAB实现拉格朗日插值、分段线性插值一、实验目的：1）学会使用MATLAB软件；2）会使用MATLAB软件进行拉格朗日插值算法和分段线性差值算法；二、实验内容：1用MATLAB实现y = 1./(x.^2+1);（-1<=x<=1）的拉格朗日插值、分段线性2.选择以下函数，在n个节点上分别用分段线性和三次样条插值的方法，计算m个插值点的函数值，通过数值和图形的输出，将插值结果与精确值进行比较，适当增加n，再作比较，由此作初步分析：(1).y=sinx;( 0≤x≤2π)(2).y=(1-x^2)（-1≤x≤1）三、实验方法与步骤：问题一用拉格朗日插值法1）定义函数：y = 1./(x.^2+1);将其保存在f.m 文件中，程序如下：function y = f1(x)y = 1./(x.^2+1);2）定义拉格朗日插值函数：将其保存在lagrange.m 文件中，具体实现程序编程如下：function y = lagrange(x0,y0,x)m = length(x); /区间长度/n = length(x0);for i = 1:nl(i) = 1;endfor i = 1:mfor j = 1:nfor k = 1:nif j == kcontinue;endl(j) = ( x(i) -x0(k))/( x0(j) - x0(k) )*l(j); endendendy = 0;for i = 1:ny = y0(i) * l(i) + y;end3）建立测试程序，保存在text.m文件中，实现画图：x=-1:0.001:1;y = 1./(x.^2+1);p=polyfit(x,y,n);py=vpa(poly2sym(p),10)plot_x=-5:0.001:5;f1=polyval(p,plot_x);figureplot(x,y,‘r',plot_x,f1)二分段线性插值：建立div_linear.m文件。
实验四、Matlab 程序设计
实验四 MATLAB 程序设计实验目的：1. 掌握建立和执行M 文件的方法2. 掌握利用if 语句实现选择结构的方法3. 掌握利用switch 语句实现多分支选择结构的方法4. 掌握利用for 语句实现循环结构的方法5. 掌握利用while 语句实现循环结构的方法6. 熟悉利用向量运算来代替循环操作的方法实验内容：1、输入一个x ，求分段函数的值，用if 语句实现。
32,1x 50,6530,6222≠⎪⎩⎪⎨⎧--≠<≤+--≠<-+=x x x x x x x x x x y 及其它且且2、输入一个百分制成绩，要求输出成绩等级A 、B 、C 、D 、E. 其中90分~100分为A ，80分~89分为B ，70分~79分为C ，60分~69分为D ，60分以下为E. 要求：分别用if 语句和switch 语句实现。
3、某公司员工的工资计算方法如下：a) 工作时数超过120小时者，超过部分加发15%；b) 工作时数低于60小时者，扣发700元；c) 其余按每小时84元计发。
试编程按输入的员工的工时数，计算应发工资。
4、根据2222213121116n ++++= π，求π的近似值。
当n 分别取100，1000，10000时，结果是多少？要求：分别用循环结构和向量运算(使用sum 函数)来实现。
5、根据12151311-++++=n y ，求： a) y<3时的最大n 值。
b) 与(1)的n 值对应的y 值。
6、考虑以下迭代公式：nn x b a x +=+1 其中a 、b 为正的常数。
(1) 编写程序求迭代的结果，迭代的终止条件为5110-+≤-n n x x ，迭代初值x 0=1.0，迭代次数不超过500次。
(2) 如果迭代过程收敛于r ，那么r 的准确值是242a b b +±-，当(a, b)的值取(1,1)、 (8,3)、(10,0.1)时，分别对迭代结果和准确值进行比较。
matlab实验四
实验四曲线绘图【实验目的】1．了解曲线的几种表示方法。
2．学习掌握MATLAB软件有关的命令。
【实验内容】绘制下列三种曲线：1．以直角坐标方程x,=表示的正、余弦线。
sin=y cosyxt=0:0.1:10;y1=sin(t);y2=cos(t);plot(t,y1,'r',t,y2,'b--');x=[1.7*pi;1.6*pi];y=[-0.3;0.8];s=['sin(t)';'cos(t)'];text(x,y,s);title('正弦和余弦曲线');legend('正弦','余弦')xlabel('时间t'),ylabel('正弦、余弦')gridaxis square2．以参数方程]y=tx表示的平面曲线（单位圆）。
∈=t t2,0[,sin,cosπ>> t=0:0.1:2*pi;>> x=cos(t);>> y=sin(t);>> plot(x,y)>> title('x=cos(t),y=sin(t)');>>axis square3． 以参数方程]20,0[,,sin 2,2cos2.02.0∈===--t t z t e y t e x t t ππ表示的空间曲线。
>> t=0:0.1:20; >> x=exp(-0.2*t).*cos(pi/2*t);>>y=pi/2*exp(-0.2*t).*sin(t);>>z=t;>>plot3(x,y,z);>> xlabel('x 轴')>> ylabel('y 轴')>> zlabel('z 轴')>> title(' x=exp(-0.2*t)*cos(pi/2*t),y=pi/2*exp(-0.2*t)*sin(t),z=t, t ∈[0：20]')4． 以极坐标方程]2,0[,1),cos 1(πϕϕ∈=+=a a r 表示的心脏线。
《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 求解状态空间模型1、实验设备MATLAB 软件2、实验目的① 学习线性定常连续系统的状态空间模型求解、掌握MATLAB 中关于求解该模型的主要函数；② 通过编程、上机调试，进行求解。
3、实验原理说明Matlab 提供了非常丰富的线性定常连续系统的状态空间模型求解(即系统运动轨迹的计算)的功能，主要的函数有：初始状态响应函数initial()、阶跃响应函数step()以及可计算任意输入的系统响应数值计算函数lsim()和符号计算函数sym_lsim()。
数值计算问题可由基本的Matlab 函数完成，符号计算问题则需要用到Matlab 的符号工具箱。
4、实验步骤① 根据所给状态空间模型，依据线性定常连续系统状态方程的解理论，采用MATLAB 编程。
② 在MATLAB 界面下调试程序，并检查是否运行正确。
习题1：试在Matlab 中计算如下系统在[0,5s]的初始状态响应，并求解初始状态响应表达式。
Matlab 程序如下：A=[0 1; -2 -3];B=[]; C=[]; D=[];x0=[1; 2];sys=ss(A,B,C,D);[y,t,x]=initial(sys,x0,0:5);plot(t,x)0011232⎡⎤⎡⎤==⎢⎥⎢⎥--⎣⎦⎣⎦x x x习题2：试在Matlab 中计算如下系统在[0,10s]内周期为3s 的单位方波输入下的状态响应。
并计算该系统的单位阶跃状态响应表达式。
Matlab 程序如下：A=[0 1; -2 -3];B=[0; 1]; C=[]; D=[];x0=[1; 2];sys=ss(A,B,C,D);[u t]=gensig('square',3,10,0.1)0011232⎡⎤⎡⎤==⎢⎥⎢⎥--⎣⎦⎣⎦x x x[y,t,x] = lsim(sys,u,t,x0)plot(t,u,t,x);（注：文档可能无法思考全面，请浏览后下载，供参考。
实验4Matlab程序设计1
实验4 Matlab程序设计1实验目的：1、掌握建立和执行M文件的方法；2、掌握实现选择结构的方法；3、掌握实现循环结构的方法。
实验内容：1. 从键盘输入一个4位整数，按如下规则加密后输出。
加密规则：每位数字都加上7,然后用和除以10的余数取代该数字；再把第一位与第三位交换，第二位与第四位交换。
2. 求分段函数的值。
2x +x-6, x <0且x式-3y = *x2—5x + 6 0Wxc 5 且x 式2及x 式3X2 _ X -1 其他用if语句实现，分别输出x=-5,-3,0,1,2,2.5,3,5时的y值。
3. 输入一个百分制成绩，要求输出成绩等级A、B、C、D、E，其中90~100分为A，80~89分为B，70~79分为C，60~69分为D，60分以下为E。
要求：（1）分别用if语句和swich语句实现。
（2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。
4. 硅谷公司员工的工资计算方法如下：（1）工作时数超过120小时者，超过部分加发15% ；（2）工作时数低于60小时者，扣发700元；（3）其余按每小时84元计发。
试编程按输入的工号和该号员工的工时数，计算应发工资。
5. 设计程序，完成两位数的加、减、乘、除四则运算。
即：输入两个两位随机整数，再输入一个运算符号，做相应的运算，并显示相应的结果。
6. 建立5X 6矩阵，要求输出矩阵的第n行元素。
当n值超过矩阵的行数时，自动转为输出矩阵的最后一行元素，并给出出错信息。
7. 产生20 个两位随机整数，输出其中小于平均数的偶数。
实验代码及实验结果1.>> a=input(' 请输入一个四位整数：');请输入一个四位整数：9988>> a1=fix(a/1000);>> a2=rem(fix(a/100),10);>> a3=rem(fix(a/10),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);>> b=a1;a1=a3;a3=b;>> b=a2;a2=a4;a4=b;>> c=a1*1000+a2*100+a3*10+a4;>> disp(c)3.>> a=input C请输入一个四位整数：一); 请输人一个四f立整數：3988 »(a/1000);>> a2=ren(fix(a/100), 10);>> a3=ren(fiK(a/10)?10);>> 血(a, ID);» al=re&(al+^ 10)；>> a^rsB(a2+7f 10)；>> a3=rem(a3+7? 10):>> 丑"“血(M+C 10);>> b=al:a1=a3:a3=b；» b=a2:a2=a4:a4=b,» c-al *100(Hai2*l 0(H-a3*l 0+ttl;>> disp(c)55662.x=input('请输入x的值：');if x<0 & x~=-3y=x92+x-6;elseif x>=0 & x<5 & x~=2 & x~=3y=x.A2+5.*x+6elsey=x.A2-x-1;end请输入孟的値：[-瓦-亠0,1, 2t2. E f 3, 5]7 =29. 0000 IL ODDO -I. 0000 -1.0000 1. OODO 2.7500 5. QOOD(1) if语句:a=i nput（'请输入成绩：'）;ig.ooooif a>=90 & a<=100disp('A');elseif a>=80 & a<90disp('B');elseif a>=70 & a<80disp('C');elseif a>=60 & a<70disp('D');elseif a<60 &a>=0disp('E');elsedisp('输入有误！'); end( 2) switch 语句：a=input(' 请输入成绩：'); switch fix(a/10)case{9}disp('A');case{8}disp('B');case{7}disp('C');case{6}disp('D');case nu m2cell(2:5)disp('E')otherwisedisp('输入有误！');end请输入咸绩：80E请输入咸绩；5E请输入成绩：1212输入有误I»4.n=i nput('请输入工号：’)；a=i nput('请输入工作小时数：'); if a>=120 y=a*84+a*84*0.15;elseif a<120 & a>=60y=a*84;elsey=a*84-700;disp(y);end请输入工号：30请输入工作小时数：231232。
matlab数学实验报告
MATLAB数学实验报告指导老师：班级：小组成员：时间：201_/_/_Matlab第二次实验报告小组成员：1题目：实验四;MATLAB选择结构与应用实验目的：掌握if选择结构与程序流程控制;重点掌握break;return;pause语句的应用..问题：问题1：验证“哥德巴赫猜想”;即：任何一个正偶数n>=6均可表示为两个质数的和..要求编制一个函数程序;输入一个正偶数;返回两个质数的和..问题分析：由用户输入一个大于6的偶数;由input语句实现..由if判断语句判断是否输入的数据符合条件..再引用质数判断函数来找出两个质数;再向屏幕输出两个质数即可..编程：function z1;z2=geden;n=input'please input n'if n<6disp'data error';returnendif modn;2==0for i=2:n/2k=0;for j=2:sqrtiif modi;j==0k=k+1;endendfor j=2:sqrtn-iif modn-i;j==0k=k+1;endendif k==0fprintf'two numbers are'fprintf'%.0f;%.0f';i;n-ibreakendendend结果分析：如上图;用户输入了大于6的偶数返回两个质数5和31;通过不断试验;即可验证哥德巴赫猜想..纪录：if判断语句与for循环语句联合嵌套使用可使程序结构更加明晰;更快的解决问题..2题目：实验四;MATLAB选择结构与应用实验目的：用matlab联系生活实际;解决一些生活中常见的实际问题..问题：问题四：在一边长为1的四个顶点上各站有一个人;他们同时开始以等速顺时针沿跑道追逐下一人;在追击过程中;每个人时刻对准目标;试模拟追击路线;并讨论.. （1）四个人能否追到一起（2）若能追到一起;每个人跑过多少路程（3）追到一起所需要的时间设速率为1问题分析：由正方形的几何对称性和四个人运动的对称性可知;只需研究2个人的运动即可解决此问题..编程：hold onaxis0 1 0 1;a=0;0;b=0;1;k=0;dt=0.001;v=1;while k<10000d=norma-b;k=k+1;plota1;a2;'r.';'markersize';15;plotb1;b2;'b.';'markersize';15;fprintf'k=%.0f b%.3f;%.3f a%.3f;%.3f d=%.3f\n';k;b1;b2;a1;a2;da=a+b1-a1/d*dt;b2-a2/d*dt;b=b+b2-a2/d*dt;-b1-a1/d*dt;if d<=0.001breakendendfprintf'每个人所走的路程为：%.3f';k*v*dtfprintf'追到一起所需要的时间为%.3f';k*dt结果分析：上图为2人的模拟运动路线;有对称性可解决所提问题..-上图为运算过程和运算结果..四个人可以追到一起;走过的路程为1.003;时间也为1.003.纪录：此题利用正方形和运动的对称性可以简便运算..3题目：实验八;河流流量估计与数据插值目的：由一些测量数据经过计算处理;解决一些生活实际问题..问题：实验八上机练习题第三题：瑞士地图如图所示;为了算出他的国土面积;做以下测量;由西向东为x轴;由南向北为y轴;从西边界点到东边界点划分为若干区域;测出每个分点的南北边界点y1和y2;得到以下数据mm..已知比例尺1:2222;计算瑞士国土面积;精确值为41288平方公里..测量数据如下：x=7.0 10.5 13.0 17.5 34 40.5 44.5 48 56 61 68.5 76.5 80.5 91 96 101 104 106 111.5 118 123.5 136.5 142 146 150 157 158 ;y1=44 45 47 50 50 38 30 30 34 36 34 41 45 46 43 37 33 28 32 65 55 54 52 50 66 66 68;y2=44 59 70 72 93 100 110 110 110 117 118 116 118 118 121 124 121 121 121 122 116 83 81 82 86 85 68;问题分析：先由题目给定的数据作出瑞士地图的草图;再根据梯形法;使用trapz语句;来估算瑞士国土的面积..编程：x=7.0 10.5 13.0 17.5 34 40.5 44.5 48 56 61 68.5 76.5 80.5 91 96 101 104 106 111.5 118 123.5 136.5 142 146 150 157 158;y1=44 45 47 50 50 38 30 30 34 36 34 41 45 46 43 37 33 28 32 65 55 54 52 50 66 66 68;y2=44 59 70 72 93 100 110 110 110 117 118 116 118 118 121 124 121 121 121 122 116 83 81 82 86 85 68;plotx;y1;'r.';'markersize';15;plotx;y2;'r.';'markersize';15;axis0 160 0 135grid;hold ont=7:158;u1=splinex;y1;t;u2=splinex;y2;t;plott;u1plott;u2s1=trapzt;u1;s2=trapzt;u2;s=s2-s1*2222*22222/10000000;fprintf'S=%.0f';s结果分析：上图为由所给数据绘制出的瑞士地图..上图为运算结果;计算出瑞士的国土面积为42472平方公里;与准确值41288较为接近..纪录：使用梯形分割的方法;trapz语句可以方便计算不规则图形面积;但存在一定误差..4题目：实验七：圆周率的计算与数值积分目的：将数值积分最基本的原理应用于matlab之中;解决一些与积分有关的问题..问题：实验七上机练习题第一题：排洪量某河床的横断面如图7.3所示;为了计算最大排洪量;需要计算其断面积;试根据所给数据m用梯形法计算其断面积..问题分析：河床断面可近似分割成若干曲边梯形;近似处理把它们当做梯形来计算面积可使问题得到简化..编程：clc;clear;x=0 4 10 12 15 22 28 34 40;y=0 1 3 6 8 9 5 3 0;y1=10-y;plotx;y1;'k.';'markersize';15;axis0 40 0 10;grid;hold ont=0:40;u=splinex;y1;t;plott;u;s=40*10-trapzt;u;fprintf's=%.2f\n';s结果分析：上图为河床的断面图..上图为计算结果面积约为180.70平方米..纪录：使用梯形法计算不规则图形面积十分简便易行..5题目：实验七：圆周率的计算与数值积分目的：使用matlab计算解决一些有关积分的问题..问题：实验七上机练习题第三题：从地面发射一枚火箭;在最初100秒内记录其加速度如下;试求火箭在100秒时的速度..Ts=0 10 20 30 40 50 60 70 80 90 100;Am/s*s=30.00 31.63 33.44 35.47 37.75 40.33 43.29 46.69 50.67 54.01 57.23;问题分析：加速度为速度的微分;已知微分求积分;类似于面积问题;可使用梯形法来计算..编程：clc;clear;x=0 10 20 30 40 50 60 70 80 90 100;y=30.00 31.63 33.44 35.47 37.75 40.33 43.29 46.69 50.6754.01 57.23;plotx;y;'k.';'markersize';15;axis0 100 20 60;grid;hold ons=0:10:100;z=splinex;y;s;plots;y;v=trapzx;y;fprintf'v=%.2f\n';v结果分析：上图为加速度变化图..上图为计算结果;求得火箭在100秒时速度约为4168.95m/s..纪录：梯形法可以推广解决许多已知微分求积分的其他问题..6题目：实验七：圆周率的计算与数值积分目的:计算曲线弧长闭曲线周长可使用微元法;ds=sqrtdx^2+dy^2;在转化微积分问题;累加即可得到结果..问题：实验七上机练习题第三题：计算椭圆想x^2/4+y^2=1的周长;使结果具有五位有效数字..问题分析：编程：s=0;dx=0.001;for x=0:0.001:1.999dy=1.-x+0.001.^2/4-1.-x.^2/4;ds=sqrtdx.^2+dy.^2;s=s+ds;ends=4*s;fprintf'the length is'fprintf'%.4f';s结果分析：上图为计算结果;给定椭圆的周长约为9.1823五位有效数字纪录：计算不规则曲线弧长;可使用微元法;划分为若干小的看做直角三角形;利用勾股定理解决..7题目：实验九人口预测与数据拟合目的：掌握一些曲线拟合的方法;了解曲线拟合常用函数..问题：用电压U=10v的电池给电容器充电;t时刻的电压Vt=U-U-V0exp-t/τ;其中V0是电容器的初始电压;τ是充电常数;由所给数据确定V0和τ..t=0.5 1 2 3 4 5 7 9;V=3.64 3.52 2.74 1.78 1.34 1.01 0.57 0.37;问题分析：题中已给出函数关系式;为指数函数曲线拟合;将所给函数式整理可得标准的exp形函数曲线;从而便于解决..编程：t=0.5 1 2 3 4 5 7 9;V=3.64 3.52 2.74 1.78 1.34 1.01 0.57 0.37;plott;V;'k.';'markersize';20;axis0 10 0 4;grid;hold onpause0.5n=8;a=sumt1:n;b=sumt1:n.*t1:n;c=sumlogV1:n;d=sumt1:n.*logV1:n;A=n a;a b;B=c;d;p=invA*Bx=0:10;y=expp1+p2*x;plotx;y;'r-';'linewidth';2结果分析：上图为电压与时间关系图..上图为计算结果;即U-V0=1.4766;所以V0=8.5234;-1/τ=-0.2835;所以τ=3.5273纪录：曲线拟合的一个重难点是选择合适的曲线函数;才能提高拟合度..8题目：实验七圆周率的计算与数值积分目的：拓展圆周率的各种计算方法;掌握其他数值的近似计算方法..问题：实验七练习2：计算ln2的近似值精确到10的-5次方（1）利用级数展开的方法来计算（2）利用梯形法计算（3）利用抛物线法问题分析：级数展开;梯形法;抛物线法是常见的近似运算方法..编程：1级数展开的方法clc;clear;n=0;r=1;p=0;k=-1;while r>=0.1e-5n=n+1;k=k*-1;p1=p+k/n;r=absp1-p;fprintf'n=%.0f;p=%.10f\n';n;p1;p=p1;end2梯形法clc;clear;f=inline'1./x';x=1:0.1:2;y=fx;p=trapzx;y;fprintf'p=%.6f\n';p3抛物线法clc;clear;f=inline'1./x';a=1;b=2;n=1;z=quadf;a;b;fprintf'z=%.10f\n';z结果分析：（1）级数展开的方法（2）梯形法3抛物线法纪录：级数展开法;梯形法;抛物线法;计算近似值时应合理利用..梯形法和抛物线法不易提高精确度;级数展开法可以提高精确度..9题目：实验八河流流量估计与数据插值目的：掌握求插值多项式的方法;并利用此计算近似值..问题：已知y=fx的函数表如下x=0.40 0.55 0.65 0.80 0.90 1.05;y=0.41075 0.57815 0.69675 0.88811 1.02652 1.25382;求四次拉格朗日插值多项式;并由此求f0.596问题分析：利用所给函数表可计算拉格朗日插值多项式..编程：function p=lagrangex;yL=lengthx;a=onesL;for j=2:La:;j=a:;j-1.*x';endx=inva*y';for i=1:Lpi=xL-i+1;endx=0.40 0.55 0.65 0.80 0.90 1.05;y=0.41075 0.57815 0.69675 0.88811 1.02652 1.25382; plotx;y;'k.';'markersize';15axis0 2 0 2grid;hold on;p=lagrangex;y;t=0:0.1:1.5;u=polyvalp;t;plott;u;'r-'a=polyvalp;0.596结果分析：上图为所求结果;估算值和插值多项式..纪录：插值多项式是一项十分实用的方法..10题目：求正整数n的阶乘：p=1*2*3*…*n=n;并求出n=20时的结果目的：练习使用循环变量解决数学问题问题：对程序：Clear;clc;n=20;p=1;for i=1:np=p*i;fprintf’i=%.0f;p=%.0f\n’;i;pend进行修改使它：利用input命令对n惊醒赋值问题分析：题中给出程序中“n=20”修改;使用input命令；讲题中的输出命令放出循环之外..编程：clear;clc;n=input'n=';p=1;for i=1:np=p*i;endfprintf'i=%.0f;p=%.0f\n';i;p结果：n=20i=20;p=2432902008176640000>>结果分析：使用input命令可以实现人机对话;使用户自由赋值；输出语句在程序中的位置对输出的结果有很大的影响；在循环内部可以在计算过城中不断输出结果；在循环之外则可以控制只输出最后结果..11题目：对于数列｛√2｝;n=1;2;…;求当其前n项和不超过1000时的n的值及合的大小..目的：运用条件循环解决文帝个项数的循环程序求解;问题：对程序：clear;clc;n=0;s=0;while s<=1000n=n+1;s=s+sqrtn;fprintf’n=%.0f;s=%.4f\n’;n;send问题分析：题中所给程序中的限制变量为上次循环之后的s;导致s超过上限后仍有一次的循环;若把循环变量改为这次的s;则可以避免这种情况的发生..编程：clear;clc;n=0;s=0;while s+sqrtn<=1000n=n+1;s=s+sqrtnfprintf'n=%.0f;s=%.4f\n';n;send结果：……s =970.8891n=128;s=970.8891s =982.2469n=129;s=982.2469s =993.6487n=130;s=993.6487>>结果分析：从结果中可以看出;最后一步为我们需要的答案;从这道题我们可以得出循环变量对一道编程的重要性..。
matlab实验报告(实验4)
学生实验报告开课学院及实验室: 机电学院2012年12月21日学院机电学院年级、专业、班姓名学号实验课程名称MATLAB程序设计成绩实验项目名称实验4: 数据和函数的可视化指导老师一、实验目的1、掌握MATLAB绘图的基本步骤和相关指令调用的先后顺序。
2、掌握MATLAB绘图指令的调用方法。
二、实验内容数学函数从形式上可以分为离散函数和连续函数。
MATLAB对这两种函数数据的可视化都提供了相应的指令。
仔细阅读教材【例5.1-1】的实现代码, 运行并保存结果；并改用stem函数, 画出【例5.1-1】的序列图。
仔细阅读教材【例5.1-2】的实现代码, 运行并保存结果；并分别使用描点和连折线方式, 画出连续函数y=xcosx的近似图形（采样点数自定, 要求画出的图尽量接近原连续函数的图）。
仔细阅读【例5.2-2】的实现代码, 理解plot指令画多条曲线的运用方法, 运行并保存结果；并使用plot函数和legend函数, 在同一个图形窗口上画出y=sint和y=sin(2t)在[0,2pi]区间上的图形, 并标出图例。
仔细阅读【例5.2-4】的实现代码, 理解图形标识选项的运用方法, 运行并保存结果；并修改代码, 把“sin(t)”字体改为正体, 大小改为20, “极大值”改为宋体。
阅读【例5.2-6】, 理解使用hold on指令画多幅图的方法, 运行并保存结果。
阅读【例5.2-8】, 理解使用subplot函数画多个子图的方法, 运行并保存结果。
（1）综合实验: 阅读以下关于通过绘制二阶系统阶跃响应综合演示图形标识的示例, 理解示例中所有图形标识指令的作用, 掌握各个图形标识指令的运用方法, 并在原指令上改动以实现以下功能:（2）把横坐标范围改为0至5pi, 纵坐标范围改为0至2；（3）把图中的横轴的刻度改为从0开始到4pi, 中间各点间隔为pi/2；纵轴刻度改为从0开始到1.5, 中间各点间隔为0.3；（4）把图中的α改为σ。
实验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）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。
MATLAB实验四：MATLAB程序设计 参考答案
function y=f(x) y=(x^3-2*x^2+x-6.3)/(x^2+0.05*x-3.14)
河南财经政法大学数学与信息科学学院
6
实验报告
f(1)*f(2)+f(3)
function f=f(a,b,c) f=g(a)*g(b)+g(c)^2; function g=g(x) g=(x^3-2*x^2+x-6.3)/(x^2+0.05*x-3.14);
a=rand(1,100); [b,c]=sort(a); index=c(1); i=1; while i<=7 if i==8 disp('You Lost!'); i=i+1; else d=input('please input the number you guess:'); if d>index disp('High'); elseif d<index disp('Low'); elseif d==index disp('You won!'); i=9;
syms i j k l n i=1;k=0; n=input('Input n:'); for j=1:n i=i*j; k=k+j; end fprintf('%d！=%d \n Sum(1,...,%d)=%d\n',n,i,n,k);
4、用 while-end 循环语句求不超过 1000 的偶数之和，并求显示出最大值。 s=0,n=0; while s<=1000 n=n+2; s=s+n; if s>1000 break end
matlab实验四函数编写与程序设计
实验四：函数编写与程序设计一、实验目的1 . 掌握M文件的创建。
2．掌握函数的编写规则。
3．掌握函数的调用。
4 . 掌握基本的输入输出函数以及显示函数的用法。
5．会用Matlab程序设计实现一些工程算法问题。
二、实验内容1 . 设计程序，完成两位数的加、减、乘、除四则运算，即产生两个两位随机整数，再输入一个运算符号，做相应的运算，显示相应的结果，并要求结果显示类似于“a=x+y=34”。
a=input('请输入一个数：')b=input('请输入一个数：');fuhao=input('请输入一个运算符号（+-*/）：','s');switch fuhaocase {'+'}he=a+b;disp(['和=',num2str(a),'+',num2str(b),'=',num2str(he)]);case {'-'}he=a-b;disp(['减=',num2str(a),'-',num2str(b),'=',num2str(he)]) case {'*'}he=a*b;disp(['乘=',num2str(a),'*',num2str(b),'=',num2str(he)]) case {'/'}he=a/b;disp(['除=',num2str(a),'/',num2str(b),'=',num2str(he)]) otherwise disp('请输入正确的符号');end2 . 求下列分段函数的值2 2 26,0 56,1,yχχχχχχχχχχχ⎧+-<≠⎪=-+≤≠≠⎨⎪--⎩且-30<10,2且3其他要求：用if语句实现，分别输出x = -5.0，-3.0，1.0，2.0，2.5，3.0，5.0时的y值。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验四 MATLAB 程序设计
1.实验目的
（1）熟练掌握MATLAB 的程序流程控制结构。
（2）掌握M 文件的结构，M 函数文件编写、使用。
（3）熟练掌握函数调用和参数传递。
2.实验仪器
（1）Matlab6.5应用软件安装版 一套
（3）PC 机 一台
3. 实验原理
依据MA TLAB 的编程的原理，编写M 函数文件，调用M 函数文件，完成曲线绘制。
4. 实验步骤
（1）利用for 、while 控制语句和sum 指令求和。
（2）使用MA TLAB 函数文件绘图。
（3）利用置换指令绘制脉冲响应曲线.
（4）利用feval 指令实现指定的函数。
5. 实验报告内容
（1）分别用for 和while 循环语句计算∑==630i i 2
K 的程序，再写出一种避免循环的计算程序。
（提示：可考虑利用MA TLAB 的sum （X,n ）函数，实现沿数组X 的第n 维求和。
）保存为M 文件.
（2）将课本例4-8子函数编程及调用演示，三个子图上的圆和多边形绘制在同一坐标系中。
并保存为M 文件。
保存输出结果图。
（3）利用置换指令subs(X,new),例4-18中的脉冲响应在t=[0,18]的曲线。
并保存为M 文件。
保存输出结果图。
（4）试利用feval （）指令计算F(x)+F 2(x)，其中F 可取‘sin ’、‘cos ’。
（提示：先编写以个M 函数function y=trif(F,x)实现F(x)+F 2(x)的计算，在编写调用函数完成F 为‘sin ’、‘cos ’的计算），并保存为M 文件。