MATLAB全部实验及答案_第1页
MATLAB全部实验及答案_第2页
MATLAB全部实验及答案_第3页
MATLAB全部实验及答案_第4页
MATLAB全部实验及答案_第5页
已阅读5页,还剩31页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、.实验项目序号1 MATLAB基本操作一、实验目的熟悉MATLAB软件环境,掌握命令窗口的使用。二、实验基本知识1、熟悉MATLAB的桌面工具和环境2、Start菜单三、主要仪器及耗材计算机四、实验内容和步骤1、 MATLAB的启动及操作界面认识说明有哪几种启动的方法;操作界面有哪些窗口;对常用的窗口作些简单的介绍。 使用Windows的“开始”菜单。 运行MATLAB系统的启动程序matlab.exe。 利用桌面上的快捷方式。 MATLAB主窗口 命令窗口(Command Window) 当前目录窗口(Current Directory) 工作空间窗口(Workspace) 命令历史窗口(C

2、ommand History)2、 下列符号中哪些是MATLAB的合法变量名?用给变量赋值的方法在机器上验证你的答案(3vars、global、help、My_exp、sin、X+Y、_input、E-4、AbCd、AB_C_。)help,My_exp,sin,AbCd,AB_C_;3、命令窗口的简单使用1)简单矩阵的输入 1 2 3A = 4 5 6 7 8 92)求12+2(7-4)32的算术运算结果4、有关向量、矩阵或数组的一些运算1)设A=15;B=20;求C=A+B与c=a+b?2)设A=1 2 3;4 5 6;7 8 9,B=9 8 7;6 5 4;3 2 1;求A*B与A.*B?

3、A*B遵循矩阵运算,A.*B遵循每个对应元素之间进行相乘3)设A=10,B=20;求C=A/B与D=AB?4)设a=1 -2 3;4 5 -4;5 -6 7请设计出程序,分别找出小于0的矩阵元素及其位置(单下标、全下标的形式),并将其单下标转换成全下标 a=1 -2 3;4 5 -4;5 -6 7;/%定义矩阵 b = a(find(a c = find(a x,y = find(a x,y = ind2sub(size(a),find(a a=1 2 3;3 4 2;5 2 3; a2 % 矩阵运算 a.2%数组与矩阵的乘法7)有一段程序设计如下,请思考并说明运行结果及其原因clear X=

4、1 2;8 9;3 6; X( : )解:转化X为列向量8)使用两种方法,创建一稀疏矩阵解:1)a = eye(5);%创建5阶单位矩阵sparse(a)%将a构造为稀疏矩阵ans = (1,1) 1 (2,2) 1 (3,3) 1 (4,4) 1 (5,5) 12) data = 1 2 2 2; ir = 1 3 2 4; jc = 1 2 4 1; s=sparse(ir,jc,data,4,4)s = (1,1) 1 (4,1) 2 (3,2) 2 (2,4) 29)写出下列指令的运行结果 A = 1 2 3 ; B = 4 5 6 ; C = 3.A C = 3 9 27 D = A

5、.B D = 1 32 7295、已知若需要计算t-1,1,取间隔为0.01,试计算出相对应的y值。 t = -1:0.01:1; y = sqrt(3)/2*exp(-4*t).*sin(4*sqrt(3)*t+pi/3) %注意要用点乘五、实验注意事项1、 注意区分大小写2、 注意数组、矩阵运算的不同之处六、思考题1、MATLAB桌面环境有几种样式?Default:缺省的界面Command Windows Only:仅包含命令行窗体(Command Window)History and Command Window: 包含两个窗体命令行窗体(Command Window)和命令历史窗体(C

6、ommand History),两个窗体并列在界面中All Tapped:启动所有MATLAB常用的图形界面工具All but Command Window Minimized:Command Window窗口显示,其他的窗口分别叠加在MATLAB窗体下2、如何使窗体浮动出来?又如何使窗体内嵌回MATLAB的界面中?在MATLAB默认的图形窗体环境下,单击命令行窗体左上角按钮,就可以将MATLAB命令行窗体弹出执行“Desktop”菜单下的“Dock Current Directory”命令,可将命令行窗口重新嵌入到MATLAB的界面中实验项目序号2 MATLAB运算基础(一)一、实验目的掌握

7、MATLAB各种表达式的书写规则及常用函数的使用二、实验基本知识1、矩阵和向量的创建方法和访问方法1)创建方法 逐个输入向量的元素 利用运算符“:”语法:X = J:INC:K 使用函数linspace和logspace语法:x = linspace(x1,x2,n),x = logspace (x1,x2,n)2)访问方法:利用向量或矩阵的索引来完成相应的操作。2、数据类型,以及操作不同数据类型数据的函数1)数据类型2)操作函数整数类型的运算函数创建逻辑类型数据的函数MATLAB的逻辑运算MATLAB的关系运算格式化字符三、主要仪器及耗材计算机四、实验内容和步骤1、已知A=12 34 -4;

8、34 7 87;3 65 7,B=1 3 -1;2 0 3;3 -2 7,求下列表达式的值:1)A+6B和A2-B+I(I为单位矩阵)2)A*B和A.*B3)A3和A.34)A,B和A(1,3,:);B22、设有矩阵A和B,A=1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25,B=3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 111)求它们的乘积C2)将矩阵C的右下角3x2子矩阵赋给D3、已知A=23 10 -0.778 0;41 -45 65 5;32 5 0 32;6 -9.54 54

9、3.14,取出其前3行构成矩阵B,其前两列构成矩阵C,其左下角3x2子矩阵构成矩阵D,B与C的乘积构成矩阵E,分别求E=10&A=10&A A=23 10 -78 0;41 -45 65 5;32 5 0 32;6 -54 92 14; a = find(A=10&A rot90(A)ans = 0 5 32 14 -78 65 0 92 10 -45 5 -5423 41 32 6右旋: rot90(A,3)ans = 6 32 41 23 -54 5 -45 10 92 0 65 -78 14 32 5 05、求S=20+21+22+23+24+210的值(提示:利用求和函数与累乘积函数。

10、) a=2*ones(1,10);%生成全为2,个数为10的行矩阵 b=cumprod(a);%得到累乘量(如21、22)等等 c=sum(b)+1c = 20476、用reshape指令生成下列矩阵,并取出方框内的数组元素。答案:a = 1:5 10:-1:6;b = 11:25;A = reshape(a,5,2); reshape(b,5,3)%生成矩阵Ad1 = A(2,2:3)d2 = A(2:4,4)d3 = A(4:5,1:3)五、实验注意事项不同的函数有不同的功能,同一函数又有不同的用法,注意正确使用各函数。六、思考题元素的单下标访问和全下标访问有何关联之处?实验项目序号3 M

11、ATLAB运算基础(二)一、实验目的掌握MATLAB中有关字符串、元胞、结构的函数的使用。二、实验基本知识有关字符串、元胞、结构的函数常用的字符串操作函数数字和字符之间的转换函数不同数值之间的转换函数格式化字符元胞数组的操作函数在cellfun函数中可用的函数结构操作函数三、主要仪器及耗材计算机四、实验内容和步骤1、求100,999之间能被21整除的数的个数(提示:先利用冒号表达式,再利用find和length函数。)答案: t = 100:999 i = find(rem(t,21)=0) length(i)ans = 432、建立一个字符串向量,删除其中的大写字母(提示:利用find函数和

12、空矩阵。) a = dadaAf b = find(a=A&a a(b) = a =dadaf3、写出使以下这段文字成为字符串的MATLAB指令。注意保持这段文字的格式。在英式用法中,引号通常是单引号,如Fire!。In GB usage quotation marks are usually single: fire!答案: disp(In GB usage quotation marks are usually single: fire!)4、用结构体数组来存储2名学生的基本情况数据,每名学生的数据包括学号、姓名、专业和3门课程的成绩。a = struct(num,1,2,name,mik

13、e,hike, major,tele,score, rand(2,1)*100,rand(2,1) *100);5、建立矩阵A并回答有关问题A1,1=1;A1,2=Brenden;A2,1=reshape(1:9,3,3);A2,2=12,34,2;54,21,3;4,23,67;1) 如何使用访问的方式将Brenden改写成BRENDEN? A1,2=upper(A1,2)A = 1 BRENDEN 3x3 double 3x3 cell2) 分别执行A(3)=和A3=后,A的值各是多少?并说明原因。 A3=A = 1 3x3 double 3x3 cell A(3)=A = 1 3x3 d

14、ouble 3x3 cell6、利用MATLAB提供的rand函数和圆整函数随机生成43整数矩阵A,进行如下操作1)A各列元素的平均值和中值平均值 A_mean=mean(A) A_mean =0.5000 0.7500 0.7500 中值 A_median=median(A)A_median = 0.5000 1.0000 1.00002)A的最大元素和最小元素及它们的位置最大值 A_max=max(max(A)A_max = 1位置: x,y=find(A=A_max) A_maxlocal=x(1) y(1)最小值 A_min=min(min(A)A_min = 0位置: xx,yy=f

15、ind(A=A_min)A_minlocal=xx(1)yy(1)3)求A的每行元素的和以及全部元素之和 A_rsum=sum(A,2)%行的和 A_sum=sum(sum(A)%全和4)分别对A的每列元素按升序、每行元素按降序排列 A_rsort=abs(sort(-A,2)五、实验注意事项字符串函数的使用细节;对元胞数组中的元胞的访问;结构数组的创建方法。六、思考题访问元胞数组的两种操作符“”和“()”之间有什么不同?使用“”直接获取元胞数组的元胞内容使用“()”直接访问元胞数组的元胞,获取的数据也是一个元胞数组。实验项目序号4 MATLAB循环结构程序设计一、实验目的1、掌握建立和执行M

16、文件的方法;2、掌握利用for语句和while语句实现循环结构的方法。二、实验基本知识MATLAB的M语言文件是纯文本格式的文件,MATLAB也提供了一个编辑器,叫作meditor,它也是系统默认的M文件编辑器。运行meditor的方法非常简单,在MATLAB命令行窗口中键入下面的指令就可以打开meditor: edit这时MATLAB将启动meditor,然后创建一个未命名的空白文件。这时用户就可以直接在编辑器中键入MATLAB指令,开发M语言文件了。M语言文件可以分为两类,其中一类是脚本文件,另外一类叫作函数文件。脚本文件是最简单的一种M语言文件。所谓脚本文件,就是由一系列的MATLAB指

17、令和命令组成的纯文本格式的M文件,执行脚本文件时,文件中的指令或者命令按照出现在脚本文件中的顺序依次执行。脚本文件没有输入参数,也没有输出参数,执行起来就像早期的DOS操作系统的批处理文件一样,而脚本文件处理的数据或者变量必须在MATLAB的公共工作空间中。 脚本文件中常用的MATLAB指令程序流程控制语句:1、for循环结构语法:for i初值:增量:终值 语句1 语句n end 说明:1)i初值:终值,则增量为1。 2)初值、增量、终值可正可负,可以是整数或小数,只须符合数学逻辑。2、while 循环结构 语法:while 逻辑表达式 循环体语句 end 说明:1)whiIe结构依据逻辑表

18、达式的值判断是否执行循环体语勾。若表达式的值为真,执行循环体语句一次、在反复执行时,每次都要进行判断。若表达式的值为假,则程序执行end之后的语句。2)为了避免因逻辑上的失误,而陷入死循环,建议在循环体语句的适当位置加break语句、以便程序能正常执行。(执行循环体的次数不确定;每一次执行循环体后,一定会改变while后面所跟关系式的值。) 3)while循环也可以嵌套、其结构如下: while逻辑表达式1 循环体语句1 while逻辑表达式2 循环体语句2 end 循环体语句3 end 三、主要仪器及耗材计算机四、实验内容和步骤1、分别使用for循环语句和while循环语句找出最小的n值,使

19、得n!10100,并求出n!。要求显示的结果为 n=70n!=1.197857e+1001e100或70!=1.197857e+1001e100程序for:s =1;for n=1:100 s=s*n; if(s10100) break; endendstr1=n=,num2str(n);str2=n!=,num2str(s);disp(str1)disp(str2)程序while:n=1;s=1;while(s3求f1 f50中:1)最大值和最小值及它们的位置,各数之和(最大值给变量MAX,其位置给变量x,最小值给变量MIN,其位置给变量y,各数之和给变量SUM)。程序:f(1)=1;f(2

20、)=0;f(3)=1;for n=4:50f(n)=f(n-1)-2*f(n-2)+f(n-3);endMAX,x=max(f)%?MIN,y=min(f)SUM=sum(f)2)正数、零、负数的个数(依次赋值给变量positive,zero,negative)。将MAX,x=max(f)和之后的语句换为:positive=length(find(f0)zero=length(find(f=0)negative=length(find(f10000 break; endend disp(i=,num2str(i) disp(a(i)=,num2str(a(i)五、实验注意事项1、for语句和w

21、hile语句的应用场合2、注意输出格式的要求六、思考题请思考下列两段程序及其运行结果,说明原因。1) clear (2)clearfor i=1:10; for i=1:10;x(i)=i.2; x(i)=i.2;end endx x(i)实验项目序号5 MATLAB选择结构程序设计一、实验目的1、掌握建立和执行M文件的方法;2、掌握利用if语句和switch语句实现选择结构的方法。二、实验基本知识MATLAB的M语言文件是纯文本格式的文件,MATLAB也提供了一个编辑器,叫作meditor,它也是系统默认的M文件编辑器。运行meditor的方法非常简单,在MATLAB命令行窗口中键入下面的指

22、令就可以打开meditor: edit这时MATLAB将启动meditor,然后创建一个未命名的空白文件。这时用户就可以直接在编辑器中键入MATLAB指令,开发M语言文件了。M语言文件可以分为两类,其中一类是脚本文件,另外一类叫作函数文件。脚本文件是最简单的一种M语言文件。所谓脚本文件,就是由一系列的MATLAB指令和命令组成的纯文本格式的M文件,执行脚本文件时,文件中的指令或者命令按照出现在脚本文件中的顺序依次执行。脚本文件没有输入参数,也没有输出参数,执行起来就像早期的DOS操作系统的批处理文件一样,而脚本文件处理的数据或者变量必须在MATLAB的公共工作空间中。脚本文件中常用的MATLA

23、B指令 程序流程控制语句: 1、if-else-end分支结构语法:if 表达式 1 语句1 else if 表达式 2(可选) 语句2 else(可选) 语句3 end end 说明:1)if结构是一个条件分支语句,若满足表达式的条件,则往下执行;若不满足,则跳出if结构。2)else if表达式2与else为可选项,这两条语句可依据具体情况取舍。3)注意:每一个if都对应一个end,即有几个if,记就应有几个end。2、switch-case结构语法:switch表达式case常量表达式1 语句组1 case常量表达式2 语句组2 otherwise 语句组n end 说明:1)switc

24、h后面的表达式可以是任何类型,如数字、字符串等。2)当表达式的值与case后面常量表达式的值相等时,就执行这个case后面的语句组如果所有的常量表达式的值都与这个表达式的值不相等时,则执行otherwise后的执行语句。三、主要仪器及耗材计算机四、实验内容和步骤1、输入一个百分制成绩,要求输出成绩等级A、B、C、D、E。其中90分100分为A,80分89分为B,70分79为C,60分69分为D,60分以下为E。要求:1)分别用if语句和switch语句实现2)对不合理的成绩应输出出错信息“输入的成绩不合理”(若成绩出现小数,则只能是“.5”)程序switch:score=input(:); i

25、f rem(score,0.5)=0 switch score case num2cell(90:0.5:100) disp( A) case num2cell(80:0.5:89) disp( B) case num2cell(70:0.5:79) disp( C) case num2cell(60:0.5:69) disp( D) case num2cell(0:0.5:59.5) disp( E) otherwise disp() end else disp(0.5)end程序if: 2、 x=input(?:); 3、 if rem(x,0.5)=04、 if (x=90) 5、 di

26、sp(A) 6、 elseif (x=80 & x=70 & x=60 & x=69) 11、 disp(D) 12、 elseif (x=0) 13、 disp(E) 14、 else15、 disp()16、 end17、 else18、 disp(?0.5)19、 end20、 设计程序,完成两位数的加、减、乘、除四则运算,即产生两个两位随机整数,再输入一个运算符号,做相应的运算,显示相应的结果,并要求结果显示类似于“a=x+y=34”。程序:x=round(rand(1)*100);y=round(rand(1)*100);disp(x=,num2str(x);disp(y=,num2

27、str(y);A=input(请输入一个运算符号:,s); %书上114页有说明switch(A)case +a=x+y; disp(a=x+y=,num2str(a);case -a=x-y;disp(a=x-y=,num2str(a);case *a=x*y;disp(a=x*y=,num2str(a);case /a=x/y;disp(a=x/y=,num2str(a);end3、求下列分段函数的值 要求:用if语句实现,分别输出x = -5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y值。其中x的值以向量的形式从键盘输入。x=-5.0,-3.0,1.0,2.0,2.5,3

28、.0,5.0ifx=0&x75b(i)=100-a(i);elseif a(i) edit这时MATLAB将启动meditor,然后创建一个未命名的空白文件。这时用户就可以直接在编辑器中键入MATLAB指令,开发M语言文件了。M语言文件可以分为两类,其中一类是脚本文件,另外一类叫作函数文件。M函数文件和脚本文件不同,函数文件不仅有自己特殊的文件格式,不同的函数还分别具有自己的工作空间。同其他高级语言类似,M函数文件也有局部变量和全局变量。函数文件的结构和编写方法:函数只能在函数体内对变量进行操作,也就是只能访问函数本身工作空间中的变量。M函数文件的结构:(1) 函数定义行(function)

29、(2) H1行(函数帮助文本的第一行) (3) 函数帮助文本(4) 函数体(5) 注释函数文件编写后,保存时,其文件名必须与函数名相同。注意:函数名不要与MATLAB自身的函数命令相同。 三、主要仪器及耗材计算机四、实验内容和步骤1、编写一个函数文件chengji.m,用于求两个矩阵的乘积和点乘,然后在命令行窗口中调用该函数(提示:使用varargout函数),要求:1)当两矩阵是尺寸相同的方阵时,返回它们的乘积和点乘。 2)当两矩阵不是尺寸相同的方阵时,则根据它们的尺寸,分别进行乘积或点乘。 3)当两矩阵的尺寸不能满足乘积或点乘的要求时,则提示“这两个矩阵不能乘”。 程序为:function

30、 varargout=chengji(a,b)x,y=size(a);xx,yy=size(b);if (xx=yy&x=y&x=xx) varargout1=a*b; varargout2=a.*b;elseif (x=xx&y=yy) varargout1=a.*b; disp() varargout2=0;elseif (x=yy&xx=y)disp()varargout1=a*b;varargout2=b*a;elseif (x=yy|xx=y)tryvarargout1=a*b;disp()varargout2=0;catchvarargout1=b*a;disp()varargou

31、t2=0;endelseerror()end2、编写一函数qiuhe.m,实现sum(A)的功能,其中A为矩阵。程序为:function y=qiuhe(a)y=;x,yy=size(a);if x=1 s=0; for k=1:yy s=s+a(1,k); end y=s;elsefor j=1:yy s=0; for i=1:x s=s+a(i,j);end y(1,j)=s;endend3、编写一个阶乘函数jiecheng.m,然后在试验三的编程题“分别使用for和while语句找出最小的n值,使得n!10100,并求出n!”程序中调用该函数。程序:function s=jiecheng

32、(n)s=1;if n=0 s=1;elseif n0error(n0)elsefor i=1:ns=s*i;endend五、实验注意事项1、注意脚本文件和函数文件的区别2、函数文件的调用六、思考题已知一函数文件如下:function y=temf(a)if length(a)2y=a;elsey=max(a)-min(a);end则在命令窗口执行有:temf(23 45 67 89)= 66 temf(78)= 78 实验项目序号7 数据插值与曲线拟合实验一、实验目的1、掌握数据统计和分析的方法;2、掌握数据插值与曲线拟合的方法及其应用。二、实验基本知识1、数据插值运算MATLAB提供的插值

33、计算函数函数说明函数说明interp1一维插值(数值查表)griddata3三维数据网格的超表面数据插值interp1q一维快速插值(数值查表)griddatan多维数据网格的超表面数据插值interp2二维插值(数值查表)mkpp产生分段多项式interp3三维插值(数值查表)pchip分段的厄密多项式interpnN维插值(数值查表)ppval计算分段多项式的数值interpft使用FFT算法的一维插值spline三次杨条插值griddata二维数据网格的表面数据插值unmkpp分段多项式的细节1)一维插值:在MATLAB中实现一维插值的函数是interp1。调用格式为:Yi=interp

34、1(X,Y,Xi,method) X和Y为已知的原始数据,分别描述采样点和采样值; Xi为需要计算的插值点,是一个向量或标量; Yi是插值结果; X和Y,Xi和Yi若都为向量,则等长,若Y和Yi为矩阵,则length(X) and size(Y,1) must be the same. method是插值方法,允许的取值有多种; 若进行插值运算时,Xi的取值超过了X的范围,则需要进行外插值运算:Yi=interp1(X,Y,Xi,method,extrap)。2)二维插值:在MATLAB中实现二维插值的函数是interp2。调用格式为:Zi=interp2(X,Y,Z,Xi,Yi,method

35、) X和Y为原始数据,是两个等长的已知向量,分别描述两个参数的采样点; Z是与参数采样点对应的函数值; Xi ,Yi是两个向量或标量,描述欲插值的点; Zi是根据相应的插值方法得到的插值结果; method是插值方法,允许的取值有多种,与一维插值函数相同; X,Y,Z也可以是矩阵的形式; 若进行插值运算时,Xi的取值超过了X的范围,则需要进行外插值运算:Zi=interp2(X,Y,Z,Xi,Yi,method,extrap)。2、曲线拟合曲线拟合的目的是用一个较简单的函数去逼近一个复杂的或未知的函数。利用MATLAB进行曲线拟合主要有两种方法:回归法拟合和多项式拟合。1)回归法拟合主要是使用

36、MATLAB的左除运算来寻找曲线拟合解析函数的系数。2)多项式拟合在MATLAB中进行多项式拟合主要使用的函数有: polyfit:用来进行拟合计算P=polyfit(x,y,n)x和y为参与曲线拟合的原始数据,n为进行拟合计算的多项式次数,函数的返回值是多项式的系数,即函数的运算结果为多项式的系数向量 polyval:用来计算多项式的数值y=polyval(p,x)p为多项式的系数,x是变量的数值,得到的结果就是函数的数值向量。三、主要仪器及耗材计算机四、实验内容和步骤1.已知原始数据,插值点为,试用四种不同的插值方法进行插值,并在一副图中绘出四种插值曲线。程序:x=0:1:4*pi;y=s

37、in(x).*exp(-x/5);xi=0:0.3:4*pi;y1=interp1(x,y,xi,nearset);y2=interp1(x,y,xi,linear);y3=interp1(x,y,xi,spline);y4=interp1(x,y,xi,cubic);plot(x,y,-m,xi,y1,*r,xi,y2,ok,xi,y3,b,xi,y4,+g)legend(,)title( ,FontName,FontSize,20)xlabel(X,FontName,FontSize,16)ylabel(Y,FontName,FontSize,16)2.考虑下面的数据,一位举重教练收集了不

38、同性别的举重动员每个年龄的最大举重磅数,他相信这三者之间存在着函数关系,请创建一个年龄数组,以便估算当前队伍中队员(要求1320岁每一岁都有)的最大举重。性别年龄重举大最1315171920男2780女2560程序:t=13 15 17 19 20;sex=0,1;y=270 290 330 350 380;250 270 310 330 360 ;x=13:20;sex_2=0,1;y_2=interp2(t,sex,y,x,sex_2)surf(x,sex_2,y_2)3.某乡镇企业2009-2015年的生产利润如下年份222015利润(万元)701221442要求用多项式拟合法预测出该企

39、业2016年和2017年的利润,并绘出拟合图。程序:x=2004:2010;y=70 122 144 152 174 196 202;k=polyfit(x,y,1);x1=2004:2013;y1=polyval(k,x1);plot(x,y,-ro,x1,y1,-b*);title(2004-2013)legend(04-10,04-13)xlabel();ylabel();text(2011,polyval(k,2011),-2011) text(2012,polyval(k,2012),-2012) grid on4. 先使用命令load加载census.mat得到一组数据(cdate

40、,pop),该数据是美国自1790到1990年(以10年为一单位)的总人口。利用回归法拟合,预测出美国在2010年的人口总数,写出拟合模型,并绘出拟合图。load(census.mat);nian=1990:2010;p=polyfit(cdate,pop,2);z=polyval(p,cdate);zz=polyval(p,nian);my=mean(pop);mz=mean(z);R=sum(pop-my).*(z-mz)./sqrt(sum(pop-my).2).*sum(z-mz).2)plot(cdate,pop,*,cdate,z,r,nian,zz,b-)xlabel( (),F

41、ontName,FontSize,16)ylabel(,FontName,FontSize,16)title( 1790-1990,FontName,FontSize,20)text(1990, 248.73,-)legend(,)五、实验注意事项1、插值运算时应变量和自变量的确定2、曲线拟合时逼近函数的确定六、思考题由离散点数据来确定函数的方法有函数插值和曲线拟合两种:若测量值是准确的,没有误差,一般用 插值 ;若测量值与真实值有误差,一般用 拟合 。实验项目序号8 MATLAB基本绘图一、实验目的掌握绘制MATLAB二维、三维和特殊图形的常用函数;掌握绘制图形的辅助操作。二、实验基本知识1

42、、基本的绘图指令 绘制函数plot函数 函数能够将向量或者矩阵中的数据绘制在图形窗体中,并且可以指定不同的线型和色彩 同一个plot函数可以一次绘制多条曲线 基本语法格式 绘制一条曲线plot(xdata,ydata,color_linestyle_marker) 绘制多条曲线plot(xdata1,ydata1,clm1,xdata2,ydata2,clm2,)2、使用子图 子图的功能 可以将图形窗体划分为多个图形显示区域,每个区域彼此独立。 使用子图的方法 用subplot函数选择绘制区域 subplot函数的基本用法 subplot(m,n,p) m图形窗体分割的行数 n图形窗体分割的列数 p选定的窗体区域的序号 (以行元素优先顺序排列)3、格式化绘图 使用函数title添加图形标题 title(string) 标题将被自动地设置在轴的正中顶部 使用函数legend添加轴的图例 Legend(string1,string2,) MATLAB可以为不同的坐标轴添加不同的坐标信息 使用函数xlabel、ylabel、zlabel分别为图形窗体的X轴、Y轴和Z轴添加轴标签 xlabel(string) ylabel(string) zlabel(string) 使用函数text向图形窗体添加文本注释 text(x,y,string) x和y是文本注释添加的坐标值4、特殊图

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论