matlab所有实验及答案李楠资料

1、第二章习题52401、矩阵丫=151，给出元素1的全下标和451002321473单下标，求出元素 100的存储位置。取出子矩阵构成的矩阵b，取出前两行构成矩阵c,转置矩阵b构成矩阵d，计算a*b、cd,c&d, c|d ,c|d解：命令为：a=-1,0,-6,8;-9,4,0,12.3;0,0,5.1,-2;0,-23,0,-7b=a(:,1 2)c=a(1 2,:)d=be=a*bf=cdg=c&dh=c|di=c|d5、使用函数，实现 A到B、C、D、E的转换1011121 4710789A=25811B=4563 69121233 21107 41654C=,D=118 529 871

2、29 6312 11 101 0，并求该矩阵的维数。3 21解：命令为：Y=5,2,4;0,15,1;45,100,23;21,47,3丫(2,3)Y(10)fin d(Y=100)sub2i nd(4 3,3,2)B=Y(2:2:4,3:-2:1) 或 B=Y(2 4,3 1)m n=size(Y)2、已知矩阵 A=1 0 -1； 241 ；-205，B=0 -10； 213； 1 1 2求 2A+B、A2-3B、A*B、B*A、A .*B，A/B、AB解：命令为：A=1 0 -1； 241；-205B=0 -1 0 ； 213 ； 1 1 2E=2*A+BF=AA2-3*BG=A*BH=B

3、*AI=A.*BJ=A/BK=AB3、利用函数产生3*4阶单位矩阵和全部兀素都为8的4*4阶矩阵，并计算两者的乘积。解：命令为：A=eye(3,4)B=8*o nes(4)C=A*B106894012.34、创建矩阵 a=，取出其前两列005.120230714 10E= 2 5113 6 12解：命令为：A=1 4 7 10;2 5 8 11;3 6 9 12B=rot90(A)C=rot90(A,3)D=fliplr(A)A(:,3)= E=A2000454012 37.矩阵A=285.1，用函数取出列向量20911004200045454008 和矩阵285.101000911004解：

4、命令为：A=2 0 0 0;45 4 0 12.3;2 8 5.2 -2;0 91 100 4B=diag(A)C=tril(A)8建立5阶魔方矩阵，求该矩阵的行列式和逆矩阵、秩以及any和all运算结果。解：命令文件为：A=magic(5)x=det(A)B=i nv(A)C=ra nk(A)D=a ny(A)E=all(A)实验三习题 abc1.创建2*2元胞数组A=100124，并求出1001 ；891 2i创建结构数组 stude nt,要求包含三个属性名：姓名，学 号，年龄，分别存放张三，10021、20李四、10022、18王五、10023、19,并写出取得李四学号的操作。 解：命

5、令为：A=abc,1001;24;89,1+2ix=A1,2(1,1)stude nt=struct( name,张三;李四;王五 , number,10021;10022;10023,age,20;18;19)stude nt(2,1). nu mber2.求多项式 p(x) 5x5 3x3 2x210,在 x 12.3时的值，和p(x)=0时候的根。解：命令为：p=5,0,-3,2,0,10;poly2sym(p)x=polyval(p,12.3)y=roots(p)11026894012.3，+3.已知矩阵A=，求A的特征值,4505202327特征向量。解：命令为：A=-1 10 26

6、 8;-9 4 0 12.3;45 0 5 -2;0 -2 32 -7X,D=eig(A)t=eig(A)4. 已知某函数如下表所示x21222324y0.355370.374610.390730.40674估算x=2 1.5时的近似值，用三次多项式拟合正弦曲线， 并求拟合多项式系数向量。解：命令为：x=21 22 23 24format longy=0.35537 0.37461 0.39073 0.40674in terp1(x,y,21.5, nearest)P,S=polyfit(x,si n(x),3)(P是一个长度为N+1的向量，代表N次多项式，S 是采样点的误差向量。)5. 将1

7、0个学生5门功课的成绩存入矩阵P中，进行(1) 求每门课的最高分，最低分，平均分，(2) 求5门课总分的最高分、最低分解：命令为：Y=90 89 91 82 75;88 82 85 89 95;79 89 93 82 76;91 89 8779 85;90 79 78 91 82;86 85 84 89 88;89 91 79 90 94;75 8595 90 87;85 78 75 90 89;98 89 87 86 80 a=max(Y)b=mi n(Y)c=mea n( Y,1)format shortf=mea n(Y ,2)*5t=max(f)h=mi n(f)或 n=sum(Y,2

8、)a=max (n)b=mi n(n)6. 用数值方法解线性方程组6x 5y 2z 5u 49x y 4z u 133x 4y 2z 2u 13x 9y 2u 11解：命令为：A=6 5 -2 5;9 -1 4 -1;3 4 2 -2;3 -9 2 0b=-4 13 1 11X=Ab7. 求方程x5x 2x 10,在x=1附近的根。解：命令为：fsolve(x*5Ax+2*x-10,1)实验四：程序设计注意：要用M文件实现，一道题，一个 M文件，命令文 件或者函数文件自选201. 求 n!n 1解：命令文件为sum=O;s=1;for n=1:20s=n *s;sum=sum+s;endsum

9、2. 求Sn a aa aaa aa a得值，其中 a是一个数字，由键盘输入，表达式中位数最多项 a的个 数，也由键盘输入。例如：输入a=3,n=4，则 s4=3+33+333+3333解：命令文件为sum=0;s=0;a=input(请输入a的值:)n=input(请输入n的值:)for i=1: ns=s+a*10A(i-1);sum=sum+s;endsumfor y=0:36-xif (x+y=36)&( 2*x+4*y=100) disp(鸡有:)disp(x)disp(兔子有)disp(y)endendend5. 产生20个两位随机整数，输出其中小于平均值的奇 数。解：命令文件为d

10、isp(随机矩阵为:)A=ra ndi nt(1,20,10,99)Aver=mea n(A) disp(小于平均值的奇数为：) for k=Aif (kaver)&rem(k,2)=1 disp(k)endend3. 计算下列分段函数的值 ：0 x 5f(x) X 3 5 x 1042x x 10解：命令文件为x=input(输入 x 的值：)if x5 &x=10f=(x-3)/4.0;elsef=2*x;endf4. 有一群鸡和兔子，加在一起头的数量是36，脚的数量是100,编程序解答鸡和兔子数量各是多少？解：命令文件为for x=0:36实验五：符号计算1. ( 1)化简表达式 f c

11、os2x 2sin2x(2)对表达式x4 x26进行因式分解。解：命令文件为syms xf仁 cos(2*x)+(si n(x)A2f2=xA4-xA2-6t1=simple(f1)t2=factor(f2)2求下列极限值(1) lim sin2x(2) lim( x2 xx2 x)x 0 sin 5xx解：命令文件为syms xf1=si n(2*x)/si n(5*x)f2=sqrt(xA2+x)-sqrt(xA2-x)t1=limit(f1,x,0)t2=limit(f2,x, inf)3.分别计算表达式4x-dx1 x2 2 0(3x2x 1)dx解：命令文件为：syms x f1=x

12、A4/(1+xA2) f2=3*xA2-x+1t1=i nt(f1,x)t2=i nt(f2,x,0,2)在x=1.723处的，精度为 20解：命令文件为:syms x n ms1= n*xA n t1=symsum(s1, n,1, n) t2=symsum(s2,n , 1,100)326. 求表达式2x 1位的代数值解：命令文件为：syms xf=(3*xA3+xA2-1)/(xA2+1) t1=subs(f,x,1.723) t2=vpa(t1,20)7. 求方程组4. 分别求xcosx得一阶导数、二阶导数和三阶导数 解：命令文件为：syms xf=x*cos(x) t1=diff(f

13、,x) t2=diff(f,x,2) t3=diff(f,x,3)5. 求下列级数之和(提示：禾U用symsum(s,v,n,m)函数或 者编程实现)(1) s1x2x2 3x3 L nxn(2) s2 1 4 9 16 L 10000x 6y z 0 x2 2xz z y 的解2 2x y z z解：命令文件为：syms x y zs1=x+6*y+zs2=xA2-2*x*z+z-ys3=xA2+y-z-zA2x,y,z=solve(s1,s2,s3)实验六绘图21.绘制行向量(1 8 6 5 4) 解：figure;A=1 8 6 5 4;plot(A,_rp)的图形。解：figure;f

14、plot(2*exp(-0.5*x).*log10(2*pi*x),1 4 -1 10,40,-.m*) grid on;2. 在区间1x4,绘制函数2e 0.5x lg(2 x)的图形，显示网格，并且显示纵轴的范围为-1 , 103. 在区间1x4，同一窗口的两个子窗口内分别绘制函数x2 1 x和sin(5x23)曲线，第一个子窗口的采用红色点线，并且要求坐标轴采用等刻度,第二个子窗口绿色点化线，并且要求坐标轴为正方形。解：figure;subplot(121);axis equal;fplot(sqrt(xA2-1)+x,1 4,r:) subplot(122);fplot(si n(5*

15、xA2+3),1 4,-.g) axis square;4. 极坐标绘制0 , 2内r=cost+sint的图形，要求用点线，颜色会红色解：figure;t=0:pi/20:2*pi;r=sin (t)+cos(t);polar(t,r,r:)10685.绘制94012 3的垂直条形图和水平条形图，条形宽为1.5。005.1202307解：B=-1 0 -6 8;-9 4 0 12.3;0 0 5.1 -2;0 -23 0 -7; subplot(121);bar(B,1.5)title Group1 subplot(122)barh(B,1.5)title Group2 6. 绘制矩阵A=1

16、 2 3 ; 0 2 7 的饼图，并让 3和2分离出来。并给图形加上标题和x轴标签，y轴标签。解：figure;A=1 2 3;0 2 7;explode=0 1 1;0 1 0;pie(A,explode)title 饼图axis onxlabel( x轴)ylabel( y 轴)Figure 1匚回反三维曲面帝等高銭的三維曲面7. x2,2 ,y 2,2,自己设定步长，产生三维网格X , Y，并在一个图形窗口的四2 2个子窗口分别绘制 z x e x y的三维网格图，带底座的三维网格图，三维曲面图 和带等高线的三维曲面图。解：figure;X,Y=meshgrid(-2:0.1:2);Z=X.*exp(-X.A2-Y .A2);subplot(2,2,1)mesh(X,Y,Z) %mesh绘制三维网格图，meshz绘制带底座的三维网格图title三维网格图subplot(2,2,2)meshz(X,Y,Z)title带底座的三维网格图subplot(2,2,3) surf(X,Y ,Z) title 三维曲面subplot(2,2,4)surfc(X,Y,Z)title带等高