物流定量技术实验报告

《物流定量技术实验报告》由会员分享，可在线阅读，更多相关《物流定量技术实验报告（39页珍藏版）》请在装配图网上搜索。

1、【最新资料,Word版，可自由编辑！】目录实验一 Matlab操作基础实验二Matlab数值计算与分析实验三Matlab程序设计实验四Matlab绘图实验六线性回归分析实验一 Matlab操作基础一、实验目的通过实验熟悉Matlab软件平台，包括：（一）熟悉Matlab软件的用户环境；（二）学会Matlab用户环境的一般操作命令。二、主要操作及命令Matlab启动与退出、主窗口及命令窗口的使用、Matlab帮助系统的用法三、实验任务（一） Matlab启动与退出1、启动 Matlab利用快捷方式或运行 Matlab.exe启动Matlab，观察Matlab的界面组成，包括主窗口、命令窗口、工作

2、空间窗口、当前目录窗口、命令历史窗口和Start按钮Hi Edit 口创聊 QioMqp-皆耳申IpEl X * KiK 冇 P Q出 Matlab在主窗口 File 菜单中选择Exit Matlab命令；在Matlab命令窗口输入Exit或Quit命令； 单击Matlab主窗口的“关闭”按钮。（二）界面操作1、主窗口的用法在 Matlab 6.5 主窗口的菜单栏，共包含 File、Edit、View、Web Window和 Help 6个菜 单项。File菜单项项实现有关文件的操作，Edit菜单项对选中的当前窗口进行操作， View菜 单项用于设置Matlab集成环境的显示方式，Wet菜单项

3、用于设置 Matlab在Web上的链接、交流等操作，Window菜单项用于关闭所有打开的编辑器窗口，Help菜单项提供帮助信息，其中DEMO选项打开演示文档。File菜单：new选项，包含4个选项M_File、Figure、Model和GUI，单击他们可打开相应的编辑窗口； Ope n选项，单击他可打开一个 Ope n对话框，用户可以在对话框中选择相应的 文件，然后Matlab将用相应的编辑器打开该文件；Close选项，该选项后面将跟随某个打开 的视窗名，单击该选项将关闭该视窗；import data 选项，单击该选项将打开一个import对 话框，用户可以选择相应的数据文件，然后将该数据文件

4、中的数据导入到Matlab工作空间；Save workspace As选项,单击该选项将打开一个 Save to MAT_File对话框，用户需要为保存 的工作空间命名；Set Path选项，单击该选项将会打开设置路径对话框，通过该对话框可以 更改Matlab执行命令时搜索的路径；Page Setuppreferences 选项,设置页面布局、页眉、 页面所用的字体；Print选项，用于打印预定义好的页画内容，也可以设置一些打印参数；Print Selection选项，当选中命令窗口内的一部分内容后，该选项将处于檄活状态，此时单击该选 项，将打开打印对话框打印选中的内容；Exit Matlab

5、选项，单击该选项将关闭 Matlab。也可以通过快捷键CTRL+Q来关闭。View菜单：Desktop Layout选项，包含6个选顶，即默认布局、只显示命令窗口、简单布局、窄的历史记录窗口、高的历史记录窗口、同时显示5个窗口； Un dock Comma nd Win dow选项，选择该菜单将使命令窗口以新的窗口形式显示，而不是作为主窗口的一个浮动窗口形式显示；几个可选的复选框菜单，包括Comma nd Win do选项、Comma nd History选项、Curre ntdirectory 选顶、Workspace选项、Launch Pad选项、Profiler 选项和Help选项，选中

6、后会 打开相应窗口，再次选择就关闭。工具栏：Matlab 6.5主窗口的工具栏共提供了 10个常用命令按钮。这些命令按钮均有对应的菜单命令，但比菜单命令使用起来更快捷、方便，体现了window风格。2、命令窗口的使用命令窗口是Matlab的主要交互窗口，用于输入命令并显示除图形以外的所有执行结果。Matlab命令窗口中的“ ”为命令提示符，表示Matlab正在处于准备状态。在命令提示符后 键入命令并按下回车键后，Matlab就会解释执行所输入的命令，并在命令后面给出计算结果。 一般来说，一个命令行输入一条命令，命令行以回车结束。但一个命令行也可以输入若干条命 令，各命令之间以逗号分隔。例如：p

7、=15,m=35如果一个命令行很长，一个物理行之内写不下，可以在第一个物理行之后加上 3个小黑点 并按下回车键，然后接着下一个物理行继续写命令的其他部分。3个小黑点称为续行符，即把下面的物理行看作该行的逻辑继续。在Matlab里，有很多的控制键和方向键可用于命令行的编辑。请作如下练习：键入以下命令绘制正弦曲线和余弦曲线：x=0:0.5:360*pi/180;plot(x,s in (x),x,cos(x);键入以下命令求方程 3x 4+7x3 +9x 2-23=0 的全部根：p=3,7,9,0,-23; x=roots(p)3、工作空间窗口、当前目录窗口、搜索路径工作空间是Matlab用于存储

8、各种变量和结果的内存空间。在该窗口中显示工作空间中所 有变量的名称、大小、字节数和变量类型说明，单击鼠标右键可对变量进行编辑、 保存和删除当前目录是指Matlab运行文件时的工作目录，只有在当前目录或搜索路径下的文件、函 数可以被运行或调用。在当前目录窗口中可以显示或改变当前目录，还可以显示当前目录下的文件并提供搜索功能。将用户目录设置成当前目录也可使用cd命令。例如，将用户目录c:mydir设置为当前目录，可在命令窗口输入命令：cd c:mydir当用户在Matlab命令窗口输入一条命令后，Matlab按照一定次序寻找相关的文件。基本的搜索过程是：(1)检查该命令是不是一个变量。(2)检查该

9、命令是不是一个内部函数。(3)检 查该命令是否当前目录下的 M文件。 检查该命令是否Matlab搜索路径中其他目录下的 M 文件。用户可以将自己的工作目录列入 Matlab搜索路径，从而将用户目录纳入 Matlab系统统 一管理。设置搜索路径的方法有：(1)用path命令设置搜索路径。例如，将用户目录c:mydir 加到搜索路径下，可在命令窗口输入命令：path(path, c: mydir ) (2)用对话框设置搜索路径，在Matlab的File菜单中选Set Path命令或在命令窗口执行 pathtool命令，将出现搜 索路径设置对话框。通过 Add Folder或Add with Sub

10、folder命令按钮将指定路径添加到搜索 路径列表中。在修改完搜索路径后，则需要保存搜索路径。4、命令历史记录窗口、启动平台窗口和Start按钮命令历史记录窗口：在默认设置下，历史记录窗口中会自动保留自安装起所有用过的命令 的历史记录，并且还标明了使用时间，从而方便用户查询。而且，通过双击命令可进行历史命 令的再运行。如果要清除这些历史记录，可以选择Edit菜单中的Clear Comma nd History命令。Matlab 6.5的启动平台窗口，即查询便捷簿窗口，用来随时向用户说明软件中包含的内 容及功能，用户可以通过该窗口查看已经安装的各种工具箱。Matlab 6.5主窗口左下角的Sta

11、rt按钮，单击该按钮会弹出一个菜单，选择其中的命令 可以执行Matlab产品的各种工具，并且可以查阅 Matlab包含的各种资源。5、Matlab帮助系统进入帮助窗口的方法：(1)单击Matlab主窗口工具栏中的Help按钮。 在命令窗口中输 入 helpwin、helpdesk 或 doc。(3)选择 Help 菜单中的“ Matlab Help” 选项。MATLA帮助命令包括help、lookfor以及模糊查询。help命令，在MATLA6.5命令窗口 中直接输入help命令将会显示当前帮助系统中所包含的所有项目，即搜索路径中所有的目录 名称，同样，可以通过help加函数名来显示该函数的帮

12、助说明；lookfor命令，help命令只 搜索出那些关键字完全匹配的结果，lookfor命令条件比较宽松；what命令，用于当前路径文 件列表；who命令，用子当前内存变量列表；whos命令，显示当前内存变量详细信息；模糊查 询：MATLAB 6.0以上的版本提供了一种类似模糊查询的命令查询方法，用户只需要输入命令 的前几个字母，然后按Tab键，系统就会列出所有以这几个字母开头的命令。演示系统：在帮助窗口中选择演示系统(Demos)选项卡，然后在其中选择相应的演示模块， 或者在命令窗口输入Demos或者选择主窗口 Help菜单中的Demos子菜单，打开演示系统。 在演示系统中，我们可以获得有

13、关 Matlab的各种介绍。请选择子文件夹，如Desktopoverview ，run this demo ，学习主窗口的有关知识。远程帮助系统：通过 MathWorks公司的主页()，可以找到很多有用的信息。四、实验心得体会体会一：学习MATLAB件的相关操作，这样可以更好的运用到以后的实验中，简化学习 中的复杂问题，直接得到答案。运用 MATLA软件可以很好的简化一些中间环节，直接运用到 实际中。很大的程度上体现了计算机的简单的运用，解决实际生活中的问题。体会二：命令窗口输入命令，得到输入命令产生的结果。如果在命令窗口输入的错误命令, 那按下回车键后，就会有红色的提示语产生。并且有红色的指

14、线指示输入的命令中存在错误的 位置，这样提示用户修改错误。这样的情况下，很大程度上减少了用户的修改操作，很大程度 上方便了用户，这样的设计更人性化。实验二MATLAB值计算与分析一、实验目的熟悉MATLABS本命令与操作；熟悉MATLAB勺矩阵运算；学会MATLAB勺多项式运算，数据统计处理；学会用MATLA求解线性方程组。二、主要操作及命令打开软件MATLA操作平台利用帮助了解 inline , size , max , sum , numel , length , prod , mean , std , corrcoef , sort , conv , fmin的功能和用法。利用MATLA

15、B令窗口，输入实验项所需的命令；观察和理解命令执行的结果，修改出现错误的命令，直到窗口得到正确的实验结果。记录主要的实验项和结果，并记下实验心得。三、实验任务变量和数据操作1数的输入及运算女口： a=5回车以后显示的值是a=5,当直接输入数值后，显示的数值就是ans=5。输入数据是复数时，如果输入b=2-5i,显示的结果是b= 2.0000 - 5.0000i, 在窗口中输入数学运算式， 显示的结果是运算后的结果。例如输入：5*9输出的结果是：ans=45,输入函数代数式：log(9), 得到的结果是ans= 2.1972。2、数组输入例如，数组的输入：a=1,3,5,7,9,11 和a=1:

16、2:11得到的最后的值是一样的。a=1 3 5 7 9 11,相比之下，要得到一样的结果通过 a=1:2:11这种方式的输入，更快捷有效。3、变量的存储和调用变量的存储通过命令：save data a b存储变量a b至U data文件中，load data是将data文件中的变量加载到工作空间，whos是显示工作空间的所有变量的相关信息。Matlab矩阵及基本运算1矩阵的输入通过直接的赋值方式进行，例如：a=2,3,5;1,5,6;8,6,9得到的矩阵值在窗口以完整显示，再输入的a=(2,3)得到的数值是矩阵的第2列第3行的数值，输入a(7)得到的结果是以第 一列开始计数到第7位的数值。2、

17、矩阵大小测试和定位先输入矩阵值后，运用合适的函数表达式计算矩阵的相应的值。例如矩阵：a=4,5,6;7,8,9;9,10,11;12,16,18;得到标准的矩阵表达式，numel(a)=12求的值为整个矩阵的元素个数,n,m=size(a)的值为：n=5,m=3表示矩阵由4行3列构成,j,k=find(a10)表示从第一列开始查询满足不小于10的元素的位置;n=length(a);表示行数和列数的最大的值。3、矩阵的块操作在矩阵的整体中取出部分元素矩阵，a(2,:)，表示在矩阵中取出第2行的整个元素组； a(1,4,:)，表示在矩阵中取出第1行和第4行两行整体兀素；a(2:4,1:2)，表示取

18、矩阵中第 2行到4行和第1列到2列同时满足的元素。a(1,3,:)=a(3,1,:)表示在矩阵中将第1行和第3行将换位置得到的新矩阵；a(:,2,3)=a(:,3,2)使矩阵的2,3列互换。a(2,:)=4是将矩阵的第2行的值全部用4代替；a(2,:)=是表示删除a矩阵中的第2行；a(4,5)=7表示 aC4, 5)=7a =45600789000101100将矩阵构造成4*5的矩阵，未定元素为0,最后的元素为7; 12161307小矩阵构造大矩阵：先定义小矩阵a的值，再编写构成大矩阵，最后运行结果。a=4,5,6;7,8,9;9,10,11;12,16,18;构建大矩阵方式：a(1:3,2:

19、3),a(2:4,1:2);a,a(:,2)，显示的结果为：4、矩阵的翻转操作矩阵的上下翻转：flipud(a)实现矩阵中每列的数值倒序排列，矩阵a的左右翻转：fliplr(a) 显示为每行的数值倒序显示结果，矩阵的逆时针旋转90： rot90(a)，将矩阵的结果逆时针旋转90显示。 flipud(a) f liplr (a) rot got a)ans -ans -ans =45g121618654S9li187勒910119B791011T391110g58ID1612is45百IS36isq7125、特殊矩阵的产生单位矩阵的产生：a=eye，产生4维的单位矩阵，产生4X 5维1矩阵：a=

20、ones(4,5);产生3X4维0矩阵：A=zeros(3,4);随机nXm维矩阵：A=rand(n,m),(产生的元素是01之间的随机元素)。产生一个区间10, 20内均匀分布的4阶随机矩阵：x=10+(20-10)*rand(4);产生mix n的正态分布随即矩阵：randn(m, n);产生1m之间整数随机排列：randperm(m);产生魔方矩阵：magic(n);6、矩阵计算求矩阵的转置： a;求矩阵的行列式： det(b),(矩阵b是行列相等的方阵)；求矩阵的秩： rank(a);求矩阵的逆： inv(b)；求矩阵的迹： trace(a),(即是求矩阵 的对角线数值的和值)；求矩阵

21、与常数的运算：矩阵除以常数的值： a/8 ,矩阵与矩阵的运 算：矩阵直接的加减必须在同阶的基础上实现，矩阵之间的乘除就不需要要求是同维矩阵。同 矩阵的幕运算：2次幕运算： a.A2相当于 a.*a的值(这是Matlab中的特殊运算，称作 点运算。)。7、结构矩阵建立学生信息数据表对数据表结构元素进行修改、删除、增加、取用。建立数据表结构和初始数据：Student=struct(xm,李红,张三,汪峰,李杰,xh,0123,0124,0125,0126,sxcj,80,76,88,86);修改数据表中的数据值：Stude nt( 3)=struct(xm,刀郎,xh,0128,sxcj,96)这

22、样将表中的第三条数据替换，得到新的数据；数据表的删除：Student=rmfield(Student,xh)表示删除数据表中的xh这一列，得到新的数据表。(三) 数据统计处理与多项式计算1、数据统计处理首先输入数据：a= a=459,362,310,851;再对输入的数据进行运算，求平均值，中位数,方差，标准差，四分位极差，平均绝对偏差数：b=a(:)将数据改改变成一纵列数;mea n( b),media n( b),var(b),std(b),iqr(b),mad(b)直接得到所要求的值。2、多项式计算对多项式p=x4+2x3-5x+6和s=x2+2x+3求和、差、积、商。P=1,2,0,-

23、5,6,S=0,0,1,2,3,P+S=1,2,1,-3,9;P-S=1,2,-1,-7,3;P.*S=0,0,0,-10,18;P/S=0.5714;多项式乘积：conv(P,S)=0,0,1,4,7,1,-4,-3,18;多项式 P 的根：(四) 方程(组)的解法多项式 w= x2+2x+6;w=1,2,6;r=roots(w),求的结果为：图(1)求线性方程组Ax=b的解：R=rank(A) ； x=Ab求的结果为：图(2)(五) 函数极值求函数f(x)x cs(x)x x*log(x)在区间1)内的最小值。e首先新建一个 Ml文件，在新文件中编写：fx=inline( (x.A3+co

24、s(x)+x*log(x)/exp(x)x,fval=fmi nbn d(fx,0,1)运行命令得到答案：x=0.5223;fval=0.3974 四、实验心得体会体会一:在针对矩阵运算时，可以任意的取矩阵的不同位置的值，进行运算，但是有些运算必须是 要满足一定的运算规则才可以。 实验二主要是对矩阵的相关运用， 矩阵实用的范围更广，更有 学习价值。对于矩阵的计算，以前线性代数学习过，但是不借助计算机，很多复杂的运算很费 时，费力，在计算机上计算很便捷，几秒钟就能得到答案。运用到实际生活中将节省很多人力 物力财力，为企业降低成本。体会二：对于这次试验的建立表结构，以前在计算机课程 VF中学习过，

25、建立表单，建立数据库进 行查询。也学习过类似的输入过程，一种就是直接输入表全部的数据，二是，一条一条的输入 表的数据最终得到完整的数据表，但是相对比较复杂。只是每个软件输入、增加、修改等命令 不一样，但是方法基本都是相似的。体会三：在解决多项式问题上，以前学习相关章节就能感到，有些函数只知道解法，但是计算过程 很复杂，不容易计算得到答案。但是借助 Matlab就可以轻松的解决问题，很简洁快速的得到 答案，节省了很多过程时间。而且还可以根据函数的表达式得到相关的图形图线，更清楚的了解函数的波动规律。实验三MATLABg序设计、实验目的（一）熟悉关系运算、逻辑运算、M函数和M文件编制方法；（二）掌

26、握MATLAB言的几种循环、条件和开关选择结构；（三）学会编写一般性程序。二、主要操作及命令（一）开启软件平台一MATLAB（二） 在D盘上以自己的学号，女口 860101、860110 （前者代表1班学号为1号的同学；后 者代表1班学号为10号的同学）建立一个文件夹，并利用文件菜单中的“ SETPATH”选项, 把自己的文件夹加入到 MATLA 的搜索路径中（注意保存），改变当前路径为自己的文件夹；（三）利用帮助了解命令disp、num2str、fix、break、int2Str、rem、var等语句的用法 及功能；（四） 利用MATLA文本编辑窗口 ,编写和调试实验项程序；注意给实验程序（

27、尤其是思考与 练习中自编的程序）取一个易记忆的名称，将程序代码及运行结果保存到自己的WOR文档, 准备平时的实验检查；（五）观察和理解程序执行的结果，若出现错误，修改、直到输出正确结果；（六）记录23个你认为重要的实验项及结果，13条学习心得体会。三、实验任务(一) 关系运算与逻辑运算数与数、矩阵与矩阵之间的关系运算，其中主要的关系运算符包括：大于、=大于等于、=7，返回的值就是：ans=0;逻辑运算包括：与(&)、或(|)、非()、异或(xor)，变量中非零数逻辑量为“真”，逻 辑运算的结果还是，如果是真就显示值为 1,逻辑假显示结果为0。例如：a=3&1得到值为： a=1，a=0&1得到值

28、为0; a=-1|5得到值是：1，而或运算中只有同时为零时才显示结果为a=0;异或运算中，只有一个是非零，一个是零显示值才是1,其余的全部为0。a=xor(3,0)和a=xor(0,5)显示值为：a=1; a=xor(3,1)和a=xor(0,0)显示的值为：a=0;非()运算是将真 变假，假变真。a=0, 0是假，但是显示的变化后显示的结果为真， a=1, 1为真变化后显示 的结果为假。如果是矩阵中运算 a=1:5 ,b=3:7 ,t1=(a-22)=0 ,0,0, 0,1 ; t2=(a3)-b=-3,-4 , -5 , -5 , -6 ; t3=t1-2=-1, -1 , -1 , -1

29、 , -2。(二) M文件1、M文件的建立创建新M文件，在M文件中输入以下命令：function v,s=ff(v0,a,t)%定义函数名和输入、输出变量v=vO+a*t;%合出输入、输出变量之间的关系s=v0*t+a*t2/2;命令窗口中输入：v,s=ff(2,5,10)得到的答案是：v=52, s=270;用函数命令inline也可以实现M文件的部分功能，例如： fu=i nli ne(2*xA2+3*x+l)显示的结果为：其他的为子函数，子函数只主函数必须在最前面，后一个M文件中可以包括一个以上的函数， 其中有个是主函数, 能被同一文件的的主函数调用，不能被不同的文件的其他函数调用 面接

30、子函数，秩序没有限制。在M文件中输入：fun ctio n out=f un c5(x)recip=reciproc(x);out=sum(recip);function output=reciproc(i nput) output=1./i nput;命令窗口输入：fun c5(2,2,5)得到的结果是：ans =1.2000(三) 程序流程控制1、if条件语句创建M文件，输入下列程序，比较两数的大小a=4;b=3;if ab;max=a;disp(最大值为：,num2str(max);elseif a=b;max=两数相等;else ;max=b;disp(最大值为：，num2str(ma

31、x);end最终得到的结果为：最大值为：42、for循环语句利用for循环求出100-200之间的所有素数。for m=101:2:200k=fix(sqrt(m);for i=2:k+1if rem(m,i)=0;break;end;end;if i=k+1disp(i nt2Str(m)end;end;得到的答案是：101，103，107，109，113，127，131，137, 139, 149，151，157, 163,167, 173，179，181，191，193，197，199；建立M文件:A=2,5,0;2,4,7 ; 3,1,9;for i=Ay=i-i(2)+i(3)end

32、得到的答案是：丫=3, y=2, y=2;3、while 循环 sum=0;i=1;while isubplot(2,2,1) fplot(cos(x),0,2*pi) subplot(2,2,2) fplot(si n(x-pi/2) ,0,2*pi) subplot(2,2,3) fplot(xA2*cos(x-pi),0,2*pi) subplot(2,2,4) fplot(exp(si n(x), 0,2*pi)(二)三维绘图1、三维曲线图格式：plot3(x,y,z,s), 其中x, y, z为均为同维向量，s为控制字符 color-li nestyle-marker。例如:t=0:

33、pi/60:10*pi;x=si n( t);y=cos(t);plot3(x,y,t,*-b)得到螺旋图为：格式：fill3(x,y,z,c);例如：x=2 1 2 ; 9 7 1;6 7 0;y=1 7 0;4,7,9;0,4,3; z=18 6;7 9 6; 1 6 1; c=1 0 0; 0 1 0; 0 0 1; fill3(x,y,z,c); grid on得到的图形为:2、三维网格图格式：mesh(x,y,z,c)%颜色由c指定的三维网格图meshc(x,y,z,c)%带有底座的三维网格图meshz(x,y,z,c)%带有等咼线的三维网格图x=1:4;y=1:5;X,Y=mesh

34、grid(x,y) ，其中X,Y是x-y平面中矩阵X，Y对应的横纵轴坐标值。练习：绘出带有底座的马鞍面,42252x=-8:8;y=-8:8;X,Y=meshgrid(x,y)Z=(X.A2/4A2-Y.A2/5A2)meshz(X,Y,Z)得到图：3、三维曲面图格式：surf(X,Y,Z,C)%画出颜色由C指定的三维曲面图surfc(X,Y,Z,C)%画出带有等高线的三维曲面图练习：绘出带有等高线的理想气体状态方程曲面，pv=nRT n=2mol:R=8.31; n=2;p=(1:20)*1e5;v=(1:20)*1e-3P,V=meshgrid(p,v) ;T=P.*V/n/Rsurfc(

35、P,V,T) 得到的图4、三维饼图格式：pie3(X,explode)x=1 3 0.5 2.5 2;ex=0 1 0 0 0;pie3(x,ex)绘出函数fx2 ey x，在-2，2，-2，2上的三维曲面图，命令:x=-2:2;y=-2:2X,Y=meshgrid(x,y)f=(X.A2+exp(Y)*abs(X)surfc(X,Y,f)得出图1:2005年4月我国豆油及其分离品进口量统计的三维饼图为图2,命令为：ex=0 0 0 0 1 0 ;pie3(x,ex)四、实验心得体会体会一：函数问题的绘图，更好的理解相关函数的一定的曲线走势，使问题解决更加简单便捷。这 样的绘图功能可以减少人工

36、绘图的工作，还能减少人工绘图出现的误差，更精确的绘制出所需要的图像。一些相关图可以在一个窗口中绘制图像，或者同个图画上操作，这样使相应比较的问题更能找出差别的地方，减少了人工比较的不准确性。本节实验主要是面对一定要求的绘图，了解和熟悉相关基本绘图的操作命令，绘制简单的 图形。按照一定的需要绘制图形，使阅读者更直观的了解图形中反应的相关信息，更形象的描 述一些问题。这样的绘图功能能让说明的人更清楚明白的表述相关信息， 减小抽象说明问题的 缺陷，更让人理解相关问题。体会二:在绘制饼形图和条形图中，可以更直观的了解相关组成部分的比例和相应的多少，绘制饼形图时，软件直接就会得出整体的构成部分各个所占比

37、例，更清楚的了解相关信息，条形图能在各个组成部分之间比较多少关系。三维绘图中可以将二维图形不能解决更进一步的加强，是 一些多维问题更形象的表现出来，立体形态的表示出来。三维坐标图将多个关系之间的数据加 以运用，体现在三维图上，更方便的解决生活和学习中的一些问题。实验六线性回归分析一、实验目的（一）学会MATLAB件中曲线拟合的方法，了解最小二乘拟合的基本原理和方法，掌握用 MATLAB作最小二乘多项式拟合的方法；（二）通过实例，学习如何用拟合方法解决实际问题。二、主要操作及命令（一）开启软件平台一MATLAB（二） 在D盘上以自己的学号建立一个文件夹，并利用文件菜单中的“SET PATH”选项

38、， 把自己的文件夹加入到 MATLA 的搜索路径中（注意保存），改变当前路径为自己的文件夹；（三）利用帮助了解命令polyfit 、polyval、polyconf、regress等语句的功能及用法。?参考实验指导书，根据你所完成的实验撰写本页；建议采用A4纸双面打印试验报告(四) 利用MATLA文本编辑窗口 ,编写和调试实验项程序；注意给实验程序(尤其是思考与练习中自编的程序)取一个易记忆的名称，将程序代码及运行结果保存到自己的WOR文档,准备平时的实验检查；(五) 观察和理解命令或程序执行的结果；(六) 记录12个重要的实验项及结果，12条学习心得体会。三、实验任务1、某市电子工业公司有1

39、4个所属企业，各企业的年设备能力与年劳动生产率统计数据如下表。试分析企业年设备能力与年劳动生产率的关系。若该公司计划新建一个设备能力为9.2千瓦/人的企业，估计劳动生产率将为多少？根据给出的数据判断两个数据之间的关系，再根据数据关系利用相关函数特性，建立数据模型。x=2.8 2.8 3.0 2.9 3.4 3.9 4.0 4.8 4.9 5.2 5.4 5.5 6.2 7.0;y=6.7 6.9 7.2 7.3 8.4 8.8 9.1 9.8 10.6 10.7 11.1 11.8 12.1 12.4scatter(x,y);%在图形中观察两变量之间的大致关系R=corr(x,y);%计算两变

40、量之间的关系度根据图形关系和变量之间的关系程度利用一元一次线性回归解决问题：Sxx=sum(x-mea n( x).A2);Syy=sum(y-mea n( y).A2);Sxy=sum(y-mea n( y).*(x-mea n( x)计算得到的线性系数：a=Sxy/Sxx,b=mea n(y)-mea n( x).*a得到的结果为a=1.4481 ; b=3.1003建立函数：Y=a*X+b若企业新建一个设备能力为9.2千瓦/人的企业，估计劳动生产率是多少？可以得 x=9.2计算 y的值。 -4 -1-将x的值代入，得到的函数值为16.4233，劳动生产率为16.4233千元/人。2、某建

41、筑材料公司的销售量因素分析下表中的数据是某建筑材料公司去年 20个地区的销售量(丫，千万)。推销开支、实际帐 目数、同类商品竞争数和地区销售潜力分别是影响建筑材料销售量的因素。试建立多元回归模型，并假设一组数据X进行销售量预测。根据数据表的数据输入相关数据，根据数据关系建立数据模型计算相关数据系数。首先建立M文件，再在文件中编写程序命令，计算得到数据结果。编写的M文件命令为：通过M文件运行得到的数据结果为：分析结果为:b。=191.9158 b o 的置信区间为(103.1071,280.7245)b1=-0.7719b1的置信区间为(-7.1445,5.6007)b2=3.1725b2的置信

42、区间为(2.0640,4.2809)b3=-19.6811b3 的置信区间为(-25.1651,-14.1972)b4=-0.4501b4的置信区间为(-3.7284,2.8283)r2=0.9034,F=35.0509,p=0.0000p0.05，回归模型y=191.9158-0.7719 x 1+3.1725x2-19.6811x 3-0.4501x 4假设数据为X=5.4,60,8,10预测得到的数据Y为：y=225.9966四、实验心得体会体会一：再根如果在构建一元一次线性回归的过程中，首先观察两变量之间在图形中呈现出来的关系,据函数关系模型特征去就算模型与实际值之间的相似度。判断是否选择此类模型去构建相似度达到要求，则去运算存在这种关系之间的系数值，计算出模型的表达关系式。体会二:虽然一元线性回归分析可以解决一定的实际问题，但是由于这样的关系之间只有两个变量相联系。但是实际生活中很多情况下都不只是一个变量去影响一个变量，很多情况下都是多个因素影响其他的变量，这就会有多个变量同时影响因变量。有些多变量与因变量之间存在关系， 但是有些不存在线性关系。