计算机仿真实验.pdf

《计算机仿真实验.pdf》由会员分享，可在线阅读，更多相关《计算机仿真实验.pdf（51页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、 计算机控制技术及仿真 实验指导姓 名：谭 富 强学 号:200810601117班 级：电 子 一班攀枝先老温机电工程学温二00、耳十月实验一 Mat lab环境语法及数学运算(验证性实验)一、实验目的1、掌 握M a t l a b软件使用的基本方法；2、熟 悉M a t l a b的数据表示、基本运算方法；3、熟 悉M a t l a b绘图命令及基本绘图控制。二、实验仪器与软件1.P C机 1台2.M A T L A B 6.X 环境三、实验原理M A T L A B环境是一种为数值计算、数据分析和图形显示服务的交互式的环境。M A T L A B 有 3 种窗口,即：命令窗口(T h

2、 e C o m m a n d Wi n d o w)、n r文件编辑窗口(T h e E d i tWi n d o w)和图形窗口(T h e F i g ur e Wi n d o w),而 S i m ul i n k 另外又有 S i m ul i n k 模型编辑窗口。1 .命令窗口(T h e C o m m a n d Wi n d o w)当M A T L A B启动后，出现的最大的窗口就是命令窗口。用户可以在提示符“”后面输入交互的命令，这些命令就立即被执行。在M A T L A B中，一连串命令可以放置在一个文件中，不必把它们直接在命令窗口内输入。在命令窗口中输入该文件

3、名，这一连串命令就被执行了。因为这样的文件都是以“.m”为 后 缀,所以称为n r文件。2.m-文件编辑窗口(T h e E d i t Wi n d o w)我们可以用m-文件编辑窗口来产生新的m-文件，或者编辑已经存在的m-文件。在M A T L A B主界面上选择菜单“F i l e/Ne w/M-f i l e”就打开了一个新的m-文件编辑窗口；选择菜单“F i l e/O p e n”就可以打开一个已经存在的m-文件，并且可以在这个窗口中编辑这个m-文件。四、实验内容:1、帮助命令使 用h e l p命令，查 找s q r t (开方)函数的使用方法;2、矩阵运算(1)矩阵的乘法已知

4、 A=l 2;3 4：B=5 5;7 8;求 A 2*BZ =1 0 5 1 1 5229 251(2)矩阵除法已知 A=l 2 3;4 5 6;7 8 9;B=l 0 0;0 2 0;0 0 3;A B,A/BC =1.0 e+0 1 6 *-0.4 50 40.9 0 0 7-0.4 50 41.80 1 4-3.6 0 291.80 1 4-1.3 51 12.7 0 22-1.3 51 1D =1.0 0 0 0 1.0 0 0 0 1.0 0 0 04.0 0 0 0 2.50 0 0 2.0 0 0 07.0 0 0 0 4.0 0 0 0 3.0 0 0 0(3)矩阵的转置及共辗

5、转置已知 A=5+i,2-i,l：6*i,4,9-i;求 A.,A C =5.0 0 0 0 +2.0 0 0 0 -1.0 0 0 01.0 0 0 0 il.O O O O i0 +4.0 0 0 09.0 0 0 0 -6.0 0 0 0 i1.0 0 0 0 iD =5.0 0 0 0 -l.O O O O i2.0 0 0 0 +1.O O O O i1.0 0 0 0(4)使用冒号选出指定元素0-4.00009.0000+6.O O O O i1.O O O O i已 知:A=l 2 3;4 5 6;7 8 9 ;求 A中 第 3 列 前 2 个元素；A中所有列第2,3 行的元素;

6、A=l 2 3;4 5 6;7 8 9 B=A(1:2,3)C=A(2:3,:)A =1 2 34 5 67 8 9B36C =4 5 67 8 9(5)方括号口用 m a g i c 函数生成一个4阶魔术矩阵，删除该矩阵的第四列A=m a g i c(4)B=A(:,1:3)A =16231351110897612414151B =16 2 35 11 109 7 64 14 153、多项式(1)求多项式p(x)=x 3+2x+4的根p=l 024 x=r oot s(p)p =10 2 4x =0.5 8 9 80.5 8 9 8-1.17 9 5(2)已知+1.7 4 4 5 i-1.7

7、 4 4 5 iA=l.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4 ,求矩阵A的特征多项式;P=1.0000-6.9 000-7 7.26 00-8 6.1300 6 04.5 5 00求特征多项式中未知数为2 0时的值；A=l.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4 p=p ol y(A)y=p ol y v a l(p,20)A =1.20003.00005.00000.9 0005.00001.7 0005.00006.00003.00009.000001.00001.00002.00003.00004.0000P=1.0000

8、-6.9 000-7 7.26 00-8 6.1300 6 04.5 5 00y =7.27 7 8 e+004把矩阵A作为未知数代入到多项式中；A=l.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4 p=p ol y(A)y=p ol y v a l(p,A)A =1.20003.00005.00000.9 0005.00001.7 0005.00006.00003.00009.000001.00001.00002.00003.00004.0000p =1.0000-6.9 000-7 7.26 00-8 6.1300 6 04.5 5 00y =1.0e+003*0

9、.38 01-0.4 5 4 5-1.9 9 5 10.4 6 01-1.9 9 5 10.209 3-1.9 9 5 1-2.8 8 8 0-0.4 5 4 5-4.8 9 7 80.6 04 60.4 35 30.4 35 30.08 4 0-0.4 5 4 5-1.16 17基本绘图命令(1)绘制余弦曲线y=c os (t),t e 0,2(2)在同一坐标系中绘制余弦曲线y=c os(t-0.25)和正弦曲线y=s i n(t-0.5),t e 0,2 JT 1-100.60.40.20-0.2-0.4-0.6-0.80.812345675、基本绘图控制t绘制 0,4 口 区间上的x l

10、=10s i n t 曲线，并要求：(1)线形为点划线、颜色为红色、数据点标记为加号;t=0:0.2:4*p ix=10*s i n(t)p l ot (t,x,，+-r，)(2)坐标轴控制：显示范围、刻度线、比例、网络线(3)标注控制：坐标轴名称、标题、相应文本；五、实验要求利用所学知识，完成上述各项实验内容，并将实验过程和实验步骤和结果写在报告中。实验二 MATLAB数值运算与绘图（验证性实验）一、实验目的1.熟悉Ma t l a b 中各类数据，尤其是矩阵的定义、赋值和运用。2.了解Ma t l a b 的矩阵分析函数以及求线性方程组的数值解；3.熟悉多项式运算函数、数值插值。二、实验仪

11、器与软件1.PC 机 1 台2.MA T LA B 6.X 环境三、实验原理1.创建矩阵的方法a.直接输入法规则：矩阵元素必须用 括住；矩阵元素必须用逗号或空格分隔；在 内矩阵的行与行之间必须用分号分隔。逗号和分号的作用：逗号和分号可作为指令间的分隔符，m a t l a b 允许多条语句在同一行出现。分号如果出现在指令后，屏幕上将不显示结果。b.用 m a t l a b 函数创建矩阵：空 阵 m a t l a b 允许输入空阵，当一项操作无结果时，返回空阵;r a n d 随机矩阵;e y e 单位矩阵;z e r o s 全部元素都为 0的矩阵;。n e s 全部元素都为1 的矩阵c.

12、矩阵的修改:可用T 键找到所要修改的矩阵，用一键移动到要修改的矩阵元素上即可修改；指令修改：可以用A（*,*）=*来修改。2 .矩阵运算a.矩阵加、减（+,一）运算规则：（1）相加、减的两矩阵必须有相同的行和列两矩阵对应元素相加减。（2）允许参与运算的两矩阵之一是标量。标量与矩阵的所有元素分别进行加减操作。b.矩阵乘（.*,./,.）运算规则：A 矩阵的列数必须等于B 矩阵的行数标量可与任何矩阵相乘。c.矩阵乘方-a n,a p,p aa p a 自乘p次累，对于p的其它值,计算将涉及特征值和特征向量，如果p是矩阵，a 是标量,a-p 使用特征值和特征向量自乘到p次毒；如 a,p 都是矩阵，a

13、 p则无意义。d.多项式运算m at l ab语言把多项式表达成一个行向量，该向量中的元素是按多项式降基排列的。f （x）=an xn+an-1xn-1+.+l o aO可用行向量p=an an-1.al +a0 表示;p o l y -产生特征多项式系数向量e.代数方程组求解m at l ab中有两种除运算左除和右除。四、实验内容1.输入下列向量（矩阵）g =1 2 3 4 ；h =L 4 3 2 1；g =1 2 3 4 h =4 3 2 1g =12 3 4h1 3 2 ：12 .分别执行以下数组点运算 s i =g +h,s 2 =g.*h,s 3=g/h,s 4 =g.2,s 5 =

14、2.hg =1 2 3 4 h =4 3 2 1s i =g +hs 2 =g.*hs 3=g.hs 4 =g.2s 5 =2.hs i =5 5 5 5s 2 =4 6 6 4s 31894s 4 =14916s 5 =16 8423.输入下列特殊矩阵 A=A=eye(10)A=ones(5,10)A=rand(10,15)A=randn(5,10)A=zeros(5,10)A =A =1000000000010000000000100000000001000000000010000000000100000000001000000000010000000000100000000001A =1

15、1111111111111111111111111111111111111111111111111A =Columns 1 through 130.4514 0.6085 0.0841 0.1210 0.2319 0.43980.1370 0.4225 0.2974 0.3759 0.1939 0.62730.0439 0.0158 0.4544 0.4508 0.2393 0.34000.8188 0.8560 0.0492 0.0099 0.9048 0.69910.0272 0.0164 0.4418 0.7159 0.0498 0.31420.4302 0.4902 0.6932 0.

16、4199 0.5692 0.39720.93420.26440.16030.312 7 0.19 0 1 0.35 33 0.8 9 2 8 0.0 7 8 4 0.365 1 0.8 7 2 90.8 9 0 3 0.8 15 9 0.65 0 1 0.7 5 37 0.6318 0.4 1360.0 12 9 0.5 8 69 0.15 36 0.2 7 31 0.64 0 8 0.39 32 0.2 37 90.7 34 9 0.4 60 8 0.9 8 30 0.7 9 39 0.2 34 4 0.65 5 20.38 4 0 0.0 5 7 6 0.67 5 6 0.2 5 4 8

17、0.19 0 9 0.5 9 15 0.64 5 80.68 7 3 0.4 5 7 4 0.5 5 2 7 0.9 2 0 0 0.5 4 8 8 0.8 37 60.68 31 0.367 6 0.69 9 2 0.8 65 6 0.8 4 39 0.119 7 0.9 6690.34 61 0.4 5 0 7 0.4 0 0 1 0.8 4 4 7 0.9 316 0.37 160.0 9 2 8 0.6315 0.7 2 7 5 0.2 32 4 0.17 39 0.0 38 1 0.664 90.1660 0.4 12 2 0.19 8 8 0.367 8 0.335 2 0.4 2

18、 5 30.0 35 3 0.7 17 6 0.4 7 8 4 0.8 0 4 9 0.17 0 8 0.4 5 8 6 0.8 7 0 40.15 5 6 0.9 0 16 0.62 5 2 0.62 0 8 0.65 5 5 0.5 9 4 70.612 4 0.69 2 7 0.5 5 4 8 0.9 0 8 4 0.9 9 4 3 0.8 69 9 0.0 0 9 90.19 11 0.0 0 5 6 0.7 334 0.7 313 0.39 19 0.5 65 7C o l u m n s 14t h r o u g h 150.7 1650.5 1130.7 7 640.4 8 9

19、 30.18 5 90.7 0 0 60.9 8 2 70.8 0 660.7 0 360.4 8 5 00.114 60.664 90.365 40.14 0 00.5 6680.8 2 300.67 390.9 9 9 40.9 6160.0 5 8 9A =-1.0 10 6-0.64 36 0.0 0 0 00.8 9 5 60.5 68 9-0.2 34 00.62 320.2 37 9 0.38 9 9 -0.9 4 9 90.614 5 0.38 0 3-0.317 90.7 310-0.2 5 5 60.118 40.7 9 9 0-1.0 0 7 8 0.0 8 8 0 0.

20、7 8 120.5 0 7 7 -1.0 0 9 1 1.0 9 5 00.5 7 7 9-0.37 7 50.314 80.9 4 0 9-0.7 4 2 0 -0.635 5 0.5 69 01.69 2 4 -0.0 19 5 -1.8 7 4 00.0 4 0 3-0.2 9 5 91.4 4 35-0.9 9 2 11.0 8 2 3-0.5 5 9 6-0.8 2 170.5 9 13-0.0 4 8 2 0.4 2 8 20.67 7 1-1.4 7 5 1-0.35 100.2 12 0-0.1315 0.4 4 37 -0.2 65 6A =M =0 00 000 00000

21、 00 000 00000 00 000 00000 00 000 00000 00 000 00004.输入下列矩阵及矩阵函数A=2 0-1；1 3 2 ;B=l 7 -1;4 2 3;2 01；M =A*B%矩阵A与B按矩阵运算相乘d e t _ B =d e t(B)%矩阵A的行列式r an k _ A =r an k (A)%矩阵A的秩i n v _ B 二i n v (B)%矩阵B的逆矩阵 V,D=e i g(B)%矩阵B的特征值矩阵V与特征向量构成的矩阵DX=A/B%A/B =A*B-1,即 XB=A,求XY =B A%B A =B-1*A,即 B Y=A,求Y0 14 -317

22、13 10d e t _ B =2 0r an k _ A =2-0.7 0 9 40.7 4 4 40.7 4 4 4-0.667 5-0.35 9 9+0.0 2 18 i-0.35 9 9 -0.0 2 18 i-0.2 2 63-0.5 5 8 7-0.0 60 7 i-0.5 5 8 7 +0.0 60 7 iD =7.2 68 0000-1.634 0+0.2 8 61i000-1.634 0 -0.2 8 61iX0.4 0 0 0 -1.4 0 0 0 3.60 0 00.0 0 0 0 1.5 0 0 0 -2.5 0 0 05.多项式运算 p=l 20-56%表示多项式px

23、)=/+2 x，-5 x+6 r r=r o o t s(p)%求多项式p的根 p p=p o l y(r r)%由根的列向量求多项式系数 s=0 0 1 2 3%表不多项式s(x)=X?+2 x+3 c=co n v (p,s)%多项式乘积 d=p o l y d e r(p)%多项式微分 x=l:0.1:2;y=p o l y v al(p,x)%计算多项式的值p=l 20-56 r r=r o o t s(p)p p=p o l y(r r)s=0 01 23 c=co n v (p,s)d=p o l y d e r(p)x=-l:0.1:2;y=p o l y v al (p,x)1

24、20-56r r =-1.8 64 7 +1.35 8 4 i-1.8 64 7 -1.35 8 4 10.8 64 7 +0.6161i0.8 64 7 -0.6161iPP=1.0 0 0 02.0 0 0 00.0 0 0 0 -5.0 0 0 0 6.0 0 0 0s =0 0 12 3-4 -30 -5C o l u m n s 1 t h r o u g h 1310.0 0 0 0 9.69 8 1 9.38 5 6 9.0 5 4 1 8.69 7 6 8.312 5 7.8 9 7 67.4 5 4 1 6.9 8 5 6 6.4 9 8 1 6.0 0 0 0 5.5 0

25、2 1 5.0 17 6C o l u m n s 14 t h r o u g h 2 64.5 62 1 4.15 36 3.8 12 5 3.5 616 3.4 2 61 3.4 336 3.614 14.0 0 0 0 4.62 61 5.5 2 9 6 6.7 5 0 1 8.32 9 6 10.312 5C o l u m n s 2 7 t h r o u g h 3112.7 4 5 6 15.67 8 1 19.1616 2 3.2 5 0 1 2 8.0 0 0 01.0 0 0 02.0 0 0 00.0 0 0 0-5.0 0 0 06.0 0 0 0d =460-5C

26、 o l u m n s 1 t h r o u g h 1310.0 0 0 0 9.69 8 1 9.38 5 6 9.0 5 4 1 8.69 7 6 8.312 5 7.8 9 7 67.4 5 4 1 6.9 8 5 6 6.4 9 8 1 6.0 0 0 0 5.5 0 2 1 5.0 17 6C o l u m n s 14 t h r o u g h 2 64.5 62 1 4.15 36 3.8 12 5 3.5 616 3.4 2 61 3.4 336 3.614 14.0 0 0 0 4.62 61 5.5 2 9 6 6.7 5 0 1 8.32 9 6 10.312

27、5C o l u m n s 2 7 t h r o u g h 3112.7 4 5 6 15.67 8 1 19.1616 2 3.2 5 0 1 2 8.0 0 0 0?1.0 0 0 0 2.0 0 0 0 0.0 0 0 0 -5.0 0 0 0 6.0 0 0 0IE r r o r:M i s s i n g o p e r at o r,co m m a,o r s e m i co l o n.-1.8 64 7-1.8 64 7 -0.8 64 7 +1.35 8 4 i1.35 8 4 i0.6161i0.8 64 7 -0.6161i1.0 0 0 02.0 0 0 0

28、0.0 0 0 0 -5.0 0 0 0 6.0 0 0 000123c0 0 1 4 7 1 -4-3 18d =4 6 0 -5y =C o l u m n s 1 t h r o u g h 1310.0 0 0 0 9.69 8 1 9.38 5 6 9.0 5 4 1 8.69 7 6 8.312 5 7.8 9 7 67.4 5 4 1 6.9 8 5 6 6.4 9 8 1 6.0 0 0 0 5.5 0 2 1 5.0 17 6C o l u m n s 14 t h r o u g h 2 64.5 62 1 4.15 36 3.8 12 5 3.5 616 3.4 2 61

29、 3.4 336 3.614 14.0 0 0 0 4.62 61 5.5 2 9 6 6.7 5 0 1 8.32 9 6 10.312 5C o l u m n s 2 7 t h r o u g h 3112.7 4 5 6 15.67 8 1 19.1616 2 3.2 5 0 1 2 8.0 0 0 06.有理多项式：G(s)10(s +3)(S+1)(S2+S+3)n=co n v(10 ,1 3)d=co n v(1 1,1 1 3)r,p,k 二r e s i d u e (n,d)p l=l-p(l),p 2=l-p 2%p 2(s)=s-p(2)d e n=co n v(p

30、 l,p 2)n u m=co n v(r l,p 2)+co n v(r 2,p l)n u m,d e n 二r e s i d u e (r,p,k)n=co n v(10 ,1 3)d=co n v(l 1,1 1 3)r,p,k =r e s i d u e (n,d)p l=l-p(l)p 2=l-p(2)d e n=co n v (p l,p 2)n u m=co n v(r,p 2)+co n v(r,p l)%定义分子多项式%定义分母多项式%进行部分分式展开定义两个极点多项式p l (s)=s-p(l),%求分母多项式d e n=p l (s)*p 2 (s)%求分子多项式%

31、根据r,p,k的值求有理多项num,den=residue(r,p,k)n=10 30d=12 4 3r=-3.3333-4.0202i-3.3333+4.0202i6.6667P=-0.5000+1.6583i-0.5000-1.6583i-1.0000k=pl=1.5000-1.6583ip2=1.5000+1.6583iden=5.0000?Un d e f i n e d f u n ct i o n o r v ar i abl e r i .E r r o r i n =C:M A T L A B 6p 5 w o r k Un t i t l e d 3.mO n l i n e

32、 7 =n u m=co n v(r l,p 2)+co n v(r 2,p l)%求分子多项式n =10 30d =12 4 3r =-3.3333-4.0 2 0 2 i-3.3333+4.0 2 0 2 i6.6667P =-0.5 0 0 0 +1.65 8 3i-0.5 0 0 0 -1.65 8 3i-1.0 0 0 0k =P l =1.5 0 0 0 -1.65 8 3ip 2 =1.5 0 0 0 +1.65 8 3id e n =5.0 0 0 0n u m-10.0 0 0 0 -12.0 60 5 i-10.0 0 0 0 +12.0 60 5 i2 0.0 0 0 0

33、n u m0.0 0 0 0 10.0 0 0 0 30.0 0 0 0d e n =1.0 0 0 0 2.0 0 0 0 4.0 0 0 0 3.0 0 0 07.函数插值运算(1)线形样条插值 x=0:10 y=s i n(x)x0=3.4 4.7 6.5 8.2 y O=i n t e r p l (x,y,xO)%线形插值 xl=0:0.1:10 y l=s i n(xl)p l o t(xl,y l,r:,x,y,b*,xO,y O,g.)%插值比较0123456789y =C o l u m n s 1 t h r o u g h 90 0.8 4 15 0.9 0 9 3 0.

34、14 11-0.7 5 68 -0.9 5 8 90.65 7 0 0.9 8 9 4-0.2 7 9 4C o l u m n s 10 t h r o u g h 110.4 12 1-0.5 4 4 0 xO =3.4 0 0 0 4.7 0 0 06.5 0 0 08.2 0 0 0y O =-0.2 18 0 -0.8 9 8 30.18 8 80.8 7 39xl =C o l u m n s 1 t h r o u g h 90 0.10 0 00.7 0 0 0 0.8 0 0 00.2 0 0 00.30 0 00.4 0 0 00.5 0 0 00.60 0 0C o l

35、u m n s 10 t h r o u g h 180.9 0 0 0 1.0 0 0 01.60 0 0 1.7 0 0 01.10 0 01.2 0 0 01.30 0 01.4 0 0 01.5 0 0 0C o l u m n s 19 t h r o u g h 2 71.8 0 0 0 1.9 0 0 02.5 0 0 0 2.60 0 02.0 0 0 02.10 0 02.2 0 0 02.30 0 02.4 0 0 0C o l u m n s 2 8 t h r o u g h 362.7 0 0 0 2.8 0 0 03.4 0 0 0 3.5 0 0 02.9 0 0

36、 03.0 0 0 03.10 0 03.2 0 0 03.30 0 0C o l u m n s 37 t h r o u g h 4 53.60 0 0 3.7 0 0 04.30 0 0 4.4 0 0 03.8 0 0 03.9 0 0 04.0 0 0 04.10 0 04.2 0 0 0C o l u m n s 4 6 t h r o u g h 5 44.5 0 0 0 4.60 0 05.2 0 0 0 5.30 0 04.7 0 0 04.8 0 0 04.9 0 0 05.0 0 0 05.10 0 0Columns 55 through 635.4000 5.50005

37、.60005.70005.80005.90006.00006.1000 6.2000Columns 64 through 726.3000 6.40007.0000 7.1000Columns 73 through 816.50006.60006.70006.80006.90007.2000 7.30007.9000 8.0000Columns 82 through 907.40007.50007.60007.70007.80008.1000 8.20008.8000 8.9000Columns 91 through 998.30008.40008.50008.60008.70009.0000

38、 9.10009.7000 9.8000Columns 100 through 1019.9000 10.0000yi=Columns 1 through 99.20009.30009.40009.50009.60000 0.09980.6442 0.7174Columns 10 through 180.19870.29550.38940.47940.56460.7833 0.84150.9996 0.99170.89120.93200.96360.98540.9975Columns 19 through 27-0.2 7 18 -0.36650.9 7 38 0.9 4 630.9 0 9

39、30.8 6320.8 0 8 50.7 4 5 70.67 5 50.5 9 8 5 0.5 15 5C o l u m n s 2 8 t h r o u g h 360.4 2 7 4 0.335 0-0.2 5 5 5 -0.35 0 80.2 39 20.14 110.0 4 16-0.0 5 8 4-0.15 7 7C o l u m n s 37 t h r o u g h 4 5-0.4 4 2 5 -0.5 2 9 8-0.9 162 -0.9 5 16-0.6119-0.68 7 8-0.7 5 68-0.8 18 3-0.8 7 16C o l u m n s 4 6 t

40、 h r o u g h 5 4-0.9 7 7 5 -0.9 9 37-0.8 8 35 -0.8 32 3-0.9 9 9 9-0.9 9 62-0.9 8 2 5-0.9 5 8 9-0.9 2 5 8C o l u m n s 5 5 t h r o u g h 63-0.7 7 2 8 -0.7 0 5 5-0.18 2 2 -0.0 8 31-0.6313-0.5 5 0 7-0.4 64 6-0.37 39-0.2 7 9 4C o l u m n s 64 t h r o u g h 7 20.0 168 0.11650.65 7 0 0.7 2 9 00.2 15 10.31

41、150.4 0 4 80.4 9 4 10.5 7 8 4C o l u m n s 7 3 t h r o u g h 8 10.7 9 37 0.8 5 0 40.9 9 8 9 0.9 8 9 40.8 9 8 70.9 38 00.9 67 90.9 8 8 20.9 9 8 5C o l u m n s 8 2 t h r o u g h 9 00.9 69 9 0.9 4 0 70.5 8 4 9 0.5 0 100.9 0 2 20.8 5 4 60.7 9 8 50.7 34 40.6630C o l u m n s 9 1 t h r o u g h 9 90.4 12 1

42、0.319 10.2 2 2 90.12 4 50.0 2 4 8-0.0 7 5 2-0.17 4 3Columns 100 through 101-0.4575-0.5440 x=01 23 45 67 89y=Columns 1 through 800.84150.90930.1411-0.7568-0.95890.6570Columns 9 through 110.98940.4121-0.5440 xO=3.40004.70006.50008.2000y0=-0.2180-0.89830.18880.8739xl=Columns 1 through 800.10000.20000.3

43、0000.40000.50000.7000Columns 9 through 160.80000.90001.00001.10001.20001.30001.5000-0.27940.60001.4000Columns 17 through 241.60002.30001.70001.80001.90002.00002.10002.2000Columns 25 through 322.40003.10002.50002.60002.70002.80002.90003.0000Columns 33 through 403.20003.90003.30003.40003.50003.60003.7

44、0003.8000Columns 41 through 484.00004.70004.10004.20004.30004.40004.50004.6000Columns 49 through 564.80005.50004.90005.00005.10005.20005.30005.4000Columns 57 through 645.60006.30005.70005.80005.90006.00006.10006.2000Columns 65 through 726.40007.10006.50006.60006.70006.80006.90007.0000Columns 73 thro

45、ugh 807.20007.90007.30007.40007.50007.60007.70007.8000Columns 81 through 888.00008.70008.10008.20008.30008.40008.50008.6000Columns 89 through 968.80008.90009.00009.10009.20009.30009.40009.5000Columns 97 through 1019.60009.70009.80009.900010.0000yi=Columns 1 through 800.09980.19870.29550.38940.47940.

46、56460.6442Columns 9 through 160.71740.78330.84150.89120.93200.96360.98540.9975Columns 17 through 240.99960.99170.97380.94630.90930.86320.80850.7457Columns 25 through 320.67550.59850.51550.42740.33500.23920.14110.0416Columns 33 through 40-0.0584-0.1577-0.2555-0.3508-0.4425-0.5298-0.6119-0.6878Columns

47、 41 through 48-0.7568-0.8183-0.8716-0.9162-0.9516-0.9775-0.9937-0.9999Columns 49 through 56-0.0 7 5 2-0.9 9 62-0.9 8 2 5-0.9 5 8 9-0.9 2 5 8-0.8 8 35-0.8 32 3-0.7 7 2 8-0.7 0 5 5C o l u m n s 5 7 t h r o u g h 64-0.6313-0.5 5 0 7-0.4 64 6-0.37 39-0.2 7 9 4-0.18 2 2-0.0 8 310.0 168C o l u m n s 65 t

48、h r o u g h 7 20.11650.2 15 10.31150.4 0 4 80.4 9 4 10.5 7 8 40.65 7 00.7 2 9 0C o l u m n s 7 3 t h r o u g h 8 00.7 9 370.8 5 0 40.8 9 8 70.9 38 00.9 67 90.9 8 8 20.9 9 8 50.9 9 8 9C o l u m n s 8 1t h r o u g h 8 80.9 8 9 40.9 69 90.9 4 0 70.9 0 2 20.8 5 4 60.7 9 8 50.7 34 40.6630C o l u m n s 8

49、9 t h r o u g h 9 60.5 8 4 90.5 0 100.4 12 10.319 10.2 2 2 90.12 4 50.0 2 4 8C o l u m n s 9 7 t h r o u g h 10 1-0.17 4 3-0.2 7 18 -0.3665 -0.4 5 7 5 -0.5 4 4 00123456789 1010.80.60.40.20-0.2-0.4-0.6-0.8-1五、实验要求利用所学知识，完成上述1至 7 项实验内容，并将实验结果写在实验报告上。六、实验思考题1.矩阵建立与有哪几种方法？a、简单矩阵的输入b、利用语句或函数产生矩阵2.矩阵的加、减、

50、乘、除运算规则是什么？只有同阶矩阵方可进行加减运算，当两个矩阵中前矩阵和后一矩阵的行数相同时，可以进行乘法运算，矩 阵 的 除 法 有 两 种 运 算 符 和“/”，分别表示左除和右除。一般地讲，x二 A B 是 A*x二 B 地解，x=B/A是 x*A二 B地解。实验三 M atlab基本编程方法（设计性实验）一、实验目的1、掌 握 M atlab软件使用的基本方法；2、熟 悉 M atlab程序设计的基本方法二、实验仪器与软件1.PC机 1 台2.MATLAB6.X 环境三、实验原理根据M A T L A B 基本数值计算、数据分析和图形函数的功能，按程序设计数学要求完成对象计算的M A