精讲多练matlab习题

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