matlab课后答案完整版

1、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

2、;>> c=97 67;-78 75;32 5;0 -12;>> d=65 5;-23 -59;54 7;>> e=b*ce =527111574-1133666419783112(2)a=50:1:100二、1、x=-74;y=-27;z=(si n(x42+y.A2)/(sqrt(ta n(abs(x+y )+pi)z =-0.09012、a=-3.0:0.1:3.0;>> b=exp(-0.3*a).*si n( a+0.3)3、x=2 4;-0.45 5;y=log(x+sqrt(1+x.A2)/2y =0.72181.0474-0.21

3、801.15624、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./BA.AB表示两个矩阵对应元素进行乘方运算A.A2表示a中的每个元素的平方AA2表示A*A例：x=1,2,3;y=4,5,6;z=x.Ayz=1 32 729指数可以是标量(如 y=2).底数也可以是 标量(如x=2)5、a=1+2i;>> b=3+4i;>> c=exp(pi*i)/6)c =0.86

4、60 + 0.5000id=c+a*b/(a+b)d =1.6353 + 1.8462i第二章二、4、( 1)y=0;k=0;>> while y<3k=k+1;y=y+1/(2*k-1);end>> display(k-1,y-1/(2*k-1)ans =56.00002.9944第三章二 1 (1)x=0:pi/10:2*pi;>> y=x-x.A3/6;>> plot(x,y)5>> x=t.*si n(t);>> y=t.*cos(t);>> plot(x,y)3020r10 .- I f 厂i0

5、.f i J_/I /xJIILX-jTJJ-10-/ J-20 . _ S«_'L-30 IIEEIIIIII-25-20-15-10-50510152025(2) x=0:pi/10:2*pi; y=(exp(-x.A2/2)/2*pi; plot(x,y)1.61.41.210.80.60.40.2001234567(3) x=-8:0.01:8; y=sqrt(64-x.A2)/2); plot(x,y)例3.4x=0:pi/100:2*pi;y1=exp(-0.5*x);y2=exp(-0.5*x).*si n( 2*x);plot(x,y1,x,y2)>&g

6、t; title('x from 0 to 2pi');>> xlabel('variable x');>> ylabel('variable y');>> text(1.5,0.5,' 曲线 y1=eA(-0.5x)');>>text(3,0.1,'曲 线y2=cos(4pix)eA-0.5x');>> lege nd('y1','y2')10.80.6y 0.4eaav 0.20-0.2-0.401234567varia

7、ble xx from 0 to 2曲线 y1=e(-0.5x)曲线 y2=cos(4x)e-0.5x2、(1)y1=2*x-0.5;t=li nspace(0,pi,100);x=si n(3*t).*cos(t);y=si n(3*t).*si n( t);>> k=fi nd(abs(y-x)<1e-2);>> t1=t(k)y1y2(4) t=0:0.1:8*pi;t1 =0.80.6fi0.410.2*0J-0.2-1.-0.4i/-0.6i/-0.8jF-1r|100.511.522.533.5-1-0.8-0.6-0.4-0.200.20.40.60

8、.814、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轴仍为线性刻度01234-0.00110-0.0181000.79331.04722.09443.1416>>

9、 z=si n(3.*(t1).*cos(t1)z =00.48410.00000.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)0.60.40.20-0.2-0.4-0.6-0.8-13、subplot(1,2,1);x=0:0.01:pi;y=si n(

10、1./x);plot(x,y)subplot(1,2,2);fplot('si n(1./x)',1,100)10.80.60.40.20-0.2-0.4-0.6-0.8-15、(1)theta=li nspace(-pi,pi,100); ro=5.*cos(theta)+4; polar(theta,ro);(2)x=li nspace(0,2*pi,100);a=1>> r=a.*(1+cos(x);polar(x,r);7、x=-30:0.1:0;>> y=0:0.1:30;>> x,y=meshgrid(x,y);>>z

11、=10.*si n( sqrt(x.A2+y.A2)./sqrt(1+x .A2+y.A2);>> meshc(x,y,z);绘制曲面图和等高线6、(1)t=0:pi/10:2*pi;>> x=exp(-t)/20).*cos(t);>> y=exp(-t)/20).*si n(t);>> z=t;>> plot3(x,y,z);(2) t=0:0.01:1; x=t;>> y=t.A2;>> z=t.A3;>> plot3(x,y,z);8、x=linspace(-3,3,100);>>

12、; y=li nspace(-3,3,100);>>xy=meshgrid(x,y);%可以将向量转化为矩阵>> fxy=-5./(1+x.A2+y.A2);>> i=fi nd(abs(x)<=0.8& abs(y)<=0.5);>> fxy(i)=NaN;>>surf(x,y,fxy)%绘制三维曲面图9、u=li nspace(1,10,100); v=li nspace(-pi,pi,100);u v=meshgrid(u,v); x=3.*u.*s in (v); y=2.*u.*cos(v); z=4*u

13、.A2; surf(x,y,z);shad ing in terp;>> ww=i nterp1(h,w,hh,'spli ne')ww =6.4903 3.5226 2.38453、x=linspace(1,10,50); y=log(x);f=polyfit(x,y,5);%求曲线的拟合>> yy=polyval(f,x);>> plot(x,y,'r-',x,yy,'g.')2.5第五章二 1、a=rand(1,30000); mea n(a)% 求平均数ans =0.5010>>b=std(

14、a)%求标准差b =0.2882>> c=max(a)c =0.9999>> d=mi n(a)d =3.5706e-00521.510.50123456789105、(1 )、(2)p1=1,2,0,7;p2=1,-2;p3=1,0,5,1;p12=co nv(p1,p2);>>p=p12+zeros(1,size(p12,2)-size(p3,2),p3;>> roots(p)ans =-3.46560.6128 + 1.6278i0.6128 - 1.6278isize(fi nd(a>0.5)/size(a)%求大于0.5的随机数个

15、数占总数的 百分比1.2400(3)a=-1,4,3;2,1,5;0,5,6;>>y1=polyval(p,a)ans =0.50322、h=466,715,950,1422,1635;>> w=7.04,4.28,3.40,2.52,2.13;%以矩阵a的每一个元素为自变量y1 =-292919519-3697-136971427篮筐-1 -0.9-0.8-0.7-0.6-0.5-0.4-0.3-0.2-0.10 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.91y=exp(-0.2*x).*cos(x)+24选 1218y2=polyvalm(

16、p,a)%以矩阵a为自变量y2 =3912084327350226934207720377558926、( 1)z=fzero('3*x-si n(x)+1',0)%求x=0时附近的根 z =-0.4903第八章二、2t=0:pi/20:2*pi;x=si n(t);y=cos(t);x1=si n( 7*t);y1=cos(7*t); h=plot(x,y,x1,y1); set(h,'marker','x','l in ewidth',2); set(gca,'xtick',-1:0.1:1);title(&#

17、39; 篮筐')2.82.62.42.221.81.61.402468101214164、t=-pi:pi/100:pi;x=cos(t);y=s in( t);z=t;h=plot(t,x,t,y,t,z);set(h,'li nestyle','-','l in ewidth',3)字符串例ch='Welcome to Beiji ng'subch=ch(12:18)个字符串(空格也算)3、x=0:pi/10:5*pi;y=exp(-0.2*x).*cos(x)+2; h=plot(x,y);set(gca,'

18、color','red','li nestyle','-','li newidth',3);text(5,2.4,'y=exp(-0.2*x).*cos(x)+2' );subch =Beiji ng>>k=fi nd(ch>='a'&ch<='z');找到所有的小写字母的位置>>ch(k)=ch(k)-('a'-'A');将小写字母变成大写字母>> char(ch)ans =WELCOM

19、E TO BEIJING>>统计小写字母的个数ans =14例：已y=1-1/2+1/3-1/4-1/100值y=0;>> n=100;>> for i=1:100; y=y+(-i)F-i)/i; end>> disp(y)0.6882绘制二维曲线图x=0:pi/100:2*pi;>>y1=0.2*exp(-0.5*x).*cos(4*pi*x);>>y2=1.5*exp(-0.5*x).*cos(pi*x); >> plotyy(x,y1,x,y2);0.2 c.clen gth(k)-0.2 tC1CL.

20、01234绘制三维图像例: x=sint+tcost y=cost-ts int z=tt=0:pi/10:10*pi;x=s in( t)+t.*cos(t);y=cos(t)-t.*si n(t);z=t;plot3(x,y,z);axis(-3030 -30 30 0 35);坐标轴的最大值与最小值title('li nein 3-D space');图形的题目>>xlabel('x');ylabel('y');zlabel('z');标注坐标>>grid加网格线三维on;line in 3-D sp

21、acezz sinx2 cosy2, x0, ,y 0, /2x,y=meshgrid(0:pi/100:pi,0:pi/100:pi/2);>> z=si n(x.A2)+cos(y42);>> mesh(x,y,z);>> axis(0 4 0 1.8 -1.5 1.5);t=0:pi/20:2*pi;subplot(1,2,1);x,y,z=cyli nder(si n(t),30);surf(x,y,z);绘制三维曲面图subplot(1,2,2);>> x,y,z=peaks(100);>> mesh(x,y,z);绘制三维

22、网格图多项式求导例：f(x)=1/xA2+5p=1;>> q=1,0,5;>> p,q=polyder(p,q)注: c=conv(a,b)表示a多项式与b多项式乘积p,r=deconv(a,b)表示 a 多项式与b多项式相除 其中p为商向量r为余数向 量p=polyder(p)表示求p的导数p=poleder(p,q)表示求 p 乘以 q的导数p,q=poleder(p,q)表示 p 除以 q的导数多项式求根例：f(x)=2xA4-12xA3+3xA2+5 p=2,-12,3,0,5;>> x=roots(p);>> p=2,-12,3,0,5;x=roots(p)求方程f(x)=0 的根x =5.72460.8997-0.3122 + 0.6229i-0.3122 - 0.6229i>> g=poly(x)已知多项式的根求多项式g =1.0000-6.00001.5000-0.0000 2.5000符号求导例7.3x=a(t-tsi nt)y=b