资源预览内容
第1页 / 共21页
第2页 / 共21页
第3页 / 共21页
第4页 / 共21页
第5页 / 共21页
第6页 / 共21页
第7页 / 共21页
第8页 / 共21页
第9页 / 共21页
第10页 / 共21页
亲,该文档总共21页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述
说明:这只是报告的雏形!说明:这只是报告的雏形! 分享智慧,分享经验,仅分享智慧,分享经验,仅 供参考,欢迎交流,切勿供参考,欢迎交流,切勿 Ctrl+C and Ctrl+V本科实验报告本科实验报告课程名称: MATLAB 及在电子信息工程中应用 实验项目: MATLAB 实验指导书二.三.四 实验地点: 中区图书馆起点机房 专业班级: 信息 1102 班 学号: 学生姓名: 李政洁 指导教师: 李鸿燕 2013 年 5 月 26 日实验二 矩阵和数组的操作一.实验环境MATLAB2012a二.实验目的1.掌握矩阵和数组的一般操作。 2.学习矩阵和数组的加减运算与乘法。 3.掌握对数组中的元素的寻访与赋值,会对数组进行一般的操作。三.实验内容% 1.创建一个 5X5 矩阵,提取主对角线以上部分A=1:5;2 4 6 8 10;3:2:11;linspace(1,6,5);11 ,12, 13, 14,15 A =1.0000 2.0000 3.0000 4.0000 5.00002.0000 4.0000 6.0000 8.0000 10.00003.0000 5.0000 7.0000 9.0000 11.00001.0000 2.2500 3.5000 4.7500 6.000011.0000 12.0000 13.0000 14.0000 15.0000%在矩阵 A 的创建中,太原理工大学信息学院李政洁一次性演示了多种创建方式的混合使用。 。 。 。B=triu(A) B =1.0000 2.0000 3.0000 4.0000 5.00000 4.0000 6.0000 8.0000 10.00000 0 7.0000 9.0000 11.00000 0 0 4.7500 6.00000 0 0 0 15.0000%矩阵 B 为提取了主对角线及以上元素,但 B 不符合题目要求,于是李政洁在下面将使用 triu(A,1)C=triu(A,1) C =0 2 3 4 50 0 6 8 100 0 0 9 110 0 0 0 60 0 0 0 0% OK,至此,李政洁同学完成了第一题,A,C 分别为所求矩阵。% 2.计算三个指定矩阵相乘D=rand(3) D =0.8147 0.9134 0.27850.9058 0.6324 0.54690.1270 0.0975 0.9575 E=magic(3) E =8 1 63 5 74 9 2 F=rand(3,4) F =0.9649 0.9572 0.1419 0.79220.1576 0.4854 0.4218 0.95950.9706 0.8003 0.9157 0.6557 G=D*E*F G =22.7418 23.2308 15.6399 23.548622.9828 23.9760 15.4310 24.78559.6742 12.0883 7.6995 15.1322% OK,至此,李政洁完成了第二题,下面我们快乐的开始第三题吧!% 3.创建一个 3X3 矩阵,并求其转置,逆矩阵H=magic(3) H =8 1 63 5 74 9 2%为了快速创建,李同学使用了函数 magic();I=H I =8 3 41 5 96 7 2% I 为所创建矩阵的转置J=inv(H) J =0.1472 -0.1444 0.0639-0.0611 0.0222 0.1056-0.0194 0.1889 -0.1028% J 为所创建矩阵 H 的逆矩阵% OK,至此,政洁完成了第三题,进入下一题% 4.用两种方法求 Ax=b 的解(A 为四阶随机矩阵,b 为四阶行向量)A1=rand(4) A1 =0.0357 0.7577 0.1712 0.04620.8491 0.7431 0.7060 0.09710.9340 0.3922 0.0318 0.82350.6787 0.6555 0.2769 0.6948 b1=1 2 3 4 b1 =1234%法一x1=inv(A1)*b1 x1 =-2.3797-0.05134.89886.1774%法二x2=A1b1 x2 =-2.3797-0.05134.89886.1774%x1=x2,这也验证了方法的正确性,李政洁同学圆满完成老师这道题要考察%的内容!% 5.创建一个随机矩阵,计算 A 的立方K=rand(4) K =0.3171 0.3816 0.4898 0.75470.9502 0.7655 0.4456 0.27600.0344 0.7952 0.6463 0.67970.4387 0.1869 0.7094 0.6551 L=K3 L =1.7085 2.1869 2.4099 2.44052.3205 2.7714 3.0098 3.05532.0843 2.6034 2.7182 2.77321.7565 2.2160 2.4594 2.4274% ok,这也是一个简单的操作,做下一个吧。 。% 6.求 100-999 之间能被 21 整除的数的个数M=100:999; %数太多了,就不显示了N=find(rem(M,21)=0); O=length(N) O =43%M:100-999 之间整数组成的一维数组或者可以说是一个行向量%N:M 中能被 21 整除的数的下标%0:这些下标的个数,亦为数组中能被 21 整除的数的个数%下面我要做第七题了% 7.设有矩阵 A,B,求 C=AXB,将右下角 3X2 子矩阵赋给DP=1:25; A3=reshape(P,5,5) A3 =1 6 11 16 212 7 12 17 223 8 13 18 234 9 14 19 245 10 15 20 25 A3=A3 A3 =1 2 3 4 56 7 8 9 1011 12 13 14 1516 17 18 19 2021 22 23 24 25%政洁同学使用维数重组和转置快速完成了题目所设定矩阵的创建B3=3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11 B3 =3 0 1617 -6 90 23 -49 7 04 13 11%A3 ,B3 分别题目中的 A,B,我这么做是想是同一报告中不要总是出现 A,BC3=A3*B3 C3 =93 150 77258 335 237423 520 397588 705 557753 890 717 D3=tril(C3,-2) D3 =0 0 00 0 0423 0 0588 705 0753 890 717%至此,李政洁完成了第七题,哈哈-%太原理工电子信息工程 1102 班李政洁-%至此,完成实验二 矩阵和数组的操作 。 。diary四.实验说明我在 MATLAB2012a 中完成实验二 矩阵和数组的操作,在这个过程中,我使用了 diary 命令将我在实验过程中的数据记录,操作记录,心得体会全部保存到指定文档并有了如上的实验。实验三 MATLAB 绘图一.实验环境MATLAB2012a二.实验目的1.掌握 MATLAB 的基本绘图命令 2.掌握运用 MATLAB 绘制一维,二维,三维图形的方法3.绘图加以修饰三.实验内容(一)书本实例1.创建一个 5X5 魔方矩阵,并画出这个矩阵的图形代码:%太原理工大学信息学院电子信息工程 1102 班李政洁%实验三实验内容和步骤A=magic(5); plot(A); title(magic(5) by LZJ)图形:11.522.533.544.550510152025magic(5) by LZJ2. 在同一坐标轴绘制 sin(x),cos(x)两条曲线代码:x=linspace(0,2*pi,50); y=sin(x); plot(x,y) title(sin and cos by LZJ) hold on z=cos(x); plot(x,z) hold off图形:01234567-1-0.8-0.6-0.4-0.200.20.40.60.81sin and cos by LZJ3.画出 y=x2,-51f(i,j)=0.5457*exp(-0.75*y(i)2-3.75*x(j)2-1.5*x(j);else if x(j)+y(i) exp_jiecheng(7) ans =50402.编写一个程序求出阶乘大于或等于 9999 的最小整数。M 文件代码:%太原理工大学信息学院信息 1102 班李政洁%计算使阶乘大于等于 9999 的 n 最小值for n=110:125;fjiecheng=exp_jiecheng(n);if fjiecheng=9999nbreak;end end四.实验说明我在 MATLAB2012a 中完成实验四 M 文件的编写,在这个过程中,我使用了 diary 命令将我在实验过程中的数据记录,操作记录,心得体会全部保存到指定文档并且使用Word2007 记录图形。同时,为了保护知识产权,我在图形的 title 上加了by LZJ ,LZJ 既李政洁的拼音首写字母简写。说明:这只是报告的雏形!说明:这只是报告的雏形! 分享智慧,分享经验,仅分享智慧,分享经验,仅 供参考,欢迎交流,切勿供参考,欢迎交流,切勿 Ctrl+C and Ctrl+V
网站客服QQ:2055934822
金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号