Matlab上机练习二答案

Matlab上机练习二

班级 学号 姓名

按要求完成题目，并写下指令和运行结果。（不需要画图） 1、 求x???4?8i3?5i2?7i1?4i7?5i??的共轭转置。 3?2i7?6i9?4i3?9i4?4i??>> 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.0000i 3.0000 - 2.0000i 3.0000 - 5.0000i 7.0000 + 6.0000i 2.0000 + 7.0000i 9.0000 - 4.0000i 1.0000 - 4.0000i 3.0000 + 9.0000i 7.0000 + 5.0000i 4.0000 - 4.0000i 2、计算a???693??241?与b???468?的数组乘积。 275????>> a=[6 9 3;2 7 5];

>> b=[2 4 1;4 6 8]; >> a.*b ans =

12 36 3 8 42 40 3、

?492??37????对于AX?B，如果A???764?，B??26?，求解

???357???28??X。

>> A=[4 9 2;7 6 4;3 5 7];

>> B=[37 26 28]’; >> X=A\\B X =

-0.5118 4.0427 1.3318

4、 a???，b??336?4???>> a=[1 2 3;4 5 6]; >> b=[8 –7 4;3 6 2]; >> a>b ans =

0 1 0 1 0 1 >> a>=b ans =

0 1 0

?125??8?74?，观察62??a与b之间的六种关系运算的结果。

1 0 1 >> a

1 0 1 0 1 0 >> a<=b ans =

1 0 1 0 1 0 >> a==b ans =

0 0 0 0 0 0 >> a~=b ans =

1 1 1 1 1 1

5、a??50.20?8?0.7?，在进行逻辑运算时，a相当于什么样的逻辑量。 相当于a=[1 1 0 1 1]。

6、 角度x??304560?，求x的正弦、余弦、正切和余切。

>> x=[30 45 60]; >> x1=x/180*pi; >> sin(x1) ans =

0.5000 0.7071 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

7、 用四舍五入的方法将数组[2.4568 6.3982 3.9375 8.5042]取整。 >> b=[2.4568 6.3982 3.9375 8.5042]; >> round(b) ans =

2 6 4 9

4813??1?543?2??，B??6?23?8?，求C1=A*B’;C2=A’*B;C3=A.*B,并求上述?36?5?98、设A??????????13?97???2?7?12?8??

所有方阵的逆阵。

>> 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]; >> C1=A*B' C1 =

19 -82 30 12 27 3 -38 54 29 >> C2=A'*B C2 =

-15 16 -24 36 63 -17 93 -105 22 6 117 -60 19 46 84 -10 >> C3=A.*B C3 =

5 16 24 -26 -18 -12 -15 72 -2 -21 108 -56 >> inv(C1) ans =

0.0062 0.0400 -0.0106 -0.0046 0.0169 0.0030 0.0168 0.0209 0.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.9553 -0.2391 -0.1997 0.2700 0.9667 0.2420 0.2021 -0.2732 -0.4473 -0.1120 -0.0935 0.1264 -1.1259 -0.2818 -0.2353 0.3182 >> inv(C3)

??? Error using ==> inv Matrix must be square.

9、设x=rcost+3t,y=rsint+3,分别令r=2,3,4，画出参数t=0~10区间生成的x~y曲线。 >> 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')

10、设f(x)=x- 4x +3x- 2x+ 6

(1) 在x=[-2,8]之间取100个点，画出曲线，看它有几个过零点。（提示：用polyval 函数） >> x=linspace(2,8,100);

>> y=polyval([1 0 -4 3 -2 6],x); >> plot(x,y,'b',x,0,'y')

542

(2) 用roots函数求此多项式的根。 t=[1 0 -4 3 -2 6] p=roots(t)

11、设x=sint, y=sin(nt+a)，

(1) 若a=1，令n =1,2,3,4,在四个子图中分别画出其曲线。

(2) 若n=2,取a=0，π/3，π/2，及π，在四个子图中分别画出其曲线。 (1) a=1; x=sin(t);y1=sin(1*t+a);y2=sin(2*t+a);y3=sin(3*t+a);y4=sin(4*t+a); subplot(2,2,1);plot(x,y1);subplot(2,2,2);plot(x,y2);subplot(2,2,3);plot(x,y3);subplot(2,2,4);plot(x,y4) (2)

n=2;x=sin(t);y1=sin(2*t+0);y2=sin(2*t+pi/3);y3=sin(2*t+pi/2);y4=sin(2*t+pi);subplot(2,2,1);plot(x,y1);subplot(2,2,2);plot(x,y2);subplot(2,2,3);plot(x,y3);subplot(2,2,4);plot(x,y4)

注：本题好像题目本身就有问题，因为“t” 取值不明，所以运行不了，画不出图来。转化为符号式来画图也是不可行的。（纯粹个人意见，可能我想错了方向，会做的同学请上传一下正确的做法）

12、 求解多项式x3-7x2+2x+40的根。 >> r=[1 -7 2 40]; >> p=roots(r); -0.2151 0.4459 0.7949 0.2707

13、符号函数绘图法绘制函数x=sin(3t)cos(t)，y=sin(3t)sin(t)的图形，t的变化范围为[0,2?]。 >> syms t

>> ezplot(sin(3*t)*cos(t),sin(3*t)*sin(t),[0,pi])

14、设

,求

x=sym('x'); y=(sin(x))^4+(cos(x))^4;diff(y,10) 15、?(9?x2)3x6dx

x=sym('x');int((sqrt((9-x^2)^3)/x^6),x)

2dx16、?

1x?x3x=sym('x');int(x+x^3,x,1,2) 17、求级数的和:

n?2 ?nn?12?sym(‘n’);symsum(n+2/2^n,n,1,inf)

?xy18、利用函数int计算二重不定积分

??( x?y)edxdyx=sym('x'); y=sym('y');z=(x+y)*exp(-x*y);a=int(z,x);int(a,y) 19、试求出如下极限。 （1）lim(3x?9x)； （2）limx??1xxyxy?1?1x?0y?0；