2023年MATLAB上机实验实验报告.pdf

《2023年MATLAB上机实验实验报告.pdf》由会员分享，可在线阅读，更多相关《2023年MATLAB上机实验实验报告.pdf（97页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、1人11_人8上机实验一一实验目的初步熟悉M A T L A B工作环境，熟悉命令窗口，学会使用帮助窗口查找帮助信息。?X命令窗口二 实验内容(1)熟悉M A T L A B平台的工作环境。(2)熟悉M ATLAB 的5个工作窗口。commandVMndowJ Command History|Current Directory|Workspace|Help|Profiler.(3)M A T L A B的优先搜索顺序。三 实验环节1.熟悉M A T L A B的5 个基本窗口 Comm and Window(命令窗口)X W orkspace(工作空间窗口)s 田田田田国TValue1 lik

2、e M A 1 2;2 1 3 5;7 1 1x13 d o w =l i n s p ac e(1 ,1 0,1 0)；t l=o n e s(3 ),t 2=o n e s(1,3),t 3=o n e s (3,1)t 4 =o n e s (3),t 4 =e y e (4)X =1zl=12 3 4 5 6 7 89 1 0t 11 1 11 1 111t2=111t 3=111t4=111111 1 110 0 0010 000 1 00001思考题:变量如何声明，变量名须遵守什么规则、是否区分大小写。答：（1）变量声明1.局部变量每个函数都有自己的局部变量，这些变量只能在定义它的

3、函数内部使用。当函数运营时，局部变量保存在函数的工作空间中，一旦函数退出，这些局部变量将不复存在。脚本(没有输入输出参数，由一系列MATLAB命令组成的M文献)没有单独的工作空间，只能共享调用者的工作空间。当从命令行调用，脚本变量存在基本工作空间中；当从函数调用，脚本变量存在函数空间中。2.全局变量在函数或基本工作空间内，用g loba I声明的变量为全局变量。例如声明a 为全局变量：g I o ba I a声明了全局变量的函数或基本工作空间，共享该全局变量，都可以给它曲赋值。假如函数的子函数也要使用全局变量，也必须用gl oba I 声明。3.永久变量永久变量用pe r s i s t en

4、t声明，只能在M文献函数中定义和使用，只允许声明它的函数存取。当声明它的函数退出时,MATLAB不会从内存中清除它，例如声明a 为永久变量：pers i st e nt a(2)变量命名规则如下：始于字母，由字母、数字或下划线组成；区分大小写；可任意长，但使用前N个字符。N与硬件有关，由函数name I e n g thm a x 返回，一般 N=63；不能使用关键字作为变量名(关键字在后面给出)；避免使用函数名作为变量名。假如变量采用函数名，该函数失效。(3)区分大小写 试说明分号、逗号、冒号的用法。答：分号：加上分号“；”其作用是将计算机结果存入内存,但不显示在屏幕上，反之，语句结尾若不加

5、“；”，则表达在语句执行后，在将计算结果存入内存的同时，还将运算结果显示出来。逗号：分割列表冒号：从什么到什么，比如1:10意思是 1到10li n space()称 为“线性等分”函数，说明它的用法。可使用help命令，格式如下：hel p I i nspace help 1inspaceLINSPACE Linearly spaced vector.LINSPACE(XI,X2)generates a row vector of 100 linearlyequally spaced points between XI and X2.LINSPACE(XI,K2,N)generates N

6、points between XI and X2.For N 2,LINSPACE returns X2.See also logspace,:.Overloaded functions or methods(ones with the same name in chelp cgualue/linspace.mReference page in Help browserdoc linspace 说 明 函 数one s()、zeros。、e y e()的用法。答；(1 )o n e s ()函数:所有元素都为1的常数矩阵；(2)ze ros()函数：所有元素都为0的矩阵；(3)e ye()函数

7、：单位矩阵；(4)lins p a c e 0函数：如a=linspa c e(nl,n 2,n3),表达在线性空 间上,行矢量的值从n l到n2(2)工作空间窗口(Workspace)。单击工作空间窗口右上角的按钮，将其从MATLAB主界面分离出来。low 在工作空间查看各个变量，或在命令窗口用who,wh os（注意大小写）查看各个变量。Woit space您 Array Editor-11.XE的 电 包 昌 St3Mt5w1 1 1;11 1;11 1 1 W0 0 0;0 0 0;01 2 3 4 5 6 7|C urre nt D ire c tory workspa c e I?

8、xy=l 2 3;4 5 67 8 9zl=l:1 0 z2=1:2:5w=linspa c e (1,1 0,1 0)tl=one s(3)t2=one s(1,3)1 3=one s 1)t4=ze ros(3)t5=e ye (4)11 I 12 13 14567_89101112132 3 41_ 1 1 1 ，X inCommand WindowX1000010000100001爸 日 期 时 昌X 在工作空间双击变量，弹出Array Edi t o r 窗口（数组编辑器窗口），即可修改变量。：s：;Array Editor-11123451111l21113111456789101

9、11213 J 使用s a v e 命令把工作空间的所有变量保存为my_v a r.m at文献。s a ve my_ v ar.matFile Edit View Graphics Debug Desktop Window Help0方 X 电 OJGC K Y Current Directory:F:皿 ATLAB7回Shortcuts How to Add 回 WLSave to MAT-File:Workspace国 彦/昌国.保存在 C D ：O M A T L A B 7Name12345w*tttt41V田田田田田田Tvalue1 1 1；1 1 11；1；10 0 04x4 1

10、2 3R e c e nt够桌面我的文档|Current Directory|Workspace iCommand Historyy=l 2 34 5 67 8 9zl=l:10z2=1:2:5w=linspace(I,10,10tl=ones(3)t2=ones(I,3)t3=ones(3,1)t4=zeros(3)t5=eye(4)我的电脑网上邻居Lxterheb声srs i m url itDozbbdehjajaknr tssst cuw0-曰00050505000w o rk文件名（S）：保存类型注）:0 10 00 0my_ va r.ma t|M A T-f ile s(*.ma

11、 t)0 01 00 1输入下列命令：c 1 ear all%清除工作空间的所有变量Work 3f，XCommand Window福 Name,|Current DirectoryWorkspace0 0 00 0 0 t5=eye(4)t5 =1 0 00 1 00 0 10 0 0 save my_var.mat clea r a ll?x7 8 9-z l=l:1 0 z2=1:2:5w=linspace(1,10,10)tl=o n e s(3)t2=ones(I,3)13=ones(3,1)t4=zeros(3)t5=eye(4)save my_var.matcle a r a ll

12、n J_III观测工作空间的变量是否被清空。使用lo a d 命令把刚才保存的变量载入工作空间。load my_var.matinuriuuib；rj nuvv LU AUU TJ vvria iy newCommand Window 福 国 施 强 昌 国,Name、Value12345t*1tttV1 2 3 4 56 7Current Directory workspace I1 1 1;1 13 1;1 1 11；1；10 0 0;0 0 0;0?xzl=l:10 L Tz2=1:2:5w=linspace(13 10,10)tl=ones(3)t2=ones(1,3)13=ones(

13、3,1)t4=zeros(3)t5=eye(4)save my_var.matclear all t4=zeros(3)t4000000000t5=eye(4)t5=I0000I0000I0000I save my_var.mat clear all load my_var.mat清除命令窗口命令：c 1 c(3)历史命令窗口(C o mma n d His t o r y)。打开历史命令窗口，可以看到每次运营MATL A B 的时间和曾在命令窗口输入过的命令，练习以下几种运用历史命令窗口反复执行输入过的命令的方法。Command Historyxz2=l:2:5 .w=linsp a c e

14、(I,10,10)tl=o nes(3)t2=o nes(I,3)t3=o nes(3 j 1)t4=zero s(3)t5=eye(4)sa ve my_ va r.ma tc lea r a lllo a d my_ va r.ma tc lc VMVI 3 5 y pIVIVUVI l-local modelsNatural logarithmaverage fitModel-B.B-splineModel-B.Syntaxgrowth modelsModel-B.newfts=log(tsobj)linearModel-B.multiple modelsModel-B.polynomi

15、alModel-B.Descriptionpolynomial splineModel-B.newfts=log(tsobj)calculates the natural logarithm(log base itransient 1 2|Model-B.in a financial time series object tsobj.It returns another time serietruncated power seriesModel-B.containing the natural logarithms.user-defined 1 2Model-B.local regressio

16、n smoothingCuive Fi.See Alsolocal weighted mean transfer Image P.LocalHost propertyInstrume.三exk,1O2.LocalPort propertyInstrume.LocalPortMode propertyInstrume.I S H o wlocation of block namesSimulinkLock blockDials&.loessCurve Fi.logFinancia.也可以通过S e a rc h选项卡查找log 2()函数的用法。Se a r c h选项卡与I ndex选项卡不同

17、，I n d e x只在专用术语表中查找，而Sear c h 搜索的是整个H T M L帮助文献。X 。昌 拗Contents|Index Search-Demos|Search for:logTitle:log(MATLAB Functions)TitleMATLAB Function ReferenceIQ ES.loqlogFunclogFinarNatural logarithmloglpMATIlogmMATISyntaxlogsigNeur：Y=log(X)Logging or Sending DataTheLogging DataGettiDescriptionLogging I

18、nformation to DiskSavirLogFileNameDataThe log function operates element-wise on arrays.Its domain includesLoggingModeDatanegative numbers,which may lead to unexpected results if used uninterLogToDiskModeDataY=log(X)returns the natural logarithm of the elements of X.For comLoggingImagnegative 2.where

19、 z-x+汽 t the complex logarithm is returned.LoggingModeImagLogging Image Data to DiskAcqulog(z)=log(abs(z)+i*atan2(y,x)Logging OPC Server DataReadLogging and Buffering FunctionsFuncExamplesLogFileNameOPC.The statement abs(log(-l)is a clever way to generate 7r.a=ones(l,3 0)forma t c om p a ct a=ones(1

20、,30);format compact aa=Columns 1 through 8111111 11Columns 9 through 16111111 11Colujnns 17 through 24111111 11Columns 25 through 30111111 a数字显示格式有shor t、Ion g、shor t e long e 等，请参照教材的列表练习一遍。format Ion g p i format long pians=3.14159265358979format sh o r t p i format short pians=3.14161 f o r mat 1

21、 ong pi format long pians=3.14159265358979form at+p i format+pians=+-p i -pians=3.变量的搜索顺序在命令窗口中输入以下指令：p i p im s =3.1 4 1 6sin(pi);s in (p i)a n s =1.2 2 4 6 e-0 1 6ex i s t C p i 0 e xis t (p i，)a n s =5pi=O;e x i s t(*p i)p i=O;e xis t (p i)a n s =1 p ip i=|0 c 1 e a r p ie x ist(*p i)p i clear pi

22、 pians=3.1416思 考 题：3次执行exist(p i)的结果同样吗？假 如 不 同 样，试 解 释 为 什 么？答：不同样，p i本来是库存函数，但 是 假 如 被 赋 值 则 系 统 默 认 被 赋 予 的 值 为p i后来 的 值，但是当执行cle a r p i之 后 所 赋 的 值 被 清 空，因此p i的 值 又 成 为3.141 6 圆 周 率p i是系统的默认常量，为 什 么 会 被 改 变 为0?答:p i本 来 是库存函数，但 是 假 如 被 赋 值 则 系 统 默 认 被 赋 予 的 值 为p i后来的值,但 是 当 执 行cl ear p i之 后 所 赋 的

23、 值 被 清 空，因 此p i的 值 又 成 为3.141 6实验二 MATLAB语言基础一、实验目的基本掌握M A T LA B向量、矩阵、数组的生成及其基本运算(区分数组运算和矩阵运算)、常用的数学函数。交接字符串的操作。二、实验内容(1)向量的生成和运算。(2)矩阵的创建、引用和运算。(3)多维数组的创建及运算。(4)字符串的操作。三、实验环节1.向量的生成和运算1)向量的生成直接输入法:A=56423 B =1;2;3;4；5B2345冒号生成发:A=l:2:10,B=l：10,C=1 0:-1：12A9357B8721 043569C874096153函数法:L ins p ac e

24、)是线性等分函数，1 ogsp a c e()是对数等分函数。A=linspace(l,1 0),B=linspac e(1,3 0,10)A=Co 1 u mn s 1 t h r ough 91.54552.36361.00 0 01.63641.090 91.727 3C ol u mns 10 th r oug h 1 81.81821.90911 .18182.00001.27272.0 9 091.36362.1 8 181.4 5 452.2 7 272.4 5452.5455C o 1 umn s1 9 th ro ugh272.63642.72732.81822.9 0 91

25、3.000 03.0 9 093.18183.27273.3 6 3 6Columns 28t h roug h 3 64.00 0 04.8 1 823.4 54 53.54 5 53.63 6 43.7 2 7 33.818 23.9 0914.09094.1818Co 1 umn s 3 7 th r oug h 454.272 74.9 0 9 1C o lumns 4 65.09096 3645.7273Co 1 u mn s5.9 0 9 14.36 3 64.45454.54554.63644.7 2735.0 000th r ou g h 545.18185.272 75.36

26、365.45455.54555.5.81 8 255 t hrou g h 636.00006.09096.1 8186.2 7 276.3 6366.4 5 4 56.545 56.6364Co 1 umns 64 th r o u gh 726.72 7 36.81826.90917.00007.090 97.18 1 87.2 727 7.36 3 67.4 5 4 5C o 1 umns 73 t h r o u g h 817.54 5 58.0 909 8.18 1 87.63648.2 7 2 77.727 37.8 18 27.9 0 9 18.0000Columns 8 2

27、t h r o u g h908.3 6 3 68.909 1 9.0 0 008.454 59.09098.5455 8.63648.72738.8 1 8 2Co 1 umn s 91through 999.18 1 89.7273 9.81 8 2Column 10010.00009.2 7279.90 9 19.3 63 69.4 54 59.54559.6 364B=C o lum n s 1 throu g h 91 .0 0 004.22 2 27.4444 10.666 71 3.88 89 1 7 J11 1 20.3 33 3 23.5556 26.7778Colum n

28、1 03 0.0000 A=logsp a c e (0,4,5)A=1 10 1 00 1 0 00 10000练习：使 用 lo g s p a c e()创 建 1-4万的有1 0 个元素的行向量。答案：A=lo g s p a c e(l,1 0,4*p i)A=1 .0e+0 10*C o lum n s 1 through 90.0 0 0 0 0.0 0 00 0.0 0 0 0 0.0000 0.00 0 0 0.0 0 0 00.0001 0.0 0 0 5 0.0 035C o lumns 10 t hrough 120.0 23 1 0.1520 1.00002)向量的运

29、算维数相同的行向量只见可以相加减，维数相同的列向量也可以相加减,标量可以与向量直接相乘除。A=l 2 3 4 5 ,B=3:7,A=1 2 3 4 5B=34567 A T=Af,B T=B;AT=12345BT=34567E1=A+B,E2=A-BEl=4 6 8 10 12E2=-2-2-2-2-2 F=AT-BT,F=-2-2-2-2 G l=3*A,G 2=B/3,Gl=3 6 9 1 2 1 5G2=1.0000 1.3 3 33 1.6667 2.0000向量的点积与叉积运算。A=one s(1,1 0);B=(1 ：10)；BT=B；El=dot(A,B)E l=5 5 E 2=

30、A*B TE2=55 c 1 ear A=1:3,B=3:5,A=1 2 3B=3 4 5 E=cr o ss(A,B)E=-2 4-22.3 3 332.矩阵的创建、引用和运算1）矩阵的创建和引用矩阵是由m x n元素构成的矩阵结构，行向量和列向量是矩阵的特殊形式。直 接 输 入 法：A=l 2 3;4 5 6A=12 34 5 6 B=1 4 72 583 69B=1472 5 83 6 9 A(l)an s=1 A(4:end)ans=5 3 6 B(:,1)ans=123 B(:)a n s =123456789 B(5)a n s =5抽取法 c l e a r A=L 1 2 3

31、4;56 7 8;9 1 0 1 1 1 2；1 3 1 4 1 5 1 6 A1234567891 01 1121 31 41 51 6 B=A (1：3,2:3)B2 367101 1C=A(L1 3,2 4)C=2 41 0 12 A (1 3;2 4)an s=195 13函数法:A=ones(3,4)A=111111111111 B=zero(3)?Und e f i ne d func ti on o r m e thod*z e ro f or inp u t argument s of typedouble1.B=z e ros(3)B0 000000 00 C=e y es(

32、3,2)?u n d ef i n ed fu notion or method e yes fo r i npu t argum e nts of t y pe,doub 1 e C=e ye(3,2)C=100 100 D=magic(3)D=8163 5 74 9 2拼接法 c 1 ear A=ones(3,4)A=111111111111 B=ze r os(3)B=00 00 0 000 0 C=eye(4)C0000000010000D 二ABD=100011000000 F=A;CF1111100100100000000拼接函数和变形函数法:cl e a rA01 A=0 1;1

33、 1 1 B=2*o n es(2)B=22 ca tans=0122022(1,A,B,A)112211 c a t(2A,B,A)?cat(2A,B,A)Er r o r:Une x p e c ted MATLAB e x p r es s ion.ca t(2,A,B,A)ans2 2 02 2 10 11111 r e p ma t(A,2,2)an s=010 101011 11 re p mat(A,2)a n s=0 10 111 1 10 1110 11 1练习：使用函数法、拼接法、拼接函数法和变形函数法，按照规定创建以下矩阵:A 为3x4的全1矩阵、B为3 x 3 的0 矩

34、阵/C 为3 x 3 的单位矩阵、D 为3 x 3 的魔方阵、E 由C 和D 纵向拼接而成，F 抽 取 E 的 2-5 行元素生成、G 由F 经变形为3 x 4 的矩阵而得、以G 为子矩阵用复制函数(repmat)生成6 x 8 的大矩阵H。答案：A=ones(3,4)A=11 1 11 1 11 1 11 B=z e ros(3,3)B=0 0 00 0 00 0 0 C=ey e(3)C=10 00 1 00 0 1 D=ma g i c(3)D=8 1 63 5 749 2E=C;D 1E=10 00 1 00 0 18 1 63 5 74 9 2 F二(2:5,:)?F=(2:5,:)

35、E rror：E x pression or statement i s inco r rec t-p o ssi b ly unbal a nee d(,o r.F=E(2:5,：)F=0100 0 18163 5 7 G=re s pace(E,3,4)?Unde f ined c o mmand/func t io n*r es p a c e,.G=r e sp a c e(F,3,4)?U n define d comma n d/fu nre s p ace.G=re s hap e(F,3,4)H=r epmat(G,2,2)c t i on 008315600703031015

36、60568800780070003560700315680072）矩阵的运算矩阵的加减、数乘与乘法已知矩阵：A=1 23-1,A123-1 B=-l 01 21B=-1 01 2 A+Ban s=0 24 1 2*Aan s=2 46-2 2*A-3*Ban s=5 43-8 A*Ba ns=1 4-4-2 矩 阵 的 逆矩阵 fb r m a tr a t ;A=1 0 1;2 1 2；0 4 6Ad =101212046 A 1=i n v(A)A l =-1/32/3-1/6-2104/3-2/31/6A*A 1a n s =100010001矩阵的除法 a =1 2 1;3 1 4；2

37、 2 1 ,b=l 1 2 ,d=b*a =121314221b =11212 c 1=b*inv(a),c 2=b/acl=6/7 3/7-4/7c2=6/73/7-4/7 c 3=inv(a)*d,c4=a bc 3=12/7-4/7?Er r o r usi n g=mldivi d eMa t rix dime n s ions mu s t a g r ee.c 3=inv(a)c 4=ad2/7-4/7c4=12/7-4/7练习：按下列规定求出各种的矩阵运算的值求矩阵的秩、特性值和特性向量、矩阵的乘幕与开方；矩阵的指数与对数矩阵的提取与翻转答案：A=6 3 43-2 5 7-48-

38、1-3 7A=6 3 4 3-2 5 7-48-1-3-7 B=rank(A)B=3 r b=rank(A)rb=3 X,Lam da=eig(A)?Erro r u s in g =e i gM a t r ix m u st b e squar e.X,Lamda=e igs(A)?E r r o r using=ei g sA m u st b e a square matr i x or a fu n ction which com p utes A*x.C=6 34-2 5 78-1-3C=6 3 4-2 5 78-1-3 X,Lamda=eigs(C)X=0.8 0 1 3-0.1

39、094-0.1 6060.3638-0.6 5 64 0.86 6 90.4749 0.7464-0.47 1 9Lamd a=9.7326000一 3.29280001.560 2 L X,L a md a =eig(C)X=0.8013-0.1094-0.16060.3 6 3 8-0.65 60.8 66 90.474 90.74 6 4-0.4 7 1 94Lamda9.7326000-3.29 2 80001.5602 X,Lamda=eig(C)X0.801 3-0.109 4-0.16 0 60.3 6 3 8-0.6 5 6 40.8 6690.4 7490.7 464-0.4

40、7 1 9L a mda=9.732 6000-3.292 80001 .5 6 02 D=AA2?Err o r us i ng=mpowe rMatr i x mu s t be sq u a r e.D=CA2D=62 29 3334 12 626 22 34 E=s q rtm(C)E=2.2 4 47+0.2 706i 0.6 9 74-0.14 0 0i 0.9422-0.3 494 i-0.5815+1.6 2 44i 2.1005-0.8 405i 1.7620-2.0970i1.97 1 9-1.8471i-0.3017+0.9 557i 0.0 2 36+2.3845i F

41、=e xpm(C)F=1.0e+004*1.065 3 0.5415 0.63230.4 8 300.2 4650.28760.6 3160.32 0 60.37 4 5 G=logm(C)Warn i n g:P rin c i p a 1 ma t r i x lo g arithm is not d e fined for A withn o n p o s i t ive real eig e n v alues.A n o n-principal m a tri xlo g a r it h m is re t urn e d.In funm at 153I n 1 o gm at

42、27G=1.7 1 29+0.4 68 6i 0.5305-0.2425i 0.5 429-0.6049 i1.1 9 3 8+2.8123i 0.3 6 5 8 一 1.4552i 0.55 1 4-3.6305i-0.0 7 4 8-3.1978i 0.7 4 19+1.6546 i 1.8333+4.1282 i H=f 1 ip 1 r(C)H=4 3 67 5-2-3-1 8 1=t riu(C)I=6 3 40 5 700 J=t r il(C)J=6 0 0-2 5 08-1-3 K=di a g(C)K=653.多维数组的创建及运算1)多维数组的创建A 1=1,2,3;4 5

43、6;7,8,9;A2=re s hape(10:18,3,3)A2=1 0 1 3 1611 14 171 2 1 5 18 T 1(：,:,l)=ones(3)；T 1(:,：,2)=zeros(3)T11 1 11 1 1T 1 (:,:,2)=0 0 00 0 00 0 0T2=ones(3,3,2)T 2(：,:,1)=11111 1111T 2(:,：,2)=1 1 111 11 1 1 T 3 =c a t(3,A 1,A 2),T 4=r e p m a t(A 1,1,1,2 J)T3(：,：,1 )=1234 5 67 8 9T 3（:,:,2）=1 0 1 3 1 61 1

44、 1 4 1 71 2 1 5 1 8T 4（:,：,1）=1 2 34 5 67 8 9T4（：,：,2）=1 2 34 5 67 8 92)多维数组的创建数组运算用小圆点加在运算符的前面表达，以区分矩阵的运算。特点是两个数组相相应的元素进行运算。A=1 :6 ；B=o n e s (1 ,6);C 1 =A+B,C 2=A-BC 12 3 4 5 6 7C2=0 1 2 3 4 5 C 3=A.*B,C 4=B /A,C 5=A.BC3=1 2 3 4 5 6C4=1.OOOO O.5OOO 0.3333 0.2500 0.2023C 5=0.16671.000 0 0.5 0 OO0.3

45、 3 330.2500 0.2 0230.1667关系运算或逻辑运算的结果都是逻辑值。I=A3,C 6=A(I)1 =0 0 0 1C6=4 5 6 A 1=A-3,I2=AI&AA 1 =-2 -1 0 1 2 312=1 1 0 111 I 3=I13=1 1 10 0 04,字符串的操作1)字符串的创建 S 1=n i ke MATLAB1S 1 =Hike MATL AB S2=I z m a s t u e n t.7S2=Pm a stuent.S 3=S 2/and；SIS3=Fm a stuent.andl 1 ike MATLA B2)求字符串长度 1 ength(S 1)a

46、 n s=1 2 s i z e(S 1)ans=1 123)字符串与一维数值数组的互相转换 CSl=ab s(S1)C S l=7 3 10 8 1 05 107 1 0 1 32 7 7 6 5 84 7665 66 CS2=doubl e(S 1)C S2=73 108 1 0 5 107 10 1 3 2 7 7 6 5 84 7 665 66 char(CS2)ans=H i k e MATLAB setst r(CS 2)ans=Hike MATLAB练习:用ch ar()和向量生成的方法创建如下字符串AaB b C c XxY yZz.S 1=6 5：90;S2=9 7：122；

47、C=SI;S2;c=c(：y；S 3=dou b 1 e(C);c h ar(S 3)an s=A a BbC c DdEeFfGgHhl i J jKkLIMmNnO o PpQqR r S sTtUuVvWwXxYyZ z实验三 MALTA B数值运算一、实验目的掌握M A T L A B 的数值运算及其运算中所用到的函数,掌握结构数组和细胞数组的操作。二、实验内容1）多项式运算。2）多项式插值和拟合。3）数值为积分。4）结构数组和细胞数组。三、实验环节1、多项式运算1）多项式表达。在 M A T L A B 中，多项式表达成向量的形式。$Q:S4+3 S3-5 S2+9 在 M A T

48、L A B 中表达为 s =1 3 -5 0 9 2）多项式的加减法相称于向量的加减法，但必须注意阶次要相同。如不同，低阶次的要补0。如多项式2 s?+3 s +9与多项式s 4+3 S 3 _ 5 s 2+4s +7相加。s l=0 0 2 3 11 s 2=1 2-5 4 7 s3=sl+s2答；s 1 =00 2 3 11s2=12-547s3=1 2-37 183）多项式的乘、除法分别用函数con v 和的de con v 实现。s 1=23 1 1 s 2=l 3-5 4 7 s3=conv(s 1 ,s 2)s4=d ec o nv(s 3,s 1)答;s 1 =2 3 11s 2

49、=13-54 7s 3=2 9 1 0 26-2 9 65 77s413-5474)多项式求根用函数roots。sl=24 2 roots(sl)答;sl=242ans=-1-15)多项式求值用函数P olyval s 1 =2 4 1 -3 pol y va 1 (s 1,3)x=l：1 0 y=polyva 1 (s 1 ,x)答；si=2 4 1-3an s=901 2 3 4 5 6 7 8 9 1 0y=Co 1 umns 1 t hroug h 84 31 90 193 3 52 579 886 1 2 8 5C o 1 umns 9 t h r ough 10练习：求(s-+?(

50、s +3)(s +l)的“商,及余数。s+2s+1 s 1=1 0 I；s2=1 3；s3=1 1;s4=l 0 2 1;q,r=d e c onv(conv(conv(s 1,s2),s3),s4)1 4答;q0 0 2-5 -1多项式插值和拟合有一组实验数据如附表1 1 所示。附 表 1-1X1234567891 0Y1 632701422604 3668 210 101 4 321960分别用拟合(二阶至三阶)和插值(线性和三次样条)的方法来估算X=9.5时Y的值。以下是实现一阶拟合的语句。x=l:1 0 y =|16 3 2 70 1 4 2 26 0 4 36 682 1010 13