资源预览内容
第1页 / 共11页
第2页 / 共11页
第3页 / 共11页
第4页 / 共11页
第5页 / 共11页
第6页 / 共11页
第7页 / 共11页
第8页 / 共11页
第9页 / 共11页
第10页 / 共11页
亲,该文档总共11页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述
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=ab= ?(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+i5,6;7,8,结果又如何?(6) 请写出完成下列计算的指令:a=1 2 3;3 4 2;5 2 3,求a2=?,a.2=?(7) 有一段指令如下,请思考并说明运行结果及其原因clear X=1 2;8 9;3 6; X( : ) %转化为列向量(8) 写出下列指令的运行结果 A = 1 2 3 ; B = 4 5 6 ; C = 3.A D = A.B 2、设有矩阵A和B,A=1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25,B=3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 111) 求它们的乘积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函数和圆整函数随机生成43整数矩阵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的值以向量的形式从键盘输入)。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.e+1001e100 或 70!=1.e+1001e100 2、已知若需要计算t-1,1,取间隔为0.01,试计算出相对应的y值。用步长为0.01的for循环 3、已知 f1=1,n=1f2=0,n=2f3=1,n=3fn=fn-1-2fn-2+fn-3,n3求f1 f50中:1)最大值和最小值及它们的位置,各数之和(最大值给变量MAX,其位置给变量x,最小值给变量MIN,其位置给变量y,各数之和给变量SUM)。2)正数、零、负数的个数(依次赋值给变量positive,zero,negative)。4、Fibonacci数组的元素满足Fibonacci规则:ak+2=ak+ak+1,(k=1,2,3,);且a1= a2=1。请设计一段程序,求出该数组中第一个大于10000的元素要求显示的结果为:i=21a(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的初值为x0=-0,6,y0=-1.2,z0=0.135,求函数f(x,y,z)=x2+2.5siny-z2x2y2的极值。 6、 求积分 7、求解泊松方程 并与精确解比较,区域为单位圆盘,在圆盘的边界=0实验五 数据插值与曲线拟合实验实验目的1掌握函数的编写规则。2掌握函数的调用。3会用Matlab程序设计实现一些工程算法问题实验内容1. 已知原始数据,插值点为,试用四种不同的插值方法进行插值,并在一副图中绘出四种插值曲线。2. 考虑下面的数据,一位举重教练收集了不同性别的举重动员每个年龄的最大举重磅数,他相信这三者之间存在着函数关系,请创建一个年龄数组,以便估算当前队伍中队员(要求1320岁每一岁都有)的最大举重,并绘制图形。重举大最1315171920男270290330350380女2502703103303603.某乡镇企业2006-2012年的生产利润如下年份2006200720082009201020112012利润(万元)70122144152174196202要求用多项式拟合法预测出该企业2013年和2014年的利润,并绘出拟合图。4、对矩阵A的列向量、行向量分别进行离散傅立叶变换、并对变换结果进行逆变换。实验六MATLAB绘图实验目的1掌握二维图形的绘制。2掌握图形的标注3了解三维曲线和曲面图形的绘制。实验内容1、 二维图形的绘制1)设计程序,在一个图形窗体中画出函数y=ex和y=logx的特性曲线。其中y=ex中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=2sin3(), r=cos3 (), r=sin2(), r=5cos3(3.5)四种图。2、三维图形的绘制分别用线框图和曲面图表现函数z=cos(x)sin(y)/y,其中x的取值为-1.5pi,1.5pi,y=x,要求:要有标题、坐标轴标签3、特殊图形的绘制求解热传导方程 ,区域为单位方形,边界上u=0,生成动态图。附加实验 MATLAB 自选项目不作全面要求,根据自己的学习能力,完成问题分析、建模、编程。1、Koch分形曲线分形原理这是一类复杂的平面曲线,可用算法描述。从一条直线段开始,将线段中间三分之一部分用等边三角形的两条边代替,形成具有5个结点的图形(图1);在新的图形中,又将图中每一直线段中间的三分之一部分都用一等边三角形的两条边代替,再次形成新的图形(图2),这时,图形中共有17个结点。图2第二次迭代图1第一次迭代这种迭代继续进行下去可以形成Koch分形曲线。在迭代过程中,图形中的点将越来越多,而曲线最终显示细节的多少将取决于迭代次数和显示系统的分辩率。模型分析考虑由直线段(2个点)产生第一个图形(5个点)的过程。设和分别为原始直线段的两个端点。现在需要在直线段的中间依次插入三个点 产生第一次迭代的图形(图1)。显然,位于点右端直线段的三分之一处, 位于点右端直线段的三分之二处;而点的位置可以看成是由点绕旋转60度(逆时针方向)而得到的,故可以处理为向量经正交变换而得到向量。算法如下:(1) ;(2) ;(3) ;在(3)中, A为正交矩阵: 算法根据初始数据(和点的坐标),产生图1中5个结点的坐标。结点的坐标数组形成一个52矩阵,矩阵的第一行为的坐标,第二行为的坐标,第五行为的坐标。矩阵的第一列元素分别为5个结点的X坐标,第二列元素分别为5个结点的Y坐标。进一步考虑Koch曲线形成过程中结点数目的变化规律。设第k次迭代产生结点数为,第k+1次迭代产生结点数为,则和之间的递推关系式为。要求:根据题目分析,编程实现算法2、双层玻璃的功效北方城镇的有些建筑物的窗户是双层的,即窗户上装两层厚度为的玻璃夹着一层厚度为的空气,如左图所示,据说这样做是为了保暖,即减少室内向室外的热量流失。我们要建立一个模型来描述热量通过窗户的热传导(即流失)过程,并将双层玻璃窗与用同样多材料做成的单层玻璃窗(如右图,玻璃厚度为
网站客服QQ:2055934822
金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号