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 =

2、a=::;

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

3、x=[2 4; 5];

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

y =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 =

+

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

d =

+

第二章

二、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 =

第三章

二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::8;

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

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

>> y=t.*cos(t);

>> plot(x,y)

例

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

y1=exp*x);

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

plot(x,y1,x,y2)

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

>> xlabel('variable x');

>> ylabel('variable y');

>> text,,'曲线y1=e^');

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

variable x

v

a

r

i

a

b

l

e

y

2、(1)

y1=2*;

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 =

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

z =

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

（2）subplot(1,2,1); >> scatter(x1,y1,10)； >> title('y=')； >> subplot(1,2,2); >> scatter(x,y,10)

3、subplot(1,2,1); x=0::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); a=1

>> r=a.*(1+cos(x)); polar(x,r);