高数之MATLAB例题解析【全】

1、数学,MATLAB,2,常用数学函数表,第一讲 MATLAB入门（一）,第二讲 MATLAB入门（二）,第三讲 一元函数的图形,第四讲 极限,第五讲 一元函数微分学,第六讲 一元函数积分学,第七讲 矩阵及其运算,目 录,3,返回,返回目录,常用数学函数表,注：函数的使用格式是函数名（自变量）,4,返回,返回目录,常用数学函数表,注：函数的使用格式是函数名（自变量）,2014-10-7-17:43,实验目的 1.了解MATLAB软件。 2.熟悉MATLAB软件的基本操作。 3.了解MATLAB软件的功能和使用方法。 实验内容 1.Desktop操作桌面的启动 2.MATLAB的启动 3.双击桌面

2、上的MATLAB图标。 4.依次点击开始，程序，MATLAB。 5.Desktop操作桌面简介 6.操作桌面的缺省外貌,返回,返回目录,第一讲 MATLAB入门（一）,5,6,通用操作界面Command Window运行入门Command Window指令窗简介,7,最简单的计算器使用法,【例1】求的算术运算结果。 （1）用键盘在MATLAB指令窗中输入以下内容 (12+2*(7-4)/32 （2）在上述表达式输入完成后，按【Enter】键， 该就指令被执行。 （3）在指令执行后，MATLAB指令窗中将显示以下结果。 ans = 2,（1）在键盘上输入下列内容 A = 1,2,3; 4,5,6

3、; 7,8,9 （2）按【Enter】键，指令被执行。 （3）在指令执行后，MATLAB指令窗中将显示以下结果： A = 1 2 3 4 5 6 7 8 9,8,【例题】指令的续行输入S=1-1/2+1/3-1/4+ .1/5-1/6+1/7-1/8 S = 0.6345,交互界面分类目录窗,该窗口显示MATLAB总包含和已安装的工具箱的帮助， 演示和产品等内容，分别用HELP,DEMOS,MA 来表示，双击相应的图标，则可以看到相关的信息。,9,工作空间浏览器简介,10,Command HistoryCommand History历史指令窗简介,11,历史指令行的再运行,【例4】再运行图1.

4、5-2所示历史指令窗中的三行指令。,12,2014-10-7-17:43,13,14,15,16,练习作业 1.熟悉MATLAB的各窗口功能，并且做简单的使用。 2.熟悉各类运算符和常量与变量的命名与使用。,返回,返回目录,17,实验目的 1.使用MATLAB软件的基本语言解决简单的问题。 2.了解M文件和M函数文件的建立方式。 3.了解程序设计的基本结构，初步了解如何对实际问题建立模型。 4.学习一维数值数组的创建、访问。 实验内容 1.M文件 M文件的建立 打开程序编辑器FILE-NEW-M-FILE。 点击新建NEW M-FILE 按钮。 在命令窗口键入EDIT命令。,返回,返回目录,第

5、二讲 MATLAB入门（二）,18,M文件的编写,【例1】计算y=sin(x)在x=0，x=pi/2，pi, 3pi/2，2pi处的函数值。 操作步骤： 1.打开程序编辑器，编写程序如下： x=0:pi/2:2*pi y=sin(x) 2.保存文件名lx1.。然后再命令窗口键入 lx1 3.命令窗口得到结果，如下 x= 0 1.5708 3.1416 4.7124 6.2832 y = 0 1.0000 0.0000 -1.0000 -0.0000,19,M函数文件,1.M函数文件是一种特殊的文件，格式为function【输出变量列表】=函数名（输入变量列表） 【例2】已知，试建立M函数文件。

6、 操作步骤： 1.打开程序编辑器，创建lx2.m文件如下： function y=lx2(x） y=(cos(x)3+2)0.5 2.命令窗口中输入 y=lx2(2) 3.得到结果为： y = 1.3885 程序结构 顺序结构 顺序结构是指依次的逐条执行程序的结构。,20,分支结构（IF-ELSE-END）,21,【例题】比较两个数a1，a2的大小，将其按从小到大的顺序进行排列。 操作步骤：1.打开程序编辑器，编辑函数文件lx3.m： function a=lx3(a1,a2) if a1 a=lx3(4,5) 3.命令窗口得到结果： a = 4 5,22,for循环格式： for 循环变量=

7、初值：步长：终值 循环体 end,【例4】计算。 1.打开程序编辑器，编辑函数文件 s=0 for n=1:1:100 s=s+n end s 在保存文件，点击运行。 命令窗口得到结果： s = 5050 WHILE循环 格式： while 表达式 循环体 end,循环结构,23,一维数组一维数组的创建,【例题】创建向量x=(1 3 5 7 9) x=1:2:9 x = 1 3 5 7 9,24,格式：x（n ） 功能：访问数组x的第n个元素,【例6】访问向量x=(1 3 5 7 9)的第3，5个元素。 x=1:2:9 x = 1 3 5 7 9 x(3) ans = 5 x(5) ans =

8、 9,一维数组的访问,25,1.编写一个函数文件，求任意给定的3个数的最小值。 打开程序编辑器，编辑函数文件 function a=lx5(a1,a2,a3) if a1a2 a=a2 elseif a1a3 a=a else a=a3 end 在保存文件为lx5,并且命名。 命令窗口得到结果： a=lx5(1,2,3) a = 1,练习作业,26,(1)打开程序编辑器，编辑函数文件 s=0 for n=0:1:50 s=s+sin(n*(pi/50) end S (2)在保存文件,并且点击运行。 （3）命令窗口得到结果： s = 31.8205,27,（1）计算并列示出小于100的Fibon

9、naci数列的各项。 打开程序编辑器，编辑函数文件 a=1; b=1; c=0; while c100 c=a+b a=b; b=c; end 在保存文件,并且点击运行。 命令窗口得到结果：, a = 1 b = 1 c = 0 c = 2 a = 1 b = 2 c = 3 a = 2 b = 3 c = 5,a = 55 b = 89 c = 144 a = 89 b = 144,a = 3 b = 5 c = 8 a = 5 b = 8 c = 13 a = 8 b = 13 c = 21,a = 13 b = 21 c = 34 a = 21 b = 34 c = 55 a = 34

10、b = 55 c = 89,3.Finbonaci 数列是这样一个数列：它的前两项都为1， 第三项是前两项之和，以后各项都是前两项之和。,(2)寻找Finbonaci数列中第一个大于10000的项及项数。 打开程序编辑器，编辑函数文件 a=1 b=1 c=0 cnt =3 while c10000 c=a+b a=b b=c end 在保存文件,并且点击运行。,返回,命令窗口得到结果： a = 1 b = 1 cnt = 3 c = 0,返回目录,2014-10-7-17:43,28,29,返回,返回目录,实验目的1.学习使用MATALAB绘制一元函数图形的方法。2.通过观察图形特征来分析函数

11、的有关性质。,预备知识 一.显函数 设一元函数y=f(x)的定义域为D，动点（x,f(x)）的运动轨迹成为 函数y=f(x)的图形,，它是一条平面曲线。,第三讲 一元函数的图形,二.参数方程,当参数t从a变化到b时，对应点(x,y)的图形是一条平面曲线。,实验内容,二维平面图形的描绘方法:,30,(一) 离散数据作图,31,(二)绘图参数,32,例1绘制散点图。其中y=1 2 5 6 3 0 4. 在编辑区中编辑文件 y=1 2 5 6 3 0 4; figure(1) plot(y,r*) 点击保存，再点击运行。 在另一窗口出现：,（三）作图步骤： （1）取自变量数据向量。 （2）计算因变量

12、向量。 （3）调用命令作图。,33,例2在同一窗口中绘制正，余弦函数的图形。,在编辑区中编辑文件 x=0:0.01:2*pi y1=sin(x) y2=cos(x) plot(x,y1,r+,x,y2,g*) 点击保存，再点击运行。 在另一窗口出现：,34,二维数值函数作图,35,在编辑区中编辑文件 x,y=fplot(sin(1/x),-0.1,0.1); n=length(x); x=linspace(-0.1,0.1,n); y=sin(1./x); subplot(1,2,1) plot(X,Y) subplot(1,2,2) plot(x,y),点击保存，再点击运行。 在另一窗口出现

13、：,36,二维符号函数作图,37,例题1用ezplot命令绘制函数y=xsin(x)的图形。 在编辑区中编辑文件 ezplot(x*sin(x)，-2,2)点击保存，再点击运行。 在另一窗口出现：,解：在编辑区中编辑文件 ezplot(t-sin(t),1-cos(t)，0,2*pi) 点击保存，再点击运行。 在另一窗口出现：,38,39,图形的标准与控制,（一）图形的标注,（二） 图形的控制,40,在编辑区中编辑文件 x1=0:0.01:1; x2=1:0.05:3; y1=2*sqrt(x1); y2=1+x2; plot(x1,y1,r+,x2,y2,b*) text(0.25,1.5,

14、y=x1/2) text(1.45,2.8,y=1+x) xlabel(自变量) ylabel(因变量) title(分段函数作图) 点击保存，再点击运行。 在另一窗口出现：,41,1.试验数据如下表所示，请将三条曲线绘制在同一图形窗口。其中氮肥用红色（+），磷肥用紫色（.），钾肥用绿色（*）。,练习作业,42,解：在编辑区中编辑文件 x=1 2 3 4 5 6 7 8 y1=179 178 179 179 339 339 339 340 y2=146 144 246 246 146 146 246 256 y3=272 472 272 472 272 472 272 470 figure(3

15、) plot(x,y1,r+,x,y2,m.,x,y3,g*) 点击保存，再点击运行。 在另一窗口出现：,43,解：在编辑区中编辑文件 ezplot(x/(1+x2),-6,6) 点击保存，再点击运行。 在另一窗口出现：,解：在编辑区中编辑文件 ezplot(log10(x),0,5) 点击保存，再点击运行。 在另一窗口出现：,2014-10-7-17:43,44,解：在编辑区中编辑文件 ezplot(t*(1-(sin(t),t*cos(t),-2,2) 点击保存，再点击运行。 在另一窗口出现：,2014-10-7-17:43,45,46,返回,返回目录,解：在编辑区中编辑文件 ezplot

16、(t2)/2,1-t,-3,3) 点击保存，再点击运行。 在另一窗口出现：,47,返回,返回目录,第四讲 极限,实验目的: 学习使用MATALAB求极限的方法。 通过图形演示、观察极限过程，加深对极限、无穷小、无穷大等基本概念的认识。,48,49,（四）无穷大,50,（1）作前100项的散点图并观察趋势；（2）求极限。,（1）解：在编辑区里： n=1:1:100; xn=1+1./n; yn=(1+1./n).n; zn=1+(-1).n./n; figure(1) plot(n,xn,r+) figure(2) plot(n,yn,m+) figure(3) plot(n,zn,g+） 在保

17、存，后运行， 在另一窗口出现：,51,52,2求极限： 解：在编辑区里： syms n a1=limit(1+1/n,n,inf) a2=limit(1+1/n)n,n,inf) a3=limit(1+(-1)n/n,n,inf) 命令窗口中的结果： a1 = 1 a2 = exp(1) a3 = 1,注：观察图形得知，当n无限增大的时候，数列无限接近1， 数列无限接近e。,53,解：在编辑区里： syms x t a1=limit(tan(x)/x,x,0); a2=limit(x/(x-1)-1/log(x),x,1); a3=limit(cos(x+t)-cos(t)/t,t,0); a

18、1 a2 a3 命令窗口中的结果： a1 = 1 a2 = 1/2 a3 = -sin(x),54,解：在编辑区里： syms x b1=limit(log(1+x)/x,x,0) x=-0.5:0.05:0.5; y1=log(1+x); y2=x; plot(x,y1,r-,x,y2,g*) 在保存，后运行， 在另一窗口出现：,55,解：在编辑区里： syms x b2=limit(1-cos(x)/x2,x,0) ezplot(1-cos(x),-0.5,0.5) hold on ezplot(x2,-0.5,0.5) hold off 在保存，后运行， 在另一窗口出现：,56,解：在编

19、辑区里： syms x b3=limit(tan(x)-sin(x)/x2,x,0) ezplot(x2,-0.5,0.5) hold on ezplot(tan(x)-sin(x),-0.5,0.5) hold off 在保存，后运行， 在另一窗口出现：,57,58,59,结果说明： （1）i1 =inf，i2 =inf，i3 =inf，i4 =0，i5 =0。 （2）观察图形可知，指数函数比幂函数趋于无穷大的速度快，幂函数比对数函数趋于无穷大的速度快。,60,（1）求极限； 解：在编辑区里： syms n a1=limit(n-1)/(n+1),n,inf) a2=limit(1/n)*(

20、sin(pi/n),n,inf) a3=limit(1/n)*(cos(n*pi)/2),n,inf) 在保存，后运行， 在另一窗口出现： 由于软件的原因，此处的答案省略,练习作业,61,（2）做3幅图，包括数列前100项的散点图。 在编辑区里： n=1:1:100 xn=(n-1)./(n+1) yn=(1./n).*(sin(pi./n) zn=(1./n).*(cos(n.*pi)./2) figure(1) plot(n,xn,r+) figure(2) plot(n,yn,m+) figure(3) plot(n,zn,g+) figure(4) plot(n,xn,r+,n,yn,

21、m+,n,zn,g+) 在保存，后运行，在另一窗口出现：,62,63,64,前三个图合在一起：,65,66,解：在编辑区里： syms x a m n i1=limit(exp(x)-exp(-x)/(sin*(x),x,0) i2=limit(xm-am)/(xn-an),x,a) i3=limit(1+x)/x)(2*x),x,inf) i4=limit(exp(1/(x-1),x,1) i1 i2 i3 i4 在保存，后运行，在另一窗口出现： i1 i2 i3 i4 由于软件的原因，此处的答案省略,67,解：在编辑区里： ezplot(1-x)/(1+x),-2,2) hold on e

22、zplot(1-(x(1/2),-2,2) hold off 在保存，后运行， 在另一窗口出现：,3.比较无穷小:,68,再在编辑区里输入：syms xb1=limit(1-x)/(1+x)/(1-(x(1/2),x,1)在保存，后运行，在另一窗口出现：由于软件的原因，此处的答案省略,解：在编辑区里： ezplot(1-x,-2,2) hold on ezplot(1-(x(1/8),-2,2) hold off,69,在保存，后运行，在另一窗口出现：,再在编辑区里输入： syms x b2=limit(1-x)/(1-(x(1/8),x,1) 在保存，后运行，在另一窗口出现： 由于软件的原因

23、，此处的答案省略,70,解：在编辑区里： ezplot(x,-2,2) hold on ezplot(2*(sinx)2),-2,2) hold off 在保存，后运行， 在另一窗口出现：,再在编辑区里输入： syms x b3=limit(x/(2*(sinx)2),x,0) 在保存，后运行，在另一窗口出现： 由于软件的原因，此处的答案省略,返回目录,返回,71,第五讲 一元函数微分学,返回,返回目录,一.导数,实验目的,1.学习使用MATALAB求一元函数微分的方法。 2.会使用软件描绘函数的曲线，切线，并从几何图形上观察其特点。 3.掌握用软件求一元函数极值的方法。,预备知识,1.显函数

24、的导数,几何意义：曲线在点处的斜率。,2.参数方程的导数,72,微分,极限,73,实验内容,1.符号导数,74,程序编写如下： syms a x f=a2*sin(a*x); a1=diff(f,x,1) a2=diff(f,x,2) a3=diff(f,a,1) 运行结果如下： a1 = a3*cos(a*x) a2 = -a4*sin(a*x) a3 = 2*sin(a*x)*a+a2*cos(a*x)*x 说明：,75,程序编写如下： syms a b t x=a*cos(t); y=b*sin(t); a1=diff(y,t,1) a2=diff(x,t,2) yx=a1a2 运行结果

25、如下： yx = -b*cos(t)/a/sin(t),76,2.一元函数的极值,77,程序编写： ezplot(2*x3-9*x2+12*x-3,-10,10) 保存后，在运行，之后在另一窗口出现：,78,在再编辑区中编辑程序： x1,y1=fminbnd(2*x3-9*x2+12*x-3,1.5,3) x1= 2.0000 y1 = 1.0000 x2,y2=fminbnd(-2*x3+9*x2-12*x+3,0.5,1.5) x2= 1.0000 y2 = -2,79,练习作业,在编辑区里输入： syms x y=log(exp(x)+(1+exp(2*x)1/2) a1=diff(f,

26、x,1) 再点击保存运行。 在命令窗口会出现: (由于软件的问题，此处答案省略。),1.求下列函数的一阶导数。,80,在编辑区里输入： syms x y=exp(-(sin(1/x) a1=diff(f,x,1) 再点击保存运行。 在命令窗口会出现: (由于软件的问题，此处答案省略),81,在编辑区里输入： syms x y=(x(1/x) a1=diff(f,x,1) 再点击保存运行。 在命令窗口会出现: (由于软件的问题，此处答案省略),82,在编辑区里输入： syms t x=2*exp(t) y=exp(-t) a1=diff(y,t,1) a2=diff(x,t,1) yx=a1/a

27、2 再点击保存运行。 在命令窗口会出现: (由于软件的问题，此处答案省略。),83,在编辑区里输入： ezplot(2*x)/(1+x2),-10,10) 再点击保存运行。 在另一窗口会出现:,2.求下列函数的极值,84,再点击保存运行。 在另一窗口会出现:,前面是一个整体函数图像，而这里是为了找到最小值而得到的一部分图像。（要根据整体函数图像来判断最小值的函数图像邻域）在编辑区里输入：ezplot(2*x)/(1+x2),-4,4),85,再在编辑区里输入：x1,y1=fminbnd(2*x)/(1+x2),-2,0)之后在命令窗口出现：x1 = -1.0000y1 = -1.0000,前面

28、介绍的是最小值，而这里是最大值。（为了找到最大值，需先将整体函数图像翻折） 在编辑区里输入： ezplot(2*x)/(1+x2),-4,4) 再点击保存运行。 在另一窗口会出现:,86,再在编辑区里输入： x1,y1=fminbnd(-(2*x)/(1+x2),0,2) 再点击保存运行。在另一窗口会出现:,87,x2 = 1.0000 y2 = -1.0000 所以：极大值f(1)=1, 极小值f(-1)=-1,88,在编辑区里输入： ezplot(x3)/3-2*x2+3*x+1,-10,10) 再点击保存运行。 在另一窗口会出现:,89,前面是一个整体函数图像，而这里是为了找到最小值而得

29、到的一部分图像。（要根据整体函数图像来判断最小值的函数图像邻域）在编辑区里输入：ezplot(x3)/3)-(2*x2)+(3*x)+1,-1,4),再点击保存运行。 在另一窗口会出现:,90,再在编辑区里输入：x1,y1=fminbnd(x3)/3)-(2*x2)+(3*x)+1,2.5,3.5)之后在命令窗口出现：x1 = 3.0000y1 = 1.0000,前面介绍的是最小值，而这里是最大值。（为了找到最大值，需先将整体函数图像翻折） 在编辑区里输入： ezplot(x3)/3)-(2*x2)+(3*x)+1),-1,4) 再点击保存运行。 在另一窗口会出现:,91,再在编辑区里输入：

30、x1,y1=fminbnd(-(x3)/3)-(2*x2)+(3*x)+1),0.5,1.5) 再点击保存运行。在另一窗口会出现:,92,x1 = 1.0000 y1 = -2.3333 极大值f（1）=2.333；极小值f（3）=1,返回,返回目录,93,第六讲 一元函数积分学,返回,返回目录,一.原函数与不定积分 1.原函数存在定理 如果原函数f(x)在区间a,b上连续，则在区间a,b上存在原函数f(x)，使得对所有的a,b上的x都有 。,实验目的 1.学习使用MATALAB求一元函数积分的方法。 2.从几何图形上直观理解定积分的定义。 3.学习使用软件解决定积分的应用问题,2.不定积分的

31、定义 函数f(x)的所有原函数的全体称为不定积分，记作,94,二.定积分,1.定积分的定义：,2.定积分的几何意义：曲边梯形的面积。,95,三.定积分的应用,1.求平面图形的面积：,3.牛顿-莱布尼茨公式 如果函数f(x)是连续函数f(x)在区间a,b上的一个原函数，则,96,一.符号积分,97,程序编写如下： syms x t f=cos(3*x+t); i1=int(f,x); i2=int(f,t); i1 i2 运行结果如下： i1 = 1/3*sin(3*x+t) i2 = sin(3*x+t),98,程序设计如下： syms x t f=cos(t2); I=limit(int(f

32、,0,x)/x,x,0) 运行结果如下： I = 1,99,二.定积分的应用,解：首先，要求出两条曲线的交点；然后，作图观察，确定积分的上下限；最后，积分求面积。 程序设计如下： xc,yc=solve(y=x2/8,y=64/(x2+16) ezplot(x2/8,-10,10) hold on ezplot(64/(x2+16),-10,10) hold off,100,运行结果如下： xc = -4 4 4*i*2(1/2) -4*i*2(1/2) yc = 2 2 -4 -4 s=int(-x2/8+64/(x2+16),x,-4,4) s = -16/3+8*pi,101,练习作业,

33、解：在编辑区里输入： syms x f=log(1+x2) i1=int(f,x) 在命令窗口中出现答案。 （但是由于软件的问题，此处答案省略）,102,解：在编辑区里输入： syms x f=(sin(x)*cos(x)/(1+(sin(x)2) i1=int(f,x) 在命令窗口中出现答案。 （但是由于软件的问题，此处答案省略）,103,解：在编辑区里输入： syms x f=sin(log(x) i1=int(f,x,1,exp(1) 在命令窗口中出现答案。 （但是由于软件的问题，此处答案省略）,104,解：在编辑区里输入： syms x f=(x2)*exp(-2*(x2) i1=in

34、t(f,x,0,inf) 在命令窗口中出现答案。 （但是由于软件的问题，此处答案省略）,105,解：在编辑区里输入： syms x t f=(sin(t)2 a=int(f,t,0,x) b=a/(x-(pi/2) c=limit(b,x,0) 在命令窗口中出现答案。 （但是由于软件的问题，此处答案省略）,106,解：在编辑区里输入： syms x t f=(1+t2)(1/2) a=int(f,t,0,x) b=a/(x2) c=limit(b,x,0) 在命令窗口中出现答案 （但是由于软件的问题，此处答案省略）,107,2.求下列平面区域的面积。,解：在编辑区里输入： ezplot(y=x

35、2,-2,2) hold on ezplot(y=x,-2,2) hold off 在另一窗口出现；,返回,108,再在编辑区内输入： x1,y1=solve(y=x2,y=x) s=int(x-x2,x,0,1) 在命令窗口中出现答案。 （但是由于软件的问题，此处答案省略）,返回,返回目录,109,第七讲 矩阵及其运算,返回,返回目录,称为一个m行n列的矩阵。 二、矩阵的线性运算,一、矩阵的定义,实验目的,1.学习矩阵的创建，使用和保存。 2.熟练掌握矩阵和数组的运算。 3.了解并初步掌握矩阵元素的标识和提取方法。,预备知识,110,二、矩阵的线性运算,1、矩阵的加减法,2、数乘矩阵,3.矩

36、阵的乘法,三.矩阵的逆：,实验内容,一、矩阵的创建,、直接输入法:,111,矩阵以“ ”为首尾标识；行与行之间用“；”来间隔； 每行的元素用“，”来间隔； 矩阵的元素可以是数字或者表达式,基本规则,解：在命令窗口中直接 输入下面的指令即可： A=1,2,3;8,0,5;2,4,6 显示的结果是： A = 1 2 3 8 0 5 2 4 6,、创建特殊的矩阵（见下一页所示）：,112,113,114,115,116,n=4 b1=eye(n) b3=eye(size(b2)-表示生成与 矩阵b1相同列，但是 要多出一行的单位矩阵。,【例题】： 先在编辑区中输入： n=4 b1=eye(n) b3

37、=eye(size(b2), n = 4 b1 = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 b3 = 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0,再保存后点击运行，则在命令窗口中出现：,补充知识,117,二、矩阵的访问与赋值,118,解： （1）取A的第1行，第2,3列元素构成子矩阵A1； 在命令窗口输入： A=2,3,4;5 6 7 A = 2 3 4 5 6 7 A1=A(1,2 3) A1 = 3 4,119,(2)取A的第2行元素构成子矩阵A2； 在命令窗口输入： A2=A(2,:) A2 = 5 6 7,(3)取A的第3列元

38、素构成子矩阵A3； 在命令窗口输入： A3=A(:,3) A3 = 4 7,(4)将A的1,2,3列顺序组成列矩阵A4； 在命令窗口输入： A4 = 2 5 3 6 4 7,120,(5)设行向量s=2,3,5,6,取A的由s确定位置的元素组成的行向量A5； 在命令窗口输入： s=2 3 5 6; A5=A(s) A5 = 5 3 4 7,将A的由s确定位置的元素赋值为B的元素，构成 矩阵A6； 在命令窗口输入： B=zeros(2); A6=A;A6(s)=B;,(7)将A的第1,3列的元素重新赋值为B的元素构成矩阵A7； 在命令窗口输入： A7=A;A7(:,1 3)=B A7 = 0 3

39、 0 0 6 0,(8)构造矩阵A8,其左上角是矩阵A;(3,3)元为1,(3,1),(3,2)元为0； 在命令窗口输入： A8=A;A8(3,3)=1 A8 = 2 3 4 5 6 7 0 0 1,121,(9)删除矩阵A的第2列，组成矩阵A9。 在命令窗口输入：, A9=A;A9(:,2)= A9 = 2 4 5 7,122,例3已知向量A=1 2 3 4 5 6 7 8 9,解：（1）将A的元素排成3行3列的矩阵B； 在命令窗口输入： A=1:9 A = 1 2 3 4 5 6 7 8 9 （1）B=reshape(A,3,3) B = 1 4 7 2 5 8 3 6 9,（2）将B的第一行与第三行互换生成矩阵C； 在命令窗口输入：