2015MATLAB原理及应用实验报告第五章
MATLAB实验报告
实验一MATLAB操作基础实验目的和要求：1、熟悉MATLAB的操作环境及基本操作方法。
2、掌握MATLAB的搜索路径及设置方法。
3、熟悉MATLAB帮助信息的查阅方法实验内容：1、建立自己的工作目录，再设置自己的工作目录设置到MA TLAB搜索路径下，再试验用help命令能否查询到自己的工作目录。
2、在MA TLAB的操作环境下验证课本；例1-1至例1-4，总结MATLAB的特点。
例1-1例1-2例1-3例1-43、利用帮助功能查询inv、plot、max、round等函数的功能。
4、完成下列操作：（1）在matlab命令窗口输入以下命令：x=0:pi/10:2*pi;y=sin(x);（2）在工作空间窗口选择变量y，再在工作空间窗口选择回绘图菜单命令或在工具栏中单击绘图命令按钮，绘制变量y的图形，并分析图形的含义。
5、访问mathworks公司的主页，查询有关MATLAB的产品信息。
主要教学环节的组织：教师讲授实验目的、开发环境界面、演示实验过程，然后同学上机练习。
思考题：1、如何启动与退出MA TLAB集成环境？启动：（1）在windows桌面，单击任务栏上的开始按钮，选择‘所有程序’菜单项，然后选择MA TLAB程序组中的MA TLABR2008b程序选项，即可启动MATLAB系统。
（2）在MA TLAB的安装路径中找到MA TLAB系统启动程序matlab.exe，然后运行它。
（3）在桌面上建立快捷方式后。
双击快捷方式图标，启动MA TLAB。
退出：（1）在MA TLAB主窗口file菜单中选择exitMATLAB命令。
（2）在MA TLAB命令窗口中输入exit或quit命令。
（3）单击MATLAB主窗口的关闭按钮。
2、简述MATLAB的主要功能。
MATLAB是一种应用于科学计算领域的数学软件，它主要包括数值计算和符号计算功能、绘图功能、编程语言功能以及应用工具箱的扩展功能。
3、如果一个MATLAB命令包含的字符很多，需要分成多行输入，该如何处理？使用‘；’隔开4、Help命令和lookfor命令有何区别？help是显示matlab内置的帮助信息一般是help 命令；而lookfor 关键词是通过关键词查找，“关键词”就是你要找的词语5、在MA TLAB环境下，建立了一个变量fac，同时又在当前目录下建立了一个m文件fac.m，如果需要运行fac.m文件，该如何处理？实验二MATLAB矩阵及运算实验目的和要求：1、掌握MATLAB数据对象的特点及运算规则2、掌握MATLAB建立矩阵的方法及矩阵处理的方法3、掌握MATLAB分析的方法实验内容：1.求下列表达式的值（1）（2）（3）（4）2.已知A=[-1,5,-4;0,7,8;3,61,7],B=[8,3,-1;2,5,3;-3,2,0] 求下列表达式的值：（1）A+6B和A^2-B+I(其中I为单位矩阵)。
MATLAB原理及应用实验报告第五章答案
《MATLAB 原理及应用》实验报告三．练习题1、编写程序绘制如下图形（时间间隔自己设定）图中曲线标注:sin(t)标注在(1.7*pi,-0.3)，cos(t)标注在(1.6*pi,0.8) 2、使用subplot指令在不同区域画出下列曲线y=sinAx;y=cosAx;y=Ae x ,y=x^2。
x 的取值范围是]15,0[，A 为自己的学号。
3、画出衰减振荡曲线t e y t 3sin 3-=及其它的包络线30t e y -=。
t 的取值范围是]4,0[π，图的标题字体大小(fontsize)为10号，字体(fontname)宋体，字体风格（fontangle ）为斜体字体，并选择合适的位置在图上标注自己的学号，姓名，以及作图日期，字体大小、字体风格自己指定。
四．实验程序及结程序为1：t=0:0.1:4*pi;y=sin(t);y1=cos(t);plot(t,y,'r',t,y1,'--b') text(1.7*pi,-0.3,'sint') text(1.6*pi,0.8,'cost')2468101214-1-0.8-0.6-0.4-0.200.20.40.60.812：x=0:0.1:15; y1=sin(1*x); subplot(2,2,1); plot(x,y1); y2=cos(1*x); subplot(2,2,2); plot(x,y2);y3=(1*exp(x)); subplot(2,2,3); plot(x,y3); y4=x.^2; subplot(2,2,4); plot(x,y4);051015-1-0.500.51051015-1-0.500.510510151234x 1060510151002003003 程序为：t=0:0.01:4*pi;y0=exp(-t/3);y1=-exp(-t/3);y=exp(-t/3).*sin(3*t);plot(t,y0,'--b',t,y,'r',t,y1,'--b');title('y=e^(-t/3)sin(3t)及它的包络线','fontsize',10','fontweigh','bold','fontangle','italic');text(9,-0.4,['学号:'],'fontsize',10','fontweigh','bold','fontangle','italic');text(9.5,-0.6,['姓名:某某某'],'fontsize',10','fontweigh','bold','fontangle','italic');text(10,-0.8,['2013.12.01']);。
matlab实验报告总结精选
matlab实验报告总结电气工程学院自动化102班 2012年12月21日实验一 MATLAB环境的熟悉与基本运算一、实验目的1．熟悉MATLAB开发环境2．掌握矩阵、变量、表达式的各种基本运算二、实验基本知识1.熟悉MATLAB环境MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器、文件和搜索路径浏览器。
2.掌握MATLAB常用命令变量与运算符变量命名规则如下：变量名可以由英语字母、数字和下划线组成变量名应以英文字母开头长度不大于31个区分大小写MATLAB中设置了一些特殊的变量与常量，列于下表。
MATLAB运算符，通过下面几个表来说明MATLAB的各种常用运算符表2 MATLAB算术运算符表3 MATLAB关系运算符表4 MATLAB逻辑运算符表5 MATLAB特殊运算的一维、二维数组的寻访表6 子数组访问与赋值常用的相关指令格式的基本运算表7 两种运算指令形式和实质内涵的异同表的常用函数表8 标准数组生成函数表9 数组操作函数三、实验内容1、新建一个文件夹2、启动，将该文件夹添加到MATLAB路径管理器中。
3、保存，关闭对话框4、学习使用help命令，例如在命令窗口输入help eye，然后根据帮助说明，学习使用指令eye5、学习使用clc、clear，观察command window、command history和workspace等窗口的变化结果。
6、初步程序的编写练习，新建M-file，保存，学习使用MATLAB的基本运算符、数组寻访指令、标准数组生成函数和数组操作函数。
注意：每一次M-file的修改后，都要存盘。
练习A：help rand，然后随机生成一个2×6的数组，观察command window、command history和workspace等窗口的变化结果。
学习使用clc、clear，了解其功能和作用。
答：clc是清除命令窗体内容 clear是清除工作区间输入C=1:2:20，则C表示什么？其中i=1,2,3,?,10。
MATLAB实验报告5
实验报告（5）：插值实验（Interpolation）1．(1)Write a MATLAB function to perform the Lagrangian form Interpolation.function yy=Lagrange(x,y,xi)m=length(x);n=length(y);if m ~= n , error('向量x与y长度须一致'); end;s=0;for i=1 : nz=ones(1,length(xi));for j=1 : nif j ~= iz=z.*(xi-x(j))./(x(i)-x(j));endends=s+z*y(i);endyy=s;保存为Lagrange2(2)Write a call back m-file to check the function you write.x=[0.5610,0.56280,0.56401,0.56521];y=[0.82741,0.82659,0.82577,0.82495];xi=[0.5625,0.5635,0.5645];yi=Lagrange2(x,y,xi)yi=2．In the following table(saved in the data file named Interpolation_Problem1_data.xls), X is the Fahrenheit（华氏的）temperature, and Y is the number of times a cricket（蟋蟀）chirps（鸣叫）in one minute.(1)Make a scatterplot(散点图) of the data.x=[46 49 51 52 54 56 57 58 59 60 61 62 63 64 66 67 68 69 70 71];y=[40 50 55 63 72 70 77 73 90 93 96 88 99 110 113 120 127 137 132 137];plot(x,y,'ro');xlabel('Fahrenheit temperature');ylabel('number of times a cricket chirps in one minute ');grid on实验报告（5）：插值实验（Interpolation ）Fahrenheit temperaturen u m b e r o f t i m e s a c r i c k e t c h i r p s i n o n e m i n u t e(2) Interpolate the data with some low order piecewise interpolating polynomials. x=[46 49 51 52 54 56 57 58 59 60 61 62 63 64 66 67 68 69 70 71];y=[40 50 55 63 72 70 77 73 90 93 96 88 99 110 113 120 127 137 132 137]; xi=45:0.05:71; yi=interp1(x,y,xi); plot(x,y,'r*',xi,yi)xlabel('Fahrenheit temperature');ylabel('number of times a cricket chirps in one minute '); grid onFahrenheit temperaturen u m b e r o f t i m e s a c r i c k e t c h i r p s i n o n e m i n u t e实验报告（5）：插值实验（Interpolation ）splineFahrenheit temperaturen u m b e r o f t i m e s a c r i c k e t c h i r p s i n o n e m i n u t ecubicFahrenheit temperaturen u m b e r o f t i m e s a c r i c k e t c h i r p s i n o n e m i n u t e(3) Fit the data by using an 18th -degree polynomial that passing through （通过） the data points as an empirical model. Analyze the appropriateness of doing so.x=[46 49 51 52 54 56 57 58 59 60 61 62 63 64 66 67 68 69 70 71];y=[40 50 55 63 72 70 77 73 90 93 96 88 99 110 113 120 127 137 132 137]; xi=46:0.05:71; p=polyfit(x,y,18); yi=polyval(p,xi); plot(x,y,'o',xi,yi); title('polyinter')xlabel('Fahrenheit temperature');ylabel('number of times a cricket chirps in one minute ');grid on实验报告（5）：插值实验（Interpolation ）Fahrenheit temperaturen u m b e r o f t i m e s a c r i c k e t c h i r p s i n o n e m i n u t e后期图像放大：Fahrenheit temperaturen u m b e r o f t i m e s a c r i c k e t c h i r p s i n o n e m i n u t e运用高次多项式是插值，虽然得出的结果光滑有较好的分析性质，但图形的震荡的频率极大，不具收敛性。
matlab实验报告
目录一、基础题 (2)二、绘图题 (3)三、句柄图形和图形用户界面制作 (9)四、Simulink仿真题 (13)五、专题实验总结 (17)六、参考文献 (17)一．基础题实验目的：1、掌握数组的创建和寻访；2、掌握MA TLAB 数组的运算；3、熟悉MA TLAB 关系操作和逻辑操作；4、掌握函数的调用。
实验要求：1、熟练掌握数组运算；2、熟悉MA TLAB 的各种数据类型及函数调用；3、初步了解M 文件及其编写、调试与运行。
实验内容：某公司投资2000万元建成一条生产线。
投产后，在时刻t 的追加成本和追加收益分别为G(t)= (百万元/年), H(t)= (百万元/年)。
试确定该生产线在合适何时停产可获最大利润？最大利润是多少？解：构造函数f(t)=H(t)-G(t)=13-t-3t 2/3=0 ;令t 1/3=x,则f(t)=-t 3-3t 2+13 可得矩阵P=[-1,-3,0,13]求最佳生产时间的源程序如下：p=[-1,-3,0,13]; x=roots(p); t=x.^3运行结果如下： t =3.6768 +21.4316i 3.6768 -21.4316i4.6465考虑到实际情况，显然两个虚数根应该舍掉。
所以将t=4.6425带入，求积分。
代码：t=4.6465; x=0:0.01:t;y=13-x-3*x.^(2/3); trapz(x,y)运行结果： ans =26.32083/225tt ++3/218t-结论：比较以上三组数据，可知最佳生产时间t=4.6465年，可获得的最大收益为 26.3208（百万元/年）。
减去20（百万元）投资，可得最终利润为6.3208（百万元）。
二．作图题实验目的：1、进一步熟悉M 文件调试过程；2、熟练掌握MA TLAB 二维、三维图形的绘制；3、掌握图形的修饰； 实验要求：1、进一步熟悉和掌握MA TLAB 的编程及调试；2、掌握二维、三维图形的绘制；3、掌握图形交互指令的使用；实验内容：1、二维绘图： （1）函数)sin(2x ey x-=,求：a. 绘制[0 ,8]区间内的图形，加注x,y 轴及图形名称；b. 最大值和最小值；c. 零值；求解：a、命令：ezplot('2*exp(-x)*sin(x)',[0,8]) %绘图title(f) , xlabel('x') ,ylabel('y') %加轴名及图名图像:b、先求最小值：从图像上可以看出，在[3,4]之间有最小值，所以命令：[xmin,fmin]=fminbnd(f,3,4)Hold onPlot(xmin,fmin)运行结果：xmin =3.9270 fmin =-0.0279增加一句标注： text(xmin,fmin,'(3.9270,-0.0279)')求最小值求最大值分析：要求最大值，即为求y=-f(x)的最小值，从图像上可以看出，在[0,2]之间有最大值命令：f=inline('-2*exp(-x)*sin(x)') ezplot(f,[0,8][xmin,fmin]=fminbnd(f,0,2) hold onplot(xmin,fmin,'*') Hold onPlot(xmin,fmin)运行结果： xmin = 0.7854 fmin = -0.6448增加标注：text(xmin,fmin,'(0.7854,-0.6448)')（如上图）综上可知：fmax= 0.6448 fmin= -0.0279c 、为了便于观察，首先增加一条直线y=0 命令：hold onezplot('0',[0,8])图像：0123456780.10.20.30.40.50.60.7x2 exp(-x) sin(x)由图像可以看出，在x=0,x=3,x=6附近有零点，所以命令及运行结果如下：fzero(f,0) %求在x=0附近的零点 ans = 0fzero(f,3) %求在x=3附近的零点 ans =3.1416fzero(f,6) %求在x=6附近的零点 ans =6.2832所以零点有三个分别为：x 1=0 x 2=3.1416 x 3=6.2832（2）在同一个画面上建立几个坐标系, 用subplot(m,n,p)命令；把一个画面分成m×n 个图形区域, p 代表当前的区域号，在每个区域中分别画一个图（函数为：;cos sin 2;cos ;sin x x u x z x y ===xxv cos sin =）； 命令：>> subplot(2,2,1)>> ezplot('sin(x)',[-2*pi,2*pi]) >> subplot(2,2,2)>> ezplot('cos(x)',[-2*pi,2*pi]) >> subplot(2,2,3)>> ezplot('2*sin(x)*cos(x)',[-2*pi,2*pi]) >> subplot(2,2,4)>> ezplot('sin(x)/cos(x)',[-2*pi,2*pi])图像：2、三维绘图：（1）绘制[-8 8]区间内函数2222)sin(),(yx y x y x f ++=的三维网格曲面图，给该图加上标题；（2）绘制[-8 8]区间内单叶双曲面 125169222=-+z y x 的三维网格曲面图，给该图加上标题；(3) 求函数 f(x,y) = 3x 2+10y 2+3xy-3x +2y 在原点附近的一个极小值点和极小值；作函数 f(x,y)在|x|<2, |y|<1内的图；解：（1）命令：ezsurf('sin(sqrt(x^2+y^2))/sqrt(x^2+y^2)',[-8,8])title('f=sin(sqrt(x^2+y^2))/sqrt(x^2+y^2)')图像：（2）命令：xa=-8:0.2:8; ya=xa;[x,y]=meshgrid(xa,ya); a=3;b=4;c=5;z1=sqrt(c.^2*(x.^2/(a.^2)+y.^2/(b.^2)-1)); z2=-sqrt(c.^2*(x.^2/(a.^2)+y.^2/(b.^2)-1)); surf(x,y,real(z1)) hold on ;surf(x,y,real(z2)) shading interp图像：(3)求极值点命令：>> f=inline('3*x(1)^2+10*x(2)^2+3*x(1)*x(2)-3*x(1)+2*x(2)')f =Inline function:f(x) = 3*x(1)^2+10*x(2)^2+3*x(1)*x(2)-3*x(1)+2*x(2)>> [x,fval]=fminsearch(f,[0,0])运行结果：x = 0.5946 -0.1892fval = -1.0811所以原点附近的极小值点为（0.5946，-0.1892），极小值为：-1.0811 绘图命令：>> x=-2:0.1:2;>> y=-1:0.1:1;>> [xx,yy]=meshgrid(x,y);>> zz=3.*xx.^2+10.*yy.^2+3.*xx.*yy-3.*xx+2.*yy>> surf(xx,yy,zz)图像：三、句柄图形和图形用户界面制作实验目的：1、熟悉句柄图形体系的对象树结构；2、熟练掌握句柄图形体系的对象属性的创建、设置、查询；3、熟练掌握句柄的获取；4、熟练掌握图形用户界面（GUI）的制作；实验要求：1、会获取和显示图形对象的句柄；2、会设置菜单和子菜单；3、会设置用户控件；实验内容：利用底层绘图指令绘制一条余弦曲线。
（最新版）MATLAB实验报告
（最新版）MATLAB实验报告实验一典型环节的MATLAB仿真一、实验目的1．熟悉MATLAB桌面和命令窗口，初步了解SIMULINK功能模块的使用方法。
2．通过观察典型环节在单位阶跃信号作用下的动态特性，加深对各典型环节响应曲线的理解。
3．定性了解各参数变化对典型环节动态特性的影响。
二、SIMULINK的使用MATLAB中SIMULINK是一个用来对动态系统进行建模、仿真和分析的软件包。
利用SIMULINK功能模块可以快速的建立控制系统的模型，进行仿真和调试。
1．运行MATLAB软件，在命令窗口栏“>>”提示符下键入simulink命令，按Enter 键或在工具栏单击按钮，即可进入如图1-1所示的SIMULINK仿真环境下。
2．选择File菜单下New下的Model命令，新建一个simulink 仿真环境常规模板。
3．在simulink仿真环境下，创建所需要的系统。
以图1-2所示的系统为例，说明基本设计步骤如下：1）进入线性系统模块库，构建传递函数。
点击simulink下的“Continuous”，再将右边窗口中“Transfer Fen”的图标用左键拖至新建的“untitled”窗口。
2）改变模块参数。
在simulink仿真环境“untitled”窗口中双击该图标，即可改变传递函数。
其中方括号内的数字分别为传递函数的分子、分母各次幂由高到低的系数，数字之间用空格隔开；设置完成后，选择OK，即完成该模块的设置。
3）建立其它传递函数模块。
按照上述方法，在不同的simulink 的模块库中，建立系统所需的传递函数模块。
例：比例环节用“Math”右边窗口“Gain”的图标。
4）选取阶跃信号输入函数。
用鼠标点击simulink下的“Source”，将右边窗口中“Step”图标用左键拖至新建的“untitled”窗口，形成一个阶跃函数输入模块。
5）选择输出方式。
用鼠标点击simulink下的“Sinks”，就进入输出方式模块库，通常选用“Scope”的示波器图标，将其用左键拖至新建的“untitled”窗口。
第5章卷积的原理及应用
一个周期取16个采样点，求该系统的响应。
解 编写MATLAB程序如下：
N＝16； n＝0：N－1； x＝［ones(1，N/4)，zeros(1，3*N/4)］； %产生输入 信号序列
15 subplot(3，1，1)；stem(n，x，filled)；
34
五、实验预习
(1)认真阅读实验原理部分，了解用MATLAB进行离散
时间系统卷积的基本原理、方法和步骤。 (2)读懂实验原理部分的有关例题，根据实验任务编写
实验程序。
(3)预习思考题：MATLAB中提供的conv卷积子函数， 使用中需满足什么条件？如果条件不满足，应如何处理？
35
六、实验报告
(1)列写已调试通过的实验任务程序，打印或描绘实验
果
程序执行的结果如图5-3所示。
22
图5-3 例5-3f1(n)、f2(n)、y(n)的波形
23
例5-4 已知一个IIR数字低通滤波器的系统函数公式为
0.13 21 0.39 63 1 0.39 63 2 0.13 21 3 z z z H(z) 1 0.34 31 9 1 0.60 43 9 2 0.20 40 7 3 z z z
%若f1与f2同长，不补0
%取长者为补0长度基础
%先将f2补得与f1同长，再将两边补最大长度的0
u＝［zeros(1，lt)，f2，zeros(1，nf2)，zeros(1，lt)］； t1＝(－lt＋1：2*lt)；
30 %先将f1补得与f2同长，再将左边补2倍最大长度的0
f1＝［zeros(1，2*lt)，f1，zeros(1，nf1)］；
MATLAB第五章实验报告
MATLAB数值计算班级：11电信一班姓名：何得中学号：20111060108实验目的：1.了解MATLAB提供的标准多项式的常用函数；2.掌握用MATLAB创建多项式；3.学会用MATLAB对数据的分析；实验仪器：MATLAB软件，电脑实验数据：>> poly2sym([1 2 6 -5 3 9])ans =x^5+2*x^4+6*x^3-5*x^2+3*x+9 >> A=[7 8 9;4 5 6;1 2 3];>> B=poly(A)B =1.0000 -15.0000 18.00000.0000>> poly2sym(B)ans=x^3-15*x^2+18*x+422583449577050 7/633825300114114700748351602688 >> v=[0.5 0.6i -0.4 -0.6i]v =0.5000 0 + 0.6000i -0.4000 0 - 0.6000i>> p=poly(v)p =1.0000 -0.1000 0.1600 -0.0360 -0.0720>> pr=poly2sym(p)pr=x^4-1/10*x^3+4/25*x^2-9/250*x-9/12 5>> p=poly([1 2 3 4 5]);>> polyval(p,6)ans =120>> p=[1 3 1 2 1];>> polyvalm(p,[3 6;8 4])ans = 7186 66668888 8297>> polyvalm(p,7)ans =3494>> p=[1 4 2];>> x=0:0.8:4;>> y=polyval(p,x)y = 2.0000 5.8400 10.9600 17.3600 25.0400 34.0000>> a=[1 -1 4 10];>> r=roots(a);r = 1.1879 + 2.4202i1.1879 -2.4202i-1.3758>> pp=poly(r)pp =1.0000 -1.0000 4.0000 10.0000>> p1=[1,0,1,1];>> p2=[0,1,-1,2];>> p1+p2ans =1 1 0 3>> p1-p2ans =1 -1 2 -1>> c=conv([1 1 2 2],[1 2 5 4])c = 1 3 9 15 18 18 8>> p1=[1,0,1,1]; >> p2=[1,-1,2];>> [q,r]=deconv(p1,p2) q = 1 1 r = 0 0 0 -1 >> g=[1 5 4 6 20 48 27 72 39 ]; >> h=polyder(g) h = 8 35 24 30 80 144 54 72 >> q=polyint(h) q =1 5 4 6 20 48 27 72 0 >> b=[4,9]; >> a=[1,7,8]; >> [r,p,k]=residue(b,a) r = 3.2127 0.7873 p = -5.5616 -1.4384 k = [] >> [b,a]=residue(r,p,k) b = 4.0000 9.0000 a = 1 7 8 x=linspace(0,2*pi,50); >> y=sin(x); >> [p,s]=polyfit(x,y,3) p = 0.0912 -0.8596 1.8527 -0.1649 s = R: [4x4 double] df: 46 normr: 0.5055>> plot(x,y,'r:*',x,polyval(p,x),'-o')>> a=ceil(rand(4,4)*12) a = 12 11 10 12 3 10 6 9 8 6 8 3 6 1 10 5 >> b=ceil(rand(4,4)*10) b = 10 1 2 3 10 4 3 2 5 9 2 1 9 1 7 8 >> x=a\b x = -0.6400 -0.3800 -0.5800 -0.3600 0.0814 1.4697 -0.2959 -0.7021 0.9883 1.0679 0.9748 0.7876 0.5752 -1.7738 0.2055 0.5972 >> x=[1 1.5 2 2.5 3 3.5 4 4.5 5]'; >> y=[-1.4 2.7 3 5.9 8.4 12.2 16.6 18.8 26.2]'; >> e=[ones(size(x)) x.^2]; >> c=e\y c = -1.0685 1.0627>> x1=[1:0.1:5]'; >> y2=[ones(size(x1)),x1.^2]*c;>> plot(x,y,'ro',x1,y2,'k')>> a=[2 4 7 4;9 3 5 6]; >> b=[8 5]'; >> x=pinv(a)*bx = -0.21510.44590.79490.2707>> x=0:0.02*pi:2*pi;>> y=sin(x);>> ymax=max(y)ymax = 1>> ymin=min(y)ymin = -1>> ymean=mean(y)ymean =2.2995e-017>> ystd=std(y)ystd =0.7071>> score=[90 80 78 98 79;100 90 68 78 69;88 77 99 87 100];>> score_sum=sum(score)score_sum = 278 247 245 263 248>> score_avg=score_sum/3.score_avg =92.6667 82.3333 81.6667 87.6667 82.6667>> score_max=max(score)score_max=100 90 99 98 100>>[score_max,score_max_student]=max(s core)score_max =100 90 99 98 100score_max_student = 2 2 3 1 3>> score_min=min(score)score_min = 88 77 68 78 69>>[score_min,score_min_student]=min(sc ore)score_min = 88 77 68 78 69score_min_studen t=3 3 2 2 2>> score_student_sum=sum(score) score_student_sum=278 247 245 263 248>>score_student_savg=score_student_sum. /3score_student_savg =92.6667 82.3333 81.6667 87.6667 82.6667>> std(score)ans =6.4291 6.8069 15.8219 10.0167 15.8219>> var(score)ans =41.3333 46.3333 250.3333 100.3333 250.3333>> x=linspace(0,2*pi,10)x =0 0.6981 1.3963 2.0944 2.7925 3.4907 4.1888 4.8869 5.5851 6.2832>> y=sin(x)y=0 0.6428 0.9848 0.8660 0.3420 -0.3420 -0.8660 -0.9848 -0.6428 -0.0000>> dy=diff(y)dy=0.6428 0.3420 -0.1188 -0.5240 -0.6840 -0.5240 -0.1188 0.3420 0.6428>> d2y=diff(y,2)d2y= -0.3008 -0.4608 -0.4052 -0.1600 0.1600 0.4052 0.4608 0.3008>> d3y=diff(y,3)d3y= -0.1600 0.0556 0.2452 0.3201 0.2452 0.0556 -0.1600>> x=[1 2 3 4 5];>> y=[4 5 6 7 8];>> cx=cov(x)cx = 2.5000>> cy=cov(y)cy = 2.5000>> cxy=cov(x,y)cxy =2.5000 2.50002.5000 2.5000>> x=-1:0.2:1;>> z=exp(x.^2);>> fx=gradient(z)fx= -0.8218 -0.6425 -0.3615 -0.1963 -0.0868 0 0.0868 0.1963 0.3615 0.6425 0.8218function ex=ex(x)ex=exp(-x.^2)return>> format long>> quad('ex',0,1,1e-6)ex = 1.00000000000000 0.981730032806080.928898599395100.778800783071400.588254444230630.473852624646280.36787944117144ex = 0.99540087747802 0.95936125127772ex =0.86171221264734 0.68574192538547ex = 0.53040332078406 0.41944631668514ans = 0.74682418072642>> quadl('ex',0,1,1e-6)ex=1.000000000000000.999184727905840.991616959101330.968441445328330.926451851978820.864349866175970.778800783071400.682385279905670.592380095006350.509706301048790.438272262810980.389185953389380.36787944117144ex = 0.99161695910133 0.926451851978820.778800783071400.592380095006350.43827226281098ans = 0.74682413398845>> quad8('ex',0,1,1e-6) ex =1.00000000000000 0.999184727905840.991616959101330.968441445328330.92645185197882 0.86434986617597 0.778800783071400.682385279905670.592380095006350.509706301048790.438272262810980.389185953389380.36787944117144ex= 0.99161695910133 0.926451851978820.778800783071400.592380095006350.43827226281098ans = 0.74682413398845x0=0:pi/5:4*pi;y0=sin(x0).*exp(-x0/10);xi=0:pi/20:4*pi;yi=interp1(x0,y0,xi,'spline');plot(x0,y0,'or',xi,yi,'--b')；Untitled2>> t=0:5:65;>> x=2:5:57;>>f=[3.1025,2.256,879.5,1835.9,2968.8,41 36.2,5237.9,6152.7,6725.3,6848.3,6403. 5,6824.7,7328.5,7857.6];>> f1=interp1(t,f,x)f1 = 1.0e+003 *0.002763900000000.353153600000001.262060000000002.289060000000003.435760000000004.576880000000005.603820000000006.381740000000006.774500000000006.670380000000006.571980000000007.02622000000000>> f2=interp1(t,f,x,'nearext') f2 = 1.0e+003 *0.003102500000000.002256000000000.879500000000001.835900000000002.968800000000004.136200000000005.237900000000006.152700000000006.725300000000006.848300000000006.403500000000006.82470000000000>> f3=interp1(t,f,x,'spline') f3 = 1.0e+003 *-0.170259056484550.306978324423981.256000890788642.269771808421483.439603875525464.589619889476705.637043766567756.422853844252316.859291256423016.653528330055636.481654623354457.04407557652658 >> f4=interp1(t,f,x,'cubic')f4 = 1.0e+003 *0.002438844000000.223194994900141.248358437730492.273647289639383.436483654858934.591276133775775.636234122067276.436239461848896.797756124209326.691730400000006.507716445924327.01855707680041[x,y,z]=peaks(6);mesh(x,y,z)[xi,yi]=meshgrid(-3:0.2:3,-3:0.2:3); z1=interp2(x,y,z,xi,yi,'nearest');z2=interp2(x,y,z,xi,yi,'linear');z3=interp2(x,y,z,xi,yi,'spline');z4=interp2(x,y,z,xi,yi,'cubic'); subplot(2,2,1)mesh(xi,yi,z1)title('nearest插值的网络图') subplot(2,2,2)mesh(xi,yi,z2)title('linear插值的网络图') subplot(2,2,3)mesh(xi,yi,z3)title('spline插值的网络图') subplot(2,2,4)mesh(xi,yi,z4)title('cubic插值的网络图')。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
《MATLAB原理及应用》实验报告
一．实验目的
1、掌握二维图形的绘制方法
2、掌握二维图形的标注
二．实验内容
1.绘制一条二维曲线
①MATLAB中最常用的就是plot指令，其具体调用格式为
plot(x) %绘制以x为纵坐标的二维曲线
plot(x,y) %绘制以x为横坐标y为纵坐标的二维曲线
说明：x和y可以是向量或矩阵。
【实验5-1】用plot(x)命令画直线。
>> x1=[1 2 3]
x1 =
1 2 3
>> plot(x1)
2.二维图形的标注
【实验5-2】在0≤x≤2πy1=2e-0.5x和y2=cos(4πx)
>> x=0:pi/100:2*pi;
y1=2*exp(-0.5*x);
y2=cos(4*pi*x);
plot(x,y1,x,y2),title('x from 0 to 2{\pi}'); %加图形标题。
xlabel('Variable X'); ylabel('Variable Y'); % 加X轴，Y轴说明。
text(0.8,1.5,'曲线y1=2e^(-0.5*x)'); %在指定位置添加图形说明。
text(2.5,1.1,'曲线y2=cos(4{\pi}x)');
legend('y1','y2') %加图例。
三．练习题
1、编写程序绘制如下图形（时间间隔自己设定） 024
6810-1-0.8
-0.6
-0.4
-0.2
00.20.40.6
0.8
1
sin(t)
cos(t)ÕýÏÒºÍÓàÏÒÇúÏß
ʱ¼ät ÕýÏÒ¡
¢ÓàÏ
Ò
ÕýÏÒÓàÏÒ
图中曲线标注:sin(t)标注在(1.7*pi,-0.3)，cos(t)标注在
(1.6*pi,0.8)
3、使用subplot 指令在不同区域画出下列曲线y=sinAx;y=cosAx;y=Ae x ,y=x^2。
x 的取值范围是]15,0[，A 为自己的学号。
程序如下：
4.画出衰减振荡曲线t e y t 3sin 3-=及其它的包络线30t e y -=。
t 的取值范围是]4,0[π，图的标题字体大小(fontsize)为10号，字体(fontname)宋体，字体风格为斜体1（it ）字体，并选择合适的位置在图上标注自己的学号，姓名，以及作图日期，字体大小、字体风格自己指定。