MTALAB实验一 Matlab使用方法和程序设计.doc

实验一 Matlab使用方法和程序设计一、 实验目的1、掌握Matlab软件使用的基本方法；2、熟悉Matlab的数据表示、基本运算和程序控制语句3、熟悉Matlab绘图命令及基本绘图控制4、熟悉Matlab程序设计的基本方法二、 实验内容：1、帮助命令使用help命令，查找 sqrt（开方）函数的使用方法； help sqrt SQRT Square root. SQRT(X) is the square root of the elements of X. Complex results are produced if X is not positive. See also sqrtm, realsqrt, hypot. Overloaded methods: darray/sqrt sym/sqrt Reference page in Help browser doc sqrt2、矩阵运算（1） 矩阵的乘法已知A=1 2;3 4; B=5 5;7 8;求A2*B A=1 2;3 4; B=5 5;7 8; S=A2*BS = 105 115 229 251（2） 矩阵除法已知 A=1 2 3;4 5 6;7 8 9;B=1 0 0;0 2 0;0 0 3;求AB,A/B A=1 2 3;4 5 6;7 8 9; B=1 0 0;0 2 0;0 0 3; s1=AB s2=A/Bs1 = 1.0e+016 * -0.4504 1.8014 -1.3511 0.9007 -3.6029 2.7022 -0.4504 1.8014 -1.3511s2 = 1.0000 1.0000 1.0000 4.0000 2.5000 2.0000 7.0000 4.0000 3.0000（3） 矩阵的转置及共轭转置已知A=5+i,2-i,1;6*i,4,9-i;求A., A A=5+i,2-i,1;6*i,4,9-i; s1=A s2=A.s1 = 5.0000 - 1.0000i 0 - 6.0000i 2.0000 + 1.0000i 4.0000 1.0000 9.0000 + 1.0000is2 = 5.0000 + 1.0000i 0 + 6.0000i 2.0000 - 1.0000i 4.0000 1.0000 9.0000 - 1.0000i（4） 使用冒号表达式选出指定元素已知： A=1 2 3;4 5 6;7 8 9;求A中第3列前2个元素；A中所有列第2，3行的元素；方括号用magic函数生成一个4阶魔术矩阵，删除该矩阵的第四列 A=1 2 3;4 5 6;7 8 9; B1=A(1,2,3) B2=A(2,3,:) C=magic(4); C(:,4)=B1 = 3 6B2 = 4 5 6 7 8 9C = 16 2 3 5 11 10 9 7 6 4 14 153、多项式（1）求多项式 的根 P=1 0 -2 -4; roots(P)ans = 2.0000 -1.0000 + 1.0000i -1.0000 - 1.0000i（2）已知A=1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4 ，求矩阵A的特征多项式;把矩阵A作为未知数代入到多项式中； A=1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4; p=poly(A) c=polyvalm(p,A)p = 1.0000 -6.9000 -77.2600 -86.1300 604.5500c = 1.0e-011 * 0.1364 -0.3350 0.0080 0.1421 -0.0455 -0.0114 -0.1627 -0.0909 -0.3084 0.1705 -0.0227 -0.2217 0.0171 -0.1265 -0.0327 0.02274、基本绘图命令（1）绘制余弦曲线 y=cos(t)，t0，2 t=0:0.05:2*pi; y=cos(t); plot(t,y)（2）在同一坐标系中绘制余弦曲线y=cos(t-0.25)和正弦曲线y=sin(t-0.5)，t0，2 t=0:0.05:2*pi;y1=cos(t-0.25);y2=sin(t-0.5);plot(t,y1),hold onplot(t,y2)5、基本绘图控制绘制0，4区间上的x1=10sint曲线，并要求：（1）线形为点划线、颜色为红色、数据点标记为加号；（2）坐标轴控制：显示范围、刻度线、比例、网络线（3）标注控制：坐标轴名称、标题、相应文本； t=0:0.1:4*pi; x1=10*sin(t); plot(t,x1,r-.x) grid on xlabel(t); ylabel(x1); title(x1=10sin(t)函数图象)6、基本程序设计（1）编写命令文件：计算1+2+n s=0;m=0; while (s2000),m=m+1; s=s+m; end s=s-m;m=m-1; s,mans = 1953 62（2）编写函数文件：分别用for和while循环结构编写程序，求2的0到n次幂的和。For循环结构的函数：function s=shiyan162a(k)s=0;for i=1:k, s=s+2i;endWh