Matlab上机练习二答案剖析

1、1、4 8i|(3 2i3 5i7 _6i2 _7i 1 4i9 4i3 _9i7的共轭转置。4 亠 4i J x=4+8i 3+5i 2-7i 1+4i 7-5i;3+2i 7-6i 9+4i 3-9i 4+4i; x ans =4.0000 - 8.0000i3.0000 - 5.0000i2.0000 + 7.0000i3.0000 - 2.0000i7.0000 + 6.0000i9.0000 - 4.0000i1.0000 - 4.0000i3.0000 + 9.0000i7.0000 + 5.0000i4.0000 - 4.0000i2、计算a = ： 7 3与b = 2 6 8的

2、数组乘积 a=6 9 3;2 7 5; b=2 4 1;4 6 8; a.*bans =12363842403、对于AX =B，如果A =492764，B = p63 5 7 一也8,求解X。 A=4 9 2;7 6 4;3 5 7; B=37 26 28； X=ABX =-0.51184.04271.3318,254、a = ,|3 6-4 a=1 2 3;4 5 6; b=8 - 4;3 6 2; abans =010101 a=bans =010101 a a a=b0110ans =000000 a=bans =1111115、a = 50.2 0-8 -0.7 ,在进行逻辑运算时，a

3、相当于什么样的逻辑量相当于 a=1 1 0 1 1。6、角度x = 3o 45 60，求X的正弦、余弦、正切和余切 x=30 45 60; x1=x/180*pi; sin (x1)ans =0.50000.70710.8660 cos(x1)ans =0.86600.70710.5000 tan( x1)ans =0.57741.00001.7321 cot(x1)ans =1.73211.00000.57747、用四舍五入的方法将数组2.4568 6.3982 3.9375 8.5042取整 b=2.4568 6.3982 3.9375 8.5042; roun d(b)ans =2649

4、1j4813 1-5 43 -2 -&设A 36-5-9，B =6-23-8，求 3= A*B ;C2=A *B;C3=A.*B并求上述2_7-12_ 8-13-97所有方阵的逆阵。 A=1 4 8 13;-3 6 -5 -9;2 -7 -12 -8; B=5 4 3 -2;6 -2 3 -8;-1 3 -9 7; 3= A*BC1 =19-823012 273-385429 C2=A*BC2 =-15 16-243663-1793-10522 6117-60194684-10 C3=A.*BC3 =5 1624-26-18 -12-1572-2 -21108-56 inv(C1)ans =0

5、.00620.0400-0.0106-0.00460.01690.00300.01680.02090.0150 inv(C2)Warning: Matrix is close to singular or badly scaled.Results may be inaccurate. RCOND = 8.997019e-019.ans =1.0e+015-0.1997 0.27000.2021 -0.2732-0.0935 0.1264-0.2353 0.3182-0.9553 -0.23910.9667 0.2420-0.4473 -0.1120-1.1259 -0.2818 inv(C3)

6、 ? Error using = inv Matrix must be square.9、设x=rcost+3t,y=rsint+3,分别令r=2,3,4，画出参数t=010区间生成的xy曲线。 t=linspace(0,10); r1=2; x1=(r1*cos(t)+3*t); y1=r1*sin(t)+3; r2=3; x2=(r2*cos(t)+3*t); y2=r2*sin(t)+3; r3=4; x3=(r3*cos(t)+3*t); y3=r3*sin(t)+3; plot(x1,y1,r,x2,y2,b,x3,y3,m)、r54210、设 f(x)=x 5- 4x4 +3x2-

7、 2x+ 6(1) 在x=-2,8之间取100个点，画出曲线，看它有几个过零点。(提示：用polyval函数) x=li nspace(2,8,100); y=polyval(1 0 -4 3 -2 6,x); plot(x,y,b,x,0,y)(2) 用roots函数求此多项式的根。t=1 0 -4 3 -2 6p=roots(t)12、求解多项式x3-7x2+2x+40的根。 r=1 -7 2 40; p=roots(r);-0.21510.44590.79490.270713、 符号函数绘图法绘制函数 x=sin(3t)cos(t)，y=sin(3t)sin(t)的图形，t的变化范围为0

8、,2二。 syms t设 F:OJ ezplot(si n(3*t)*cos(t),si n(3*t)*si n( t),0,pi)14、 x=sym(x); y=(sin(x)八4+(cos(x)八4;diff(y,10)15、xx=sym(x);int(sqrt(9-xA2)A3)/xA6),x)16、2 dx1 x x3x=sym(x);int(x+xA3,x,1,2)门n亠217、求级数的和： 、呼心2sym( n;symsum(n+2/2An,n,1,inf)18、利用函数int计算二重不定积分.(x，y)eydxdyx=sym(x); y=sym(y);z=(x+y)*exp(-x

9、*y);a=int(z,x);int(a,y) 佃、试求出如下极限。1(1) lim (3x - 9x)x ；xipdxy(1) syms x;f=(3Ax+9Ax)A(1/x);l=limit(f,x,inf)ans=9(2) syms x y;f=x*y/(sqrt(x*y+1)-1);limit(limit(f,x,0),y,0)ans =2(3) syms x y;f=(1-cos(xA2+yA2)*exp(xA2+yA2)/(xA2+yA2);limit(limit(f,x,0),y,0)ans =020、已知参数方程丿x incost，试求出和必yL/3= costtsi ntdx

10、 dx syms t;x=log(cos(t);y=cos(t)-t*si n( t);diff(y,t)/diff(x,t) ans =-(-2*si n( t)-t*cos(t)/si n( t)*cos(t) f=diff(y,t,2)/diff(x,t,2);subs(f,t,sym(pi)/3) ans =3/8-1/24*pi*3A(1/2)21、假设 f (x, y)=xy 42edt，试求-2y ex2 syms x y t s=i nt(exp(-tA2),t,0,x*y); x/y*diff(f,x,2)-2*diff(diff(f,x),y)+diff(f,y,2)ans

11、 =2*xA2*yA2*exp(-xA2*yA2)-2*exp(-xA2*yA2)-2*xA3*y*exp(-xA2*yA2)22试求出下面的极限(1)Hm2121二22 -142 -1 62 -1(2n)21111 1(2) lim n( + + + + ) syms k n;symsum(1/(2*kF2-1),k,1,inf)n +兀 n 十2兀n +3兀nans =1/2 limit(symsum(1/(2*k)A2-1),k,1, n),n,inf) ans =1/2(2) limit (n *symsum(1/( nA2+k*pi),k,1, n),n,inf)ans =123、假

12、设一曲线数据点为 x=0:2:4*pi y=sin(x).*exp(-x/5) 试将x的间距调成0.1,并用下列方法进行内插：(1)(2)(3)(4)线性内插法(method= linea)样条内插法(method= spline 二次多项式内插法(method= cubic)多项式拟合法：直接利用6次多项式去通过7个数据点请将这些内插法的结果及原先的数据点画在同一个图上 25、阅读以下程序并在指定位置上写注释(6分) t=(0:0.02:2)*pi;x=si n( t);y=cos(t);z=cos(2*t);plot3(x,y,z,r-,x,y,z,bd) view(-82,58) box

13、 onlege nd(图形,视角)hold on%(1)绘制一条红色连续曲线和监色菱形离散点%(2)以方位角-82度，俯角58度观察视冬%显示坐标轴的矩形框%(3)标注图例%保持图形x=-4:4;y=x; X,Y=meshgrid(x,y);Z=X.a2+Y.a2；surf(X,Y,Z);%(5)生成平面网格坐标矩阵%(6)绘制三维曲面26、 编写一个M函数文件fun_es(x)，计算如下函数：y = 0.5ex/3 - x2 sin(x)，其中参数可以 为标量，也可以为向量。在 Matlab指令窗口键入如下命令检验该函数：fun _es(3)fun cti on fun _es(x)y=0.

14、5.*exp(x./3)-x.A2.*si n(x)27、 编写m函数文件求半径为r的圆的面积和周长。fun cti on f=circle(r)S=pi*rA2L=2*pi*rx,yMatlab上机实验四28、编制m文件，等待键盘输入，输入密码123,密码正确，显示输入密码正确，程序结束; 否则提示，重新输入。fun cti on f=zy2()t=input(Input key(key is 123)if t=123goodelseYou are fool, i nput aga inzy2; /从新输入；end1 1 1 14.3已知y= 2222。当n=100时，求y的值123n步骤：

15、建立函数，命名为 Cbb.m求平方倒数的和。其程序如下function F=Cbb( n)F=0;k=1;while k=nF=F+(1/k)/k;k=k+1;endend在工作空间中求y值y= Cbb(100)y=1.63504.5编制程序，计算1+2+*2000时的最大n值步骤：建立函数文件 Dbb.m,求累加和小于n的最大数，使n=2000程序如下function k=Dbb( n)k=1;f=0;while fnf=f+k;k=k+1;endk=k-2;end将2000作为输入值，求n.在工作空间中输入如下：Dbb(2000) ans =634.6分别用for和while 循环结构编写

16、程序，求出K=2i =0分别用for和while循环结构编写命令文件Ebb,Fbb.并在工作空间中求KEbb程序如下f=1;J(1)=1；for i=1:63J(i+1)=J(i)*2;f=f+J(i+1);End求K程序：Ebb.m1.8447e+019Fbb程序如下1=1;H=1;i=1;while i=63H=H*2;l=l+H;i=i+1;end在工作空间中求KFbb.mI1.8447e+0191. 倘若要是生成一个4行5列的矩阵a=ra nd(4,5)2. 倘若要是对x进行赋值，从5到25，间隔是0.1 x=5:0.1:25;倘若要是对x进行赋值，从2到20，中间间隔100个点 x=

17、li nspace(2,20,100);4.用3种方法将图画在同一个窗口中 x=5:0.1:25;y1=si n(x).*cos(2*x);y2=si n(2*x).*cos(x);plot(x,y1,x,y2) 有一点必须注意：在式子与式子之间，必须要用点乘 要想对所画的图形进行修饰，可以用以下的语句，要点是必须是先画出图形，再写有关对图 形修饰的语句plot(x,y1,x,y2);plot(x,y1,r-,x,y2,b-.);xlabel(横坐标);title(标题);lege nd(si nx*cosx,si n2x*cosx);(” 加示例) 怎样改变线宽？以及怎样给一个中间空着的正方

18、形的边缘及内部填充颜色 plot(x,y,markeredge,g);plot(x,y,markerfacecolor,b);5关于三维图的画法 饼图 pie 火柴杆图 stem(x,y) 阶梯图6. 关于排序问题 a=rand(4,5);随机生成 4*5 的矩阵sort(a);按列排序，且每列由上到下是由小到大sort(a,1);按列排序，且每列由上到下是由小到大sort(a,2);按行排序，且每行由左到右是由小到大倘若想要改变由大到小还是想改变由小到大，需要用到 增长 ascend减小 descend例如：按行进行排序，且每一行自左至右是递减的,即：按照行从高到低排序a=1 2 3;4 5

19、 6;7 8 9sort(a,2,descend)a =123456789ans=3216549877. 关于矩阵问题求矩阵中每一列的平均值mean(a)mean(a,1) 求矩阵中每一行的平均值mean(a,2)从矩阵a中找出最大值max(max(a)内层的 max 相当于找出每一列的最大值，外层的 max 相当于找出每一列最大值的最大 值或者a=rand(3,4) max(a(:) 即 a 的所有行，所有列 倘若对矩阵a前两行的2, 3 ,4列求最大值 max(a(1:2,2:4) 求矩阵a的行列式值|a|b=det(a)det语句是求n*n的行列式的值求矩阵a的逆矩阵b=inv(a)求矩阵a的转置 b=ainv 是矩阵求逆的意思对于矩阵 a:求出矩阵a中每一列中的最大值，最小值a=1 2 3;4 5 6;7 8 9 b=max(a) c=min(a)求出矩阵a中每一行中的最大值，最小值a=1 2 3;4 5 6;7 8 9 b=max(a,2) c=min(a,2)8. 关于取整问题a=5/2 b=ceil(a)向大取c=floor(a)向小取9. 关于多形式的乘法和除法例如：求 x2+2x+3 与 4x2+5x+6 的乘积 a=1 2 3b=4 5 6 c=conv(a,b) 求 x2+2x+3 除