matlab上机实验题

“MATLAB及其在通信中的应用”上机实验一1

一、MATLAB环境与帮助

1．简述command window的作用。

（5分）答案:输入命令，输出结果。具体有：

（1）输入语句（2）输入变量名（3）输入系统命令/函数（4）输入.m 文件名

（5）输出上述的运算结果，包括数值或图形

2．（1）简述command history的作用；

（2）列举出四个在command window中快速操作command history内容的快捷键。

（5分）答案：（1）存储运行过的命令，便于以后调用或修改

（2）快捷键：方向键上、下、左、右。

3．简述workspace的作用，并说明workspace的菜单上的五个按钮的功能。

（5分）

答案：（1）列出程序创建的所有变量

（2）加入新变量，打开所有变量，装载数据文件，保存空间或变量，删除变量

4．简述清除command window、command history、workspace中内容的方法

（5分）答案：清除command window内容：在command window键入clc

清除command history内容：command history窗口上点右键，选择clear command history

清除工作空间的内容：在command window键入clear

5．（1）简述current directory的作用；

（2）在桌面新建一个文件夹，命名为自己的学号，如（1110920101），将current directory

1完成后以附件形式发送到邮箱commu_matlab2014@https://www.wendangku.net/doc/5b19040732.html,。邮件主题为“班级”加“下划线”加

“姓名”加“下划线”加“实验一”，如“通信1101_姓名_实验一”，word文件名类似，如“通信1101_姓名_实验一”。

改为新建文件夹的路径，简述操作方法,并对更改后的current directory截图，填入到下面空白处。

（5分）答案：（1）保存变量、m文件的默认路径。

（2）matlab快捷方式-→右键→快捷方式→起始位置

6．关闭command window，command history，current directory，workspace后，如何再显示这些窗口？如何同时显示这4个窗口并dock，这里dock是什么意思？undock呢？

（5分）答案：DeskTop->Layout->Default

dock的意思是停靠，undock的意思是非停靠

7．对含有command window，command history，current directory，workspace和help窗口（dock 后）的MATLAB界面进行截图，将截图，填入到下面空白处。（注意截图时注意大小比例，5个窗口标题栏的字迹要清楚）

（5分）

8．查看MATLAB的帮助，写出一、两个专业相关的matlab toolbox的名称。

（5分）答案：Control System Image Processing Neural Network

Communications Signal Processing

9．在MATLAB帮助里搜索sum函数和size函数，简述两个函数的功能。

（5分）答案：sum求矩阵内元素的和，size求矩阵的尺寸。

10．简述fix、floor、ceil、round、mod、rem函数的作用，并比较mod和rem函数的差别

（5分）

答案：fix ：朝零方向取整 floor ： 朝负无穷大方向取整

ceil: 朝正无穷大方向取整 round:

朝最近的整数取整（四舍五入） mod ：模数（即有符号的除后取余） rem ：除后余数

二、 MATLAB 里的矩阵和数组

1． 编辑.m 文件：以字符串形式输入自己的学号，命名为No ，并默认No 为16进制格式，将No 转化为十进制格式后，命名为No_dec 。 令n=mod(No_dec,5)+3，生成n 乘n 全1矩阵，利用whos 语句查看变量类型，所编程序并对输出结果截图，填入到下面空白处。

（5分）

2． 编辑.m 文件：生成函数3

2x 4y x =++的曲线图，定义域将所编程序及绘制的曲线图截图，填入到下面空白处。

（5分）

3．编辑.m

（1）用15位浮点格式显示y（8）的数值；

（2）用5

将所编程序和输出结果截屏，填入到到下面空白处。

（5分）答案：y(8)=1.393418503875663e+002

sum(y)=2.1562e+003 +2.0000e+000i

4. 运行如图程序。

（5分）

（1）说明= = 的运算符类型以及其与= 的差别

（2）简要解释eps的意义与功能。

（3）写出A、B、C的值，并对其运算结果进行说明

答案：（1）== 是关系运算符，=是赋值运算

（2）eps 浮点相对误差限,eps(a)是|a|与大于|a|的最小的浮点数之间的距离，距离越小表示精度越高。

（3）A=[0 1 0]; B=[1 1 1];C=[0 1 0]

计算机无法区分出1和1+eps*0.4, 可以区分出1和1+eps

计算机无法区分出2和2+eps*0.4, 无法区分出2和2+eps

计算机无法区分出2和2+eps（2）*0.4, 可以区分出2和2+eps（2）

5.一句话编程：在命令窗口，只编写一条语句，实现以下功能，将代码和结果截屏到下面

空白处

（1）生成3*4

（5分）答案：B=105+3*rand(3,4)

（2）生成2*3阶矩阵，矩阵元素服是从期望为5，方差为9的高斯分布的随机数。

（5分）答案：C=5+3*randn(2,3)

（3）矩阵

，对矩阵A进行简单操作，生

成矩阵

（5分）答案：B A=（1:9）

（4）

，对矩阵A进行简单操作，生成

矩阵

（5分）

答案：B=[A A(:,1)*4]

（5）调用A中元素及使用ones函数，生成矩阵

（5分）

答案：B=5*ones(2,2)+A(2:3,2:3)

（6）生成矩阵1*6阶的矩阵A，A中元素服从等比分布，A(1)=10，A(6)=100

（5分）答案：A=logspace(1,2,6)

Matlab上机实验答案

Matlab上机实验答案 实验一 MATLAB运算基础 1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。 >> z1=2*sin(85*pi/180)/(1+exp(2)) z1 = >> x=[2 1+2i; 5]; >> z2=1/2*log(x+sqrt(1+x^2)) z2 = - + + -

>> a=::; >> z3=(exp.*a)-exp.*a))./2.*sin(a++log(+a)./2) (>> z33=(exp*a)-exp*a))/2.*sin(a++log(+a)/2)可以验证z3==z33，是否都为1） z3 = Columns 1 through 5 + + + + + Columns 6 through 10 + + + + + Columns 11 through 15 + + + + + Columns 16 through 20 + + + + +

Columns 21 through 25 + + + + + Columns 26 through 30 + + + + + Columns 31 through 35 + + + + + Columns 36 through 40 + + + + + Columns 41 through 45 + + + + + Columns 46 through 50

+ + + + + Columns 51 through 55 + + + + + Columns 56 through 60 + + + + + Column 61 + (4) 2 2 4 2 01 112 2123 t t z t t t t t ?≤< ? =-≤< ? ?-+≤< ? ，其中t=0:: >> t=0::; >> z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3).*(t.^ 2-2.*t+1) z4 =

matlab上机习题2

实验二 基本矩阵操作 实验目的： ① 掌握matlab 变量和数据操作； ② 掌握matlab 矩阵的创立、拆分及特殊矩阵； ③ 掌握matlab 运算，掌握matlab 在矩阵分析中的应用，掌握稀疏矩阵的存储方式和创建方法； ④ 掌握矩阵的基本运算和结构体数组的创建。 ⑤ 了解字符串处理函数。 实验要求:给出程序和实验结果。 实验内容： 1、利用列向量()1,2,3,,10T 建立一个范得蒙矩阵A ，并利用位于矩阵A 的奇数行偶数列的元素建立一个新的矩阵B ，须保持这些元素的相对位置不变。 2、矩阵的基本运算与点运算的区别。 3、给出矩阵的两种存储方式的联系和区别，这两种存储方式在实际应用中主要应用于具有什么特点的矩阵？ 4*．将字符串'very good'转换为等值的整数。 5．按水平和竖直方向分别合并下述两个矩阵： 6． 分别删除第5题两个结果的第2行。 7． 分别将第5题两个结果的第2行最后3列的数值改为[11 12 13]。 8． 分别查看第5题两个结果的各方向长度。 9*． 分别判断pi 是否为字符串和浮点数。 10．分别将第5题两个结果均转换为2?9的矩阵。 11．计算第5题矩阵A 的转置。 12．分别计算第5题矩阵A 和B 的A +B 、A .* B 和A \B 。 13．判断第5题矩阵A 和B 中哪些元素值不小于4。 14*．分别用函数strcat()和矩阵合并符合并如下字符串：' The picture is '和' very good '。 15*．创建字符串数组，其中元素分别为‘Picture ’和'Pitch '。 16*．在第14题结果中查找字符串'e'。 17*．在第15题结果中匹配字符串'Pi'。

matlab作业

上机实践练习五 实验五：MATLAB 数据可视化 实验目的： 掌握MATLAB 二维、三维图形绘制，掌握图形属性的设置和图形修饰；掌握图像文件的读取和显示。 学时：2学时 实验内容： (1)二维图形绘制。 (2)隐函数绘图 (2) 三维曲线和三维曲面绘制。 (3) 图像文件的读取和显示。 1.写出实现下列图形的M文件。 t=[0:0.1:pi]; y1=sin(2*t); y2=cos(3*t); plot(t,y1,'b-o'); hold on plot(t,y2,'k:o'); xlabel('ê±??') ylabel('·ù?μ') text(1,sin(2*1),'\fontsize{16}\leftarrowsin(2t) ');

text(2.4,cos(3*(2.4)),'\fontsize{16}cos(3t)\rightarrow ',... 'HorizontalAlignment','right') 2. 将窗口分割成4个区域，并且在[] 0,2π区间上绘制1s i n 2 , 6y x π? ? =+ ?? ? ()224c o s 3,y x =()3c o s y x =，()35sin 1y x =+四条曲线，并且给每一个图形添加标 题和标注。 clear x=0:pi/10:2*pi; y1=sin(2.*x+pi/6); y2=4.*cos(3.*x.^2); y3=cos(x); y4=5.*sin(x+1); subplot(2,2,1); plot(x,y1); title('y1'); text(1,sin(2*1+pi/6),'sin(2*x+pi/6)') subplot(2,2,2); plot(x,y2) title('y2');text(2,4.*cos(3*2^2),'4*cos(3*x^2)') subplot(2,2,3); plot(x,y3)

MATLAB上机实验(答案)

MATLAB工具软件实验(1) （1）生成一个4×4的随机矩阵，求该矩阵的特征值和特征向量。程序： A=rand(4) [L,D]=eig(A) 结果： A = 0.9501 0.8913 0.8214 0.9218 0.2311 0.7621 0.4447 0.7382 0.6068 0.4565 0.6154 0.1763 0.4860 0.0185 0.7919 0.4057 L = -0.7412 -0.2729 - 0.1338i -0.2729 + 0.1338i -0.5413 -0.3955 -0.2609 - 0.4421i -0.2609 + 0.4421i 0.5416 -0.4062 -0.0833 + 0.4672i -0.0833 - 0.4672i 0.4276 -0.3595 0.6472 0.6472 -0.4804 D = 2.3230 0 0 0 0 0.0914 + 0.4586i 0 0 0 0 0.0914 - 0.4586i 0 0 0 0 0.2275 （2）给出一系列的a值，采用函数 22 22 1 25 x y a a += - 画一组椭圆。 程序： a=0.5:0.5:4.5; % a的绝对值不能大于5 t=[0:pi/50:2*pi]'; % 用参数t表示椭圆方程 X=cos(t)*a; Y=sin(t)*sqrt(25-a.^2); plot(X,Y) 结果： （3）X=[9,2,-3,-6,7,-2,1,7,4,-6,8,4,0,-2], （a）写出计算其负元素个数的程序。程序： X=[9,2,-3,-6,7,-2,1,7,4,-6,8,4,0,-2]; L=X<0; A=sum(L) 结果： A =

Matlab上机实验

Matlab 上机实验 一、 实验目的 1、 掌握绘制MATLAB 二维、三维和特殊图形的常用函数； 2、 熟悉并掌握图像输入、输出及其常用处理的函数。 二、 实验内容 1 绘制函数的网格图和等高线图。42 2cos cos y x ye x z +-= 其中x 的21个值均匀分布在[-5,5]范围，y 的31个值均匀分布在[0,10]，要求将产生的网格图和等高线图画在同一个图形窗口上。 2 绘制三维曲面图，使用纯铜色调色图阵进行着色，并进行插值着色 处理。?????===s z t s y t s x sin sin cos cos cos 230,20ππ≤≤≤≤t s 3 已知 ???????>++≤+=0),1ln(210,22x x x x e x y π 在-5<=x<=5区间绘制函数曲线。 4 已知y1=x2，y2=cos(2x)，y3＝y1*y2，其中x 为取值-2π~2π的等差数列(每次增加0.02π)，完成下列操作： a) 在同一坐标系下用不同的颜色和线型绘制三条曲线，给三条曲线添加图例； b) 以子图形式，分别用条形图、阶梯图、杆图绘制三条曲线，并分别给三个图形添加标题“y1=x^2”，“y2=cos(2x)”和“y3＝

y1*y2”。 5 在xy 平面内选择区域[][],,-?-8888 ，绘制函数 z =的三 种三维曲面图。 6 在[0,4pi]画sin(x),cos(x)(在同一个图象中); 其中cos(x)图象用红色小圆圈画.并在函数图上标注 “y=sin(x)”, “y=cos(x)” ,x 轴,y 轴,标题为“正弦余弦函数图象”. 7 分别用线框图和曲面图表现函数z=cos(x)sin(y)/y ，其中x 的取值为 [-1.5pi,1.5pi]，y=x ，要求：要有标题、坐标轴标签 8 有一组测量数据满足-at e =y ，t 的变化范围为0~10，用不同的线型和标记点画出a=0.1、a=0.2和a=0.5三种情况下的曲线，并加入标题和图列框（用代码形式生成） 9 2 2y x xe z --=，当x 和y 的取值范围均为-2到2时，用建立子窗口 的方法在同一个图形窗口中绘制出三维线图、网线图、表面图和带渲染效果的表面图 10 x= [66 49 71 56 38]，绘制饼图，并将第五个切块分离出来。 11 用sphere 函数产生球表面坐标，绘制不通明网线图、透明网线图、表面图和带剪孔的表面图。 12 以自己的个人画像或照片(JPG)为对象，读入该图像并了解图像的信息，同时利用所学函数对其进行灰度、二值、旋转及缩放等处理，并以PNG 形式输出。

matlab上机实验指导书

MATLAB应用基础实验指导书

第一章 MATLAB及其工作环境介绍 (1) 1.1 MATLAB简介 (1) 1.2 MATLAB的工作环境介绍 (1) 1.3 MATLAB的基本管理命令 (4) 第二章 MATLAB的数值计算功能 (5) 2.1 变量与赋值语句 (5) 2.2 MATLAB矩阵 (5) 2.3 MATLAB表达式 (10) 2.4 MATLAB常用数学函数 (11) 2.5 矩阵的基本运算 (12) 2.6 数组运算 (16) 2.7 多项式及其运算 (17) 第三章 MATLAB程序设计入门 (19) 3.1 M文件 (19) 3.2 数据的输入输出 (21) 3.3 全局变量和局部变量 (23) 3.4 程序流程控制 (23) 第四章 MATLAB的符号运算功能 (28) 4.1 建立符号对象 (28) 4.2 符号算术运算 (29) 4.3 符号微积分运算 (32) 4.4 符号函数的可视化 (34) 第五章 MATLAB的可视化功能 (37) 5.1 二维图形 (37) 5.2绘制三维图形 (42) 5.3 特殊坐标图形 (44) 5.4 图形句柄 (45)

第一章 MATLAB及其工作环境介绍 1.1 MATLAB简介 MATLAB是matrix和laboratory前三个字母的缩写，意思是实验室矩阵。MATLAB 语言是一种广泛应用于工程计算及数值分析领域的新型高级语言，自1984年由美国MathWorks公司推向市场以来，经过十多年的发展与完善，MATLAB已发展成为由MATLAB语言、MATLAB工作环境、MATLAB图象处理系统、MATLAB数学函数库和MATLAB 应用程序接口五大部分组成的集数值计算、图形处理、程序开发为一体的功能强大的体系。MATLAB由“主包”和三十多个扩展功能和应用学科性的工具箱组成。 MATLAB具有以下基本功能： ●数值计算功能 ●符号计算功能 ●图形处理及可视化功能 ●可视化建模及动态仿真功能 MATLAB语言是以矩阵计算为基础的程序设计语言，语法规则简单易学。其指令格式与数学表达式非常相近，用MATLAB编写程序犹如在便笺上列写公式和求解，因而被称为“便笺式”的编程语言。另外，MATLAB还具有功能丰富和完备的数学函数库及工具箱，大量繁杂的数学运算和分析可通过调用MATLAB函数直接求解，大大提高效率，其程序编译和执行速度远远超过了传统的C和FORTRAN语言，因而用MATLAB 编写程序，往往可以达到事半功倍的效果。在图形处理方面，MATLAB可以给数据以二维、三维乃至四维的直观表现，并在图形色彩、视角、品性等方面具有较强的渲染和控制能力，使技术人员对大量原始数据的分析变得轻松和得心应手。 MATLAB的上述特点，使它深受工程技术人员及科技专家的欢迎，并成为应用学科计算机辅助分析、设计、仿真、教学等领域不可缺少的基础软件。目前MATLAB已成为国际上公认的最优秀的科技应用软件。 1.2 MATLAB的工作环境介绍 一、MATLAB的工作环境

河南城建学院MATLAB上机实验答案

一熟悉Matlab工作环境 1、熟悉Matlab的5个基本窗口 思考题： （1）变量如何声明，变量名须遵守什么规则、是否区分大小写。 答：变量一般不需事先对变量的数据类型进行声明，系统会依据变量被赋值的类型自动进行类型识别，也就是说变量可以直接赋值而不用提前声明。变量名要遵守以下几条规则：?变量名必须以字母开头，只能由字母、数字或下划线组成。 ?变量名区分大小写。 ?变量名不能超过63个字符。 ?关键字不能作为变量名。 ?最好不要用特殊常量作为变量名。 （2）试说明分号、逗号、冒号的用法。 分号：分隔不想显示计算结果的各语句；矩阵行与行的分隔符。 逗号：分隔欲显示计算结果的各语句；变量分隔符；矩阵一行中各元素间的分隔符。 冒号：用于生成一维数值数组；表示一维数组的全部元素或多维数组某一维的全部元素。 （3）linspace（）称为“线性等分”函数，说明它的用法。 LINSPACE Linearly spaced vector. 线性等分函数 LINSPACE(X1, X2) generates a row vector of 100 linearly equally spaced points between X1 and X2. 以X1为首元素，X2为末元素平均生成100个元素的行向量。 LINSPACE(X1, X2, N) generates N points between X1 and X2. For N < 2, LINSPACE returns X2. 以X1为首元素，X2为末元素平均生成n个元素的行向量。如果n＜2，返回X2。 Class support for inputs X1,X2: float: double, single 数据类型：单精度、双精度浮点型。 （4）说明函数ones（）、zeros（）、eye（）的用法。 ones（）生成全1矩阵。 zeros（）生成全0矩阵。 eye（）生成单位矩阵。 2、Matlab的数值显示格式

matlab上机实验五

上机五 1、 众所周知，水是地球上所有生命赖以生存的基础。没有水，一切生命创造的精彩都将不复存在。当今世界，经济在高速发展，我们对于水需求更大，然而我们却在面临前所未有的水危机，水污染的恶化更使水短缺雪上加霜。我们的水资源正在遭受各种污染的侵袭，水污染严重破坏生态环境、影响人类生存，要想实现人类社会的可持续发展，首先要解决水污染问题。 由有害化学物质造成水的使用价值降低或丧失称之为水污染。水的污染有两类：一类是自然污染；另一类是人为污染。而后者是主要的。水污染可根据污染杂质的不同而主要分为化学性污染、物理性污染和生物性污染三大类。水中杂质按尺寸分，可分为溶解物、胶体颗粒和悬浮物3种。有些杂质可以用基于高浓度、外加计量反应试剂为基础的传统的物化方法(如沉降、吸附、湿式氧化等) 以及生化技术等进行处理。而对于天然水体和饮用水中低浓度、高毒性、难降解污染物 (如多溴联苯醚、全氟辛酸(磺酸)、消毒副产物、内分泌干扰物、PPCPs(抗生素)等) 很难用前述传统的物化方法和生化技术等技术进行处理，迫切需要提出建立新型的高效选择性检测和消除的原理和方法。 问题：附件中给定的数据是利用动态光反射仪器测量出水中某污染物粒径随时间的变化值，请就给定的数据拟合出粒径随时间变化的曲线和分布。尝试拟合出相应的函数？ Time (s) Aggregation Size (nm纳米) 0 50.72 21 67.96 42 73.13 63 82.6 84 83.6 105 89.09 126 93.16 147 101.1 168 106.3 189 104.6

210 108.6 231 115.1 252 112 273 122.1 294 132 315 132.9 336 131.9 357 128.7 378 142.8 399 142.6 420 152.7 441 152.6 462 147.8 483 149.5 504 156.2 525 170.1 546 159.4 567 167.3 588 171.5 609 165.1 630 178.6 651 174.5 672 169.8 693 174.4 714 179 735 176.1 756 165.6 777 166 798 181.4 819 190.2 840 185.2 861 188.4 882 187.4 903 197.9 924 193.7 945 210.6 966 212.5 987 196.5 1008 195.3 1029 215.3 1050 208.2 1071 217.5 1092 206.7 1113 232.3

MATLAB上机实验练习题答案

数学建模 MATLAB上机实验练习题 1、给出一个系数矩阵A[234;541;132]，U=[123]，求出线性方程组的一个精确解。 2、给出两组数据x=[00.30.81.11.62.3]’y=[0.820.720.630.600.550.50]’，我们可以简单的 认为这组数据在一条衰减的指数函数曲线上，y=C1+C2e-t通过曲线拟合求出这条衰减曲线的表达式，并且在图形窗口画出这条曲线，已知的点用*表示。 3、解线性方程 4、通过测量得到一组数据： 5、已知一组测量值 6、从某一个过程中通过测量得到： 分别采用多项式和指数函数进行曲线拟合。 7、将一个窗口分成四个子窗口，分别用四种方法做出多峰函数的表面图（原始数据法，临近 插值法，双线性插值法，二重三次方插值法） 8、在同一窗口使用函数作图的方法绘出正弦、余弦、双曲正弦、双曲余弦。分别使用不同的 颜色，线形和标识符。 9、下面的矩阵X表示三种产品五年内的销售额，用函数pie显示每种产品在五年内的销售额

占总销售额的比例，并分离第三种产品的切片。 X=19.322.151.6 34.270.382.4 61.482.990.8 50.554.959.1 29.436.347.0 10、对应时间矢量t，测得一组矢量y t00.30.8 1.1 1.6 2.3 y0.50.82 1.14 1.25 1.35 1.40 采用一个带有线性参数的指数函数进行拟合，y=a0+a1e-t+a2te-t，利用回归方法求出拟合函数，并画出拟合曲线，已知点用圆点表示。 11、请创建如图所示的结构数组（9分） 姓名编号指标 江明顺071023身高：176，体重：82 于越忠060134身高：168，体重：74 邓拓050839身高：182，体重：77 12、创建如图所示的元胞数组。（9分） 13、某钢材厂从1990年到2010年的产量如下表所示，请利用三次样条插值的方法计算1999年该钢材厂的产量，并画出曲线，已知数据用‘*’表示。要求写出达到题目要求的MATLAB 操作过程，不要求计算结果。 年份19901992199419961998200020022004200620082010 产量（万吨）75.99591.972105.711123.203131.669150.697179.323203.212226.505249.633256.344 14、在一次化学动力学实验中，在某温度下乙醇溶液中，两种化合物反应的产物浓度与反应时间关系的原始数据如下，请对这组数据进行三次多项式拟合，并画出拟合曲线，已知数据如下。 time=[2.55.07.510.013.017.020.030.040.050.060.070.0] res=[0.290.560.771.051.361.522.002.272.813.053.253.56]

Matlab上机实验答案 (1)

Matlab上机实验答案 实验一MATLAB运算基础 1. 先求下列表达式的值,然后显示MATLAB工作空间的使用情况并保存全部变量。 >> z1=2*sin(85*pi/200)/(1+exp(2)) z1 = 0.2375 >> x=[2 1+2i;-0.45 5]; >> z2=1/2*log(x+sqrt(1+x^2)) z2 = 0.7120 - 0.0253i 0.8968 + 0.3658i 0.2209 + 0.9343i 1.2041 - 0.0044i 2.9,,2.9, 3.0

>> 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) (>> z33=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2)可以验证z3==z33,是否都为1） z3 = Columns 1 through 5 0.7388 + 3.2020i 0.7696 + 3.2020i 0.7871 + 3.2020i 0.7920 + 3.2020i 0.7822 + 3.2020i Columns 6 through 10 0.7602 + 3.2020i 0.7254 + 3.2020i 0.6784 + 3.2020i 0.6206 + 3.2020i 0.5496 + 3.2020i Columns 11 through 20 0.4688 + 3.2020i 0.3780 + 3.2020i 0.2775 + 3.2020i 0.2080 + 3.2020i 0.0497 + 3.2020i

电子科技大学matlab-上机实验五

Matlab 与数值分析实验报告（五） 沈汉男 一、实验的性质和任务 通过上机实验，使学生对数值积分、微分方程求解方法有一个初步的理解。 二、实验内容 1.选用复合Simpson 公式，计算 并用Matlab 的符号运算工具箱计算其精确值。比较结果，找出问题原理，提出解决问题的方法。 2.求积分方程 t t e ds s y e e t y --=?10 )(12)(的数值解和精确解，分析二者的差异。 3.利用Euler 法对不同的步长求下面初值问题的数值解： ???=-=10 )0(20)()('y t y t y 并通过绘图，与方程的解析解进行比较。 三、实验过程 实验1 a.实验过程： （1）、创建文件simpr1.m 如下： function s=simpr1(f,a,b,n)h=(b-a)/(2*n);s1=0;s2=0;for k=1:n x=a+h*(2*k-1);s1=s1+feval('f',x);end for k=1:(n-1)x=a+h*2*k;s2=s2+feval('f',x);end s=h*(feval('f',a)+feval('f',b)+4*s1+2*s2)/3 dx x x ?--+1 1)1ln()1ln(

创建文件f.m 如下： function y=f(x)y=(log(1+x))*(log(1-x)); 在主界面中输入 simpr1('f',-1,1,10) 所得结果为： s = -Inf ans = -Inf （2）在主界面中输入程序 clear all syms x y=int(log(1+x)*log(1-x),'x','-1','1') eval(y) 所得结果为 y = 2*log(2)^2-pi^2/3-log(16)+4 ans = -1.1016 b.实验结果分析： 由上面的实验结果可以看出有时候复合Simpson 公式并不能求出所需解实验2： a.实验过程：原积分方程为 t t e ds s y e e t y --=?10 )(12)(对其进行化简和变形，得： ?? ???+==21)0(' e y y y 原问题转化为一个常微分初值问题。

Matlab编程与应用习题和一些参考答案

Matlab 上机实验一、二 3.求下列联立方程的解???????=+-+-=-+=++-=--+4 1025695842475412743w z y x w z x w z y x w z y x >> a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10]; >> b=[4;4;9;4]; >> c=a\b 4．设???? ??????------=81272956313841A ，??????????-----=793183262345B ，求C1=A*B’;C2=A’*B;C3=A.*B,并求上述所有方阵的逆阵。 >> A=[1 4 8 13;-3 6 -5 -9;2 -7 -12 -8]; >> B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7]; >> C1=A*B' >> C2=A'*B >> C3=A.*B >> inv(C1) >> inv(C2) >> inv(C3) 5．设 ?? ????++=)1(sin 35.0cos 2x x x y ，把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线。 >> x=linspace(0,2*pi,101); >> y=cos(x)*(0.5+(1+x.^2)\3*sin(x)); >> plot(x,y,'r') 6．产生8×6阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。并求该矩阵全体数的平均值和均方差。 （mean var ） a=randn(8,6) mean(a) var(a) k=mean(a) k1=mean(k) i=ones(8,6) i1=i*k1 i2=a-i1 i3=i2.*i2 g=mean(i3) g2=mean(g)

MATLAB上机习题一

MATLAB上机习题一 请按以下步骤完成上机实验： 1）在FTP上下载“MATLAB上机习题一.doc”文件，所有习题列在该文件内； 2）在MATLAB中完成所有习题，并将屏幕截图粘贴到相应习题后面； 3）如果习题是问答题，请将答案写在题目后； 4）如果有的习题要求提供脚本文件，请将脚本文件内容拷贝到相应习题后； 5）将文件保存并重命名为“自己的学号-姓名”，例如“20110771-张三.doc”； 6）上传该文件到FTP的相关目录。 1. 运行MATLAB软件，拖放、关闭界面上的子窗口，并恢复到原始试图。 2. 采用鼠标及命令两种方式将桌面添加到MATLAB搜索路径列表的起始及最后位置。 3. 采用鼠标及命令方式将当前工作目录设置为桌面指向的文件夹。 4. 通过使用帮助确定内置变量ispc 的含义。 6. 观察MATLAB中关键字、字符串、注释的字体显示颜色。 5. 创建double类型的变量，并进行计算 1）a=87，b=190，计算a+b、a-b、a*b；

2）创建uint8类型的变量，数值与（1）中相同，进行相同的计算，观察计算结果与预想的是否一致。 6. 计算如下表达式： 1）()sin 60 2）3e 3）3cos 4π?? ??? 4）2 7562323336 +?-?? 7. 设u=2，v=3，计算： 1）4log uv v 2）() 2 2u e v v u +- 3 8. 计算如下表达式： 1）()()3542i i -+ 2）()sin 28i - 9. 判断下面语句的运算结果，并与MATLAB 计算结果做比较：

1）420< 2）420≤ 3）420== 4）4~20= 5）''''b B < 10. 设a=39，b=58，c=3，d=7，判断下面表达式的值与MATLAB 计算结果做比较： 1）a b > 2）a b < 3）&&a b b c >> 4）a d == 5）|a b c > 6）~~d 11. 判断下列变量名哪些是合法的MATLAB 变量： 1）fred 7）fred! 2）book_1 8）book-1 3）2ndplace 9）Second_Place 4）#1 10）No_1 5）vel_5 11）vel.5 6）tan 12）while 提示，可以使用isvarname 对上述变量名进行检验。 12. 编写脚本程序，命名为barbell.m ，完成如下计算： 如图1所示的杠铃，

南邮Matlab上机实验报告五

Matlab上机实验报告 实验名称：MATLAB语言平台与SIMULINK工具箱班级：自动化二班 学号： B11050216 姓名：李鹏飞 南京邮电大学

2013年5月13日 一、目的与任务 1．学习了解MATLAB语言环境； 2．练习MATLAB命令； 4．掌握SIMULINK工具箱的应用。 二、内容与要求 1．实验内容 （1）学习了解MATLAB语言环境； （2）练习MATLAB命令的基本操作； （3）掌握SIMULINK仿真环境的使用方法; (4) 掌握线性系统仿真常用基本模块的用法 2．要求 按照实验文档的要求与步骤完成实验，撰写实验报告。 三、实验基本知识： 1.熟悉MATLAB环境： MATLAB桌面和命令窗口、命令历史窗口、帮助信息浏览器、工作空间浏览器、文件和搜索路径浏览器。 2.掌握MATLAB常用命令 3.了解SIMULINK模块库中各子模块基本功能 四、上机练习 1.学习使用help命令，例如在命令窗口输入help eye，然后根据帮助说明，学习使用指令eye（其它不会用的指令，依照此方法类推）>> help eye eye Identity matrix. eye(N) is the N-by-N identity matrix. eye(M,N) or eye([M,N]) is an M-by-N matrix with 1's on the diagonal and zeros elsewhere. eye(SIZE(A)) is the same size as A. eye with no arguments is the scalar 1. eye(..., CLASSNAME) is a matrix with ones of class specified by CLASSNAME on the diagonal and zeros elsewhere. Note: The size inputs M and N should be nonnegative integers. Negative integers are treated as 0. Example: x = eye(2,3,'int8');

Matlab上机实验二

1. 对Runge 函数2()1/(125)R x x =+在区间[-1,1]作下列插值逼近，并和()R x 的图 像进行比较，并对结果进行分析。 （1） 用等距节点1,0.1,020i x ih h i =-+=≤≤，绘出它的20次Newton 插值 多项式的图像。 （2） 用节点21cos()(0,1,2...,20)42 i i x i π+==，绘出它的20次Lagrange 插值多项式的图像。

（3） 用等距节点1,0.1,020i x ih h i =-+=≤≤，绘出它的分段线性插值函数 的图像。 （4） 用等距节点1,0.1,020i x ih h i =-+=≤≤，绘出它的三次自然样条插值 函数的图像。

结果分析： 等距节点的Newton插值在边界附近呈现不收敛的状态，也就是所谓的Runge 现象。这是因为被插函数2 ()1/(125) R x x =+及其任意阶导数并不是一致有界的。 由Lagrange插值的余项 (n1) 1 () ()()()() (1)! n n n f R x f x L x x n ξ ω + + =-= + ，由于高阶导数值 非常大，在等距节点的高次插值中，余项的值很大，亦即误差特别大，表现为Runge现象。在实际中，除非被插函数的光滑性十分好，否则等距节点的高次插值一般不使用。 切比雪夫多项式零点插值没有出现Runge现象，这是因为对于Lagrange插 值的余项 (n1) 1 () ()()()() (1)! n n n f R x f x L x x n ξ ω + + =-= + ，切比雪夫多项式零点插值可以 将 1() n x ω +项控制在 1 2n 以下（ 1 [1,1] max() n x x ω + ∈- = 1 2n ），由此可得切比雪夫多项式零点 插值的误差范围 (n1) [1,1] () max() 2(1)! n n x f R x n ξ + ∞ ∈- ≤ + ，比等距节点的高次插值小很多，因此 对于被插函数的光滑性要求要比等距节点的插值也低很多。 分段低次线性插值是常用的一种插值方式。线性插值的优点在于计算量小，但是从图中可以看出，其光滑性不够。三次样条插值要求其一阶导数是连续的，因此其准确性好，光滑性也较优。

MATLAB上机题汇总(基本题,抽3题)合肥学院

1设置matlab 的工作环境，将工作目录设置为d:\work ，添加搜索目录d:\example 2在matlab 的命令窗口里完成如下计算，其中t 的值分别取-1,0,1，表达式如下： 4/3)2 t y e π-= 3自行产生一个5行5列的数组，得到最中间的三行三行矩阵。 4用magic 产生一个5*5的矩阵，将这个矩阵的第二行与第三行互换位置 5求方程组的根 x 1+4x 2-3x 3=2 2x 1+5x 2-x 3=11 x 1+6x 2+x 3=12 6已知：一个多项式的系数向量是p=[1 -6-72 -27]，求这个多项式的根。 7已经两个多项式的系数分别是：[1 2 3 4]和[1 4 9 16]，请求这两个多项式的乘积，及商和余数。 8给定一个多项式的根是[-5 -3+4i -3-4i],求原来的多项式 9 A=[2 3 4； 1 5 7； 6 2 5] 用什么函数，保证第一列排序的时候，其他列跟着变化。最后的结果是： [1 5 7； 2 3 4； 6 2 5] 10补充题：电路分析 电路分析常常涉及到对方程组的求解，这些方程常常是利用描述电流进入和离开节点的电流方程，或者描述电路中网络回路上的电压的电压方程得出的。下图描述了3个回路电压的方程

方程式如下： 11121222132423432 532()0()()0()0V R I R I I R I I R I R I I R I I R I V -++-=??-++-=??-++=? 假设5个电阻值为已知，2个电压值也为已知，求3个电流值。 11自行产生一个5行5列的数组，用两种方法得到最中间的三行三行矩阵。 12根据a=reshape(-4:4,3,3)做一个矩阵，然后（1）取出所有大于0的元素构成一个向量（可推广到从一个矩阵里筛选出符合条件的元素组成一个向量）（2）将原矩阵中大于0的元素正常显示，而小于等于0的元素全部用0来表示(可推广到将原矩阵中不符合条件的全用0来表示，符合条件的值不变)。(思考：大于0的正常显示，小于等于0的用-10来表示) 13建立如下一个元胞数组，现在要求计算第一个元胞第4行第2列加上第二个元胞+第三个元胞里的第二个元素+最后一个元胞的第二个元素。 a={pascal(4),'hello';17.3500,7:2:100} 14建立一个结构体的数组，包括3个人，字段有姓名，年龄，分数，其中分数由随机函数产生一个3行10列的数据表示了有10门课程，每门课程有三个阶段的分数。问题是： 问题1，如何找到第2个人的分数并显示出来 问题2，如何找到第2个人的每门课程3个阶段的平均分数并显示出来 问题3，全班同学（指这3个学生）的10门课程的每门课程的平均分如何计算出来？要求放到一个数组里。 问题4，找到这个班所有同学的姓名放到一个元胞数组里 15给定一个图像文件，格式是jpg ，通过inportdata 引入这个文件，查看引入后数据保存是一个数组，是100*100*3的一个数组，这是一个三维的数组，表明有100*100个点，每个点有RGB 三个方面决定的。现在，要求将这些点的中间部分30个象素宽度的正方形用黑点表示。 16手动构造一个长宽页是1000*1000*3的图像，每个图像的点有三个0-255之间的随机值构成。然后使用image 函数显示这个图像，再用imwrite 函数保存这个图像，使这个图像用图片浏览器也能够查看。 17给定一个矩阵 a = 1 2 3 4 5 6 编写一个M 函数，要求输入是a ，输出有三个：平均数，标准差，秩。程序运行

matlab上机习题5matlab7.0二维绘图

实验五二维绘图 实验目的： ①掌握绘制数据曲线图的方法； ②掌握绘制其他坐标系下的数据曲线图和统计分析图的方法； ③掌握绘制隐函数图形的方法。 ④掌握图形修饰处理方法； 实验要求:给出程序和实验结果。 实验内容： 8. 编制程序，该程序绘制两条曲线，x的取值在[0,2pi]，易pi/10为步长，一条是正弦曲线，一条是余弦曲线，线宽为6个象素，正弦曲线为绿色，余弦曲线为红色，线型分别为实线和虚线。给所绘的两条曲线增添图例，分别为“正弦曲线”和“余弦曲线”。 9. 在同一坐标内，分别用不同线型和颜色绘制曲线y1= 和y2=(πx)，标记两曲线交叉点。 10. 在0≤x≤2区间内，绘制曲线y1=和y2=cos(4πx)，并给图形添加图形标注。 11.重新绘制第一题所描述的曲线，将正弦曲线和余弦曲线分别画在两个子图中，子图竖向排列。

12、绘制r=sin(t)cos(t)的极坐标图； 13、分别以条形图、阶梯图、杆图和填充图形式绘制曲线y=2sin(x)。 实验程序与结果： 1 x=-2::2; y=sin(x).*cos(x); plot(x,y,'-r') -2 -1.5 -1 -0.5 0.5 1 1.5 2 -0.5-0.4-0.3-0.2-0.100.10.20.30.4 0.5 2 ezplot('x^2/9+y^2/16-1',[-5,5,-5,5]);

x y x 2/9+y 2/16-1 = 0 -5 -4 -3 -2 -1 01 2 3 4 5 -5-4-3-2-101234 5 3 x1=-2::2; x2=-2::2; y1=sin(x2).*x1; y2=cos(x1).*x2; plot3(x1,x2,y1,'d',x1,x2,y2,'d')

D实验五 M文件和MATLAB程序设计

实验五M文件和MATLAB程序设计 一、实验目的 matlab作为一种高级计算机语言，不仅可以命令行方式完成操作，也具有数据结构、控制流、输入输出等能力，本次实验通过熟悉和掌握m文件的建立与使用方法，以及函数与控制程序流程语句的使用，使学生具备一定的编程和程序调试能力。 1．掌握M文件的使用方法。 2．掌握if语句和switch语句的使用 3. 掌握循环语句的使用 4. 通过练习理解MATLAB编程方法。 二、实验原理 1．m文件 用matlab语言编写的程序，称为m文件。M文件根据调用方式的不同分为两类，命令文件（Script file）和函数文件（Function file）。区别? 2．程序控制结构 1）顺序结构 2）选择结构 （1）if语句a) 单分支if语句b) 双分支if语句c) 多分支if语句 （2）switch 语句 （3）try语句 3）循环结构 （1）for 语句 （2）while语句 （3）break语句、continue语句、return使用,区别? 3．函数文件 function 输出形参表＝函数名（输入形参表） 注释说明部分 函数体语句 注意事项？ 三、实验要求 1．首先上机练习PPT中各种流程控制语句的有关实例。 2．然后上机练习下面的实验习题。 四、实验习题

1．数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。重复此过程，最终得到的结果为1。如： 2→1 3→10→5→16→8→4→2→1 6→3→10→5→16→8→4→2→1 运行下面的程序，按程序提示输入n=1,2,3,5,7，8,9等数来验证这一结论。 %classic "3n+1" problem from number theory. while 1 n=input('Enter n,negative quits:'); if n<=0 break end a=n; while n>1 if rem(n,2)==0 n=n/2; else n=3*n+1; end a=[a,n]; end a end Enter n,negative quits: 2． 编程求满足∑=>m i i 11000020的最小m 值。 a=0; i=1; while (a<100000) a=a+pow2(i); i=i+1; end m=i-1 3． 编写一个函数，计算下面函数的值，给出x 的值，调用该函数后，返回y 的值。 function [y]=myfun1(x)

MATLAB实验题答案

1、求以下变量的值，并在MATLAB中验证。( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ；a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 8 2、下列运算是否合法，为什么？如合法， 结果是多少？ >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 2 5 8 11 >> result4=b*d result4 = 31 22 22 40 49 13 >> result5=[b;c']*d result5 = 31 22 22 40 49 13 -5 -8 7 >> result6=a.*b result6 = 2 8 -3 4 1 5 30 >> result7=a./b result7 = 0.5000 0.5000 -3.0000 4.0000 1.6667 1.2000 >> result8=a.c Attempt to reference field of non-structure array. >> result9=a.\b result9 = 2.0000 2.0000 -0.3333 0.2500 0.6000 0.8333 >> result10=a.^2 result10 = 1 4 9 16 25 36 >> result11=2.^a result11 = 2 4 8 16 32 64 3、用MATLAB求解下面的的方程组。 (1) ? ? ? ? ? ? ? ? ? ? ? ? - = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? - - - - 1 7 4 13 2 3 1 5 11 2 2 2 3 15 9 2 1 2 7 4 3 2 1 x x x x >> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0] >> B=B' >> x=inv(A)*B (2) ? ? ? ? ? ? ? = - + + = - - = - + + = + + 5 6 5 3 3 3 3 2 8 2 1 w z y x w y x w z y x z y x >> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6] >> B2=[1;8;3;5] >> x2=inv(A1)*B2 4、已知 ? ? ? ? ? ? ? ? ? ? ? ? - - - - = 13 2 3 1 5 11 2 2 2 3 15 9 2 1 2 7 A