实验四 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的绘图

西安邮电大学《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基本操作一、实验目的2．学习使用图形函数计算器命令funtool及其环境。

3. 学习使用help命令进行帮助4. 掌握向量与矩阵的创建以及矩阵的基本操作5. 掌握数组与矩阵的概念二、实验内容熟悉Matlab操作环境，认识命令窗口、内存工作区窗口、历史命令窗口；学会使用format命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；掌握数组与矩阵的概念；学会使用help命令进行帮助；学会使用who和whos命令查看内存变量信息；学会使用图形函数计算器funtool；1.命令窗口的简单使用（1）简单矩阵的输入（自由创建）x=[1 3 5;2 4 6]x =1 3 52 4 6（2）求[12+2×(7-4)]÷32的算术运算结果，总结算术运算符先级[12+2*(7-4)]/3^2ans =22.有关向量、矩阵或数组的一些运算（1）设A=15；B=20；求C=A+B与c=a+b的区别A=15;B=20;>> C=A+BC =35>> c=a+bUndefined function or variable 'a'.（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是对应元素相乘)A=[1 2 3;4 5 6;7 8 9];B=[9 8 7;6 5 4;3 2 1];>> A*Bans =30 24 1884 69 54138 114 90>> A.*Bans =9 16 2124 25 2421 16 9（3）设a=10，b=20；求i=a/b与j=a\ba=10;>> b=20;>> i=a/bi =0.5000>> j=a\bj =2（4）设a=[1 -2 3;4 5 -4;5 -6 7]请设计出程序，分别找出小于0的矩阵元素的线性索引以及行列索引（sub2ind/ind2sub）。

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实验四

MATLAB实验四MATLAB程序设计试验报告说明：1 做试验前请先预习，并独立完成试验和试验报告。

2 报告解答方式：将MATLAB执行命令和最后运行结果从命令窗口拷贝到每题的题目下面，请将报告解答部分的底纹设置为灰色，以便于批阅。

3 在页眉上写清报告名称，学生姓名，学号，专业以及班级。

3 报告以Word文档书写。

文档命名方式: 学号+姓名+_（下划线）+试验几.doc 如：110400220张三_试验1.doc4 试验报告doc文档以附件形式发送到******************。

凡文档命名不符合规范，或者发送方式不正确，不予登记。

一目的和要求（1）熟练掌握MATLAB的程序设计流程控制结构（2）熟练掌握M文件的结构（3）熟练掌握函数调用和参数传递（4）熟练掌握内联函数的使用（5）了解程序性能剖析窗口二内容和步骤1 （必做）编写M脚本文件，实现：1）从键盘逐个输入学生的某科成绩，计算学生的平均分并输出显示。

2）统计学生成绩分布，其中90分以上为A类，80~89分为B类，70~79分为C类，60~69分为D类，60分以下为E类，统计各类的学生人数，并打印显示。

3）将学生成绩保存到结构数组student. score域和统计成绩分布信息保存数组level中并存作为student. level域，最后把结构数组student保存到score. mat文件中。

（提示：1）成绩输入采用input函数，可采用while循环，提示用户连续输入学生成绩，直到用户输入负数成绩，终止成绩录入。

2）成绩保存可采用save函数）score=0;n=0;sum=0;a=0；b=0;c=0;d=0;e=0;while score>=0if score>=90a=a+1;elseif score>=80elseif score>=70c=c+1;elseif score>=60d=d+1;elsee=e+1;ends um=sum+score;n=n+1;score=input('score:');endaverage=sum/(n-1)fprintf(' A类人数：%d',a)fprintf(' B类人数:%d',c)fprintf(' C类人数:%d',c)fprintf(' D类人数:%d',d)fprintf(' E类人数:%d',e-1)2（必做）有一分数序列：2/1，3/2，5/3，8/5，13/8，21/13... 采用循环迭代求出这个数列的前n项之和。

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 输出形参表=函数名（输入形参表）注释说明部分函数体语句注意：函数名的命名规则与变量名相同。当输出形参多于一个时，应用方括号括起来。

实验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实验报告(实验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）把图中的α改为σ。

实验四用Matlab求解方程

均可以得到实根 2.0946.
温州大学城市学院 1. fzero 的基本用法 [x, fv, ef, out, jac]=fsolve(@F, x0, opt, P1, P2,…)
2 2 x1 + x2 = 4 例 11 求方程组 2 的解. 2 x1 − x2 = 1
Matlab 编程代码: 先建立文件 fangchz.m function y=fangchz(x,a,b,c,d) y(1)=x(1)^2+a*x(2)^2-b; %当a=1,b=4时为第1个方程 y(2)= x(1)^2+c*x(2)^2-d; %当c=-1,d=1时为第2个方程输入程序: x0=[2,2];
3 2 p 1 = x + 21 x + 20 x 例 12 计算多项式
的根以及由多项式的根得出系数. >> p1=[1 21 20 0]; >> roots(p1) >> poly([0;-20;-1]) % 计算多项式的根 % 计算多项式的系数
温州大学城市学院例 13 求下列方程的根 f = function y=fun(x)
例 4 求函数 f=3x2+4x+5 对变量 x 的一阶导数.
温州大学城市学院 (3) 积分格式: int (f, t, a, b) 功能: 求函数 f 对变量 t 从 a 到 b 的定积分. 例 4 求函数 f=ax2+bx+c 对变量 x 的不定积分. >> syms a b c x >> f=a*x^2+b*x+c >> int(f) 求函数 f 对变量 x 从 1 到 5 的定积分, 命令为 >> int(f,1,5)

实验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 实验四

实验四姓名：XX 学号：XXXXXXXXXX7.1创建一个正弦信号的仿真模型。

创建模型如下：运行结果：7.2在工作空间计算变量t 和y，将其运算的结果作为系统的输入。

在Command Window中输入：>> t=0:0.1:10;y=sin(t);t=t';y=y';并创建如下模型：运行结果：7.3建立二阶系统的仿真模型。

输入信号源使用阶跃信号，系统使用开环传递函数，接受模块使用示波器来构成模型。

方法一：构建模型如下所示：运行结果：方法二：构建模型如下所示：仿真结果如下：7.4控制部分为离散环节，被控对象为两个连续环节，其中一个有反馈环，反馈环引入了零阶保持器，输入为阶跃信号。

构建模型如下：调整T、Tk d(t)、y(t)的波形变化：T=Tk=0.1 (a) d(k)示波器显示y(t)示波器显示T=0.6 Tk=0.1 (a) d(k)示波器显示y(t)示波器显示T=0.6 Tk=0.6 (a) d(k)示波器显示y(t)示波器显示7.5将例7.4 中的模块结构参数用变量表示，结构图如图7.35所示。

构建模型如下：在Command Window中输入如下指令：T=0.1;Tk=0.6;k=0.03;zt1=1.44;zt2=-1.26;zt3=1;zt4=-1;tf11=6.7;tf12=0.1;tf13=1;tf21=1;tf22=3;tf23=1;运行结果：7.6打开【例7.4】建立的模型，将控制对象中的第一个连续环节中的反馈环建立为一个子系统。

创建子函数后：打开子函数如下：7.7 在【例7.6】的基础上建立新子系统，将【例7.6】模型的控制对象中的第一个对象环节整个作为一个子系统。

构建系统如下：打开子系如下：7.8建立一个用使能子系统控制正弦信号为半波整流信号的模型。

构建模型如下：运行结果如下：7.9建立一个用触发子系统控制正弦信号输出阶梯波形的模型。

构建如下模型：运行结果如下：创建模型如下：创建子系统后如下：封装之后系统为：打开二阶系统其参数输入对话框如下：。

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]的曲线。