MATLAB实习一答案.doc

《MATLAB实习一答案.doc》由会员分享，可在线阅读，更多相关《MATLAB实习一答案.doc（9页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、如有侵权，请联系网站删除，仅供学习与交流MATLAB实习一答案【精品文档】第 9 页MATLAB实习一答案实习一 数值计算和符号计算1、“左除”与“右除”有什么区别？ 在通常情况下，左除x=ab是a*x=b的解，右除x=b/a是x*a=b的解，一般情况下，abb/a。 2、，在进行逻辑运算时，a相当于什么样的逻辑量。 =相当于3、在sin(x)运算中，x是角度还是弧度？ 在sin(x)运算中，x是弧度，MATLAB规定所有的三角函数运算都是按弧度进行运算。 4、角度，求x的正弦、余弦、正切和余切 x=30 45 60; x1=x/180*pi; sin(x1) ans = 0.5000 0.7

2、071 0.8660 cos(x1) ans = 0.8660 0.7071 0.5000 tan(x1) ans = 0.5774 1.0000 1.7321 cot(x1) ans = 1.7321 1.0000 0.5774 5、用四舍五入的方法将数组2.4568 6.3982 3.9375 8.5042取整。 b=2.4568 6.3982 3.9375 8.5042; round(b) ans = 2 6 4 9 6、将矩阵、和组合成两个新矩阵： =（1）组合成一个4x3的矩阵，第一列为按列顺序排列的a矩阵元素，第二列为按列顺序排列的b矩阵元素，第三列为按列顺序排列的c矩阵元素，即

3、（2）按照a、b、c的列顺序组合成一个行矢量，即 a=4 2;5 7; b=7 1;8 3; c=5 9;6 2; % (1) d=a(:) b(:) c(:) d = 4 7 5 5 8 6 2 1 9 7 3 2 % (2) e=a(:);b(:);c(:) e = 4 5 2 7 7 8 1 3 5 6 9 2 或利用（1）中产生的d e=reshape(d,1,12) ans = 4 5 2 7 7 8 1 3 5 6 9 2 7 将(x-6)(x-3)(x-8)展开为系数多项式的形式。 a=6 3 8; pa=poly(a); ppa=poly2sym(pa) ppa = x3-17

4、*x2+90*x-144 8求解多项式x3-7x2+2x+40的根。 r=1 -7 2 40; p=roots(r); 5.0000 4.0000 -2.0000 9求解在x=8时多项式(x-1)(x-2) (x-3)(x-4)的值。 p=poly(1 2 3 4); polyvalm(p,8) ans = 840 10计算多项式乘法(x2+2x+2)(x2+5x+4)。 c=conv(1 2 2,1 5 4) c = 1 7 16 18 8 11计算多项式除法(3x3+13x2+6x+8)/(x+4)。 d=deconv(3 13 6 8,1 4) d = 3 1 2 12对下式进行部分分式

5、展开： a=1 3 4 2 7 2; b=3 2 5 4 6; r,s,k=residue(b,a) r = 1.1274 + 1.1513i 1.1274 - 1.1513i -0.0232 - 0.0722i -0.0232 + 0.0722i 0.7916 s = -1.7680 + 1.2673i -1.7680 - 1.2673i 0.4176 + 1.1130i 0.4176 - 1.1130i -0.2991 k = 13计算多项式的微分和积分。 p=4 12 14 5; pder=polyder(p); pders=poly2sym(pder) pint=polyint(p);

6、 pints=poly2sym(pint) pders = 12*x2-24*x-14 pints = x4-4*x3-7*x2+5*x 14解方程组。 a=2 9 0;3 4 11;2 2 6; b=13 6 6; x=ab x = 7.4000 -0.2000 -1.4000 15矩阵，计算a的行列式和逆矩阵。 a=4 2 -6;7 5 4 ;3 4 9; ad=det(a) ai=inv(a) ad = -64 ai = -0.4531 0.6562 -0.5937 0.7969 -0.8437 0.9062 -0.2031 0.1562 -0.0937 16 y=sin(x)，x从0到

7、2p，Dx=0.02p，求y的最大值、最小值、均值和标准差。 x=0:0.02*pi:2*pi; y=sin(x); ymax=max(y) ymin=min(y) ymean=mean(y) ystd=std(y) ymax = 1 ymin = -1 ymean = 2.2995e-017 ystd = 0.7071 17在不加专门指定的情况下，以下符号表达式中的哪一个变量是独立自由变量？sym(sin(w*t) wsym(a*exp(-X) asym(z*exp(j*theta) z18求以下符号矩阵的行列式值和逆，所得结果应用“子表达式替换”简洁化。A=sym(a11,a12,a13;

8、a21,a22,a23;a31,a32,a33);det(A)ans =a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22inv(A)ans =(a22*a33-a23*a32)/(a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22), -(a12*a33-a13*a32)/(a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a

9、22), (a12*a23-a13*a22)/(a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22) -(a21*a33-a23*a31)/(a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22), (a11*a33-a13*a31)/(a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22), -(a11*a23-a13*a2

10、1)/(a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22)(a21*a32-a22*a31)/(a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22), -(a11*a32-a12*a31)/(a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22), (a11*a22-a12*a21)/(a11*a22*a33-a11*a2

11、3*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22)subexpr(ans,w)w = a11*a22*a33-a11*a23*a32-a21*a12*a33+a21*a13*a32+a31*a12*a23-a31*a13*a22ans = (a22*a33-a23*a32)/w, -(a12*a33-a13*a32)/w, (a12*a23-a13*a22)/w -(a21*a33-a23*a31)/w, (a11*a33-a13*a31)/w, -(a11*a23-a13*a21)/w (a21*a32-a22*a31)/w, -(a

12、11*a32-a12*a31)/w, (a11*a22-a12*a21)/w19(1)通过符号计算求的导数；(2)然后根据此结果，求和。（1）syms t y = abs(sin(t);d=diff(y,t)d =abs(1,sin(t)*cos(t)(2)subs(d,t,0)ans = 1subs(d,t,pi/2)ans = 020在区间，画出曲线，并计算。syms t x zy = sin(t)/t;z = int(y,t,0,x);ezplot(z,x,0,2*pi)subs(z,x,4.5)ans = 1.6541404143792421求方程的解。x,y = solve(x2+y

13、2=1,x*y=2,x,y)x = -1/2*(1/2*5(1/2)+1/2*i*3(1/2)3+1/4*5(1/2)+1/4*i*3(1/2) -1/2*(1/2*5(1/2)-1/2*i*3(1/2)3+1/4*5(1/2)-1/4*i*3(1/2) -1/2*(-1/2*5(1/2)+1/2*i*3(1/2)3-1/4*5(1/2)+1/4*i*3(1/2) -1/2*(-1/2*5(1/2)-1/2*i*3(1/2)3-1/4*5(1/2)-1/4*i*3(1/2)y = 1/2*5(1/2)+1/2*i*3(1/2) 1/2*5(1/2)-1/2*i*3(1/2) -1/2*5(1/

14、2)+1/2*i*3(1/2) -1/2*5(1/2)-1/2*i*3(1/2)22求一阶微分方程，的解。s = dsolve(Dx=a*t2+b*t,x(0)=2)s =1/3*a*t3+1/2*b*t2+223求边值问题，的解。f,g = dsolve(Df=3*f+4*g,Dg=-4*f+3*g,f(0)=0,g(0)=1,x)f =exp(3*x)*sin(4*x)g =exp(3*x)*cos(4*x)24求微分方程的通解，并求任意常数为1时的解。s = dsolve(4*y*Dy+5*x=0,x)s = 1/2*(-5*x2+4*C1)(1/2) -1/2*(-5*x2+4*C1)(1/2)syms C1f = subs(s,C1,1)f = 1/2*(-5*x2+4)(1/2) -1/2*(-5*x2+4)(1/2)