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表示自然指數(shù)函數(shù)log自然對數(shù)函數(shù)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指數(shù)可以是標量(如 y=2).底數(shù)也可以是 標量(如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);%函數(shù)使用全對數(shù)坐標，x,y均采用常用對數(shù)刻度subplot(2,2,4);semilogy(t,y);%函數(shù)使用半對數(shù)坐標，y軸為常用對數(shù)刻度,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)% 求平均數(shù)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的隨機數(shù)個

15、數(shù)占總數(shù)的 百分比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>>統(tǒng)計小寫字母的個數(shù)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加網(wǎng)格線三維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、網(wǎng)格圖多項式求導例：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為余數(shù)向 量p=polyder(p)表示求p的導數(shù)p=poleder(p,q)表示求 p 乘以 q的導數(shù)p,q=poleder(p,q)表示 p 除以 q的導數(shù)多項式求根例：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