1、第 1 章練習題答案1- 2.A=1.2 ;B=-4.6 ;C=8.0 ;D=3.5 ;E=-4.0; T=atan(2*pi*A+E/(2*pi*B*C)/D)1- 3 clearx=45;(si n(x* pi/180)+sqrt(35)/(72A(1/5)1- 4clear a=5.67;b=7.811;exp(a+b)/log10(a+b)1- 5clearx=3;y=sqrt(x)-6*(x+1/x)+(x-3.2)A2/(x+7.7)A31- 6 clear r=15;d=2*rc=2*pi*rs=pi*r*r1- 7cleara=8.5;b=14.6;c=18.4;s=(a+b+

2、c)/2; area=sqrt(s*(s-a)*(s-b)*(s-c)第 2 章練習題答案2- 1 c learA= 3 1 1;2 1 2;1 2 3;B=1 1 -1;2 -1 0;1 -1 1;2*A+B,4*(AA2)-3*(BA2),A*B,B*A,A*B-B*A2- 2 clearA=1/3 0 0;0 1/4 0;0 0 1/7;B=inv(inv(A)-eye(3)*62- 3 clearB=1 2 -3 -2;0 1 2 -3;0 0 1 2;0 0 0 1; C=1 2 0 1;0 1 2 0;0 0 1 2;0 0 0 1; A=(inv(2*eye(4)-inv(C)*

3、B)*inv(C)'2- 4 clearA=2 -1;-1 2;B=0 -2;-2 0;X=(1/2)*(B+2*A)2- 5 clearA=2 -3 O 2;1 5 2 1;3 -1 1 -1;4 1 2 2;B=8;2;7;12;X=AB2- 6A=3 12 4 7 O 8 1;p=poly2str(A,'x')x=roots(A)2-7A=3 12 4 7 O 8 1;B=1 -3 5 -15 O;div,rest=deconv(A,B)第三章3- 1. x=sym('x');limit(cos(sqrt(x)A( pi /x),x,O,'

4、right')ans =exp(-1/2*pi)3- 2.x=sym('x');f=sym('(3*si n(x)+xA2*cos(1/x)/(1+cos(x)*log(1+x)');limit(f,x,O)ans =3/23- 3 .x=sym('x') ;f=sym('(sqrt(4*x2+x-1)+x+1)/sqrt(x2+si n(x)');limit(f,x,-inf)ans =13- 4.s yms x y;f=sym('(x2+y2)(x2*y2)');limit(limit(f,x,O),y

5、,O)ans =13-5.x=sym('x');y=sym('(ta n( sqrt(x+sqrt(x+sqrt(2*x)2');diff(y,x)ans=ta n(x+(x+2A(1/2)*xA(1/2)A(1/2)A(1/2)*(1+ta n(x+(x+2A(1/2)*xA(1/2 )A(1/2)A(1/2)A2)/(x+(x+2A(1/2)*xA(1/2)A(1/2)A(1/2)*(1+1/2/(x+2A( 1/2)*xA(1/2)A(1/2)*(1+1/2*2A(1/2)/xA(1/2)3-6.x=sym('x');y=sym('

6、cos(x2)*si n(1/x)A2');diff(y,x)ans =-2*si n(xA2)*x*si n(1/x)A2-2*cos(xA2)*si n(1/x)*cos(1/x)/xA23-7.x=sym('x');y=sym('sqrt(s in( x)-s in( x)3)');int(y,x,O,pi)ans =4/33- 8 .x=sym('x'); y=sym('sqrt(x+1)/(x-1)/x');int(y,x)ans=-(1+x)/(x-1)A(1/2)*(x-1)*(ata n(1/(-1+xA2

7、)A(1/2)-log(x+(-1+xA2)A(1 /2)/(1+x)*(x-1)A(1/2)3- 9. dsolve('D2y+4*Dy+4*y=exp(-2*x)','x')ans =1/2*ex p(-2*x)*(2*C2+2*x*C1+xA2)3-10.dsolve('xA2*Dy+x*y=yA2','y(1)=1')ans =x*exp(1/x)/(exp(1/x)+exp(1/x*t)*x-exp(1/x*t)第四章4- 1.plot(1 2 4 0 5 10 11 21 3 1)4- 2.t=0:pi/100:2*p

8、i;x=cos(t);y=sin(t);plot(x,y);axis('square');4- 3.R=1 5 10 20;I=0:0.1:10;U=I'*R;plot(I,U);4- 4.mon=1:12;temp=0.2 2.3 8.7 18.5 24.6 32.1 36.8 37.1 28.3 17.8 6.4 -3.2; rainf=4.6 3.6 2.1 2.9 3.0 2.7 2.2 2.5 4.3 3.4 2.1 3.7;plot(mon,temp,'r-p',mon,rainf,'b:s');xlabel('mon

9、th');ylabel('temperature and rainfall');%title('Temperature and Rainfall Curve');%text(3,1,'rainfall');%text(9,28.3,'temperature');flag=12;while(flag) text(mon(flag),temp(flag),num2str(mon(flag),num2str(temp(flag);text(mon(flag),rainf(flag),num2str(mon(flag),num2s

10、tr(rainf(flag)flag=flag-1;end4- 5.A=1 1 1 1 1 1 1;1 2 2 2 2 2 1;1 2 2 3 2 2 1;1 2 2 2 2 2 1;1 1 1 1 1 11;plot(A);4- 6.x=-10:0.5:10;y=x;X,Y=meshgrid(x,y);Z=XA2+YA2;surf(X,Y,Z);4-7.clccleartheta=(0:0.1:2)*pi;phi=(0:0.1:2)*pi;x=zeros(length(theta),length(phi);y=x;z=x;for i=1:length(theta)for j=1:length

11、(phi) x(i,j)=2*sin(theta(i)*cos(phi(j); y(i,j)=3*sin(theta(i)*sin(phi(j);endz(i,:)=4*cos(theta(i);endaxis square subplot(2,2,1); mesh(x,y,z) title(' 三維 ') subplot(2,2,2); mesh(x,y,z) view(90,0) title('X 軸視圖 ') subplot(2,2,3); mesh(x,y,z) view(0,0) title('Y 軸視圖 ') subplot(2,2,

12、4); mesh(x,y,z) view(0,90) title('Z 軸視圖 ')第五章5- 1：clear;x=-2:0.01:2;for ii=1:401if(x(ii)<0) y(ii)=x(ii)+1;elseif (x(ii)<1)&(x(ii)>=0) y(ii)=1;elseif(x(ii)>=1) y(ii)=x(ii)A3;endendplot(x,y);5- 2：function y=fenduan(x)if(x<0)y=x+1;elseif (x<1)&(x>=0)y=1;elseif(x>

13、=1)y=x3;end5- 3：建立文件 arcsin.m 如下： function y=arcsin(x) if abs(x)>1disp(' 輸入參數的絕對值必須小于 1');return;endn=0;y=0;u=x;while u>epsu=(jieche ng(2* n)*xA(2* 門+1)/(2八(2* n)*(jieche ng(門)八2*(2* n+1); y=y+u;n=n+1;end建立文件 jiecheng.m 如下：function y=jiecheng(n) if(n<0)disp(' 輸入參數必須為大于或等于 0 的整數

14、'); return;endif n=0 | n=1y=1;elsey=n*jiecheng(n-1);end5- 4：建立 shuliang.m 文件如下：function tu,ji=shuliang(tou,jiao) tu=jiao/2-tou;ji=2*tou-jiao/2;在命令窗口中執行： >>a,b=shuliang(36,100)5- 5：for i=2:999 ge=rem(i,10); shi=rem(fix(i/10),10); bai=fix(i/100);if rem(ge+shi+bai),2) & isprime(i) disp(i)

15、;endend5- 6：function y=issushu(n) y=1;for ii=2:fix(n/2) if(rem(n,ii)=0) y=0;endend習題 1(程序段后面的是其運行的結果)例題 1.2 function lt12 a=1+2*i;b=3-4*i;disp('a*b=disp('a+b=',num2str(a+b) disp('a-b=',num2str(a-b)',num2str(a*b) disp('a/b=',num2str(a/b) b=-2+6i a*b=11+2ia/b=-O.2+O.4i>>例題 1.3 function z=lt13 x=(-3.5)*pi/18O;y=6.7*pi/18O; z=sin(abs(x)+abs(y)/sqrt(cos(abs(x+y);>> lt12 a+b=4-2i a->> lt13ans =O.1772 例題 1.4 function p1=lt14 r=O.O1O7;n=2O15-2OOO;pO=12.9533E8; ans =1.5195e+OO9>>例題 1.5 fun