精講多練matlab習(xí)題

1、第一章-4.0計(jì)算:2.設(shè) 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 .設(shè)x 45 ,計(jì)算sin(x) 、355 72 x=pi/180*45; (sin(x)+sqrt(35)/72八(1/5) ans = 2.81584 .設(shè) a 5.67,b 7.811,計(jì)算: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.計(jì)算 y Jx 6(x -) (x 3.2) /(x 7.7)在 x 3 時(shí)的值。 x x=3; y=sqrt(x)-6*(x+1/x)+(x-3.2)八2/(x+7.7)八3y =-18.26796.已知圓的半徑為15,求其直徑，周長(zhǎng)及面積。 r=15;d=2*r; c=pi*2*rc =94.2478 s=pi*r*rs =706.85837已知某三角形的三個(gè)邊的邊長(zhǎng)為8.5,14.6和18.4,求該三角形的面積。提示：area Js(s a)(s b)(s c)其中：s

3、 (a b c)/2, a、b、C分別為三角形三邊邊長(zhǎng)。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 .設(shè)矩陣求：(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.設(shè)三階矩陣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.設(shè)(2E_ 1 TC 1B)AT其中E是4階單位矩陣,AT是4階矩陣A的轉(zhuǎn)置。其中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.設(shè)二階矩陣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 .求多項(xiàng)式3x6 12x5 4x4 7x3 8x 1被(x 3)(x3 5x)除后的結(jié)果。 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，請(qǐng)繪圖表示。 A=1 2

11、4 0 5 10 11 21 3 1; plot(A)123456789102.繪制單位圓。提示：使用axis( square)命令保證圖形的縱橫坐標(biāo)刻度比例 相同。 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 ,假設(shè)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、某地區(qū)一年中每月的平均氣溫和平均降雨量如下表 4 - 5所示，請(qǐng) 畫出其圖形，要求標(biāo)注出坐標(biāo)軸、數(shù)據(jù)點(diǎn)位置、數(shù)據(jù)點(diǎn)大小等。表4 - 5 溫度-降雨量數(shù)據(jù)月份濕度降雨量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.已知矩陣請(qǐng)繪圖表示。丁 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.繪制由函數(shù)x2/4 y2/9 z2/16 1形成的立

15、體圖，并通過改變觀察點(diǎn)來獲得該圖形在各個(gè)坐標(biāo)平面上的平面投影。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.使用命令文件，畫出下列分段函數(shù)所表示的曲線。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.計(jì)算上述分段函數(shù)的

17、值，要求能夠根據(jù)用戶對(duì)x值的不同輸入，程序給出相對(duì)應(yīng)的結(jié)果。function y=fenduan(x)%fenduan用于計(jì)算分段函數(shù)%x為任意實(shí)數(shù) 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。提示：結(jié)束條件可用|u|e,其中u為通項(xiàng)公式，e為滿足所求精度的極小值。function y=arcsin(x) if abs(x)1disp（輸入的參數(shù)的絕對(duì)值必須小于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 .編寫一個(gè)判斷任意輸入的正整數(shù)是否為素?cái)?shù)的函數(shù)文件。function y=issushu(n)y=1;for ii=2:fix(n/2)if (re