1、第一章-4.0計算:2.設 A=1.2,B =-4.6,C =8.0,D = 3.5,E2 A ET arctan 2 BCD 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)/D0.4112 x=pi/180*45; (sin(x)+sqrt(35)/72八(1/5) ans = 2.81583 .設x 45 ,計算sin(x) 、355 72 x=pi/180*45; (sin(x)+sqrt(35)/72八(1/5) ans = 2.81584 .設 a 5.67,b 7.811,計算:e(a b)e(a b)lOg

2、10 a=5.67;b=7.811; exp(a+b)/log10(a+b)ans =6.3351e+005一1235.計算 y Jx 6(x -) (x 3.2) /(x 7.7)在 x 3 時的值。 x x=3; y=sqrt(x)-6*(x+1/x)+(x-3.2)八2/(x+7.7)八3y =-18.26796.已知圓的半徑為15,求其直徑，周長及面積。 r=15;d=2*r; c=pi*2*rc =94.2478 s=pi*r*rs =706.85837已知某三角形的三個邊的邊長為8.5,14.6和18.4,求該三角形的面積。提示：area Js(s a)(s b)(s c)其中：s

3、 (a b c)/2, a、b、C分別為三角形三邊邊長。a=8.5;b=14.6;c=18.4; s=(a+b+c)/2; area=sqrt(s*(s-a)*(s-b)*(s-c) area =第二章60.61061 .設矩陣求：(1)2A B2_ 2(2)4 A 3B(3)ABBA(5)AB BA A=3 1 1;2 1 2;1 2 3；B=1 1 -1;2 -1 0;1 -1 1;(1) 2*A+B ans =73161433 7(2) 4*A*A-3*B*B ans =42 21 3840 19 4640 33 56(3) A*Bans =61 -26 -108 -42(4) B*Aa

4、ns =40041022 2(5) A*B-B*Aans =21 -22 -206 -602.設三階矩陣A B,滿足A1BA 6A BA,其中求矩陣B 。 A=1/3 0 0; 0 1/4 0;0 0 1/7; B=6*inv(inv(A)-eye(3)3.設(2E_ 1 TC 1B)AT其中E是4階單位矩陣,AT是4階矩陣A的轉置。其中10002100321023211000210002101021 B=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)*B

5、)*inv(C)-2100歡迎下載1.000001 -21 4.設二階矩陣A、B、X ,滿足X其中2102AB1 22 0求矩陣X。 A=2 -1;-1 2;B=0 2;-2 0; X=(2*A+B)/2X =20-222x1 x1 3x1 4x15 .求解線性方程組 3x2 2x485x2 2x3 x42X2 X3 X47x2 2x3 2x412 A=2 -3 0 2;1 5 2 1;3 -1 1 -1;4 1 2 2; B=8;2;7;12; x=AB x =3.00000.0000-1.00006 .求解一元六次方程3x6 12x5 4x4 7x3 8x 1 0的根。 A=3 12 4

6、7 0 8 1; x=roots(A) x =-3.8230-0.5275 + 0.8497i-0.5275 - 0.8497i0.5007 + 0.6749i0.5007 - 0.6749i-0.12347 .求多項式3x6 12x5 4x4 7x3 8x 1被(x 3)(x3 5x)除后的結果。 a=3 12 4 7 0 8 1;A=3; b=poly(A); c=1 0 5 0; d=conv(b,c); div,rest=deconv(a,d)div =3 2152rest =000 103 55 7881第三章歡迎下載1.求極限歡迎下載lim cos x xx 0x=sym(x);l

7、imit(cos(sqrt(x)八(pi/x),x,0,right)ans =exp(-1/2*pi)2.求極限3sin xlimx21x cos一 x1 cosx In 1 xx=sym(x);f=sym(3*sin(x)+xA2*cos(1/x)/(1+cos(x)*log(1+x);limit(f,x,0)ans =3/23.求極限4x2 x 1 x 1lim,二x, x2 sin xx=sym(x);f=sym(sqrt(4*xA2+x-1)+x+1)/sqrt(xA2+sin(x);limit(f,x,-inf)ans =14.求極限22、2 2嗎(xy )x yy osyms x

8、y;f=sym(xA2+yA2)A(xA2*yA2);limit(limit(f,x,0),y,0) ans =15 .已知2 I /y tg x x 2x求yx=sym(x);y=sym(tan(sqrt(x+sqrt(x+sqrt(2*x)A2);diff(y,x) ans =tan(x+(x+2A(1/2)*xA(1/2)A(1/2)A(1/2)*(1+tan(x+(x+2A(1/ 2)*xA(1)八(1)A(2)A2)/(x+(x+2A(1*xA(2)A(i)a (1/2)*(1 + 1/2/(x+2A(1/2)*xA(1/2)A(1/2)*(1 + 1/2*2A(1/2)/xA(1/

9、 2)6 .已知 y cos(x2)sin21 x求yx=sym(x);y=sym(cos(xA2)*sin(1/x)A2);歡迎下載diff(y,x) ans =-2*sin(xA2)*x*sin(1/x)A2-2*cos(xA2)*sin(1/x)*cos(1/x)/xA27.求積分 sin x sin3 xdx0x=sym(x);y=sym(sqrt(sin(x)-sin(x)A3);int(y,x,0,pi)ans =4/38.求積分x=sym(x);y=sym(sqrt(x+1)/(x-1)/x);int(y,x)ans =-(1+x)/(x-1)A(1/2)*(x-1)*(atan

10、(1/(-1+xA2)A(1/2)-log(x+(-1+xA2)A(1/2)/(1+x)*(x-1)A(1/2)9.求下列微分方程的通解y 4y 4y e 2xdsolve(D2y+4*Dy+4*y=exp(-2*x),x)ans =1/2*exp(-2*x)*(2*C2+2*x*C1+xA2)10.求微分方程22X y xy y , y |x 1 1dsolve(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)第四章1.已知向量1 2 4 0 5 10 11 21 3 1，請繪圖表示。 A=1 2

11、4 0 5 10 11 21 3 1; plot(A)123456789102.繪制單位圓。提示：使用axis( square)命令保證圖形的縱橫坐標刻度比例 相同。 ezplot(xA2+yA2=1) axis squarex2+y2=11.5廠 1 -_J*、0.5.F.fiiy 0-III-0.5 工/J-1 -L-1.5,一-1.5-1-0.500.511.5x3.繪制伏安特性曲線：U RI ,假設R分別為1,5,10和20。R=1 5 10 20;1=0:0.1:10;U=I*R;plot(I,U)legend R=1 R=5 R=10 R=20grid on伏安特性曲線電流/A4.

12、某地區一年中每月的平均氣溫和平均降雨量如下表 4 - 5所示，請 畫出其圖形，要求標注出坐標軸、數據點位置、數據點大小等。表4 - 5 溫度-降雨量數據月份濕度降雨量10.24.622.33.638.72.1418.52.9524.63.0632.12.7736.82.2837.12.5928.34.31017.83.4116.42.112-3.23.7month=1;12;degree=0.2 2.3 8.7 18.5 24.6 32.1 36.8 37.1 28.3 17.8 6.4-3.2；rain=4.6 3.6 2.1 2.9 3.0 2.7 2.2 2.5 4.3 3.4 2.1

13、3.7;plot(month,degree,+,month,rain,p) month=1:12;degree=0.2 2.3 8.7 18.5 24.6 32.1 36.8 37.1 28.3 17.8 6.4-3.2；rain=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(month,degree,+,month,rain,p)grid onxlabel(month)ylabel(degree&rainfall)text(5,24.6,leftarrowdegree);text(3,2.1,leftarrowrainfall);t

14、itle(Temperature&rainfall);歡迎下載40-5 05.已知矩陣請繪圖表示。丁 degreerainfallr_d卜1( |11Temperature&rainfall35302520151050248106month12A=ones(5,7);A(2:4,2:6)=2;plot(A)1111112221122211 12 22 22 21 112221111116.繪制z x2 y2的三維立體圖 x=-5:0.1:5; y=x; X,Y=meshgrid(x,y); Z=X.A2+Y.A2; surf(X,Y,Z)7.繪制由函數x2/4 y2/9 z2/16 1形成的立

15、體圖，并通過改變觀察點來獲得該圖形在各個坐標平面上的平面投影。theta=(0:0.01:2)*pi;phi=theta;x=zeros(length(theta),length(phi);y=x;z=x;for i=1:length(theta)for j=1:length(phi)x(i,j)=2*sin(theta(i)*cos(phi(j);y(i,j)=3*sin(theta(i)*sin(phi(j);z(i,j)=4*cos(theta(i);endendaxis squaresubplot(2,2,1)mesh(x,y,z)title(三維)subplot(2,2,2)mesh

16、(x,y,z)view(90,0)title(x 軸視圖)subplot(2,2,3)mesh(x,y,z)view(0,0)title(y 軸視圖)subplot(2,2,4)mesh(x,y,z)view(0,90)title(z 軸視圖)第五章1.使用命令文件，畫出下列分段函數所表示的曲線。x 1,x 0y 1,0 x 13x ,1 xx=-2:0.01:2;for ii=1:401if(x(ii)0)y(ii)=x(ii)+1;elseif (x(ii)=0)y(ii)=1；elseif (x(ii)=1)y(ii)=x(ii)A3;endendplot(x,y)2.計算上述分段函數的

17、值，要求能夠根據用戶對x值的不同輸入，程序給出相對應的結果。function y=fenduan(x)%fenduan用于計算分段函數%x為任意實數 if(x0)y=x+1;elseif(x=0)y=1;elsey=xA3;end3.編程求arcsinx1痢52g4g5(2n)!x2n 122n (n!)2(2n 1)K ,其中 |x| 1。提示：結束條件可用|u|e,其中u為通項公式，e為滿足所求精度的極小值。function y=arcsin(x) if abs(x)1disp（輸入的參數的絕對值必須小于1）;return;end n=0;y=0;u=x;while uepsu=(fac(

18、2*n)*xA(2*n+1)/(2A(2*n)*(fac(n)A2*(2*n+1);y=y+u;n=n+1;endfunction y=fac(n)if(n fo門=2:999a=rem(i,10);b=rem(fix(i/10),10);c=fix(i/100);if rem(a+b+c),2)&isprime(i) disp(i)；endend357 23 29 41 43 47 61 6783 89 113 131 137 139 151 157 173 179191193197199223227229241 263 269281 283 311 313 317 331 337 353 359 373379397401409421443449461 463 467487557571577593599601607 641 643647 661 683 719 733 739 751 757 773 797809821 823827829 863881883 887 911919937 953971977 9919976 .編寫一個判斷任意輸入的正整數是否為素數的函數文件。function y=issushu(n)y=1;for ii=2:fix(n/2)if (re