MATLAB实验报告(二) 西安邮电大学
MATLAB实验报告
MATLAB实验报告一、实验目的本次 MATLAB 实验旨在深入了解和掌握 MATLAB 软件的基本操作和应用，通过实际编程和数据处理，提高解决问题的能力，培养编程思维和逻辑分析能力。
二、实验环境本次实验使用的是 MATLAB R2020a 版本，运行在 Windows 10 操作系统上。
计算机配置为英特尔酷睿 i5 处理器，8GB 内存。
三、实验内容（一）矩阵运算1、矩阵的创建使用直接输入、函数生成和从外部文件导入等方式创建矩阵。
例如，通过｀1 2 3; 4 5 6; 7 8 9` 直接输入创建一个 3 行 3 列的矩阵；使用｀ones(3,3)｀函数创建一个 3 行 3 列元素全为 1 的矩阵。
2、矩阵的基本运算包括矩阵的加减乘除、求逆、转置等。
例如，对于两个相同维度的矩阵｀A` 和｀B` ，可以进行加法运算｀C ＝ A ＋ B` 。
3、矩阵的特征值和特征向量计算通过｀eig` 函数计算矩阵的特征值和特征向量，加深对线性代数知识的理解和应用。
（二）函数编写1、自定义函数使用｀function` 关键字定义自己的函数，例如编写一个计算两个数之和的函数｀function s ＝ add(a,b) s ＝ a ＋ b; end` 。
2、函数的调用在主程序中调用自定义函数，并传递参数进行计算。
3、函数的参数传递了解值传递和引用传递的区别，以及如何根据实际需求选择合适的参数传递方式。
（三）绘图功能1、二维图形绘制使用｀plot` 函数绘制简单的折线图、曲线等，如｀x ＝ 0:01:2pi; y ＝ sin(x)； plot(x,y)｀绘制正弦曲线。
2、图形的修饰通过设置坐标轴范围、标题、标签、线条颜色和样式等属性，使图形更加清晰和美观。
3、三维图形绘制尝试使用｀mesh` 、｀surf` 等函数绘制三维图形，如绘制一个球面｀x,y,z ＝ sphere(50)； surf(x,y,z)｀。
（四）数据处理与分析1、数据的读取和写入使用｀load` 和｀save` 函数从外部文件读取数据和将数据保存到文件中。
MATLAB实验报告
MATLAB实验报告实验报告课程名称 MATLAB基础及应用专业班级电子xxxx姓名学号电气与信息学院实验二 MATLAB 数值计算（一）一实验目的：1．掌握数组的创建与运算方法；2. 掌握矩阵的创建与运算方法；3. 掌握数组的运算方法和矩阵运算方法的区别；4．掌握线性方程的求解方法二实验装置：计算机三实验内容：1．数组的创建和运算创建两个含5个元素的一维数组，并求这两个数组的四则运算。
2．矩阵的创建和运算（1）创建两个3×3的矩阵，并求这两个矩阵的四则运算。
（2）创建一个4×4的矩阵，并求这个矩阵的行列式值。
（3）线性方程的求解解方程组=6613753467294x 。
四实验要求:写出实验程序与仿真结果。
1. >> a=linspace(2,10,5) a =2 4 6 8 10>> b=linspace(1,9,5) b =1 3 5 7 9>> a+bans =3 7 11 15 19>> a-bans =1 1 1 1 1>> a.*bans =2 12 30 56 90>> a./bans =2.0000 1.3333 1.2000 1.1429 1.1111 >> a.\b ans =0.5000 0.7500 0.8333 0.8750 0.9000 2.(1) >> a=[1 1 1;2 2 2;3 3 3];>> b=[4 4 4;5 5 5;6 6 6];>> a+bans =5 5 57 7 79 9 9>> a-bans =-3 -3 -3-3 -3 -3-3 -3 -3>> a*bans =15 15 1530 30 3045 45 45>> a/bWarning: Matrix is singular to working precision. ans =NaN NaN NaNNaN NaN NaNNaN NaN NaN>> a\bWarning: Matrix is singular to working precision. ans =NaN NaN NaNNaN NaN NaNInf Inf Inf(2).>> a=[1 1 1 1;2 2 2 2;3 3 3 3;4 4 4 4];>> det(a)ans =(3).>> a=[4 9 2;7 6 4;3 5 7];>> b=[13;6;6];>> x=inv(a)*bx =-0.53081.7109-0.1374实验三 MATLAB数值计算（二）一实验目的：1．掌握多项式的创建与运算方法；2. 掌握基本的数据分析方法；二实验装置：计算机三实验内容：1．多项式创建输入系数矢量，创建多项式x^3-2*x^2+5*x+3。
matlab实验报告
Matlab实验报告实验二图像处理一、实验目的（1）通过应用MA TLAB语言编程实现对图像的处理，进一步熟悉MATLAB软件的编程及应用；（2）通过实验进一步掌握图像处理的基本技术和方法。
二、实验内容及代码㈠.应用MA TLAB语言编写显示一幅灰度图像、二值图像、索引图像及彩色图像的程序，并进行相互之间的转换首先，在matlab页面中的current directory下打开存放图像的文件夹。
1.显示各种图像⑴显示彩色图像：①代码：>> mousetif=imread('tif.TIF');>> image(mousetif)显示截图：②代码：>> mousetif=imread('tif.TIF');>> imshow(mousetif)显示截图：③代码：mousetif=imread('tif.TIF');subimage(mousetif)显示截图：显示截图：⑵显示二值图像①代码：>> I=imread('单色bmp.bmp');>> imagesc(I,[0 2])显示截图：②代码：>> I=imread('单色bmp.bmp');>> imshow(I,2)显示截图：③代码：>> I=imread('单色bmp.bmp');>> subimage(I)显示截图：⑶显示灰度图像①代码：>> I1=imread('256bmp.bmp');>> imagesc(I1,[0,256])显示截图：代码：>> I1=imread('256bmp.bmp');>> colormap(gray);>> subplot(1,2,1);>> imagesc(I1,[0,256]);>> title('灰度级为[0 256]的mouse.bmp图');>> subplot(1,2,2);>> imagesc(I1,[0,64]);>> colormap(gray);>> title('灰度级为[0 64]的mouse.bmp图');显示截图：⑷显示索引图像代码：>> mousetif=imread('tif.TIF');>> [x,map]=rgb2ind(mousetif,0.7);>> subplot(1,2,1);>> imshow(mousetif);>> title('彩色图像');>> subplot(1,2,2);>> imshow(x,map);>> title('索引图');截图显示：彩色图像索引图2.图像之间的转换⑴索引图像转灰度图像代码：>> II1=ind2gray(x,map);>> imshow(x,map);>> title('索引图像');>> imshow(II1);>> title('灰度图像');截图显示：索引图像灰度图像⑵ 索引图像转彩色图像⑶ 索引图像转二值图像⑷ 灰度图像转索引图像代码：>> I1=imread('256bmp.bmp');>> subplot(1,2,1);>> imagesc(I1,[0,256]); >> colormap(gray); >> title('灰度图');>> [x,map]=gray2ind(I1,256); >> subplot(1,2,2); >> image(x);>> colormap(map); >> title('索引图')；显示截图：灰度图2004006008001000100200300400500600索引图2004006008001000100200300400500600⑸灰度图像转彩色图像⑹灰度图像转二值图像⑺彩色图像转灰度图像⑻彩色图像转索引图像代码：>> mousetif=imread('tif.TIF');>> [x,map]=rgb2ind(mousetif,0.7);>> subplot(1,2,1);>> imshow(mousetif);>> title('彩色图像');>> subplot(1,2,2);>> imshow(x,map);>> title('索引图');显示截图：彩色图像索引图⑼彩色图像转二值图像⑽二值图像转索引图像⑾二值图像转彩色图像⑿二值图像转灰度图像㈡应用MA TLAB工具箱演示一幅图像的傅里叶变换、离散余弦变换，观察其频谱图。
Matlab实验报告
实验一：Matlab操作环境熟悉一、实验目的1．初步了解Matlab操作环境。
2．学习使用图形函数计算器命令funtool及其环境。
二、实验内容熟悉Matlab操作环境，认识命令窗口、内存工作区窗口、历史命令窗口；学会使用format命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；学会使用who和whos命令查看内存变量信息；学会使用图形函数计算器funtool，并进行下列计算：1．单函数运算操作。
求下列函数的符号导数(1)y=sin(x); (2) y=(1+x)^3*(2-x);求下列函数的符号积分(1)y=cos(x);(2)y=1/(1+x^2);(3)y=1/sqrt(1-x^2);(4)y=(x1)/(x+1)/(x+2)求反函数(1)y=(x-1)/(2*x+3); (2) y=exp(x); (3) y=log(x+sqrt(1+x^2));代数式的化简(1)(x+1)*(x-1)*(x-2)/(x-3)/(x-4);(2)sin(x)^2+cos(x)^2;(3)x+sin(x)+2*x-3*cos(x)+4*x*sin(x);2．函数与参数的运算操作。
从y=x^2通过参数的选择去观察下列函数的图形变化(1)y1=(x+1)^2(2) y2=(x+2)^2(3) y3=2*x^2 (4) y4=x^2+2 (5) y5=x^4 (6)y6=x^2/23．两个函数之间的操作求和(1)sin(x)+cos(x) (2) 1+x+x^2+x^3+x^4+x^5乘积(1)exp(-x)*sin(x) (2) sin(x)*x商(1)sin(x)/cos(x); (2) x/(1+x^2); (3) 1/(x-1)/(x-2);求复合函数(1)y=exp(u) u=sin(x) (2) y=sqrt(u) u=1+exp(x^2)(3) y=sin(u) u=asin(x) (4) y=sinh(u) u=-x实验二：MATLAB基本操作与用法一、实验目的1.掌握用MATLAB命令窗口进行简单数学运算。
西安邮电大学光学实验matlab仿真结果分析与程序
光学实验实验报告课程名称：光学实验*名：***学院：电子工程学院系部：光电子技术系专业：电子科学与技术年级：科技1201学号：********指导教师：**2014年12 月24 日光波在介质中界面上的反射及透射特性一．实验目的：1．掌握反射系数及透射系数的概念；2．掌握反射光与透射光振幅和相位的变化规律； 3．掌握布儒斯特角和全反射临界角的概念。
二．实验原理：1 反射定律和折射定律光由一种介质入射到另一种介质时，在界面上将产生反射和折射。
现假设二介质为均匀、透明、各向同性介质，分界面为无穷大的平面，入射、反射和折射光均为平面光波，其电场表示式为)(0r k t i l l l l e E E ⋅--=ω l =i, r, t式中，脚标i 、r 、t 分别代表入射光、反射光和折射光；r 是界面上任意点的矢径，在图2-1所示的坐标情况下，有r=ix+jy图2-1 平面光波在界面上的反射和折射 图2-2 k i 、k r 、k t 三波矢关系根据电磁场的边界条件，可以得到如下关系)(0)(t i r i tr i =⋅-=⋅-==r k k r k k ωωω 这些关系表明：①入射光、反射光和折射光具有相同的频率；②入射光、反射光和折射光均在入射面内，k i 、k r 和k t 波矢关系如图2-2所示。
进一步可得tt i i r r i i sin sin sin sin θθθθk k k k == 或tt i i r r i i sin sin sin sin θθθθn n n n ==即介质界面上的反射定律和折射定律，它们给出了反射光、折射光的方向。
折射定律又称为斯涅耳(Snell)定律。
2 菲涅耳公式 s 分量和p 分量通常把垂直于入射面振动的分量称做s 分量，把平行于入射面振动的分量称做p 分量。
为讨论方便起见，规定s 分量和p 分量的正方向如图2-3所示。
图2-3 s 分量和p 分量的正方向反射系数和透射系数 假设介质中的电场矢量为)(i 0e r k t l l l E E ⋅--=ω l =i, r, t其s 分量和p 分量表示式为)(i 0e r k t lm lm l E E ⋅--=ω m =s,p则定义s 分量、p 分量的反射系数、透射系数分别为tmtm m im rmm E E t E E r 0000==菲涅耳公式假设界面上的入射光、反射光和折射光同相位，根据电磁场的边界条件及s 分量、p 分量的正方向规定，可得ts rs s E E E i =+和2tp 1rp 1ip cos cos cos θθθH H H =-利用E H εμ=，上式变为22ts 11rs is cos cos )(θθn E n E E =-再利用折射定律，消去E ts ，经整理可得)sin()sin(1212is rs θθθθ+-=E E 根据反射系数定义，得到)sin()sin(2121θθθθ+--=s r221111cos cos cos 2θθθn n n t s +=将所得到的表示式写成一个方程组，就是著名的菲涅耳公式：212122112*********tan tan tan tan cos cos cos cos )sin()sin(θθθθθθθθθθθθ+--=+-=+--==n n n n E E r is rs s 2121211221122121002sin 2sin 2sin 2sin cos cos cos cos )tan()tan(θθθθθθθθθθθθ+-=+-=+-==n n n n E E r iprp p 21121121112100221111212100cos cos cos 2)cos()sin(sin cos 2cos cos cos 2)sin(sin cos 2θθθθθθθθθθθθθθθθn n n E E t n n n E E t iptp p is ts s +=-+==+=+==这些系数首先是由菲涅耳用弹性波理论得到的，所以又叫做菲涅耳系数。
MATLAB实验报告2
实验二用MATLAB实现线性系统的时域分析[实验目的]1．研究线性系统在典型输入信号作用下的暂态响应；2．熟悉线性系统的暂态性能指标；3．研究二阶系统重要参数阻尼比ξ对系统动态性能的影响；4．熟悉在MATLAB下判断系统稳定性的方法；5．熟悉在MATLAB下求取稳态误差的方法。
[实验原理]MATLAB中有两类用于求解系统时域响应的方法。
一、用MATLAB函数（命令）进行暂态响应分析1 求取线性连续系统的单位阶跃响应的函数——step基本格式为：step(sys)、step(num,den)、step(A,B,C,D)、step(sys,t) 、step(sys1,sys2,…,t) 、y=step(sys,t)、[y,t]=step(sys)、[y,t,x]=step(sys)其中模型对象的类型如下：sys = tf(num,den) 多项式模型、sys = zpk(z,p,k) 零点极点模型、sys = ss(a,b,c,d) 状态空间模型参数无t，表示时间向量t的范围自动设定。
参数有t，表示给定时间向量t，应该有初值，时间增量，末值，如t=0:0.01:2。
前5种函数可以绘出阶跃响应曲线；后3种函数不绘阶跃响应曲线，而是返回响应变量y，时间向量t，以及状态变量x。
2 求取线性连续系统的单位脉冲响应的函数——impulse基本格式为：impulse(sys)、impulse(num,den) 、impulse (sys,tf)、impulse (sys,t) 、impulse (sys1,sys2,…,t) 、y=impulse(sys,t) 、[y,t]=impulse(sys) 、[y,t,x]=impulse (sys)3 求取线性连续系统的单位斜坡响应MATLAB没有直接求系统斜坡响应的功能函数。
在求取控制系统的斜坡响应时，通常用阶跃响应函数step()求取传递函数为G (s)/s的系统的阶跃响应，则其结果就是原系统G (s)的斜坡响应。
matlab实验报告
（2）将上述信号左平移或右平移 得到两个新信号y2和y3，利用hold命令在同一图中显示三条曲线并通过线型和颜色加以区分。
（3）利用subplot命令将y，y2和y3以子图的形式显示。
3.
实验结果及分析
1．（1）
t1 =
0.0780
实验结果及分析
（2）
t2 =
0.0780
实验结果及分析
（3）
t3 =
0
实验结果及分析
2.
实验
名称
实验四Matlab在通信系统中的应用
实验目的
1、练习通信原理中各种调制方法及简单通信系统的Matlab编程实现；
2、学习Matlab与Simulink的混合编程。
（2）y变成虚数后，重复上述运算。
（3）创建矩阵 ，并进行如下操作，取出矩阵A的前两行元素，生成一个新矩阵B。生成一个3*3的单位矩阵B，与矩阵A进行加、减、乘、除、相等、不相等运算。
（4）随机生成2个复数x1和x2，计算x3=x1÷x2，x3的实部real_x3、虚部image_x3、模abs_x3、幅角angle_x3，并把x1，x2，x3，real_x3，image_x3，abs_x3，angle_x3保存至文件complex_学号.mat。利用save命令保存数据为mat格式文件后，用clear命令清除工作空间中的变量，再利用load命令从文件读入数据。
2、流程控制语句相关操作。
（1）使用for语句及while语句求解1至100整数中奇数的和。
（2）求连续自然数的和，当和大于等于1000时，求最后一个自然数以及自然数的和。
国家开放大学《Matlab语言及其应用》实验报告(第二章--Matlab基本程序设计)
国家开放大学《Matlab语言及其应用》实验报告
——Matlab基本程序设计
姓名：学号：
1．M文件
M文件是一个文本文件，可以用任何文本编辑器来建立和编辑，通常使用 Matlab 自带的 M 文件编辑器。
例：编写一个脚本文件将华氏温度转化为摄氏温度。
c=5
9
(f−32)
●新建一个 M 文件 f2cs.m，内容如下：
clear; % 清除当前工作空间中的变量
f=input('Please input Fahrenheit temperature:');
c=5*(f-32)/9;
fprintf('The centigrade temperature is %g\n',c);
●在命令窗口中输入 f2cs，即可执行该 M 文件。
脚本文件运行后，文件中所产生的所有变量都驻留在 Matl ab 当前工作空间，可以用 whos 查看。
2.程序控制结构
例：计算函数y=cos(x+1)+√x2+1 x≥10d的值。
clear
x=input('请输入x的值:');
if x>=10
y=cos(x+1)+sqrt(x*x+1);
end
y
请输入x的值:12
y =12.9490
3.程序调试
例：拼写错误
clear
t=0:0.1:10;
y=son(t);
例：函数错误
调试器的使用：MatLab的调试器和M文件编辑器集成一体。
实验二MATLAB程序设计含实验报告
实验二MATLAB 程序设计一、 实验目的1.掌握利用if 语句实现选择结构的方法。
2.掌握利用switch 语句实现多分支选择结构的方法。
3.掌握利用for 语句实现循环结构的方法。
4.掌握利用while 语句实现循环结构的方法。
5.掌握MATLAB 函数的编写及调试方法。
二、 实验的设备及条件计算机一台（带有MATLAB7.0以上的软件环境）。
M 文件的编写：启动MATLAB 后，点击File|New|M-File ，启动MATLAB 的程序编辑及调试器（Editor/Debugger ），编辑以下程序，点击File|Save 保存程序，注意文件名最好用英文字符。
点击Debug|Run 运行程序，在命令窗口查看运行结果，程序如有错误则改正三、 实验内容1.编写求解方程02=++c bx ax 的根的函数（这个方程不一定为一元二次方程，因c b a 、、的不同取值而定），这里应根据c b a 、、的不同取值分别处理，有输入参数提示，当0~,0,0===c b a 时应提示“为恒不等式!”。
并输入几组典型值加以检验。
（提示：提示输入使用input 函数）2.输入一个百分制成绩，要求输出成绩等级A+、A 、B 、C 、D 、E 。
其中100分为A+，90分～99分为A ，80分～89分为B ，70分～79分为C ，60分～69分为D ，60分以下为E 。
要求：（1）用switch 语句实现。
（2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。
（提示：注意单元矩阵的用法）3.数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。
重复此过程，最终得到的结果为1。
如：2?13?10?5?16?8?4?2?16?3?10?5?16?8?4?2?1运行下面的程序，按程序提示输入n=1,2,3,5,7等数来验证这一结论。
请为关键的Matlab 语句填写上相关注释，说明其含义或功能。
西安邮电大学matlab实验二
实验二 Matlab 基本操作与绘图一、实验目的：1、 熟悉并掌握矩阵的基本操作与基本运算。
2、 熟练使用常用的二维和三维绘图函数。
二、实验内容：1、 创建下列矩阵134********,122438,[1,0,8,3,6,2,4,23,46,6]78276853a b c -⎡⎤⎡⎤⎢⎥⎢⎥===-⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦，对上述矩阵做如下运算：1）求a+b ，a*b ，a.*b ，a/b ，a./b ，a^2，a.^2的结果。
2）求c 中所有元素的平均值、最大值、中值；将c 中元素按从小到大顺序排列。
3）求d=b(2:3, [1,3])。
2、 已知矩阵A=[1 2 3; 8 9 10; 11 20 38]，分别求此矩阵的共轭矩阵、转置矩阵、逆矩阵、行列式的值、矩阵的迹、矩阵的秩。
3、 求解下面的方程组。
⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⨯⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----017413231511222315921274321x x x x 4、 使用Subplot 函数，把图形窗口分成两个部分，分别绘制sin2t 和3cos3t 曲线，t 范围: 0-5。
要求前者用红色实线，数据点形状为上三角，有网格线，x 轴加标注“x”，y 轴加标注“y=sin2t”，加图题目“sin2t 的曲线”；后者用蓝色虚线，数据点形状为圆形，无网格线，x 轴加标注“x”，y 轴加标注“y=3cos3t”，加图题目“3cos3t 的曲线”。
5、绘制22()2,3,3x y Z xye x y +-=-≤≤的三维网线图和三维曲面图。
三、实验报告要求：针对每一个实验内容，在实验过程记录里写出相应程序语句，并附上Matlab 执行结果。
四、实验过程记录：1.（1）.>> a=[1 3 4;5 12 44;7 8 27]; >> b=[-7 8 4;12 24 38;68 -5 3];>> d=a+bd =-6 11 817 36 8275 3 30>> d=a*bd =301 60 1303101 108 6081883 113 413 >> d=a.*bd =-7 24 1660 288 1672476 -40 81直接做a/b>> d=a/bd =0.0966 0.0945 0.0080-3.6125 1.5838 -0.5778-1.9917 0.9414 -0.2682另B可逆方阵，计算a/b，如下>> d=a*inv(b)d =0.0966 0.0945 0.0080-3.6125 1.5838 -0.5778-1.9917 0.9414 -0.2682>> d=a./bd =-0.1429 0.3750 1.00000.4167 0.5000 1.15790.1029 -1.6000 9.0000>> d=a^2d =44 71 244373 511 1736236 333 1109 >> d=a.^2d =1 9 1625 144 193649 64 729（2）>> c=[1,0,8,3,6,2,-4,23,46,6];>> mean(c)ans =9.1000>> max(c)ans =46>> median(c)ans =4.5000>> sort(c)ans =-4 0 1 2 3 6 6 8 23 46 （3）>> d=b(2:3, [1,3])d =12 3868 32.>> A=[1 2 3; 8 9 10; 11 20 38]A =1 2 38 9 1011 20 38共轭矩阵>> conj(A)A =1 2 38 9 1011 20 38转置矩阵>> A'ans =1 8 112 9 203 10 38逆矩阵>> inv(A)ans =-2.2540 0.2540 0.11113.0794 -0.0794 -0.2222-0.9683 -0.0317 0.1111行列式的值>> det(A)ans =-63矩阵的迹>> trace(A)ans =48矩阵的秩>> rank(A)ans =3.>> a=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13]; >> b=[4;7;-1;0];>> x=inv(a)*bx =0.49790.14450.0629-0.08134．t=0:pi/100:5;y1=sin(2*t);y2=3*cos(3*t);subplot(2,1,1);plot(t,y1,'r-^');xlabel('x');ylabel('y=sin(2*t)');title('sin2t的图像');grid onsubplot(2,1,2);plot(t,y2,'b:o');xlabel('x');ylabel('y=cos(3*t)');title('cos3t的图像');grid off5.x=linspace(-3,3,25);y=linspace(-3,3,25);[x,y]=meshgrid(x,y);z=x.*y.*exp(-(x.^2+y.^2)/2); mesh(x,y,z)xlabel('x axes');ylabel('y axes');zlabel('z axes');x=linspace(-3,3,25);y=linspace(-3,3,25);[x,y]=meshgrid(x,y);z=x.*y.*exp(-(x.^2+y.^2)/2); surf(x,y,z)五实验小结：通过本次Matlab实验，增强了动手的能力，知道了使用一些简单的函数，明白了.m文件的建立，以及将数据文件保存到本地的磁盘中，以.mat文件为后缀。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
西安邮电学院
《Matlab》
实验报告
（二）
2011- 2012 学年第 1 学期
自动化
专业：
自动0903 班级：
学号：
姓名：
2011 年10月15日
实验二 MATLAB 的基本计算
一、实验目的
1．掌握建立矩阵的方法。
2．掌握MATLAB 各种表达式的书写规则以及常用函数的使用。
3．能用MATLAB 进行基本的数组、矩阵及符号运算。
4．掌握矩阵分析的方法以及能用矩阵求逆法解线性方程组。
二、实验设备及条件
计算机一台（带有MATLAB6.5或以上版本的软件环境）。
三、实验内容
1．利用diag 等函数产生下列矩阵。
⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=032570800a ⎥⎥⎥⎦
⎤
⎢⎢⎢⎣⎡=8040
507
2
b
2．利用reshape 函数将1题中的a 和b 变换成行向量。
3．产生一个均匀分布在（-5，5）之间的随即矩阵（10×2），要求精确到小数点后一位。
4．已知：
⎥⎥⎥⎦⎤⎢⎢⎢
⎣⎡-=765
3877
34434
12A ⎥⎥⎥⎦
⎤
⎢⎢
⎢⎣⎡--=7312
033
2
1
B 求下列表达式的值：
(1) B A K *611+=和I B A K +-=12（其中I 为单位矩阵）
(2) B A K *21=和B A K *.22=
(3) 331^A K =和3.32^A K =
(4) B A K /41=和A B K \42=
(5) ],[51B A K =和]2:);],3,1([[52^B A K =
5．下面是一个线性方程组：
⎥⎥⎥⎦
⎤
⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡52.067.095.03216/15
/14
/15/14/13/14/13/12/1x x x (1) 求方程的解
(2) 将方程右边向量元素3b 改为0.53，再求解，并比较3b 的变化和解的相对变
6．利用randn 函数产生均值为0，方差为1的6×6正态分布随机矩阵C ，然后统计C 中大于-0.3，小于0.3的元素个数t
7．建立一个矩阵A ，找出矩阵中值大于1的元素，并将它们重新排列成行向
量B ，例如⎥⎥⎥⎥⎦
⎤⎢⎢⎢⎢⎣⎡-----=21
55
126308
394382510A 。
8．读以下程序，解释指令sprintf、fprintf、sscanf的作用？
>>rand('state',0);
>>a=rand(2,2);
>>s1=num2str(a)
>>s_s=sprintf('%.10e\n',a)
>>fprintf('%.5g\\',a)
>>s_sscan=sscanf(s_s,'%f',[3,2])
sprintf :某个东西输入到内存
fprintf :保存数据
Sscanf：用它可以从字符串中取出整数、浮点数和字符串
实验报告小结：通过本次Matlab实验，增强了动手的能力，知道了使用一些简单的函数，明白了.m文件的建立，以及将数据文件保存到本地的磁盘中，以.mat 文件为后缀。
对Matlab解释执行的方式有了一定的了解，Matlab的高封装性使得它简单易用，成为流行的软件，在多个领域都有广泛的应用，特别是仿真，科学计算的一些领域，Matlab值得深入的学习，在生活中会给我们带来很大的方便。
同时在学习方法上也有了一定的提高，有一些的东西必须自学，在没有老师帮助的情况下，去积极主动地解决问题，发现一些新方法，对未来人生的规划也增强的明确的认识。