搜索
第一章实验内容1.2.例1-1x=[0:0.5:360]*pi/180;plot(x,sin(x),x,cos(x));例1-2p=[3,7,9,0,-23];x=roots(p)x=-1.8857-0.7604+1.7916i-0.7604-1.7916i1.0732例1-3quad('x.*log(1+x)',0,1)ans=0.2500例1-4a=[2,-3,1;8,3,2;45,1,-9];b=[4;2;17];x=inv(a)*bx=0.4784-0.87930.40543.（1）helpinvinvMatrixinverse.inv(X)istheinverseofthesquarematrixX.AwarningmessageisprintedifXisbadlyscaledornearlysingular…（2）helpplotplotLinearplot.plot(X,Y)plotsvectorYversusvectorX.IfXorYisamatrix,thenthevectorisplottedversustherowsorcolumnsofthematrix,whicheverlineup.IfXisascalarandYisavector,disconnectedlineobjectsarecreatedandplottedasdiscretepointsverticallyatX…（3）helpmaxmaxLargestcomponent.Forvectors,max(X)isthelargestelementinX.Formatrices,max(X)isarowvectorcontainingthemaximumelementfromeachcolumn.ForN-Darrays,max(X)operatesalongthefirstnon-singletondimension...（4）helproundroundRoundtowardsnearestinteger.round(X)roundstheelementsofXtothenearestintegers.4.（1）x=0:pi/10:2*pi;y=sin(x)y=Columns1through1000.30900.58780.80900.95111.00000.95110.80900.58780.3090Columns11through200.0000-0.3090-0.5878-0.8090-0.9511-1.0000-0.9511-0.8090-0.5878-0.3090Column21-0.0000（2）根据提示操作，选择y变量，并绘制y图形，生成了如下图形，图形是一个正弦函数图像。思考练习1.与其他窗口程序的启动与退出类似，双击图标打开，单击X图标关闭2.MATLAB是矩阵实验室（MatrixLaboratory）的简称用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案。3.如果一个命令行很长，需要分成多行输入，可以在第一个物理行之后加上三个小黑点…并按下enter键，然后接着在下一个物理航继续输入命令的其他部分。三个小黑点称为续行符，即把下面的物理行看作该行的逻辑继续。4.help命令只搜索出那些与关键字完全匹配的结果，lookfor命令对搜索范围内的小文件进行关键字搜索，条件比较宽松。5.假如建立了一个变量fac，同时在当前目录下建立了一个M文件fac.m文件。如果在命令窗口输入fac,按照MATLAB的搜索顺序，应该在屏幕上显示fac的值。如果没有建立fac变量，则执行fac.m文件。第二章实验指导1.（1）w=sqrt(2)*(1+0.34245*10^-6)w=1.4142（2）a=3.5;b=5;c=-9.8x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2))/(tan(b+c)+a)x=0.9829(3)a=3.32;b=-7.9;y=2*pi*a^2*((1-pi/4)*b-(0.8333-pi/4)*a)y=-128.4271(4)t=[2,1-3*i;5,-0.65]t=2.00001.0000-3.0000i5.0000-0.6500z=0.5*exp(2*t)*log(t+sqrt(1+t^2))z=1.0e+04*0.0057-0.0007i0.0049-0.0027i1.9884-0.3696i1.7706-1.0539i2.A=[-1,5,-4;0,7,8;3,61,7];B=[8,3,-1;2,5,3;-3,2,0];(1)A+6*Bans=4723-10123726-15737A^2-B+eyeans=-18-216182353311022868526(2)A*Bans=141416-105121125328180A.*Bans=-815403524-91220B*Aans=-110-157228533-128(3)A/Bans=1.2234-0.92552.9787-0.94682.3511-0.95744.61703.872313.8936B\Aans=-0.5106-8.6170-1.12770.734017.57451.8085-0.8830-21.21280.4043(4)[A,B]ans=-15-483-10782533617-320[A([1,3],:);B^2]ans=-15-4361773371173713-20193.A=[23,10,-0.778,0;41,-45,65,5;32,5,0,32;6,-9,54,3.14]23.000010.0000-0.7780041.0000-45.000065.00005.000032.00005.0000032.00006.0000-9.000054.00003.1400(1)k=find(10=A&A=25)k=15A(k)ans=23.000010.0000(2)B=A(1:3,:)B=23.000010.0000-0.7780041.0000-45.000065.00005.000032.00005.0000032.0000C=A(:,1:2)C=231041-453256-9D=A(2:4,3:4)D=65.00005.0000032.000054.00003.1400E=B*CE=1.0e+03*0.9141-0.22391.20802.71501.1330-0.1930(3)EDans=010001E&Dans=110111E|Dans=111111~Eans=000000~Dans=0010004.H=hilb(5);P=pascal(5);Hh=det(H)Hh=3.7493e-012Hp=det(P)Hp=1Th=cond(H)Th=4.7661e+005Tp=cond(P)Tp=8.5175e+003条件数越趋近于1，矩阵的性能越好，所以帕斯卡矩阵性能更好。5.A=[-29,6,18;20,5,12;-8,8,5]A=-2961820512-885[V,D]=eig(A)V=0.71300.28030.2733-0.6084-0.78670.87250.34870.55010.4050D=-25.3169000-10.518200016.8351V为A的特征向量，D为A的特征值。它们之间满足A*V=V*D思考练习1.在MATLAB中6+7i是一个复数常量，6+7*i则是一个表达式。i是虚数单位，而I是单位向量。2.A.*B表示A和B单个元素之间对应相乘，A*B是按矩阵乘法的规则乘。A./B表示A矩阵除以B矩阵的对应元素，B.\A等价于A./B。A/B表示A*inv（B）而B\A表示inv(B)*A，对于矩阵运算一般两式不等。3.(1)A(7)=[](2)abs（x）(3)reshape(x,3,4)(4)k=find(t==0);t(k)=0(5)ones(size(A))(6)diag(A)B=diag(diag(A))4.3+randn(25,20)5.(1)A=[1,-1,2,-3;5,1,-4,2;3,0,5,2;11,15,0,9]A=1-12-351-423052111509diag(A)ans=1159triu(A)ans=1-12-301-4200520009tril(A)ans=100051003050111509inv(A)ans=0.13080.13950.0593-0.00060.0785-0.1163-0.12440.07970.0378-0.09300.11050.0087-0.29070.02330.1349-0.0209det(A)ans=-1.7200e+03rank(A)ans=4norm(A,1)ans=20norm(A,2)ans=21.3196norm(A,inf)ans=35cond(A,1)ans=10.7558cond(A)ans=7.6905cond(A,inf)ans=16.4419trace(A)ans=16(2)B=[0.43,43,2;-8.9,4,21]B=0.430043.00002.0000-8.90004.000021.0000diag(B)ans=0.43004.0000triu(B)ans=0.430043.00002.000004.000021.0000tril(B)ans=0.430000-8.90004.00000inv(B)ErrorusinginvMatrixmustbesquare.det(B)ErrorusingdetMatrixmustbesquare.rank(B)ans=2norm(B,1)ans=47norm(B,2)ans=43.4271norm(B,inf)ans=45.4300cond(B,1)Errorusingcond(line35)Aisrectangular.Usethe2norm.cond(B,2)ans=1.9354cond(B,inf)Errorusingcond(line35)Aisrectangular.Usethe2norm.lo9trace(B)Errorusingtrace(line13)Matrixmustbesquare.6.A=[34,NaN,Inf,-Inf,-pi,eps,0];all(A)ans=0any(A)ans=1isnan(A)ans=0100000isinf(A)ans=0011000isfinite(A)ans=10001117.结构矩阵，结构矩阵名.成员名=表达式，其中表达式应理解为矩阵表达式建立单元矩阵和一般矩阵相似，只是矩阵元素用大括号括起来。8．用稀疏存储方式可以可以节省内存空间;运算规则上，稀疏矩阵与普通矩阵一样，只是矩阵的存储方式不同。在运算过程中，稀疏存储矩阵可以直接参与运算。当参与运算的对象不全是稀疏存储矩阵时，所得结果一般是完全存储形式。第三章实验指导1.n=input('\n请输入一个任意三位数：');ifn100||n=1000||fix(n)~=nerror('输入不合要求，请输入三位整数。')enddig=@(k)rem(fix(n/10^k),10);n1=arrayfun(dig,0:2);fprintf('反向结果是：%i%i%i\n',n1)请输入一个任意三位数：369反向结果是：9632.(1)用if语句实现x=input('输入一个百分制成绩');ifx=90;disp('A');e