大连东软信息学院Matlab实验报告-第一次作业
MATLAB第一次实验报告
MATLAB第一次实验报告1. 编程实现以下数列的图像，用户能输入不同的初始值以及系数。
并以x，y为坐标显示图像x(n+1) = a*x(n)-b*(y(n)-x(n)^2);y(n+1) = b*x(n)+a*(y(n)-x(n)^2)程序：function [y1,y2]=draw(a,b,c,x1,y1)x=[];y=[];x(1)=x1;y(1)=y1;n=c;for i=1:nx(i+1)=a*x(i)-b*(y(i)-x(i)^2);y(i+1)=b*x(i)+a*(y(i)-x(i)^2);endy1=x;y2=y;plot(y1,y2)输入：draw（5,1,1,1,2）2. 编程实现奥运5环图，允许用户输入环的直径。
程序：clcclearr=input('input a mnmber:')color1=['b','k','r'];color2=['y','g'];t=linspace(0,2*pi);a1=r+0.2*r;a2=3*r+0.4*r;a3=5*r+0.6*r;a=[a1,a2,a3];b1=2*r+0.3*r;b2=4*r+0.5*r;b=[b1,b2];for i=1:3x=r*cos(t)+a(i);y=r*sin(t);plot(x,y,color1(i),'linewidth',5) axis([0,7*r,-2*r,1.5*r]);hold onaxis equalendfor n=1:2x=r*cos(t)+b(n);y=r*sin(t)-r;plot(x,y,color2(n),'linewidth',5) hold onendaxis off输入：6（圆的半径）3. 实现对输入任意长度向量元素的冒泡排序的升序排列。
不允许使用sort函数。
matlab大作业实验报告,《Matlab程序设计》期末实验报告-大作业2015.doc
matlab⼤作业实验报告,《Matlab程序设计》期末实验报告-⼤作业2015.doc《MATLAB程序设计》实验报告学院： 学号： 姓名：⼀、题⽬：1、(10分)已知矩阵，⽤Matlab代码实现以下要求：(1)将矩阵赋给变量A，并在屏幕上显⽰A；(2)将A按列进列逆序重排，重排后的矩阵赋给变量B，并在屏幕上显⽰B；(3)⽤reshape命令将A重排为⼀个2⾏6列矩阵并赋给变量C；(4)将A重排为⼀个列向量，将其赋给变量D，求D的平均值，在屏幕上显⽰D和它的平均值；(5)⽤命令查看变量A的维数，并显⽰运⾏结果。
2、(10分)写代码实现以下要求：构造菜单项‘Plot’，菜单项Plot有两个⼦菜单项Plot sin(选择此项后执⾏画出曲线，线型为虚线，线条颜⾊为红⾊)和Plot cos(选择此项后执⾏画出曲线 ，线型为实线，线条宽度为2)。
3、(20分)已知，实现下列操作：(1)在同⼀个图形窗⼝，同⼀坐标系下⽤不同的颜⾊和线型绘制三条曲线，并添加图例来区分三条曲线(5分)。
(2)⽤subplot命令，以⼦图的⽅式绘制三条曲线，图形排列⽅式为三⾏⼀列(5分)。
(3) 分别⽤直⽅图(bar)、棒状图(stem)和填充图(fill)绘制三条曲线，以⼦图⽅式绘制，排列⽅式为3⾏3列，共9幅⼦图(10分)。
4、(10分)⽤surf命令绘制曲⾯图形，⽤shading interp命令进⾏插值着⾊处理并添加垂直颜⾊棒。
5、(15分)⾃2011年9⽉1⽇起，我国实⾏新的个⼈所得税征收办法，起征点为3500元，请⽤If-else if-else-end结构实现⼈⼯输⼊⽉收⼊后能计算出个⼈所得税的缴纳额并显⽰⽉收⼊10000元时应缴纳的税款。
级数应纳税所得额x(元)税率备注1x<=15003%x指⽉收⼊扣除起征点3500元之后的余额；215008000045%同上6. (10分)⽤while-end循环结构计算级数和的值，输⼊n值，能计算出f的值，并显⽰结果。
第一次MatlAB实验报告
实验一 数学模型的MatlAB 描述一、MatlAB 中数学模型的表示MatlAB 中数学模型的表示主要有三种基本形式：传递函数分子/分母多项式模型、传递函数零极点增益模型和状态空间模型。
它们各有特点，有时需在各种模型之间进行转换。
1. 传递函数分子/分母多项式模型当传递函数为nn n nm m m m a s a sa s ab s b sb sb s G ++++++++=----11101110)(时，在MatlAB 中，直接用分子、分母的系数表示，即],,,[10m b b b num = ],,,[10n a a a d e n = ),()(d e n num tf s G =例1：已知一系统的传递函数2432328()3842ss G s s s s s ++=++++，提取其分子和分母多项式，并绘制零极点图*。
>> num=[3 2 8]; den=[1 3 8 4 2]; G=tf(num,den) [tt,ff]=tfdata(G ,'v') pzmap(G)grid on 结果如下：2. 传递函数的零极点增益模型当传递函数为)())(()())(()(1010n m p s p s p s z s z s z s Ks G ------=时，在MatlAB 中，用[z,p,k]矢量组表示，即[][][]),,()(,,,,,1010k p z zpk s G K k p p p p z z z z n m ====例2:求一传递函数222328()(61)(24)s s G s s s s s ++=++++的零极点及其增益，并绘制零极点图。
>> num=[3 2 8];den=conv([1 6 12],[1 2 4]); G=tf(num,den) GG=zpk(G)[z,p,k]=zpkdata(G ,'v') pzmap(G) grid结果如下：3.状态空间模型当系统的数学模型为状态空间表达式{BuAXX+='+=DuCXY时，在MatlAB中，用[A、B、C、D]矩阵组表示，即系统表示为()DCBAss、、、实验3:已知两系统串联，其中2132265452s sGs s s++=+++，222712432s sGs s++=++，求此串联系统的总的传递函数。
Matlab实验报告
实验结果及分析实验1：程序如下x=1:10y=2*x;plot(x,y)仿真结果：实验结果分析：仿真结果是条很规则的直线，X轴和Y轴一一对应，清楚明了，而序又特别简单。
所以用Maltab 软件很方便地画出规则的直线，方便研究。
实验结果及分析1、A=2、A=1A=实验结果及分析实验三 Matlab在信号与系统中的应用实验名称实验1、掌握信号与系统课程中基本知识的Matlab编程、仿真方法目的实验原理实验1程序：b=[1];a=[1 1];p=;t=0:p:5;x=exp(-3*t);subplot(1,2,1);impulse(b,a,0:p:5);title('冲激响应');subplot(1,2,2);step(b,a,0:p:5);title('阶跃响应');实验内容<设计性实验>1、用MATLAB在时域中编程求解y′(t)+y(t)=f(t), f(t)= exp(-3t)ε(t)的冲激响应、阶跃响应。
在simulink仿真环境下，设计系统框图，分析系统的冲激响应、阶跃响应。
<设计性实验>（选做）2、用MATLAB在时域中编程求解y′(t)+y(t)=f(t), f(t)=(1+exp(-3t))ε(t)的冲激响应、阶跃响应，要求用conv编程实现系统响应。
在simulink仿真环境下，设计系统框图，分析系统的冲激响应、阶跃响应。
实验结果及分析实验1仿真结果：simulink仿真环境下冲激响应阶跃响应实验名称实验四 Matlab在数字信号处理中的应用实验结果及分析实验1仿真结果:6khz12kHZ。
MATLAB实验报告
MATLAB实验报告MATLAB实验报告姓名：专业：学号：实验⼀MATLAB环境的熟悉与基本运算⼀、实验⽬的：1．熟悉MATLAB开发环境2．掌握矩阵、变量、表达式的各种基本运算⼆、实验基本知识：1.熟悉MATLAB环境:MATLAB桌⾯和命令窗⼝、命令历史窗⼝、帮助信息浏览器、⼯作空间浏览器⽂件和搜索路径浏览器。
2.掌握MATLAB常⽤命令3.MATLAB变量与运算符变量命名规则如下：（1）变量名可以由英语字母、数字和下划线组成（2）变量名应以英⽂字母开头（3）长度不⼤于31个（4）区分⼤⼩写MATLAB中设置了⼀些特殊的变量与常量，列于下表。
MATLAB运算符，通过下⾯⼏个表来说明MATLAB的各种常⽤运算符表2MATLAB算术运算符表3MATLAB关系运算符表4MATLAB逻辑运算符表5MATLAB特殊运算4.MATLAB的⼀维、⼆维数组的寻访表6⼦数组访问与赋值常⽤的相关指令格式5.MATLAB的基本运算表7两种运算指令形式和实质内涵的异同表6.MATLAB的常⽤函数表8标准数组⽣成函数表9数组操作函数三、实验内容1、学习安装MATLAB软件。
2、学习使⽤help命令，例如在命令窗⼝输⼊helpeye，然后根据帮助说明，学习使⽤指令eye（其它不会⽤的指令，依照此⽅法类推）3、学习使⽤clc、clear，观察commandwindow、commandhistory和workspace等窗⼝的变化结果。
4、初步程序的编写练习，新建M-file，保存（⾃⼰设定⽂件名，例如exerc1、exerc2、exerc3……），学习使⽤MATLAB的基本运算符、数组寻访指令、标准数组⽣成函数和数组操作函数。
注意：每⼀次M-file的修改后，都要存盘。
四、实验结果练习A：（1）helprand，然后随机⽣成⼀个2×6的数组，观察commandwindow、commandhistory和workspace等窗⼝的变化结果。
matlab实验报告(实验1)
开课学院及实验室：机电学院 2012年11 月16日
（2）在E盘新建一个文件夹并把其加入搜索路径。
截屏相关过程结果并保存。
（3）逐步按照以下步骤进行并记录各个步骤的结果：
1）在指令窗中输入clear指令；
5）再次在指令窗中输入clear指令后，输入
8）把cow.m暂时移动到另一个无关的文件夹，如“我的文档”。
（2）保存全部变量为数据文件Mydata.mat；在对应文件夹中找到该文件，截屏后保存结果；
（4）删除全部内存变量；
Editor/Debugger和脚本编写初步
把以下指令编写为M脚本文件，并运行：
t=0:pi/50:4*pi;
y0=exp(-t/3); y=exp(-t/3).*sin(3*t); plot(t,y,'-r',t,y0,':b',t,-y0,':b') grid
什么？此方法也是查询函数注释的有效方法，请牢记。
实验一MATLAB操作基础实验报告
闽 江 学 院 电 子 系 实 验 报 告学生姓名：班级：学 号：课程：MATLAB 程序设计一、 实验题目：MATLAB 操作基础 二、 实验地点：实验楼A210 实验目的：1、熟悉MATLAB 的操作环境及基本操作方法。
2、掌握MATLAB 的搜索路径及其设置方法。
3、熟悉MATLAB 帮助信息的查阅方法。
三、 实验内容：1、先建立自己的工作目录，再将自己的工作目录设置到MATLAB 搜索路径下，再试验用help 命令能否查询到自己的工作目录。
2、在MATLAB 环境下验证例1-1至1-4，并完成以下题目： （1）绘制右图所示图形 （2）求383、利用MATLAB 的帮助功能分别查询inv ，plot 、max 、round 等函数的功能及用法。
4、在工作空间建立一个变量a ,同时在当前目录下建立一个M 文件：a.m ，试在命令窗口输入a ，观察结果，并解释原因。
四、 实验环境（使用的软硬件）： MATLAB 7.0 五、 实验结果：实验程序： （1）求38。
程序： x=8;plot(x, x^(1/3));0.10.20.30.40.50.60.70.80.91-1-0.8-0.6-0.4-0.200.20.40.60.81得到的结果： ans= 2（2）绘制图像 程序：x=[0:0.001:1]；plot(x, sin(2*pi*x),x, 2*x-1,x,0); 得到结果如下图：（3）查询plot 、inv 、round 、max 等函数的功能及用法 分别输入：输入：help plot 得到：PLOT Linear plot.PLOT(X,Y) plots vector Y versus vector X. If X or Y is a matrix, then the vector is plotted versus the rows or columns of the matrix, whichever line up. If X is a scalar and Y is a vector, length(Y) disconnected points are plotted.PLOT(Y) plots the columns of Y versus their index.If Y is complex, PLOT(Y) is equivalent to PLOT(real(Y),imag(Y)). In all other uses of PLOT, the imaginary part is ignored.Various line types, plot symbols and colors may be obtained with0.10.20.30.40.50.60.70.80.91-1-0.8-0.6-0.4-0.200.20.40.60.81PLOT(X,Y,S) where S is a character string made from one elementfrom any or all the following 3 columns:b blue . point - solidg green o circle : dottedr red x x-mark -. dashdot c cyan + plus -- dashed m magenta * stary yellow s squarek black d diamondv triangle (down)^ triangle (up)< triangle (left)> triangle (right)p pentagramh hexagramFor example, PLOT(X,Y,'c+:') plots a cyan dotted line with a plus at each data point; PLOT(X,Y,'bd') plots blue diamond at each data point but does not draw any line.PLOT(X1,Y1,S1,X2,Y2,S2,X3,Y3,S3,...) combines the plots defined by the (X,Y,S) triples, where the X's and Y's are vectors or matrices and the S's are strings.For example, PLOT(X,Y,'y-',X,Y,'go') plots the data twice, with a solid yellow line interpolating green circles at the data points.The PLOT command, if no color is specified, makes automatic use of the colors specified by the axes ColorOrder property. The default ColorOrder is listed in the table above for color systems where the default is blue for one line, and for multiple lines, to cyclethrough the first six colors in the table. For monochrome systems, PLOT cycles over the axes LineStyleOrder property.PLOT returns a column vector of handles to LINE objects, onehandle per line.The X,Y pairs, or X,Y,S triples, can be followed byparameteralue pairs to specify additional propertiesof the lines.输入：help inv得到：INV Matrix inverse.INV(X) is the inverse of the square matrix X.A warning message is printed if X is badly scaled ornearly singular.输入：help round得到：ROUND Round towards nearest integer.ROUND(X) rounds the elements of X to the nearest integers.输入：help max得到：MAX Largest component.For vectors, MAX(X) is the largest element in X. For matrices,MAX(X) is a row vector containing the maximum element from eachcolumn. For N-D arrays, MAX(X) operates along the firstnon-singleton dimension.[Y,I] = MAX(X) returns the indices of the maximum values in vector I. If the values along the first non-singleton dimension contain morethan one maximal element, the index of the first one is returned.MAX(X,Y) returns an array the same size as X and Y with thelargest elements taken from X or Y. Either one can be a scalar.[Y,I] = MAX(X,[],DIM) operates along the dimension DIM.When complex, the magnitude MAX(ABS(X)) is used, and the angleANGLE(X) is ignored. NaN's are ignored when computing the maximum.Example: If X = [2 8 4 then max(X,[],1) is [7 8 9],7 3 9]max(X,[],2) is [8 and max(X,5) is [5 8 59], 7 5 9].inv 求矩阵的逆plot 画图max 求最大值round 向0取整六、思考练习：1、help命令和look for命令有什么区别？lookfor 关键词在所有M文件中找“关键词”，比如：lookfor inv（即寻找关键词“inv”）其实就和我们平时用CTRL+F来查找“关键词”是一样的而help是显示matlab内置的帮助信息用法：help 命令，比如 help inv ，作用就是调用inv这个命令的帮助2、什么是工作空间？假定有变量A与B存在于工作空间中，如何用命令保存这两个变量？下次重新进入MATLAB后，又如何装载这两个变量？①工作空间：工作空间是由系统所提供的特殊变量和用户自己使用过程生成的所有变量组成的一个概念上的空间；②保存变量：保存当前工作空间，比如要保存的路径为f:\matlab\work\matlab.mat,输入命令：save f:\matlab\work\matlab.mat ；③载入变量：输入命令：load f:\matlab\work\matlab.mat 。
matlab第一次实验报告
一、心得体会在此之前只接触过Matlab一段时间，通过本次上机实验初次体验，体会到Mathematica语法规则简单，操作使用方便，和其他语言交互，编程灵活方便，功能十分强大。
Mathematica和Matlab都可以实现某个方案，不同之处在于对不同的问题，处理的方法繁琐程度不一样。
软件掌握多少都不算多，类似的软件还有MathCAD、Maple等，关键在处理具体问题的时候那种软件更适合，或者那种软件处理问题普遍性强，而Mathematica和Matlab分别代表了符号计算和数据处理（或者数值计算）的两个方面普遍适用的软件。
个人感觉Mathematica符号计算强，便于公式推导和函数式制图。
Mathematica界面的窗口和菜单栏是分离式的，便于边计算边对照查看文档，其公式可以复制为公式格式，可直接粘贴于Word等文本中，也便于以公式的形式输入，但该功能在大量符号计算中可能作用不大，编写代码的执行效率肯定比手动输入公式要高的多。
输入和输出都是在同一个界面，便于查看每条代码的结果。
函数使用方式丰富，调用函数可以编写专门的文件，也可直接在输入界面中调用。
可把想到的代码都写成函数的形式，包括各种形式的代码，便于调用。
将所有的算法最终转化为表达式，以函数的形式体现。
例如绘图将表达式的变量范围指定即可，便于函数制图。
Matlab很强大，处理什么都比较得心应手，是款很好的数学式编程工具与计算工具，修改编辑很方便，但体积比较大。
Mathematica 个人感觉侧重于图像处理，尤其那些多边形，多维空间制图，曲面等很实用，而且软件体积比较少，运行速度快，是款和MATLAB 相辅相成的好工具。
通过以上，我自认为较好的解决方案是数据处理用Matlab ，然后将数据导入Mathematica 可视化处理。
符号计算就是Mathematica 了。
以后还需多进行实验，体会Mathematica 的功能。
二、指令运行关于Mathematica 绘图的基本功能以不同输出格式画曲面2||,2||,)(22≤≤=+-y x e z y x Plot3D[Exp[-(x^2+y^2)],{x,-2,2},{y,-2,2}]螺旋线：x=sint,y=cost,z=t/3的图形ParametricPlot3D[{Sin[t],Cos[t],t/8},{t,0,15},Axes →False]圆环面x=cost(3+cosu),y=sint(3+cosu),z=sinu,u≤t≤2π,0≤u≤2π的图形ParametricPlot3D[{Cos[t](3+Cos[u]),Sin[t](3+Cos[u]),Sin[u]},{t,0,2Pi},{u,0, 2Pi},Axes False]三、问题求解（Help）。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
MATLAB教程练习题
1．利用Matlab中的函数生成3×3单位阵，2×3全0阵和4×3全1阵；
2．利用Matlab中的函数生成10个在区间[14,99]上均匀分布的随机整数。
3. 计算向量（2，4，6，8）的最大值，最小值，平均值，中值【查一下这个函数】；
4. 在行向量（2，4，6，8）和（1，3，5，7）之间实施加减乘除及幂运算；
5. 在向量（1，3，5）与标量2之间实施加减乘除及幂运算。
6. 求 121211130⎛⎫ ⎪ ⎪ ⎪⎝⎭－的转置、行列式【查一下这个函数】、逆矩阵【查一下这个函数】。
7.要求利用diag 函数和其他函数生成数组A =。
8.将上面生成的数组A = ，转换成行向量B。
9. 产生一均匀分布在（-5，5）之间的5*3阶随机矩阵C，并要求所有的数值为整数。
10. 利用函数产生一个（0,1）间均匀分布的6*6随机矩阵D，然后统计D中大于等于0.6的元素个数。
11.产生一个4*5的随机矩阵E，要求为整数，且矩阵的元素范围为[0,100]，计算出矩阵的最大值所在的位置。
12. 产生100到1的等差数列，数据点为20个，赋值给F
（1）取F的第10个单元的值，赋值给F1；
（2）取F的前10个单元的值，赋值给F2；
（3）取F的后5个单元的值，赋值给F3；
（4）取F的偶数单元的值，赋值给F4；
（5）将F的数据倒排，赋值给F5。
13. 定义矩阵G=[1 2 3;4 5 6;7 8 9]；
（1）抽取A中的第1、2行，第2、3列元素，赋值给G1；
（2）抽取A中的第1、4个元素作为新矩阵G2的第1行，抽取第2、3个元素作为新矩阵G2的第2行；
resleep(1 4 2 3 , 2,2)
（3）对矩阵G1和G2进行求和、求差、求乘和点乘。