matlab课后答案完整版

ones表示1矩阵

zeros表示0矩阵

ones(4)表示4x4的1矩阵

zeros(4)表示4x4的0矩阵

zeros(4,5)表示4x5的矩阵

eye(10,10)表示10x10的单位矩阵

rand(4,5)表示4x5的伴随矩阵

det(a)表示计算a的行列式

inv(a)表示计算a的逆矩阵

Jordan(a)表示求a矩阵的约当标准块rank(a)表示求矩阵a的秩

[v,d]=eig(a)对角矩阵

b=a’表示求a矩阵的转置矩阵

sqrt表示求平方根

exp表示自然指数函数

log自然对数函数

abs绝对值

第一章

一、5（1）

b=[97 67 34 10;-78 75 65 5;32 5 -23 -59]; >> c=[97 67;-78 75;32 5;0 -12];

>> d=[65 5;-23 -59;54 7];

>> e=b*c

e =

5271 11574

-11336 664

1978 3112

（2）a=50:1:100

二、1 、x=-74;

y=-27;

z=(sin(x.^2+y.^2))/(sqrt(tan(abs(x+y)))+pi) z =

-0.0901

2、a=-3.0:0.1:3.0;

>> b=exp(-0.3*a).*sin(a+0.3)

3、x=[2 4;-0.45 5];

y=log(x+sqrt(1+x.^2))/2

y =

0.7218 1.0474

-0.2180 1.1562

4、a*b表示a矩阵和b矩阵相乘

a.*b表示a矩阵和b矩阵单个元素相乘

A(m,n)表示取a矩阵第m行，第n列

A(m,:)表示取a矩阵第m行的全部元素

A(:,n)表示取a矩阵的第n列全部元素

A./B表示a矩阵除以b矩阵的对应元素，

B.\A等价于A./B

A.^B表示两个矩阵对应元素进行乘方运算

A.^2表示a中的每个元素的平方

A^2表示A*A

例：x=[1,2,3];

y=[4,5,6];

z=x.^y

z=

1 3

2 729

指数可以是标量（如y=2）.底数也可以是标量（如x=2）

5、a=1+2i;

>> b=3+4i;

>> c=exp((pi*i)/6)

c =

0.8660 + 0.5000i

d=c+a*b/(a+b)

d =

1.6353 + 1.8462i

第二章

二、4、（1）

y=0;k=0;

>> while y<3

k=k+1;

y=y+1/(2*k-1);

end

>> display([k-1,y-1/(2*k-1)])

ans =

56.0000 2.9944

第三章

二1（1）

x=0:pi/10:2*pi;

>> y=x-x.^3/6;

>> plot(x,y)

（2）x=0:pi/10:2*pi; y=(exp(-x.^2/2))/2*pi;

plot(x,y)

（3）x=-8:0.01:8;

y=sqrt((64-x.^2)/2);

plot(x,y)

（4）t=0:0.1:8*pi; >> x=t.*sin(t);

>> y=t.*cos(t); >> plot(x,y)

例3.4

x=0:pi/100:2*pi;

y1=exp(-0.5*x);

y2=exp(-0.5*x).*sin(2*x);

plot(x,y1,x,y2)

>> title('x from 0 to 2{\pi} ');

>> xlabel('variable x');

>> ylabel('variable y');

>> text(1.5,0.5,'曲线y1=e^(-0.5x)');

>> text(3,0.1,'曲线y2=cos(4{\pi}x)e^{-0.5x}');

>> legend('y1','y2')

variable x

v

a

r

i

a

b

l

e

y

2、(1)

y1=2*x-0.5;

t=linspace(0,pi,100);

x=sin(3*t).*cos(t);

y=sin(3*t).*sin(t);

>> k=find(abs(y-x)<1e-2);

>> t1=t(k)

t1 =

0 0.7933 1.0472

2.0944

3.1416

>> z=sin(3.*(t1)).*cos(t1)

z =

0 0.4841 0.0000 0.0000 -0.0000

>> plot(t,x,t,y,'k:',t1,z,'bp');

（2）subplot(1,2,1);

>> scatter(x1,y1,10)；

>> title('y=2x-0.5')；

>> subplot(1,2,2);

>> scatter(x,y,10)

3、subplot(1,2,1);

x=0:0.01:pi;

y=sin(1./x);

plot(x,y)

subplot(1,2,2);

fplot('sin(1./x)',[1,100])

4、t=0:pi:2*pi;

y=1./(1+exp(-t));

subplot(2,2,1);

%图形窗口的分割

bar(t,'group'); %绘制柱形图（分组）

subplot(2,2,2);

barh(t,'stack'); %绘制柱形图（堆积）

subplot(2,2,3);

loglog(t,y); %函数使用全对数坐标，

x,y均采用常用对数刻度subplot(2,2,4);

semilogy(t,y); %函数使用半对数坐标，y轴为常用对数刻度，x轴仍为线性刻度

10

10

10

10

10

10

10

5、（1）

theta=linspace(-pi,pi,100);

ro=5.*cos(theta)+4;

polar(theta,ro);

（2）

x=linspace(0,2*pi,100);