实验二MATLAB数值计算 (2).doc

实验二 MATLAB数值计算1.实验目的 （1）掌握MATLAB变量的使用 （2）掌握MATLAB数组的创建，（3）掌握MATLAB数组和矩阵的运算。 （4）熟悉MATLAB多项式的运用2.实验仪器 （1）Matlab6.5应用软件安装版 一套 （3）PC机 一台3. 实验原理 矩阵运算和数组运算在MATLAB中属于两种不同类型的运算，数组的运算是从数组元素出发，针对每个元素进行运算，矩阵的运算是从矩阵的整体出发，依照线性代数的运算规则进行。4. 实验步骤 （1）使用冒号生成法和定数线性采样法生成一维数组。（2）使用MATLAB提供的库函数reshape，将一维数组转换为二维和三维数组。（3）使用逐个元素输入法生成给定变量，并对变量进行指定的算术运算、关系运算、逻辑运算。（4）使用MATLAB绘制指定函数的曲线图，将所有输入的指令保存为M文件。5. 实验报告内容 （1）在0,2*pi上产生50个等距采样数据的一维数组，用两种不同的指令实现。a=0:2*pi/49:2*pi %b=linspace(0,2*pi,50)a = Columns 1 through 4 0 0.1282 0.2565 0.3847 Columns 5 through 8 0.5129 0.6411 0.7694 0.8976 Columns 9 through 12 1.0258 1.1541 1.2823 1.4105 Columns 13 through 16 1.5387 1.6670 1.7952 1.9234 Columns 17 through 20 2.0517 2.1799 2.3081 2.4363 Columns 21 through 24 2.5646 2.6928 2.8210 2.9493 Columns 25 through 28 3.0775 3.2057 3.3339 3.4622 Columns 29 through 32 3.5904 3.7186 3.8468 3.9751 Columns 33 through 36 4.1033 4.2315 4.3598 4.4880 Columns 37 through 40 4.6162 4.7444 4.8727 5.0009 Columns 41 through 44 5.1291 5.2574 5.3856 5.5138 Columns 45 through 48 5.6420 5.7703 5.8985 6.0267 Columns 49 through 50 6.1550 6.2832 （2）将一维数组A=1:18，转换为29数组和233数组。 a=1:18a = Columns 1 through 7 1 2 3 4 5 6 7 Columns 8 through 14 8 9 10 11 12 13 14 Columns 15 through 18 15 16 17 18 c=reshape(1:18,2,9)c = Columns 1 through 7 1 3 5 7 9 11 13 2 4 6 8 10 12 14 Columns 8 through 9 15 17 16 18d=reshape(1:18,2,3,3)d(:,:,1) = 1 3 5 2 4 6d(:,:,2) = 7 9 11 8 10 12d(:,:,3) = 13 15 17 14 16 18 （3）A=0 2 3 4 ;1 3 5 0,B=1 0 5 3;1 5 0 5，计算矩阵A转置与B乘积,数组A、B乘积，计算A&B,A|B,A,A= =B,AB。a=0 2 3 4;1 3 5 0%M文件b=1 0 5 3;1 5 0 5d=a*bc=a.*bm=(a&b)n=(a|b)p=(a=b)k=(ab)a = 0 2 3 4 1 3 5 0b = 1 0 5 3 1 5 0 5d = 1 5 0 5 5 15 10 21 8 25 15 34 4 0 20 12c = 0 0 15 12 1 15 0 0m = 0 0 1 1 1 1 0 0n = 1 1 1 1 1 1 1 1p = 0 0 0 0 1 0 0 0k = 0 1 0 1 0 0 1 0 （4）绘制y= 0.5-t*t*sin(t),t=0,pi并标注峰值和峰值时间,添加标题y= 0.5-t*t*sint，将所有输入的指令保存为M文件a=0.5b=1/3t=0:0.01:piy=a*exp(b*t)-t.*t.*sin(t)y_max,t_max=max(y)t_text=t=,num2str(t(t_max)y_text=y=,num2str(y_max)max_text=char(maximum,t_text,y_text)tit=y=0.5*exp(,num2str(b),t)-t*t*sin(t)hold onplo