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.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);
a=1
>> r=a.*(1+cos(x));
polar(x,r);
270
6、（1）
t=0:pi/10:2*pi;
>> x=exp((-t)/20).*cos(t); >> y=exp((-t)/20).*sin(t); >> z=t;
>> plot3(x,y,z);
（2）t=0:0.01:1;
x=t;
>> y=t.^2;
>> z=t.^3;
>> plot3(x,y,z);
7、x=-30:0.1:0;
>> y=0:0.1:30;
>> [x,y]=meshgrid(x,y);
>>
z=10.*sin(sqrt(x.^2+y.^2))./sqrt(1+x. ^2+y.^2);
>> meshc(x,y,z);
绘制曲面图和等高线
8、x=linspace(-3,3,100);
>> y=linspace(-3,3,100);
>> [x y]=meshgrid(x,y); %可以将向量转化为矩阵
>> fxy=-5./(1+x.^2+y.^2);
>> i=find(abs(x)<=0.8 &
abs(y)<=0.5);
>> fxy(i)=NaN; >>
surf(x,y,fxy) %绘制三维曲面图
9、
u=linspace(1,10,100);
v=linspace(-pi,pi,100);
[u v]=meshgrid(u,v);
x=3.*u.*sin(v);
y=2.*u.*cos(v);
z=4*u.^2;
surf(x,y,z);
shading interp;
第五章
二1、a=rand(1,30000);
mean(a) %求平均数 ans =
0.5010
>>
b=std(a)
%求标准差
b =
0.2882
>> c=max(a)
c = 0.9999 >> d=min(a)
d =。