精講多練matlab習(xí)題

1、第一章2.設(shè)計(jì)算： 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)/DT =0.4112 x=pi/180*45; (sin(x)+sqrt(35)/72(1/5)ans =2.81583.設(shè)，計(jì)算 x=pi/180*45; (sin(x)+sqrt(35)/72(1/5)ans =2.81584.設(shè)計(jì)算： a=5.67;b=7.811; exp(a+b)/log10(a+b)ans = 6.3351e+0055.計(jì)算在時的值。 x=3; y=sqrt(x)-6*(x+1/x)+(x-3.2)2/(x+7.7)3y =-

2、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，求該三角形的面積。提示： 其中：、分別為三角形三邊邊長。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.6106第二章1.設(shè)矩陣 求： A=3 1 1;2 1 2;1 2 3;B=1 1 -1;2 -1 0;1 -1 1;（1） 2*A+Bans =7 3 1 6 1 4 3 3 7

3、（2） 4*A*A-3*B*Bans =42 21 38 40 19 46 40 33 56 （3） A*Bans = 6 1 -2 6 -1 0 8 -4 2（4） B*Aans = 4 0 0 4 1 0 2 2 2（5） A*B-B*Aans = 2 1 -2 2 -2 0 6 -6 02.設(shè)三階矩陣，滿足其中 求矩陣。 A=1/3 0 0; 0 1/4 0;0 0 1/7; B=6*inv(inv(A)-eye(3)B =3 0 0 0 2 0 0 0 13.設(shè)，其中是4階單位矩陣，是4階矩陣的轉(zhuǎn)置。其中 求矩陣。 B=1 2 -3 -2;0 1 2 -3;0 0 1 2;0 0 0

4、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)*inv(C)A =1 0 0 0 -2 1 0 0 1 -2 1 0 0 1 -2 14.設(shè)二階矩陣、，滿足，其中 求矩陣。 A=2 -1;-1 2;B=0 2;-2 0; X=(2*A+B)/2X = 2 0-2 25.求解線性方程組 A=2 -3 0 2;1 5 2 1;3 -1 1 -1;4 1 2 2; B=8;2;7;12; x=ABx =3.0000 0.0000 -1.00001.00006.求解一元六次方程的根。 A=3 12 4 7 0 8 1; x

5、=roots(A)x = -3.8230 -0.5275 + 0.8497i -0.5275 - 0.8497i 0.5007 + 0.6749i 0.5007 - 0.6749i -0.1234 7.求多項(xiàng)式被除后的結(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 21 52rest = 0 0 0 103 55 788 1第三章1. 求極限x=sym(x); limit(cos(sqrt(x)(pi/x),x,0,right) ans =exp(-1/2*pi)2

6、. 求極限x=sym(x); f=sym(3*sin(x)+x2*cos(1/x)/(1+cos(x)*log(1+x); limit(f,x,0) ans =3/23. 求極限x=sym(x); f=sym(sqrt(4*x2+x-1)+x+1)/sqrt(x2+sin(x); limit(f,x,-inf) ans =14. 求極限syms x y; f=sym(x2+y2)(x2*y2); limit(limit(f,x,0),y,0) ans =15.已知求yx=sym(x); y=sym(tan(sqrt(x+sqrt(x+sqrt(2*x)2); diff(y,x) ans =ta

7、n(x+(x+2(1/2)*x(1/2)(1/2)(1/2)*(1+tan(x+(x+2(1/2)*x(1/2)(1/2)(1/2)2)/(x+(x+2(1/2)*x(1/2)(1/2)(1/2)*(1+1/2/(x+2(1/2)*x(1/2)(1/2)*(1+1/2*2(1/2)/x(1/2)6. 已知求yx=sym(x); y=sym(cos(x2)*sin(1/x)2); diff(y,x) ans =-2*sin(x2)*x*sin(1/x)2-2*cos(x2)*sin(1/x)*cos(1/x)/x27.求積分x=sym(x); y=sym(sqrt(sin(x)-sin(x)3)

8、; 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)(1/2)*(x-1)*(atan(1/(-1+x2)(1/2)-log(x+(-1+x2)(1/2)/(1+x)*(x-1)(1/2)9.求下列微分方程的通解dsolve(D2y+4*Dy+4*y=exp(-2*x),x) ans = 1/2*exp(-2*x)*(2*C2+2*x*C1+x2)10. 求微分方程dsolve(x2*Dy+x*y=y2,y(1)=1) ans =x*exp(1/x)/(exp(1

9、/x)+exp(1/x*t)*x-exp(1/x*t)第四章1.已知向量,請繪圖表示。 A=1 2 4 0 5 10 11 21 3 1; plot(A)2.繪制單位圓。 提示：使用axis(square)命令保證圖形的縱橫坐標(biāo)刻度比例相同。 ezplot(x2+y2=1) axis square3.繪制伏安特性曲線：，假設(shè)分別為1,5,10和20。R=1 5 10 20;I=0:0.1:10;U=I*R;plot(I,U)legend R=1 R=5 R=10 R=20grid on4.某地區(qū)一年中每月的平均氣溫和平均降雨量如下表4 - 5所示，請畫出其圖形，要求標(biāo)注出坐標(biāo)軸、數(shù)據(jù)點(diǎn)位置、數(shù)

10、據(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 3.7;plot(month,degree,+,month,rain,p) month=1:12;degree

11、=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);title(Temperature&rainfall);5.已知矩陣 請繪圖表示。A=ones(5,7);A(

12、2:4,2:6)=2;plot(A)6.繪制的三維立體圖。 x=-5:0.1:5; y=x; X,Y=meshgrid(x,y); Z=X.2+Y.2; surf(X,Y,Z)7.繪制由函數(shù)形成的立體圖，并通過改變觀察點(diǎn)來獲得該圖形在各個坐標(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)*

13、sin(phi(j); z(i,j)=4*cos(theta(i); end endaxis squaresubplot(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,4)mesh(x,y,z)view(0,90)title(z軸視圖) 第五章1.使用命令文件，畫出下列分段函數(shù)所表示的曲線。x=-2:0.01:2;for ii=1:401 if(x(ii)0) y(ii)=x(ii)+

14、1; elseif (x(ii)=0) y(ii)=1; elseif (x(ii)=1) y(ii)=x(ii)3; end endplot(x,y)2.計(jì)算上述分段函數(shù)的值，要求能夠根據(jù)用戶對值的不同輸入，程序給出相對應(yīng)的結(jié)果。function y=fenduan(x)%fenduan用于計(jì)算分段函數(shù)%x為任意實(shí)數(shù)if(x0) y=x+1;elseif(x=0) y=1;else y=x3;end3.編程求，其中。 提示：結(jié)束條件可用，其中為通項(xiàng)公式，為滿足所求精度的極小值。function y=arcsin(x)if abs(x)1 disp(輸入的參數(shù)的絕對值必須小于1); retur

15、n;endn=0;y=0;u=x;while ueps u=(fac(2*n)*x(2*n+1)/(2(2*n)*(fac(n)2*(2*n+1); y=y+u; n=n+1;endfunction y=fac(n)if(n for i=2:999 a=rem(i,10); b=rem(fix(i/10),10); c=fix(i/100); if rem(a+b+c),2)&isprime(i) disp(i); end end3 5 7 23 29 41 43 47 61 67 83 89 113 131 137 139 151 157 173 179 191 193 197 199 22

16、3 227 229 241 263 269 281 283 311 313 317 331 337 353 359 373 379 397 401 409 421 443 449 461 463 467 487 557 571 577 593 599 601 607 641 643 647 661 683 719 733 739 751 757 773 797 809 821 823 827 829 863 881 883 887 911 919 937 953 971 977 991 9976.編寫一個判斷任意輸入的正整數(shù)是否為素數(shù)的函數(shù)文件。function y=issushu(n)y=1

17、;for ii=2:fix(n/2) if (rem(n,ii)=0) y=0; endend7.編寫階乘函數(shù)，并調(diào)用該函數(shù)來生成一個的階乘表，如下圖所示：1122364245120672075040840320936288010function y=fac(n)if(nif n=0|n=1 y=1; else y=n*fac(n); endfor a=1:10 b=fac(a); A(a)=b; endB=1:10; a=A; b=B; c=b ac =11223642451206720750408403209362880103628800第六章3.利用Simulink仿真兩個正弦信號相乘，即計(jì)算。