数学软件复习-matlab

复习资一、基本概念、命令、语句格式、命令表达式与数学、工程中常用的表达形式十分相似；结构、输入输出及面象等特点的高级语言； 分号隔开，逗号告诉显示结果；分号说明除了这一条命令外还有下一条命令等待输入，这时将不会给出运行的中间结果； 用path命令扩展的搜索路径仅在当 义n个等距离数据点数x1=a,x2,…,xn=b，而不是数据点之间的增量；的元素指定了不的元素；另：的数组，也表现在绝对索引上，是按从点乘符号.*对数组cd施加元素对元素的乘法运算，不带点的乘号表示矩阵乘B.\A表示对应元素相除，得到矩阵(aij/bij)；B./A得到矩阵 数值数组的预定维一般采用zeros函数，元胞数组预定维采用cell函数，结构型数组用struct函数预定维；元胞数组：元胞数组是的一种特殊数据类型，可以将元胞数组看做一种无开，例如a={' ',20;ones(2,3),1:10}，b={[109],'book'}。元胞数组中数据的，可通胞。大括号用于表示元胞的内容，小括号表示指定的元胞，例如a{2,1}(2,3)表示a中第2行第1列的元胞中第2行第3列元胞。中的下列函数用于元胞数组操作： 利用图形方式显示元胞数组cell2mat：将元胞转变成为元胞数组num2cell：将数值数组转变成为元胞数组deal：将输入参数赋值给输出cell2struct：将元胞数组转变成为结构struct2cell： 信息即可设置为结构体类型，其可以包含、学号、、籍贯等字段。='';student.numble=20141898;student.sex='男';student.nationality='中国';studentname:''sex男' 结合使用冒号运算和大括号可以方便地元胞数组中的多个元素，并把的A(:,j)A 函数cond(A)[V，D]=eig(A)det(A)计算行列式inv(A)求矩阵的逆A为方阵，poly(A)A的特征多项式；若A为向量，poly(A)A的元素为全部根的多项式；roots(P)求多项式P的根；函数polyfit 上的M函数文件所调用clear命令从内存中删除变量；clearclearall删除所有的内存变量；clearx1 另一种是函数句柄演算函数feval； 求数aa=[12],b=[234;678],请比numel(a),length(a),numel(b),length(b),size(a),size(b)之间的异同。构成Y，产生一个矩形网格的坐标(X,Y)。矩阵运算基础的第1行和第4行。 中如何建立矩阵 73，并将其赋予变量419 419 >>a=[573;49 运算规则，如矩阵a与b相乘（a*b）时必须满足a的列数等于b的行数。算的运算符前加一个点即为数组运算，如a*b为矩阵乘，a.*b为数组乘。37

67

9

3>>a=[535;374;79>>b=[242;679;83ans= 9

37

29

13

75i4>>x=[4+8i3+5i2-7i1+4i7-5i;3+2i7-6i9+4i3-9i>>ans4.0000- 3.0000-3.0000- 7.0000+2.0000+ 9.0000-1.0000- 3.0000+7.0000 4.0000-7a

9 41 与b 27 46>>a=[693;27>>b=[241;46ans= “7676 35 >>A=[492;764;35>>B=[3726X=10a

256a 78>>a=[123;456;78ans=>>ans11a1 5 874ab ，b 3

2>>a=[123;45>>b=[8–74;36ans=010101>>ans010101>>ans101010>>ans101010>>ans000000>>ans11111112a0a=[11011]在sin(x)运算中，x是弧度 14x30>>x=[3045>>ans=

ans= ans= ans= 15用四舍五入的方法将数组 >> ans= 616将矩阵a

2、b71和c

9 7 8 6组一个43的矩阵，第一列为按列顺序排列的a矩阵元素，第二列为按4756 56 21 3245277813569>>a=[42;5>>b=[71;8>>c=[59;6 >>d=[a(:)b(:)d475586219732 >>e >>e=resh ans= 三、重 close logspacerepmat铺放模块数组信息显示(Messagedisplay) 控制语句(Controlflow) 所用函数的输出变量数目realmin最小可用正实数 最大可用正实数时间和日期(Timeanddates) 时钟 占用CPU时 CLOCK 符串操作(Stringoperations) 取整和求余函数(Roundingandremainder) f 方程求解(SolutionofEquation) 四、M文件编globalM文件中定义全局变量时，如如果x直接键入x以观察其中所包含的及内容，也可以键入x.观察指定的内容，例x=color:data:[12.700035.90008.400012.200016.600018.8000ans=ifififif说明break语句、continuereturn体，继续执行循环体外的下一行语句。在嵌套循环中，breakforn=if

disp(['Divisibleby7:'return语句用于终止当前令序列，并返回到调用的函数或键盘，也用于终keyboad方式，在中，被调用的函数运行结束后会自动返回调用函return语句时将return被调用函数的某一位置，根据某种条件迫使被调用函数提前结束并返0.1（①0.*andn1n0.1的正态分布的随机噪声；②三点线性滑动平均就是依次取每三个相邻数的平均值作为新的数x12=(x(1)+x2)+(3))/31)=(x(2)x3)x(4))/3……）fori=2:n-ya(i)=sum(y(i-functionc=collatz(n)%%Classic“3n+1”Ploblemfromnumbertheorywhilec=[c有一组学生的考试成绩（见表根据规定，成绩在100分时为满分，成绩在90~99之80~8960~7960分王张刘李陈杨于黄郭赵成%划分区域：满分(100)，优秀(90-99)，良好(80-89)，及格(60-79)，不及格(<60)fori=1:nfor case Rank(i,:)='满分case Rank(i,:)='优秀case Rank(i,:)='良好case Rank(i,:)='及格 disp('disp(['学生 得分 等级 for 学生得 王张刘李陈杨于黄郭赵：');ifk~=1&k~=2：');ifk==1 %elseif %str=['C';' ','转换后的温度 ',keyboardkeyboard命令时，是否可以在工作空间浏览器中观察和修改函数中的局部变量？退出keyboard状态令是什么？编一个函数文件，在其中设置keyboard，观察变量情况。键入return命令，退出keyboard状态，控制权交还给系统，程序将继续运行。11.M上，这些变量常用来设置缺省输入变量，并决定用户所希望的输出变量。例如，考虑函数linspace：%%%%%ifnargin==2果用户只用两个输入参量调用linspace，例如linspace(0,10)，linspace产生1003，例如，linspace(0,10,50)，第三个参量决定数五、的符号运算功求 对通项S求和，其中k为变量。且从0变到k-1。 对通项S求和，指定其中v为变量。且v从0变到v-1。 对通项S求和，其中k为变量。且从a变到b。 对通项S求和，指定其中v为变量。且v从a变到b。例：键入k=sym('k');symsum(k)得ans=ans=ans=求导 求表达式S对变量v的n阶导数。 ans 0, ansA 1/x/y,dsdxdy=diff(diff(S,'x'),'y')dsdxdy 1/x^2/y^2,i*exp(i*x*y)-p1=diff(p,'x')p1=y=xf(x2)的导数，可键入命令：p1=diff(p,'x')p1=p='x*y(x)-p1= dy 表达式P中自变量趋于零时的极限。 表达式P中自变量趋于a时的极限。limit(P,x,a,'left')表达式P中自变量x趋于a时的左极限。limit(P,x,a,'right')表达式P中自变量x趋于a时的右极限。limit(P)得limit(P,x,0,'right')limit(P,h,0)得ans=limit(v,x,inf,'left')得ans=[exp(a), 以va,bP例如：键入int('-2*x/(1+x^2)^2')得ans=键入 ans=键入 int('2*x','sin(t)','log(t)')anslog(t)^2-int('[exp(t),exp(a*t)]')ans c=[[Warning:Systemisrankdeficient.Solutionisnotunique.X=[[[[代数 对方程Pv求解。solve(p1,P2,…，Pn,v1,v2,…，vn)P1，P2，…，Pnsolve('p+sin(x)=r')得：ans

P1='x^2+x*y+y=3';P2='x^2- x=[[y [-au2v2再例：要求解方程组uvP1='a+u^2+v^2=0';P2='u- 得u[1/2+1/2*(-1-[1/2-1/2*(-1-v[-1/2+1/2*(-1-[-1/2-1/2*(-1-

[x,y]=solve('sin(x+y)-exp(x)*y=0','x^2- 得x y Warning:3equationsin4variables.u=x1:[1x1sym]x2:[1x1sym]x3:[1x1x4:[1x1 ansans=ans=ans=6．解符号微分方 t。如果输入的初字母y表式函数，Dy表示y对t的一阶导数；Dny表示yt的n

xyx[x,y]=dsolve('Dx=y','Dy=- x=y=例如求,f3f4g,g4f3g,f(0)=0,g(0)=1的解，可输入指令f=gy(0)

y(0)

y=dsolve('D3y=- 得yans= [-六、绘二维图 x轴和y轴均为线性刻度（Linearscale） x轴和y轴均为对数刻度（Logarithmicscale） x轴为对数刻度，y轴为线性刻度 x轴为线性刻度，yplot(x,sin(x),x,plot(x,sin(x),'c',x,cos(x),plot(x,sin(x),'co',x,cos(x),'g*');绘图中用到的直线属性包括：LineStyle：线MarkerTypeMarkerSizeMarkerFaceColor：标记点的填充颜标记--：-线标记+o*.xsd点标记^v><ph点标记rgbcmykw颜红绿蓝黄黑白plot(t,sin(2*t),'-axis([0,6,-1.2,1.2]);xlabel('Input xylabel('Function ytitle('TwoTrigonometric %图形标legend('ysin(x)','y %图形注grid %显示格subplot(2,2,1);plot(x,subplot(2,2,2);plot(x,subplot(2,2,3);plot(x,subplot(2,2,4);plot(x,其他各种二维绘图函数 a=0.1、a=0.2a=0.5ye-at，并用箭头线标识出各曲线a的取值。解答>>>>y1=exp(->>y2=exp(->>y3=exp(->>plot(t,y1,'-ob',t,y2,':*r',t,y3,'->>>>>>三维图plot3——调用格式：plot3(x,y,z)——x,y,z是长度相同的向量plot3(x,y,z,s)——带开关量plot3(x1,y1,z1,'s1',x2,y2,z2,'s2',…)gridon(off) text(x,y,z,‘string’)三维图形标注 三维多边形的绘 调用格式：mesh(z)zn×m的矩阵，x与y[x,y]=meshgrid(x,y)——2维网格mesh(x,y,z)——绘三维网线图，x,y,z分别为三的坐标位 调用格式：surf(x,y,z)——绘制三维曲面图，x,y,zsurfl(X,Y,Z) shadingfaceted——shadingflat——shadinginterp——zxex2y2xy的取值范围均为-22时，用建立子窗口的方法在同一个解答>>[x,y]=meshgrid([->>z=x.*exp(-x.^2->>>> >>title('plot3>> >>title('mesh>> >>title('surf>> surf(x,y,z),shading>>title('surf(x,y,z),shading七、中的运算符和特殊字符说符号符号用途说 help 建矩阵的语句中指示一行元素的结束，例如m=[xyz;ijk] ci行的所有列元素；a(1:3,4)13行元素 创建单元矩阵(cellarray) 八、最例：Findxthatminimizessubject0≦x1,,enterthecoefficients:f=[-5;-4;-6]A [1- b=[20;42;l