实验二matlab数值计算 (2)

《实验二matlab数值计算 (2)》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验二MATLAB数值计算1.实验目的（1）掌握MATLAB变量的使用（2）掌握MATLAB数组的创建，（3）掌握MATLAB数组和矩阵的运算。（4）熟悉MATLAB多项式的运用2.实验仪器（1）Matlab6.5应用软件安装版一套（3）PC机一台3.实验原理矩阵运算和数组运算在MATLAB中属于两种不同类型的运算，数组的运算是从数组元素出发，针对每个元素进行运算，矩阵的运算是从矩阵的整体出发，依照线性代数的运算规则进行。4.实验步骤（1）使用冒号生成法和定数线性采样法生成一维数组。（2）使用MATLAB提供的库函数reshape，将一维数组转换为二维和

2、三维数组。（3）使用逐个元素输入法生成给定变量，并对变量进行指定的算术运算、关系运算、逻辑运算。（4）使用MATLAB绘制指定函数的曲线图，将所有输入的指令保存为M文件。5.实验报告内容（1）在[0,2*pi]上产生50个等距采样数据的一维数组，用两种不同的指令实现。a=0:2*pi/49:2*pi%b=linspace(0,2*pi,50)a=Columns1through400.12820.25650.3847Columns5through80.51290.64110.76940.8976Columns9through121.02581.15411.

3、28231.4105Columns13through161.53871.66701.79521.9234Columns17through202.05172.17992.30812.4363Columns21through242.56462.69282.82102.9493Columns25through283.07753.20573.33393.4622Columns29through323.59043.71863.84683.9751Columns33through364.10334.23154.35984.4880Columns37through40

4、4.61624.74444.87275.0009Columns41through445.12915.25745.38565.5138Columns45through485.64205.77035.89856.0267Columns49through506.15506.2832>>（2）将一维数组A=1:18，转换为2×9数组和2×3×3数组。>>a=1:18a=Columns1through71234567Columns8through14891011121314Columns15through1815161718>>c=reshape(1:18,2,9

5、)c=Columns1through71357911132468101214Columns8through915171618>>d=reshape(1:18,2,3,3)d(:,:,1)=135246d(:,:,2)=791181012d(:,:,3)=131517141618（3）A=[0234;1350],B=[1053;1505]，计算矩阵A转置与B乘积,数组A、B乘积，计算A&B,A

6、B,~A,A==B,A>B。a=[0234;1350]%M文件b=[1053;1505]d=a'*bc=a.*bm=(a&b)n=(a

7、b)p=(a==b)k=(a

8、>b)a=02341350b=10531505d=150551510218251534402012c=00151211500m=00111100n=11111111p=00001000k=01010010（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

9、=['y=',num2str(y_max)]max_text=char('maximum',t_text,y_text)tit=['y=0.5*exp(',num2str(b),'t)-t*t*sin(t)']holdonplot(t,y,'b')plot(t(t_max),y_max,'r.')text(t(t_max)+0.3,y_max+0.05,max_text)title(tit),xlabel('t'),ylabel('y'),holdoff