《Matlab 基础与应用》实验任务书
实验一 Matlab 基本操作
一、实验目的
1、认识Matlab 软件环境，熟悉Matlab 编程的特点；
2、学习信号的产生与画图，矩阵生成与运算，流程控制等操作。
二、实验内容
1、编写程序产生如下信号，并按要求绘制相应的曲线。利用saveas 命令把图形分别保存成.fig 格式和.jpg 格式的图像文件。
（1）信号形式sin (2)4y
x π=+，利用plot 命令绘制该信号在02x π≤≤区间内，且采样间隔为1000π
的曲线。填加网格线，标题为“正弦信号曲线”，横坐标显示“x ”，纵坐标显
示“y ”。
（2）将上述信号左平移或右平移2π
得到两个新信号y2和y3，利用hold 命令在同一图
中显示三条曲线并通过线型和颜色加以区分。
（3）利用subplot 命令将上述三个信号利用子图来显示；利用figure 命令将上述三个信号分别在不同的窗体中同时显示。
（4）z = 0:0.1:40；x = cos(z)；y = sin(z)，利用plot3命令绘制三维曲线并填加相应的坐标轴标注。（通信选做）
<思考与练习>
试一试，如何添加图注和为曲线在指定位置添加注释。
2、向量、矩阵生成及相关运算。
（1）x=1:3; y=4:6; x*y=?; x.*y=?; x’*y=?; x*y’=?; x.’*y=?; x*y.’=?
（2）创建矩阵1
234
56789A ??=??????，并进行如下操作，取出矩阵A 的前两行元素，生成一个
新矩阵B 。生成一个单位矩阵，与矩阵A 进行加减运算。练习利用save 命令保存数据为mat 格式，并利用load 命令从外部文件读入数据。
（3）利用wavread 与imread 命令读取一个wav 格式的语音文件和jpg 格式的图像文件，并利用sound 和image/imshow 命令将其播放或显示出来。利用auwrite 把语音文件写成.au 格式的文件，利用imwrite 把图像文件写成bmp 格式的文件。
（4）利用矩阵求逆和乘法运算求解方程组
<思考与练习>
?????=++-=+-=++-4436
322321321321x x x x x x x x x
想一想，进行矩阵的运算时，其维数需满足什么要求。
3、流程控制语句相关操作。
（1）使用for语句及while语句求解1至100整数中奇数的和。
（2）求连续自然数的和，当和大于等于1000时，求最后一个自然数以及自然数的和。
（3）Matlab函数是可以递归调用的，即在函数的内部调用函数自身，试用递归方式编写n!函数。（选做）
<思考与练习>
用Matlab编程建立九九乘法表。
实验二 Matlab在信号与系统中的应用
一、实验目的
1、掌握信号与系统中Matlab编程的特点。
2、学习微分方程、卷积、系统函数的Matlab编程实现。
二、实验内容
1、按照卷积原理步骤求如下两幅图中
1()
f t和2()
f t信号卷积后的波形。
2、给定x(n)=nR16(n)，h(n)=R8(n)，求x和h卷积y，并用函数stem画出x、h和y的图形。
3、用Matlab在时域中编程求解y′(t)+y(t)=f(t), f(t)= exp(-3t)ε(t)的冲激响应、阶跃响应及零状态响应。在simulink仿真环境下，设计系统框图，分析系统的冲激响应、阶跃响应。
实验三 Matlab在数字信号处理中的应用
一、实验目的
1、学习DFT和FFT的原理及Matlab编程实现方法。
2、学习FFT应用于卷积运算和谱分析时的Matlab编程实现方法。
二、实验内容
1、用三种不同的DFT程序计算x(n)=R8(n)的傅里叶变换X(e jw),并比较三种程序计算机运行时间。
（1）用for loop 语句的M函数文件dft1.m，用循环变量逐点计算X(k)；
（2）编写用MATLAB矩阵运算的M函数文件dft2.m, 完成下列矩阵运算；
（3）调用FFT库函数，直接计算X(k)；
（4）分别利用上述三种不同方式编写的DFT程序计算序列x(n)的傅立叶变换X(e jw)，并画出
相应的幅频和相频特性，再比较各个程序的计算机运行时间。
2、研究高密度频谱与高分辨率频谱
设有连续信号
Xa(t)=cos(2π×6.5×10^3t)+cos(2π×7×10^3t)+cos(2π×9×10^3t)
以采样频率fs=32kHz 对信号Xa(t)采样,分析下列三种情况的幅频特性。
（1） 采集数据长度N=16点，做N=16点的DFT ，并画出幅频特性。
（2） 采集数据长度N=16点，补零到256点，做N=256点的DFT ，并画出幅频特性。
（3） 采集数据长度N=256点，做N=256点的DFT ，并画出幅频特性。
<思考与练习>
观察三幅不同频率特性图，分析和比较它们的特点以及形成的原因。
实验四 Matlab 在通信系统中的应用
一、实验目的
1、练习通信原理中各种调制方法及简单通信系统的Matlab 编程实现；
2、学习Matlab 与Simulink 的混合编程。
二、实验内容（1，2选做一个）
1．用Matlab 产生一个频率为2Hz 、功率为1的正弦信源m(t)，设载波频率为10Hz ，A 分别等于0.5、1、2，试画出：
（1）AM 调制信号；
（2）调制信号的功率谱密度；
（3）相干解调后的信号波形。
<思考与练习>
（1）通过得到的曲线，分析调制深度与失真的关系。
（2）尝试用DSB-SC 调制方法重做上题。
（3）绘出加入高斯白噪声后的解调信号。
2．消息信号()cos 2sin m t t t ππ=+，采用频率调制中心频率为100Hz 的载波。频偏常数为10F M K =。
（1）画出该调频信号的波形。
（2）求出该调频信号的振幅谱。
（3）用鉴频器解调该调频信号，并与输入信号比较。
<思考与练习>
确定调制系数，加载白噪声后，重新绘出消息信号和解调信号的频谱及解调器输出的功率谱密度。
3．在Matlab 与Simulink 混合仿真环境下，建立简单通信系统的仿真模型，其中信源是间隔为0.0001s 的随机二进制数，由Matlab 编程封装成Simulink 模块，采用BPSK 调制与解调，经过AWGN 信道。利用Matlab 编程调用此Simulink 文件，并绘制信噪比为0、2、4、6、8、10时的误码率曲线。
<思考与练习>
试一试，观察信号眼图。
要求：
（1）通信与电子专业完成四个实验，其他专业完成实验一和实验二；（2）9月30号之前交实验报告。
（3）第4个实验的后一个小时上机考试。
实验一 Matlab使用方法和程序设计 一、实验目的 1、掌握Matlab软件使用的基本方法； 2、熟悉Matlab的数据表示、基本运算和程序控制语句； 3、熟悉Matlab绘图命令及基本绘图控制； 4、熟悉Matlab程序设计的基本方法。 二、实验内容 1、帮助命令 2、矩阵运算 （1）矩阵的乘法和乘方 已知A=[1 2;3 4]:B=[5 5;7 8]:求A^2*B ( 2 )矩阵除法 已知A=[1 2 3:4 5 6:7 8 9]:B=[1 0 0:0 2 0:0 0 3],求矩阵左除A\B,右除A/B。 ( 3 )矩阵的转置及共轭转置 已知A=[5+i,2-i,1;6*i,4,9-i],求该复数矩阵的转置A',共轭转置A.' ( 4 )使用冒号选出指定元素 已知：A=[1 2 3:4 5 6:7 8 9];求A中第3列前2个元素;A中第2、3行元素。 ( 5 )方括号[] 用magic函数生成一个4阶魔术矩阵，删除该矩阵的第四列 3、多项式 （1）求多项式P(x)=x3-2x-4的根 ( 2 )已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4].,构造多项式，并计算多项式值为20的解。 4、基本绘图命令 ( 1 ) 绘制余弦曲线 ] 2,0[ ), cos(π ∈ =t t y 。 ( 2 ) 在同一坐标系中绘制曲线 ] 2,0[ ), 5.0 sin( ), 25 .0 cos( yπ ∈ - = - =t t y t 5、基本绘图控制 绘制 ] 4,0[π 区间上的y=10sint曲线，并要求： （1）线形为点划线，颜色为红色，数据点标记为加号； （2）坐标轴控制：显示范围，刻度线，比例，网络线； （3）标注控制：坐标轴名称，标题，相应文本。 6、基本程序设计 （1）编写命令文件：计算1+2+....+n<2000时的最大n值； (2)编写函数文件：分别用n和which循环结构编写程序，求2的0到n次幂的和. n=input('输入正数n:') ji=1: for i=1:n; ji=ji+2^i: end ji input('输入正数n:') ji-1:i-1: While i<=n ji=ji+2^i; i=i+1; end ji (3)如果想对一个变量x赋值，当从键盘输入y或Y时，x自动赋为1；当从键盘输入n或N时，x自 动赋为0；输入其他字符时终止程序。 k=input('shuruX:'.'s'): if k=='y' k=='Y' x=1 ; else k=='n' k=='N' x=0; else ruturn end >> n=input('输入正数n:') 输入正数n:20 n =20
实验一 MATLAB 基本操作及运算 一、 实验目的 1、 理解Matlab 数据对象的特点； 2、 掌握基本Matlab 运算规则； 3、 掌握Matlab 帮助的使用方法； 二、 实验的设备及条件 计算机一台（带有MATLAB7.0以上的软件环境）。 三、 实验内容 要求建立一个名为experiment01.m 的，把与实验内容1-7相关的实验命令都放入该文件中，题与题之间用相应注释分割。注意对实验中出现的相关函数或变量，请使用help 或doc 查询相关帮助文档，学习函数的用法。 1、 建立以下标量： 1） a=10 2） b=2.5×1023 3） c=2+3i ，（i 为虚数单位） 4） d=3/2πj e ,（j 为虚数单位，这里要用到exp ，pi ） 2、 建立以下向量： 1） aVec=[3.14 15 9 26] 2） bVec=????? ???????18228871.2 3） cVec=[5 4.8 … -4.8 -5 ] （向量中的数值从5到-5，步长为-0.2） 4） dVec=[100 100.01 … 100.99 101] （产生1到10之间的等对数间隔向量，参考logspace ，注意向量的长度） 3、 建立以下矩阵： 1）???? ??????=2222 aMat aMat 一个9×9的矩阵，其元素全为2；（参考ones 或zeros ）
2）??????? ?????????=1000005000001 bMat bMat 是一个9×9的矩阵，除主对角上的元素为[1 2 3 4 5 4 3 2 1]外，其余元素均为0。（参考diag ）。 3）100 20109212291111 =cMat cMat 为一个10×10的矩阵，可有1:100的向量来产生（参考reshape ） 4）???? ??????=NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN dMat dMat 为3×4的NaN 矩阵，（参考nan ） 5）?? ????---=8710225113eMat 6）产生一个5×3随机整数矩阵fMat ，其值的范围在-3到3之间。（参考rand 和floor 或ceil ） 4、 使用题1中的变量计算下列等式的x,y,z 的值： 1） ) 6/)15((11--+=a e x 2） g g h h b a y /121,)(=+=提示π,参考sqrt 。 3） c c a d c d c R z ))3/sin()]))([(log(π-+= ，其中R 表示取括号内复数的实数部分，c 表示c 的共轭复数，log 是自然对数。（参考real ，conj ，log ） 5、 使用题2中的向量求解一下等式： 1）)25.2/(22 25.221 cVec e xVec -=π， 其中cVec 指的是题2 中定义的向量cVec ，一下雷同。 2）22)(bVec aVec yVec T +=,T aVec 表示aVec 的转置 3） )/1(log 10dVec zVec =，10log 表示已10为底的对数，参考log10 6、 使用题2和题3中所产生的向量和矩阵计算以下等式，注意本题的操作
实验2 MATLAB基础知识 一、实验目的 1.熟悉MATLAB的数据类型 2.熟悉MATLAB的基本矩阵操作 3.熟悉MATLAB的运算符 4.熟悉MATLAB的字符串处理 二、实验内容 1.创建结构体DataTypes,属性包含MATLAB支持的所有数据类型，并通过赋值构 造结构体二维数组。 DataTypes.char=char([65]); DataTypes.string='hello'; DataTypes.int=100; DataTypes.single=1.560 DataTypes = char: 'A' string: 'hello' int: 100 single: 1.5600 >> DataTypes(2).char=char([66]); DataTypes(2).string='kugou'; DataTypes(2).int=200; DataTypes(2).single=3.14 DataTypes = 1x2 struct array with fields: char string int single 2.用满矩阵和稀疏矩阵存储方式分别构造下述矩阵。 A=
答： A = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 B = A = 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 A = 0 0 0 0 1 1 0 0 0 0 1 0 0 0 0 1 A = 0 1 0 0 1 1 0 0 0 0 1 0 0 0 0 1 A = 0 1 0 0 1 0 0 0 0 0 1 0 0 0 0 1
实验一MATLAB 工作环境熟悉及简单命令的执行 一、 实验目的：熟悉MATLAB 的工作环境，学会使用MATLAB 进行一些简单的运算。 二、 实验内容：MATLAB 的启动和退出，熟悉MATLAB 的桌面（Desktop ），包括菜单 （Menu ）、工具条（T oolbar ）、命令窗口(Command Window)、历史命令窗口、工作 空间(Workspace)等；完成一些基本的矩阵操作；学习使用在线帮助系统。 三、实验步骤： 1、启动MATLAB ，熟悉MATLAB 的桌面。 2、在命令窗口执行命令完成以下运算，观察workspace 的变化，记录运算结果。 （1）（365-52?2-70）÷3； （2）area=pi*2.5^2； （3）已知x=3，y=4，在MATLAB 中求z ： ()23 2y x y x z -= ； （4）将下面的矩阵赋值给变量m1,在workspace 中察看m1在内存中占用的字节数。 m1=???? ? ???? ???11514412679810115133216 执行以下命令 >>m1( 2 , 3 ) >>m1( 11 ) >>m1( : , 3 ) >>m1( 2 : 3 , 1 : 3 ) >>m1( 1 ,4 ) + m1( 2 ,3 ) + m1( 3 ,2 ) + m1( 4 ,1) （5）执行命令>>help abs 查看函数abs 的用法及用途，计算abs( 3 + 4i )； （6）执行命令 >>x=0:0.1:6*pi; >>y=5*sin(x); >>plot(x,y) （7）运行MATLAB 的演示程序，>>demo ，以便对MATLAB 有一个总体了解。 四、思考题 1、以下变量名是否合法？为什么？ （1）x2； （2）3col ； （3）_row ； （4）for ；
第1次实验基于MATLAB的图像基本操作 二、实验内容和要求： 1．实现图像Baboon.bmp（MATLAB自带）的读入（可使用imread）和显示（可使用imshow）操作，代码加上足够的注释，需要建立一个M文件实现。 I=imread('F:\标准图像\Baboon.bmp');//读入图像 imshow(I);//显示图像 2．编程实现将一幅RGB图像转换为二值图像，并在一个窗口同时显示处理过程中得到的每一个图像和原图像，同时需要给图像加上标题。（原始数据可以是任意的RGB图像）。需要新建一个M文件实现。 figure,subplot(1,3,1),imshow(I(:,:,1)),title('R'); subplot(1,3,2),imshow(I(:,:,2)),title('G'); subplot(1,3,3),imshow(I(:,:,3)),title('B'); 3.计算图象统计参数： 读取图像（文件名为‘cameraman.tif’）; 最大值 最小值 均值 K=imread('cameraman.tif'); d_max=max(K(:)) d_min=min(K(:)) d_mean=mean(K(:)) 4.利用帮助系统了解im2double,imresize,image函数的作用和语法，并利用这些函数处理已知图像pout.tif（MATLAB自带）并显示处理前后效果。 J=imread('pout.tif'); J1=im2double(J); figure,subplot(1,2,1),imshow(J),title('Before') subplot(1,2,2),imshow(J1),title('After') J2=imresize(J,0.3); figure,subplot(1,2,1),imshow(J),title('Before') subplot(1,2,2),imshow(J2),title('After') figure,subplot(1,2,1),imshow(J),title('Before') subplot(1,2,2),image(J);title('After') 1
实验一 MATLAB基础知识 一、实验目的 1．熟练掌握Matlab的启动与退出； 2. 熟悉Matalb的命令窗口，常用命令和帮助系统； 3. 熟悉Matalb的数据类型； 4. 熟悉Matlab的基本矩阵操作，运算符和字符串处理 二、实验设备 1．方正电脑 2．MATLAB软件 三、实验内容 1. 已知矩阵 A=[3 4 -1 1 -9 10;6 5 0 7 4 -16;1 -4 7 -1 6 -8;2 -4 5 -6 12 -8;-3 6 -7 8 -1 1;8 -4 9 1 3 0] 写出完成下列操作的命令: (1) 将矩阵A的第2-5行中第1，3，5列元素赋给矩阵B； (2) 删除矩阵A的第7号元素； (3) 将矩阵A的每个元素加上30； (4) 求矩阵A的大小和维数； (5) 将矩阵A的右下角3*2矩阵构成矩阵C； (6) 输出[-5,5]范围内的所有元素； 程序： A=[3 4 -1 1 -9 10;6 5 0 7 4 -16;1 -4 7 -1 6 -8;2 -4 5 -6 12 -8;-3 6 -7 8 -1 1;8 -4 9 1 3 0] %(1) B=A(2:5,[1 3 5]) %(2) A2=A; A2(7)=[] %(3) A3=A+30 %(4) length(A) size(A) %(5) C=A(end-2:end,end-1:end) %(6) A6=A.*(A>=-5 & A<=5) 结果: A =
3 4 -1 1 -9 10 6 5 0 7 4 -16 1 -4 7 -1 6 -8 2 -4 5 -6 12 -8 -3 6 -7 8 -1 1 8 -4 9 1 3 0 B = 6 0 4 1 7 6 2 5 12 -3 -7 -1 A2 = Columns 1 through 7 3 6 1 2 -3 8 5 Columns 8 through 14 -4 -4 6 -4 -1 0 7 Columns 15 through 21 5 -7 9 1 7 -1 - 6 Columns 22 through 28 8 1 -9 4 6 12 -1 Columns 29 through 35 3 10 -16 -8 -8 1 0 A3 = 33 34 29 31 21 40 36 35 30 37 34 14
实验一Matlab基本操作 题目： 1．利用基本矩阵产生 3x3 和15x8 的单位阵，全1 阵，全0 阵，均匀分布的随 机阵（[-1，1]之间），正态分布随机阵（方差4，均值1） 2．利用diag()函数和rot90()产生下列矩阵： 然后求解a 阵的逆矩阵aa 及b 阵的特征值和对应特征向量，并利用reshape 将 aa 阵变换成行向量。 3．产生一均匀分布在(-5,5)随机阵(50x2)，精确到小数点后一位。 4．编程实现当α∈[-π，π]，间隔为1o 时，求解正弦和余弦的值，并利用plot() 函数绘制正弦，余弦曲线。 5．利用rand 函数产生(0,1)间均匀分布的10x10 随机矩阵a,然后统计a 中大于等于0.6 的元素个数。 6．利用randn 函数产生均值为0，方差为1 的10x10 正态分布随机阵，然后统计其中大于-0.5，小于0.5 的元素个数。 7．编程实现下表功能： 8．有一矩阵a,找出矩阵中其值大于1 的元素，并将他们重新排列成列向量b。 9．在一保定市区9 月份平均气温变化测量矩阵temp_Baoding_sep 中(48x30)，存在有奇异值（大于42o C，小于0o C）,编程实现删除奇异值所在的行。 10．在给定的100x100 矩阵中，删除整行内容全为0 的行，删除整列内容全为0 的列。 程序： 1. %3X3矩阵 a1=eye(3) a2=ones(3) a3=zeros(3) a4=1-2*rand(3) a5=2*randn(3)+1 %15X8矩阵 b1=eye(15,8) b2=ones(15,8) b3=zeros(15,8) b4=1-2*rand(15,8) b5=2*randn(15,8)+1 运行结果：
MATLAB讲义 第一章 MATLAB系统概述 1.1 MATLAB系统概述 MATLAB（MATrix LABoratory）矩阵实验室的缩写，全部用C语言编写。 特点： （1）以复数矩阵作为基本编程单元，矩阵运算如同其它高级语言中的语言变量操作一样方便，而且矩阵无需定义即可采用。 （2）语句书写简单。 （3）语句功能强大。 （4）有丰富的图形功能。如plot,plot3语句等。 （5）提供了许多面向应用问题求解的工具箱函数。目前，有20多个工具箱函数，如信号处理、图像处理、控制系统、系统识别、最优化、神经网络的模糊系统等。 （6）易扩充。 1.2 MATLAB系统组成 （1）MATLAB语言 MATLAB语言是高级的矩阵、矢量语言，具有控制流向语句、函数、数据结构、输入输出等功能。同时MATLAB又具有面向对象编程特色。MATLAB语言包括运算符和特殊字符、编程语言结构、字符串、文件输入/输出、时间和日期、数据类型和结构等部分。 （2）开发环境 MATLAB开发环境有一系列的工具和功能体，其部分具有图形用户界面，包括MATLAB桌面、命令窗口、命令历史窗口、帮助游览器、工作空间、文件和搜索路径等。 （3）图形处理 图形处理包括二维、三维数据可视化，图像处理、模拟、图形表示等图形命令。还包括低级的图形命令，供用户自由制作、控制图形特性之用。 （4）数学函数库 有求和、正弦、余弦等基本函数到矩阵求逆、求矩阵特征值和特征矢量等。 MATLAB数学函数库可分为基本矩阵和操作、基本数学函数、特殊化数学函数、线性矩阵函数、数学分析和付里叶变换、多项式和二重函数等。 （5）MATLAB应用程序接口（API） MATLAB程序可以和C/C++语言及FORTRAN程序结合起来，可将以前编写的C/C++、FORTRAN语言程序移植到MATLAB中。 1.3 MATLAB的应用围包括： MATLAB的典型应用包括： ●数学计算 ●算法开发 ●建模、仿真和演算 ●数据分析和可视化 ●科学与工程绘图 ●应用开发（包括建立图形用户界面） 以矩阵为基本对象 第二章 Matlab基础
实验一matlab基本操作 一、实验目的 熟悉matlab的安装与启动；熟悉matlab用户界面；熟悉matlab功能、建模元素；熟悉matlab优化建模过程。 二、实验设备与工具 1.计算机 2.matlab软件 三、实验步骤 1. 了解matlab的硬件和软件必备环境； 2. 启动matlab； 3. 学习优化建模过程。 四、实验报告要求 1. 写出matlab系统界面的各个构成；以及系统布局区的组成；以及每一部 分的功能； 2. 优化建模过程应用举例 五、实验内容 （一）、Matlab操作界面 1.命令窗口（command window） 2.命令历史窗口（command history） 3.工作空间管理窗口（workspace） 4.当前路径窗口（current directory） （二）、优化建模过程应用举例 1、简单矩阵 123 456 789 A ?? ?? =?? ?? ?? 的输入步骤。 （1）在键盘上输入下列内容 A = [1,1,3; 4,5,6; 7,8,9] （2）按【Enter】键，指令被执行。 （3）在指令执行后，MATLAB指令窗中将显示以下结果： A = 1 2 3 4 5 6 7 8 9
2、矩阵的分行输入。 A=[1,2,3 4,5,6 7,8,9] A = 1 2 3 4 5 6 7 8 9 3、指令的续行输入 S=1-1/2+1/3-1/4+ ... 1/5-1/6+1/7-1/8 S = 0.6345 4、画出衰减振荡曲线t e y t 3sin 3-=及其它的包络线3 0t e y -=。t 的取值范围是]4,0[π。 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 5、画出2222) sin(y x y x z ++=所表示的三维曲面。y x ,的取值范围是]8,8[-。 clear;x=-8:0.5:8; y=x';
Experiment 1. Fundamental Knowledge of Matlab (II) 【Experimental Purposes】 1、熟悉并掌握MATLAB的工作环境。 2、运行简单命令，实现数组及矩阵的输入输出，了解在MATLAB下如何绘图。【Experimental Principle】 1. Vectors Let's start off by creating something simple, like a vector. Enter each element of the vector (separated by a space) between brackets, and set it equal to a variable. For example, to create the vector a, enter into the MATLAB command window (you can "copy" and "paste" from your browser into MATLAB to make it easy): a = [1 2 3 4 5 6 9 8 7] MATLAB should return: a = 1 2 3 4 5 6 9 8 7 To generate a series that does not use the default of incrementing by 1, specify an additional value with the colon operator (first:step:last). In between the starting and ending value is a step value that tells MATLAB how much to increment (or decrement, if step is negative) between each number it generates. To generate a vector with elements between 0 and 20, incrementing by 2(this method is frequently used to create a time vector), use t = 0:2:20 t = 0 2 4 6 8 10 12 14 16 18 20 Manipulating vectors is almost as easy as creating them. First, suppose you would like to add 2 to each of the elements in vector 'a'. The equation for that looks like: b = a + 2 b = 3 4 5 6 7 8 11 10 9 Now suppose, you would like to add two vectors together. If the two vectors are the same length, it is easy. Simply add the two as shown below: c = a + b c = 4 6 8 10 12 14 20 18 16 Subtraction of vectors of the same length works exactly the same way. MATLAB sometimes stores such a list in a matrix with just one row, and other times in a matrix with just one column. In the first instance, such a 1-row matrix is called a row-vector; in the
实验一 MATLAB 运算基础 1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。 (1) 0 12 2sin 851z e =+ (2) 21ln(2 z x =，其中2 120.45 5i x +??=? ?-?? (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+= ++=--L (4) 2242011 122123t t z t t t t t ?≤
4. 完成下列操作： (1) 求[100,999]之间能被21整除的数的个数。 (2) 建立一个字符串向量，删除其中的大写字母。解：(1) 结果： (2). 建立一个字符串向量例如： ch='ABC123d4e56Fg9';则要求结果是：
实验二 MATLAB 矩阵分析与处理 1. 设有分块矩阵33322322E R A O S ?????? =? ??? ，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证2 2 E R RS A O S +?? =???? 。 解: M 文件如下； 5. 下面是一个线性方程组： 1 231 1 12340.951110.673450.521114 5 6x x x ?? ??????????????=??? ??????????????????? (1) 求方程的解。 (2) 将方程右边向量元素b 3改为0.53再求解，并比较b 3的变化和解的相对变化。 (3) 计算系数矩阵A 的条件数并分析结论。 解： M 文件如下： 123d4e56g9
MATLAB基本操作 实验报告 课程名称： 院系： 专业班级： 学号： 学生姓名： 指导教师： 开课时间：至学年第学期
一、学生撰写要求 按照实验课程培养方案的要求，每门实验课程中的每一个实验项目完成后，每位参加实验的学生均须在实验教师规定的时间内独立完成一份实验报告，不得抄袭，不得缺交。 学生撰写实验报告时应严格按照本实验报告规定的内容和要求填写。字迹工整，文字简练，数据齐全，图表规范，计算正确，分析充分、具体、定量。 二、教师评阅与装订要求 1.实验报告批改要深入细致，批改过程中要发现和纠正学生实验报告中的问题，给出评语和实验报告成绩，签名并注明批改日期。实验报告批改完成后，应采用适当的形式将学生实验报告中存在的问题及时反馈给学生。 2.实验报告成绩用百分制评定，并给出成绩评定的依据或评分标准（附于实验报告成绩登记表后）。对迟交实验报告的学生要酌情扣分，对缺交和抄袭实验报告的学生应及时批评教育，并对该次实验报告的分数以零分处理。对单独设课的实验课程，如学生抄袭或缺交实验报告达该课程全学期实验报告总次数三分之一以上，不得同意其参加本课程的考核。 3.各实验项目的实验报告成绩登记在实验报告成绩登记表中。本学期实验项目全部完成后，给定实验报告综合成绩。 4.实验报告综合成绩应按课程教学大纲规定比例（一般为10-15%）计入实验课总评成绩；实验总评成绩原则上应包括考勤、实验报告、考核（操作、理论）等多方面成绩； 5.实验教师每学期负责对拟存档的学生实验报告按课程、学生收齐并装订，按如下顺序装订成册：实验报告封面、实验报告成绩登记表、实验报告成绩评定依据、实验报告（按教学进度表规定的实验项目顺序排序）。装订时统一靠左侧按“两钉三等分”原则装订。
MATLAB基础及应用实验指导书
MATLAB基础及应用课程实验指导书 实验一软件环境和基本操作熟悉 一、实验目的 ①熟悉MATLAB 主界面，并学会简单的菜单操作。 ②学会简单的矩阵输入与信号输入。 ③掌握部分绘图函数。 二、实验内容 熟悉MATLAB操作环境，通过简单操作熟悉MATLAB的基本使用方法。 三、实验要求与步骤 1、用户第一次使用MATLAB 时，建议首先在屏幕上键入demo 命令，它将启动MATLAB 的演试程序，用户可在此演示程序中领略MATLAB 所提供的强大的运算与绘图功能。也可以键入help 进行进一步了解。 MATLAB 启动界面如图所示： MATLAB 语言最基本的赋值语句结构为：变量名列表= 表达式。表达式由操作符或其它字符，函数和变量名组成，表达式的结果为一个矩阵，显示在屏幕上，同时输送到一个变量中并存放于工作空间中以备调用。如果变量名和“= ”省略，则ans 变量将自动建立，例如键入：1900/81 ，得到输出结果：ans =23.4568 。 MATLAB 中变量命名的原则要求必须以英文字母开头，文件夹名字中可以
还包括、下划线和数字，不要使用其他的字符，更不要单纯使用数字或者中文名命名，有时在运行MATLAB 中一些莫名的错误可能就是不规范的命名引起的。这种规则包括将来为自己编写的脚本文件、函数文件命名以及为使用的变量命名也应遵循这个规则。 2、用户工作目录和当前目录的建立和设置 1）为管理方便，每个用户在使用MATLAB 前，尽量为自己建立一个专门的工作目录，即“用户目录”，用来存放自己创建应用文件。例如首先打开资源管理器，在E 驱动器下可以根据自己喜好建立一个新文件夹，但应注意：该文件夹必须以英文开头，文件夹名字中可以包括字母、下划线和数字，不要使用其他的字符，更不要单纯使用数字或者中文名命名，有时在运行MATLAB 中一些莫名的错误可能就是不规范的命名引起的，这与MATLAB 中为变量和文件名命名原则是一致的。尽管MATLAB\work 允许用户存放用户文件，但最好把它仅作为临时工作目录来使用。 2）为了用户运行方便，在开始工作前可把用户定义的目录设置为当前目录，方法是可直接在MATLAB 命令窗口桌面上方通过目录设置栏来实现。 3）把自己的工作目录加到MATLAB 搜索路径上，这样可以在任何情况下方便地调用自己的编写文件。MATLAB 工作时，基本搜索过程为：首先在工作空间，即MATLAB 内存中进行检查，看输入的指令是不是变量；如不是，则检查输入指令是不是内建函数（比如sin 函数等）；如不是，则在当前目录上，检查是否有与输入指令相同的M 文件存在；如还没有，则在MATLAB 定义的搜索路径其他目录中，检查是否有该M 文件存在。设置方法是通过File 菜单下设置路径对话框进行。 3、熟悉简单的矩阵输入 1）从屏幕上输入矩阵A=[1 2 3 ;4 5 6;7 8 9] 回车 A=[1,2,3;4,5,6;7,8,9] 回车 观察输出结果。 2）试用回车代替分号，观察输出结果。 3）输入size(A) ，观察结果。 4）输入矩阵B=[9,8,7;6,5,4;3,2,1];回车 输入矩阵C=[4,5,6;7,8,9;1,2,3];回车，然后分别键入 A B C 回车观察结果。 5）选择File | new 菜单中的M-file ，输入B=[9,8,7;6,5,4;3,2,1] ；保存为
实验一 MATLAB基础训练 一、实验目的 本次上机实验主要练习使用Matlab的基本操作和基础知识，包括数组（复数、向量、矩阵、结构体数组等）的创建和数组元素的操作和运算、矩阵的运算、Matlab的运算符（尤其是点运算‘.’）、脚本M文件和函数M文件的编写、Matlab文件的编程（基本的流程控制结构）、基本的二维和三维绘图方法以及图形的标注等。希望通过本次实验使大家尽量在短时间内（4学时）掌握Matlab的基本操作和基础知识，为后面的实验项目奠定基础。 二、实验原理 参见PPT中有关内容。 三、实验内容 1. 上机练习课件中的例子。 2. 设两个复数a=1+2i, b=3-4i, 计算a+b, a-b, a*b, a/b, a和b的模。 3. 计算下式的结果，其中x=-3.5°, y=6.7° （提示：①应将角度单位由度转换为数学函数所能处理的弧度值；②求根函数sqrt，取绝对值函数abs,具体用法用help查询） 4. 对矩阵 123 456 789 A ?? ?? =?? ?? ?? , 实现下列操作: (1) 左右翻转(fliplr命令) (2) 上下翻转(flipud命令) (3) 利用cat命令分别将A扩展成3×6和6×3的矩阵 (3) 分别提取A的第2行,第2列,对角线元素 (4) 删除A的第2行2列的元素 （提示：将矩阵元素赋空阵[]可以删除元素，注意此时元素的访问只能使用单下标的方式。观察删除元素后，A中元素的排列方式的变化） 5. 创建[0,2π]区间上拥有100个等间隔元素的列向量 ．．．x．（Matlab默认是行向量），并绘制 y=sin(x1/3)的函数图像。 6. 创建如下图所示的单结构体数组。
实验4 MATLAB 基本编程 一、实验目的 1. 熟悉MATLAB 的脚本编写 2. 熟悉MATLAB 的函数编写 3. 熟悉MATLAB 的变量使用 4. 熟悉MATLAB 的程序控制结构 二、实验内容 1. 分别选用if 或switch 结构实现下述函数表示。 -1 x ≤-a f(x)= x/a -a
************************ MATLAB上机指导书 ************************ 昆明理工大学机电学院 彭用新 2015年3月
实验三符号计算 一、操作部分：在命令窗口执行命令完成以下运算，记录运算结果。 1.findsym：帮助我们获取系统定义的自变量 f= sym('sin(a*x+b*y)'); findsym(f) 2.numden（获取分子分母）, sym2poly,（获取多项式时系数）poly2sym（根据多项式系 数获得符号表达式） [n,d]=numden(sym('x*x+y')+sym('y^2')) p=sym('2*x^3+3*x^2+4'); sym2poly(p) x=[2,3,0,4]; poly2sym(x) 3. collect ：合并同类项；expand：展开多项式；horner: 分解成嵌套形式；factor：因式 分解；simplify: 对表达式化简 syms x y; collect(x^2*y+y*x-x^2-2*x) collect((x+y)*(x^2+y^2+1), y) syms x y; expand((x-2)*(x-4)) syms x;horner(x^3-6*x^2+11*x-6) syms x;factor(x^3-6*x^2+11*x-6) syms x;simplify((x^2+5*x+6)/(x+2)) 4. finverse ：求得符号函数的反函数。 syms x y; finverse(1/tan(x)) f= x^2+y; finverse(f,y) finverse(f) https://www.doczj.com/doc/5f8990274.html,pose 求符号函数的复合函数 syms x y; f = 1/(1 + x^2); g = sin(y); compose(f,g) 6. subs :表达式替换。 syms a b;subs(a+b,a,4)
实验一MATLAB基本操作及运算 一、实验目的 二、实验的设备及条件 三、实验内容 1、建立以下标量： 1） a=3 2） b=5+ 3 j，（j为虚数单位） 3） c=e j 2 / 3 2、建立以下向量： 2.71 38 1） Vb= 28 82 2） Vc=[4 3.8-3.8 -4 ] （向量中的数值从 4 到 -4，步长为 -0.2） 3、建立以下矩阵： 1） 3 L 3 Ma M O M 3 L 3 Ma 为一个 7×7的矩阵，其元素全为 3. 2） 1 11 L91 2 12 O92 Mb M M O M 10 20 L100
Mb 为一个 10× 10的矩阵 . 3） 1 14 5 Mc 2 5 17 3 23 8 4、使用题 1 中的变量计算下列等式的x,y,z的值： 1） x1 1 1 e( (a 15)/6) 2）x2 ( a 15 x1) 3）x3 ln( R [( b c)(b c)]sin( a / 3)) ，其中R表示复数实部。 5、求解函数值 y e ct2/(2.252)，其中 c 取值见题 1， t 的取值范围为题 2 中行 向量 Vc。 6、使用题 1 和题 3 中所产生的标量和矩阵计算等式 Mx a Mc (Mc) 1 (Mc )T 其中 * 为矩阵所对应行列式的值，参考det 。 7、函数的使用和矩阵的访问。 1）计算矩阵 Mb 每一列的和，结果应为行向量形式。 2）计算整个矩阵 Mb 的平均值。 3）用向量 [1 1 1] 替换 Mb 的最上一行的值 4）将矩阵 Mb 的第 2~5 行，第 3 到 9 列的元素所构成的矩阵赋值给矩阵SubMb。 5）删除矩阵 Mb 的第一行； 6）使用函数 rand 产生一个 1× 10的向量 r ，并将 r 中值小于 0.5 的元素设置为 0。 8、已知 CellA（1, 1）＝‘中国’， CellA（ 1， 2）＝‘北京’， CellA（ 2， 1）是一个 3 乘 3 的单位阵， CellA（ 2, 2）＝[1 2 3]，试用 MATLAB创建一个 2 ×2 的细胞数组 CellA。 9、已知结构数组student 中信息包含有姓名，学号，性别，年龄和班级，试用 MATLAB创建相应的结构数组 student。该数组包含有从自己学号开始连续 5 个同学的信息（如果学号在你后面的同学不足 5 个则往前排序），创建完成后查看自己的信息。
实验一 MATLAB 基础知识 一、实验目的 1. 熟练掌握Matlab 的启动与退出； 2. 熟悉Matalb 的命令窗口，常用命令和帮助系统； 3. 熟悉Matalb 的数据类型； 4. 熟悉Matlab 的基本矩阵操作，运算符和字符串处理 、实验设备 1. 方正电脑 2. MATLAB^件 三、实验内容 1. 已知矩阵 A=[3 4 -1 1-9 10;6 5 0 7 4 -16;1 -4 7 -1 6 -8;2 -4 5 -6 12 -8;-3 6 -7 8 -1 1;8 -4 9 1 3 0］ 写出完成下列操作的命令： (1) 将矩阵A 的第2-5行中第1 , 3, 5列元素赋给矩阵 B ; (2) 删除矩阵A 的第7号元素； (3) 将矩阵A 的每个元素加上30 ； (4) 求矩阵A 的大小和维数； (5) 将矩阵A 的右下角3*2矩阵构成矩阵C ； (6) 输出［-5,5］范围内的所有元素； 4. 已知 3 5 4 2 1 2 67 A 34 45 7 ,B 2 8 74 求下列表达式的值： (1) A*B 和 A.*B (2) A A 3 和 A.A 3 (3) A/B 和 B\A (4) [A,B]和[A([1,3]);BA2] 87 90 15 9 3 0 i ab 5. 设 a=1+2i,b=-3+4i, c e 6 ,求 c 的值。 2. 设 x=-74 o ,y=27 o ,求 sin (x 2 y 2) ta n x y 的 值。 3.当a 取，，, ?,,，时，求 e 0.3a sin(a 0.3)在各点的函数 值。
a b 实验二MATLAB程序设计 、实验目的 1 ?学会编写简单函数m文件，调用及调试函数m文件； 2 .学会MATLA啲简单程序编写。 二、实验设备 1. 方正电脑 2. MATLAB^件 三、实验内容 1. 从键盘输入一个4位整数，按如下规则加密后输出。加密规则：每位数字都加上7 然后用和除以10的余数取代该数字；然后将第1位数与第3位数交换，第2位数与第4位数交换。 2. 硅谷公司员工的工资计算方法如下： (1) 工作时数超过120h者，超过部分加发15% (2) 工作时数低于60h者，扣发700元； (3) 其余按每小时84元计发。 试编程按输入的工号和该员工的工时数，计算应发工资。 1 1 1 1 3. 根据 2 2 2 L 2，求的近似值。当n分别取100, 1000, 10000时， 6 122232n2 结果是多少？ (要求：分别用循环结构和向量运算(使用sum函数)来实现) 1 1 1 4 .根据y 1 L ，求 3 5 2n 1 (1) y<3 时的最大n值； (2) 与(1)的n值对应的y值。考虑 a 5.考虑一下迭代公式x n 1，其中，a,b为常数。 b X n (1)编写程序求迭代结果，迭代的终止条件为X n 1 X n 10 5，迭代初值X0 1.0，迭 代次数不超过500次。 (2)如果迭代过程收敛于r，那么r的准确值是 b b 4a，当⑻卩的值取⑴“、 2 (8,3)、(10,时，分别对迭代结果和准确值进行比较。
实验一 MATLAB 基本操作 一、实验目的 1. 学习和掌握MA TLAB 的基本操作方法 2. 掌握命令窗口的使用 3. 熟悉MATLAB 的数据表示、基本运算 二、实验内容和要求 1. 实验内容 1) 练习MATLAB7.0或以上版本 2) 练习矩阵运算与数组运算 2. 实验要求 1) 每位学生独立完成，交实验报告 2) 禁止玩游戏！ 三、实验主要软件平台 装有MATLAB7.0或以上的PC 机一台 四、实验方法、步骤及结果测试 1. 实验方法：上机练习。 2. 实验步骤： 1) 开启PC ，进入MA TLAB 。 2) 使用帮助命令，查找sqrt 函数的使用方法 答: help sqrt 3) 矩阵、数组运算 a) 已知 ??????????=987654321A ，???? ??????=963852741B ，求)2()(A B B A -?+ 答： A=[1, 2, 3; 4, 5, 6; 7, 8, 9]; B=[1, 4, 7; 2, 5, 8; 3, 6, 9]; (A+B)*(2*B-A) b) 已知?? ????-=33.1x ，??????=π24y ，求T xy ，y x T c) 已知??????????=987654321A ，???? ??????=300020001B ，求A/B, A\B. d) 已知???? ??????=987654321A ，求：(1) A 中第三列前两个元素；(2) A 中所有第二行元素；(3) A 中四个角上的元素；(4) 交换A 的第1、3列。(5) 交换A 的第1、2行。(6) 删除A 的第3列。
e) 已知[]321=x ，[]654=y ，求：y x *.，y x /.，y x \.，y x .^， 2.^x ，x .^2。 f) 给出x=1,2,…,7时，x x sin 的值。 3）常用的数学函数 a ）随机产生一个3x3的矩阵A ，求：(1) A 每一行的最大、最小值，以及最大、最小值所在的列；(2) A 每一列的最大、最小值，以及最大、最小值所在的行；(3) 整个矩阵的最大、最小值；(4) 每行元素之和；(5) 每列元素之和；(6) 每行元素之积；(7) 每列元素之积。 b) 随机产生两个10个元素的向量x ，y 。（1） 求x 的平均值、标准方差。（2） 求x ，y 的相关系数。（3）对x 排序，并记录排序后元素在原向量中的位置。 4) 字符串操作函数 建立一个字符串向量（如‘ABc123d4e56Fg9’），然后对该向量做如下处理： (1) 取第1～5个字符组成的子字符串。 (2) 将字符串倒过来重新排列。 (3) 将字符串中的小写字母变成相应的大写字母，其余字符不变。 (4) 统计字符串中小写字母的个数。