2017年度matlab实验内容

  • 格式:doc
  • 大小:552.37 KB
  • 文档页数:14

2017年文化素质课MATLAB实验实验一、MATLAB基本操与运算基础【实验目的】(1)熟悉MATLAB基本环境,掌握MATLAB变量的使用(2)掌握MATLAB数组的创建(3)掌握MATLAB数组和矩阵的运算【实验内容及步骤】熟悉建立数组的方法:逐个元素输入法、冒号法、特殊方法(使用函数linspace建立)1、有关向量、矩阵或数组的一些运算(1)设A=15;B=20;求C=A+B与c=a+b?(2)设A=[1 2 3;4 5 6;7 8 9],B=[9 8 7;6 5 4;3 2 1];求A*B与A.*B?(3)设a=10,b=20;求i=a/b=?与j=a\b= ?(4)设a=[1 -2 3;4 5 -4;5 -6 7](5)在MATLAB命令行窗口运行A=[1,2;3,4]+i*[5,6;7,8];看结果如何?如果改成运行A=[1,2;3,4]+i[5,6;7,8],结果又如何?(6)请写出完成下列计算的指令:a=[1 2 3;3 4 2;5 2 3],求a^2=?,a.^2=?(7)有一段指令如下,请思考并说明运行结果及其原因clearX=[1 2;8 9;3 6];X( : ) %转化为列向量(8)写出下列指令的运行结果>> A = [ 1 2 3 ]; B = [ 4 5 6 ];>> C = 3.^A>> D = A.^B2、设有矩阵A和B,A=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 2223 24 25],B=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11]1)求它们的乘积C2)将矩阵C的右下角3x2子矩阵赋给D3、完成下列操作1)求[100,999]之间能被61整除的数及其个数(提示:先利用冒号表达式,再利用find和length函数。

)2)建立一个字符串向量,删除其中的大写字母(提示:利用find函数和空矩阵。

)4、求S=20+21+22+23+24+……+210的值(提示:利用求和函数与累乘积函数。

)5、熟悉mean( )、median( )、max( )、min( )、sort( )、round( )等的使用利用MATLAB提供的rand函数和圆整函数随机生成4X3整数矩阵A,进行如下操作1)A各列元素的平均值和中值2)A的最大元素和最小元素及它们的位置3)求A的每行元素的和以及全部元素之和4)分别对A的每列元素按升序、每行元素按降序排列6 、熟悉一些特殊矩阵的建立方法,如单位矩阵、零矩阵、随机矩阵、对角矩阵、全1矩阵、三角矩阵、魔方矩阵等,矩阵元素的处理(1) 高效地创建如下图矩阵A;(2) 将A的第1行乘以2, 第2行乘以4, …,第10行乘以20得到矩阵B;(3) 提取矩阵B中1,3,5,7,9行, 3,4,5,9列的交叉位置处元素, 构成新的矩阵C.(4)熟悉矩阵分解,如Cholesky分解,LU分解、QR分解、Schur分解等已知X=pascal(3);A=[1,2,3;4,5,6;7,8,0],分别求[R,P]=Chol(X);求[L,U]=lu(A),[Q,R]=qr(A)实验二、MATLAB程序设计实验目的熟悉M文件的概念、建立及运行,掌握Matlab的流程控制语句、函数及脚本文件的编程、调试方法。

实验内容及步骤1、分别用if语句和switch语句实现:输入一个百分制成绩,要求输出成绩等级A、B、C、D、E。

其中90分~100分为A,80分~89分为B,70分~79为C,60分~69分为D,60分以下为E。

(注意:对不合理的成绩应输出出错信息“输入的成绩不合理”(若成绩出现小数,则只能是“.5”)2、设计程序,完成两位数的加、减、乘、除四则运算,即随机输入两个整数,再输入一个运算符号,做相应的运算,显示相应的结果,并要求结果显示类似于“a=x+y=34”。

3、用if语句实现,求下列分段函数的值(注意:分别输出x = x1,x2,…,xn,其中xi为任意数值,X的值以向量的形式从键盘输入)。

2 2 26,0 56,1,yχχχχχχχχχχχ⎧+-<≠⎪=-+≤≠≠⎨⎪--⎩且-30<10,2且3其他4、熟悉矩阵的卷积和多项式的乘法,(1)展开多项式(s2+2s+2)(s+4)(s+1)(2)求多项式(x3+2x2+3x+4)(10x2+20x+30)的卷积实验三 MATLAB 程序设计实验内容及步骤1、 分别使用for 循环语句和while 循环语句找出最小的n 值,使得n!>10100,并求出n!。

要求显示的结果为n=70n!=1.197857e+100>1e100 或 70!=1.197857e+100>1e1002、已知⎪⎭⎫ ⎝⎛+⋅=-334sin 234πt e y t 若需要计算t ∈[-1,1],取间隔为0.01,试计算出相对应的y 值。

用步长为0.01的for 循环3、已知 f 1=1,n=1f 2=0,n=2 f 3=1,n=3f n =f n-1-2f n-2+f n-3,n>3求f 1 ~ f 50中:1)最大值和最小值及它们的位置,各数之和(最大值给变量MAX ,其位置给变量x ,最小值给变量MIN ,其位置给变量y ,各数之和给变量SUM )。

2)正数、零、负数的个数(依次赋值给变量positive ,zero ,negative )。

4、Fibonacci 数组的元素满足Fibonacci 规则:a k+2=a k +a k+1,(k=1,2,3,…);且a 1= a 2=1。

请设计一段程序,求出该数组中第一个大于10000的元素 要求显示的结果为: i=21 a(i)=10946 或者a(21)=10946实验四 athlab 文件实验目的1.定制自己的工作环境。

2.编写简单的M 文件。

3.保存内存工作区中的变量到.mat 文件。

4.学会只用Matlab 帮助。

实验内容及步骤1、 编写一个函数文件chengji.m ,用于求两个矩阵的乘积和点乘,然后在命令行窗口中调用该函数(提示:使用varargout 函数)要求:1)当两矩阵是尺寸相同的方阵时,返回它们的乘积和点乘。

2)当两矩阵不是尺寸相同的方阵时,则根据它们的尺寸,分别进行乘积或点乘 3)当两矩阵的尺寸不能满足乘积或点乘的要求时,则提示“这两个矩阵不能乘”2、设向量V=(5,-3,2),求V 的方向余弦。

建立一个函数文件f.m :3、编写一个求一元二次方程ax2 +bx+c=0的根程序。

4、计算分段函数值5、 设x,y,z 的初值为x 0=-0,6,y 0=-1.2,z 0=0.135,求函数f(x,y,z)=x 2+2.5siny-z 2x 2y 2的极值。

6、求积分30π⎰7、求解泊松方程21u -∇= 并与精确解比较,区域为单位圆盘,在圆盘的边界=0cos(1)1010x x y x ⎧+=⎪=⎨⎪≠⎩实验五数据插值与曲线拟合实验实验目的1.掌握函数的编写规则。

2.掌握函数的调用。

3.会用Matlab程序设计实现一些工程算法问题实验内容1.已知原始数据)5/piy==,,插值点为pix-exp(x*).sin(*4:1:0x=,试用四xi*4:3.0:0种不同的插值方法进行插值,并在一副图中绘出四种插值曲线。

2.考虑下面的数据,一位举重教练收集了不同性别的举重动员每个年龄的最大举重磅数,他相信这三者之间存在着函数关系,请Array创建一个年龄数组,以便估算当前队伍中队员(要求13~20岁每一岁都有)的最大举重,并绘制图形。

3.某乡镇企业2006-2012年的生产利润如下4、对矩阵A的列向量、行向量分别进行离散傅立叶变换、并对变换结果进行逆变换。

实验六MATLAB绘图实验目的1.掌握二维图形的绘制。

2.掌握图形的标注3.了解三维曲线和曲面图形的绘制。

实验内容1、二维图形的绘制1)设计程序,在一个图形窗体中画出函数y=e x和y=logx的特性曲线。

其中y=e x 中x取值为(-2:0.1:2),y=logx中x取值为(0.1:0.1:5)。

要求:要有图例、标题('二维图')、坐标轴标签('x轴数据'和'y轴数据');每条曲线的clm自己设置。

2)绘制tan图形,其中x的取值为(0:pi/30:pi/2)。

要求:要有图例、标题、坐标轴标签('x轴数据'和'y轴数据');绘制出的图形应如下右图所示。

3)绘制曲线x=tcos(3t)y=tsin2t 其中-π≤t≤π,步长取π/100。

要求:要有图例、标题、坐标轴标签。

4)编制M 文件,绘制图案:r=2sin 3(θ), r=cos 3 (θ), r=sin 2(θ), r=5cos 3(3.5θ)四种图。

2、三维图形的绘制分别用线框图和曲面图表现函数z=cos(x)sin(y)/y ,其中x 的取值为[-1.5pi,1.5pi],y=x ,要求:要有标题、坐标轴标签3、特殊图形的绘制求解热传导方程21uu t∂-∇=∂ ,区域为单位方形,边界上u=0,生成动态图。

附加实验 MATLAB 自选项目不作全面要求,根据自己的学习能力,完成问题分析、建模、编程。

1、Koch 分形曲线 分形原理这是一类复杂的平面曲线,可用算法描述。

从一条直线段开始,将线段中间三分之一部分用等边三角形的两条边代替,形成具有5个结点的图形(图1);在新的图形中,又将图中每一直线段中间的三分之一部分都用一等边三角形的两条边代替,再次形成新的图形(图2),这时,图形中共有17个结点。

这种迭代继续进行下去可以形成Koch 分形曲线。

在迭代过程中,图形中的点将越来越多,而曲线最终显示细节的多少将取决于迭代次数和显示系统的分辩率。

模型分析考虑由直线段(2个点)产生第一个图形(5个点)的过程。

设1P 和5P 分别为原始直线段的两个端点。

现在需要在直线段的中间依次插入三个点234,,P P P 产生第一次迭代的图形(图1)。

显然,2P 位于1P 点右端直线段的三分之一处, 4P 位于1P 点右端直线段的三分之二处;而3P点的位置可以看成是由4P 点绕2P 旋转60度(逆时针方向)而得到的,故可以处理为向量24P P 经正交变换而得到向量23P P 。

算法如下:(1) 2151()/3P P P P =+-; (2) 41512()/3P P P P =+-; (3)3242()TP P P P A =+-⨯;在(3)中, A 为正交矩阵:图2 第二次迭代图1 第一次迭代cos sin33sin cos33A ππππ⎡⎤-⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦算法根据初始数据(1P 和5P 点的坐标),产生图1中5个结点的坐标。