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 32 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)
01234567-40-35-30-25-20-15-10-505
（2）x=0:pi/10:2*pi; y=(exp(-x.^2/2))/2*pi; plot(x,y) 0123456700.20.40.60.811.21.41.6 （3）x=-8::8; y=sqrt((64-x.^2)/2); plot(x,y) -8-6-4-2024680123456 （4）t=0::8*pi; >> x=t.*sin(t); >> y=t.*cos(t); >> plot(x,y) -25-20-15-10-50510152025-30-20-100102030 例 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')
01234567-0.4-0.200.20.40.60.81x from 0 to 2 variable xvariable y曲线y1=e(-0.5x)曲线y2=cos(4x)e-0.5x y1y2
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 = 0 >> z=sin(3.*(t1)).*cos(t1) z = 0 >> plot(t,x,t,y,'k:',t1,z,'bp'); 00.511.522.533.5-1-0.8-0.6-0.4-0.200.20.40.60.81 （2）subplot(1,2,1); >> scatter(x1,y1,10)； >> title('y=')； >> subplot(1,2,2); >> scatter(x,y,10) -1-0.8-0.6-0.4-0.200.20.40.60.81-1-0.8-0.6-0.4-0.200.20.40.6 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]) 01234-1-0.8-0.6-0.4-0.200.20.40.60.812040608010000.10.20.30.40.50.60.70.80.9 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轴仍为线性刻度
1230246802468123
100.5100.710-0.01810-0.0010246810-0.310-0.210-0.1
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); 0.5 1 1.5 2302106024090270120300150
3301800
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);
-1-0.500.51-1-0.500.5102468
（2）t=0::1; x=t; >> y=t.^2; >> z=t.^3; >> plot3(x,y,z);
00.20.40.60.8100.5
100.20.40.60.81
7、x=-30::0;
>> y=0::30; >> [x,y]=meshgrid(x,y); >> z=10.*sin(sqrt(x.^2+y.^2))./sqrt(1+x.^2+y.^2); >> meshc(x,y,z); 绘制曲面图和等高线
-30-20-10001020
30-4-20246
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)<= & abs(y)<=; >> fxy(i)=NaN; >> surf(x,y,fxy) %绘制三维曲面图
-4-2024-4-2024-4-3-2-10
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; -40-2002040-20-10010200100200300400 第五章 二1、a=rand(1,30000); mean(a) %求平均数 ans = >> b=std(a) %求标准差 b = >> c=max(a) c = >> d=min(a) d = size(find(a>)/size(a) %求大于的随机数个数占总数的百分比 ans = 2、h=[466,715,950,1422,1635]; >> w=[,,,,]; >> hh=[500,900,1500]; >> ww=interp1(h,w,hh,'spline') ww = 3、x=linspace(1,10,50); y=log(x); f=polyfit(x,y,5); %求曲线的拟合 >> yy=polyval(f,x); >> plot(x,y,'r-',x,yy,'g.')
1234567891000.511.522.5
5、（1）、（2） p1=[1,2,0,7]; p2=[1,-2]; p3=[1,0,5,1]; p12=conv(p1,p2); >> p=p12+[zeros(1,size(p12,2)-size(p3,2)),p3]; >> roots(p) ans =
+ -
（3） a=[-1,4,3;2,1,5;0,5,6]; >> y1=polyval(p,a) %以矩阵a的每一个元素为自变量 y1 = -29 291 95 19 -3 697 -13 697 1427 >> y2=polyvalm(p,a) %以矩阵a为自变量