matlab设计作业.doc
matlab作业
上机实践练习五实验五：MATLAB 数据可视化实验目的：掌握MATLAB 二维、三维图形绘制，掌握图形属性的设置和图形修饰；掌握图像文件的读取和显示。
学时：2学时实验内容：(1)二维图形绘制。
(2)隐函数绘图(2) 三维曲线和三维曲面绘制。
(3) 图像文件的读取和显示。
1.写出实现下列图形的M文件。
t=[0:0.1:pi];y1=sin(2*t);y2=cos(3*t);plot(t,y1,'b-o');hold onplot(t,y2,'k:o');xlabel('ʱ¼ä')ylabel('·ùÖµ')text(1,sin(2*1),'\fontsize{16}\leftarrowsin(2t) ');text(2.4,cos(3*(2.4)),'\fontsize{16}cos(3t)\rightarrow ',... 'HorizontalAlignment','right')2. 将窗口分割成4个区域，并且在[]0,2π区间上绘制1sin 2,6y x π⎛⎫=+⎪⎝⎭()224cos 3,y x =()3cos y x =，()35sin 1y x =+四条曲线，并且给每一个图形添加标题和标注。
clearx=0:pi/10:2*pi; y1=sin(2.*x+pi/6); y2=4.*cos(3.*x.^2); y3=cos(x); y4=5.*sin(x+1); subplot(2,2,1); plot(x,y1); title('y1');text(1,sin(2*1+pi/6),'sin(2*x+pi/6)') subplot(2,2,2); plot(x,y2)title('y2');text(2,4.*cos(3*2^2),'4*cos(3*x^2)') subplot(2,2,3); plot(x,y3)title('y3');text(2,cos(2),'cos(x)')subplot(2,2,4);plot(x,y4)title('y4');text(2,5.*sin(3),'5*sin(x+1)')3.已知y1=x2，y2=cos(2x)，y3＝y1*y2，其中x为取值-2π~2π的等差数列(每次增加0.02π)，完成下列操作(1)在同一坐标系下用不同的颜色和线型绘制三条曲线，给三条曲线添加图例。
MATLAB课程作业
注意：11月15日是交作业的最后期限。
（交word和MATLAB程序清单，word以“班级-学号-姓名”命名，每个MATLAB程序以”zuoye_题号”命名后，放在一个文件夹下（文件夹命名为“作业程序”），如zuoye_1_1.m），将word和“作业程序”文件夹放在一个文件夹（文件夹以“班级-学号-姓名”）交到班长处，班长统一交到教学办公室。
第一部分程序设计1.1、用MA TLAB可以识别的格式输入下面两个矩阵：(1) 矩阵A的维数；(2) 矩阵A中的元素a41的值；(3) 修改矩阵A的元素，使a41 =3.0；(4) 矩阵A中最后2行和最后3列交汇形成的子矩阵的值。
(5)求出A和B的乘积矩阵C，并将C矩阵的右下角2X 3子矩阵赋给D矩阵。
1.2、已知111121111,131,111214A B⎡⎤⎡⎤⎢⎥⎢⎥=-=-⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求(1)AB-2A，(2)A*B，(3)A﹒*B，(4) AB-BA 1.3、解线性方程：1.4、解方程组：1.5、用MA TLAB语言实现下面的分段函数：1.6、已知x=[1 2 3 ],y=[4 5 6 ]，试计算z=x.*y 、x.\y 和x./y.1.7、分别用for 和while 循环语句编写程序，求出1.8、 已知在平面坐标中两点(x 1, y 1)和(x 2, y 2)之间的距离计算公式为 ()()222121y y x x L -+-=(1) 利用命令文件的形式，编写求解该距离的M 文件dis1.m ；(2) 利用函数文件的形式，编写求解该距离的M 文件dis2.m ；(3) 给定两点坐标的值(2,3)和(8, －5)，试分别调用命令文件dis1.m 和函数文件dis2.m 求解该两点间距离的值。
1.9、求解方程x 5+6x 3一3x 2＝10的5个根，并将其位置用五角星符号标记在复平面上，要求横纵坐标袖的刻度等长，注明虚轴和实轴，在title 位置上写㈩方程。
matlab作业
实验三某厂利用a、b、c三种原料生产A、B、C三种产品，已知生产每种产品在消耗原料方面的各项技术条件和单位产品的利润，以及可利用的各种原料的量，是制定适当的生产规划使得该厂的总利润最大。
对于这个问题，可先将该问题转化为数学模型。
Z=2x1+4x2+3x3X1+3x2+2x3<=802x1+x2+2x3<=403x1+4x2+2x3<=60X1,x2,x3>=0据此可以画出图形，可依靠图形来求出Z的最大值。
c=[-2 -4 -3]A=[1 3 2;2 1 2;3 4 2]b=[60 40 80]lb=[0 0 0][x,f,exitflag,output]=linprog(c,A,b,[],[],lb,[])将上述程序输入matlab。
就可得出结果。
结果如下：x =4.000012.000010.0000f =-86.0000exitflag =1output =iterations: 7cgiterations: 0algorithm: 'lipsol'由上述结果可以知道当A、B、C分别生产4,12,10 公斤时，总利润最大，可达到86。
某工厂生产甲、乙、丙三种产品，单位产品所需工时分别为2、3、1个；单位产品所需原料分别为3、1、5公斤；单位产品利润分别为2、3、5元。
工厂每天可利用的工时为12个，可供应的原料为15公斤。
为使总利润为最大，试确定日生产计划和最大利润。
对于这个问题，可先将该问题转化为数学模型。
Z=2x+3y+5z2x+3y+z<=123x+y+5z<=15据此可以画出图形，可依靠图形来求出Z的最大值。
c=[-2 -3 -5]A=[2 3 1;3 1 5]b=[12;15]lb=[0 0 0][x,Z,exitflag,output]=linprog(c,A,b,[],[],lb,[])将上述程序输入matlab。
就可得出结果。
MATLAB作业4-1
MATLAB作业4-1MATLAB 作业41、⽤2sin(103)y t =+在（0，3）区间内⽣成⼀组较稀疏的数据，并⽤⼀维数据插值的⽅法对给出的数据进⾏曲线拟合，并将结果与理论曲线相⽐较。
解：2、⽤242231(,)sin()3x y f x y e xy x y x y--=++原型函数⽣成⼀组⽹络数据或随机数据，分别拟合出曲⾯，并和原曲⾯进⾏⽐较。
解：⽹格数据：随机数据：x∈-区间内的光滑函数曲线，⽐较各种插值算法3、假设已知⼀组数据，试⽤插值⽅法绘制出(2,4.9)>> x=[-2,-1.7,-1.4,-1.1,-0.8,-0.5,-0.2,0.1,0.4,0.7,1,1.3,...1.6,1.9,2.2,2.5,2.8,3.1,3.4,3.7,4,4.3,4.6,4.9];y=[0.10289,0.11741,0.13158,0.14483,0.15656,0.16622,0.17332,...0.1775,0.17853,0.17635,0.17109,0.16302,0.15255,0.1402,...0.12655,0.11219,0.09768,0.08353,0.07019,0.05786,0.04687,...0.03729,0.02914,0.02236];解：x y在(0.1,0.1)~(1.1,1.1)区域内的点进⾏插值，并⽤三维曲⾯4、假设已知实测数据由下表给出，试对(,)>> [x,y]=meshgrid(0.1:0.1:1.1);z=[0.83041,0.82727,0.82406,0.82098,0.81824,0.8161,0.81481,0.81463,0.81579,0.81853,0.82304 ;0.83172,0.83249,0.83584,0.84201,0.85125,0.86376,0.87975,0.89935,0.92263,0.94959,0.9801; 0.83587,0.84345,0.85631,0.87466,0.89867,0.9284,0.96377,1.0045,1.0502,1.1,1.1529;0.84286,0.86013,0.88537,0.91865,0.95985,1.0086,1.0642,1.1253,1.1904,1.257,1.3222;0.85268,0.88251,0.92286,0.97346,1.0336,1.1019,1.1764,1.254,1.3308,1.4017,1.4605;0.86532,0.91049,0.96847,1.0383,1.118,1.2046,1.2937,1.3793,1.4539,1.5086,1.5335;0.88078,0.94396,1.0217,1.1118,1.2102,1.311,1.4063,1.4859,1.5377,1.5484,1.5052;0.89904,0.98276,1.082,1.1922,1.3061,1.4138,1.5021,1.5555,1.5573,1.4915,1.346;0.92006,1.0266,1.1482,1.2768,1.4005,1.5034,1.5661,1.5678,1.4889,1.3156,1.0454;0.94381,1.0752,1.2191,1.3624,1.4866,1.5684,1.5821,1.5032,1.315,1.0155,0.62477;0.97023,1.1279,1.2929,1.4448,1.5564,1.5964,1.5341,1.3473,1.0321,0.61268,0.14763];解：5、习题3和4给出的数据分别为⼀元数据和⼆元数据，试⽤分段三次样条函数和B样条函数对其进⾏拟合。
MATLAB作业
MATLAB作业⼀、必答题：1. MATLAB系统由那些部分组成？答：MATLAB系统主要由开发环境、MATLAB语⾔、MATLAB数学函数库、图形功能和应⽤程序接⼝五个部分组成。
2. 如何启动M⽂件编辑/调试器？答：在操作界⾯上选择“建⽴新⽂件”或“打开⽂件”操作时，M⽂件编辑/调试器将被启动。
在命令窗⼝中键⼊“edit”命令也可以启动M⽂件编辑/调试器。
3. 存储在⼯作空间中的数组能编辑吗？如何操作？答：存储在⼯作空间的数组可以通过数组编辑器进⾏编辑：在⼯作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输⼊修改内容即可。
4. 在MATLAB中有⼏种获得帮助的途径？答：在MATLAB中有多种获得帮助的途径：（1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器；（2）help命令：在命令窗⼝键⼊“help” 命令可以列出帮助主题，键⼊“help 函数名”可以得到指定函数的在线帮助信息；（3）lookfor命令：在命令窗⼝键⼊“lookfor 关键词”可以搜索出⼀系列与给定关键词相关的命令和函数（4）模糊查询：输⼊命令的前⼏个字母，然后按Tab键，就可以列出所有以这⼏个字母开始的命令和函数。
5. 有⼏种建⽴矩阵的⽅法？各有什么优点？答：（1）以直接列出元素的形式输⼊；（2）通过语句和函数产⽣；（3）.在m⽂件中创建矩阵；（4）从外部的数据⽂件中装⼊。
6. 命令⽂件与函数⽂件的主要区别是什么？答：命令⽂件: M⽂件中最简单的⼀种，不需输出输⼊参数，⽤M ⽂件可以控制⼯作空间的所有数据。
运⾏过程中产⽣的变量都是全局变量。
运⾏⼀个命令⽂件等价于从命令窗⼝中顺序运⾏⽂件⾥的命令，程序不需要预先定义，只要依次将命令编辑在命令⽂件中即可。
函数⽂件：如果M⽂件的第⼀个可执⾏⾏以function开始，便是函数⽂件，每⼀个函数⽂件定义⼀个函数。
电气工程软件训练(三)——Matlab 作业
D1(江苏大学《电气工程软件训练三》课程设计报告设计题目：MATLAB专业班级：J电气1401学生姓名：唐鹏学生学号：4141127007指导老师：完成日期：江苏大学京江学院一MATLAB课程设计的目的和要求1.MATLAB软件功能简介MATLAB的名称源自Matrix Laboratory，1984年由美国Mathworks公司推向市场。
它是一种科学计算软件，专门以矩阵的形式处理数据。
MATLAB将高性能的数值计算和可视化集成在一起，并提供了大量的内置函数，从而被广泛的应用于科学计算、控制系统和信息处理等领域的分析、仿真和设计工作。
MATLAB 软件包括五大通用功能，数值计算功能（Nemeric）、符号运算功能（Symbolic）、数据可视化功能（Graphic）、数字图形文字统一处理功能（Notebook）和建模仿真可视化功能（Simulink）。
其中，符号运算功能的实现是通过请求MAPLE内核计算并将结果返回到MATLAB命令窗口。
该软件有三大特点，一是功能强大；二是界面友善、语言自然；三是开放性强。
目前，Mathworks公司已推出30多个应用工具箱。
MATLAB在线性代数、矩阵分析、数值及优化、数值统计和随机信号分析、电路与系统、系统动力学、次那好和图像处理、控制理论分析和系统设计、过程控制、建模和仿真、通信系统以及财政金融等众多领域的理论研究和工程设计中得到了广泛应用。
2.MATLAB课程设计的目的本次课程设计主要是为了使学生了解MATLAB软件的基本知识，熟悉MATLAB的上机环境，掌握MATLAB数值运算、程序设计、二维/三维绘图、符号运算、Simulink仿真等相关知识，并初步具备将一般数学问题转化为对应的计算机进行处理的能力，以便为今后进一步的学习打下坚定基础。
二MATLAB课程内容1 MATLAB语言基础实验目的：基本掌握MATLAB 向量、矩阵、数组的生成及其基本运算(区分数组运算和矩阵运算)、常用的数学函数。
matlab课程设计大作业
matlab课程设计大作业一、教学目标本课程的教学目标是使学生掌握MATLAB基本语法、编程技巧以及MATLAB 在工程计算和数据分析中的应用。
通过本课程的学习，学生将能够熟练使用MATLAB进行简单数学计算、线性方程组求解、函数图像绘制等。
1.掌握MATLAB基本语法和编程结构。
2.了解MATLAB在工程计算和数据分析中的应用。
3.熟悉MATLAB的函数库和工具箱。
4.能够使用MATLAB进行简单数学计算。
5.能够使用MATLAB求解线性方程组。
6.能够使用MATLAB绘制函数图像。
7.能够利用MATLAB进行数据分析和处理。
情感态度价值观目标：1.培养学生对计算机辅助设计的兴趣和认识。
2.培养学生团队合作和自主学习的能力。
二、教学内容本课程的教学内容主要包括MATLAB基本语法、编程技巧以及MATLAB在工程计算和数据分析中的应用。
1.MATLAB基本语法：介绍MATLAB的工作环境、基本数据类型、运算符、编程结构等。
2.MATLAB编程技巧：讲解MATLAB的函数调用、脚本编写、函数文件编写等编程技巧。
3.MATLAB在工程计算中的应用：介绍MATLAB在数值计算、线性方程组求解、图像处理等方面的应用。
4.MATLAB在数据分析中的应用：讲解MATLAB在数据采集、数据分析、数据可视化等方面的应用。
三、教学方法本课程采用讲授法、案例分析法、实验法等多种教学方法相结合的方式进行教学。
1.讲授法：通过讲解MATLAB的基本语法、编程技巧以及应用案例，使学生掌握MATLAB的基本知识和技能。
2.案例分析法：通过分析实际工程案例，使学生了解MATLAB在工程计算和数据分析中的应用。
3.实验法：安排上机实验，使学生在实际操作中巩固所学知识，提高实际编程能力。
四、教学资源本课程的教学资源包括教材、实验设备、多媒体资料等。
1.教材：选用《MATLAB教程》作为主要教材，辅助以相关参考书籍。
2.实验设备：为学生提供计算机实验室，配备有MATLAB软件的计算机。
MATLAB作业6(2016)作业本
MATLAB 作业61、在图形绘制语句中，若函数值为不定式NaN ，则相应的部分不绘制出来，试利用该规律绘制sin()z xy =的表面图，并剪切下2220.5x y +≤的部分。
解：>> [x,y]=meshgrid(-1:.1:1); z=sin(x.*y);ii=find(x.^2+y.^2<=0.5^2); z(ii)=NaN; surf(x,y,z)2、试求解下面的无约束最优化问题。
222222211433222424min 100()(1)90()(1)10.1[(1)(1)]19.8(1)(1)xx x x x x x x x x x -+-+-+-+-+-+--解：>> f=inline(['100*(x(2)-x(1)^2)^2+(1-x(1))^2+',...'90*(x(4)-x(3)^2)+(1-x(3)^2)^2+',...'10.1*((x(2)-1)^2+(x(4)-1)^2)+',...'19.8*(x(2)-1)*(x(4)-1)'],'x');x=fminunc(f,ones(7,1))x =1.0e+002 *0.105464467987131.112320667672340.06782323911149-1.115047464577260.010000000000000.010000000000000.010000000000003、 试用图解法求解下面的非线性规划问题，并用数值求解算法验证结果。
321211221212min (44)20..100,0x x x x x s t x x x x +-+-+≥⎧⎪-+-≥⎨⎪≥≥⎩ 解：>> [x1,x2]=meshgrid(0:0.02:1,1:0.02:2);z=x1.^3+x2.^2+4*x1+4;ii=find(x1-x2+2<0); z(ii)=NaN;ii=find(-x1.^2+x2-1<0); z(ii)=NaN;ii=find(x1<0); z(ii)=NaN; ii=find(x2<0); z(ii)=NaN;surf(x1,x2,z)function [c,ce]=exc6f4(x)ce=[];c=[x(1)^2-x(2)+1];>> f_opt=inline('x(1)^3+x(2)^2+4*x(1)+4','x');A=[-1 1]; B=2; Aeq=[]; Beq=[]; xm=[0;0];x=fmincon(f_opt,[0;1],A,B,Aeq,Beq,xm,[],'exc6f4');4、 试求解此线性规划问题：6712341236723571,2,,7min 421..390x x x x x x x x x x x s t x x x x x ++++=⎧⎪-+--+=⎪⎨+++=⎪⎪≥⎩解：>> [x1,x2]=meshgrid(0:0.02:1,1:0.02:2);z=x1.^3+x2.^2+4*x1+4;ii=find(x1-x2+2<0); z(ii)=NaN;ii=find(-x1.^2+x2-1<0); z(ii)=NaN;ii=find(x1<0); z(ii)=NaN; ii=find(x2<0); z(ii)=NaN;surf(x1,x2,z)function [c,ce]=exc6f4(x)ce=[];c=[x(1)^2-x(2)+1];>> f_opt=inline('x(1)^3+x(2)^2+4*x(1)+4','x');A=[-1 1]; B=2; Aeq=[]; Beq=[]; xm=[0;0];x=fmincon(f_opt,[0;1],A,B,Aeq,Beq,xm,[],'exc6f4');5、 试求解下面的二次型规划问题，并用图示的形式解释结果。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1.在同一个图形窗口内画出衰减震荡曲线及其包络线,取值范围是[0,6pi]. 代码：
t=(0:pi/100:6*pi);
y1=exp(-3*t).*cos(t/2);
y2=exp(-3*t);
t3=pi*(0:9)/9;
y3=exp(-3*t3);
plot(t,y1,'r--',t,y2,'b',t3,y3,'bo')
2.画出所表示的三维曲面。
X,y 的取值范围是[-9,9]。
x=linspace(-9,1,9);
y=linspace(-9,1,9);
[x,y]=meshgrid(x,y);
z=cos((sqrt(2*x.^2+2*y.^2))/sqrt(x.^2+y.^2));
surf(x,y,z);
3.求和当n =100 时的值。
sum=0;
for i=1:100
sum=sum+1/factorial(i);
end
disp(sum);
1.7183
4. 求1000 个元素的随机数向量A 中大于0.5 的元素个数。
A=rand(1,1000);
count=0;
for i=1:length(A)
if A(i)>0.5
count=count+1;
end
end
count
count =
514
5.画正态分布2 的概率密度函数曲线，产生个相应的随机数，N (1,4 m 10000 画出直方图和带正态密度曲线的直方图。
将随机数的频率曲线与概率密数曲线画在一起进行比对。
x=-20:20;
y=normpdf(x,1,4);
a=1+4*randn(1,10000);
[N,h]=hist(a);
y=h(1)-h(2);
subplot(2,2,1);plot(x,y);
subplot(2,2,2);bar(h,N/(10000*y),1);
subplot(2,2,3);plot(x,y);
hold on;
bar(h,N/(10000*y),1);
6.对不同的参数n ，画出2
()n χ分布的概率密度函数曲线，讨论n 的不同变化对曲线的影响。
x=0:0.1:30; y1=chi2pdf(x,2) plot(x,y1,'r') hold on
y2=chi2pdf(x,4); plot(x,y2,'b'); y3=chi2pdf(x,6); plot(x,y3,'y'); axis([0,30,0,0.2])
7.求2
2
()(9)n χχ=分布的双侧临界值并画图（变换不同的α）。
alpha=0.5:0.01:0.975; lambda=icdf('chi2',alpha,9); plot(alpha,lambda)
设计题II ：国内产值与人力资本的回归分析
改革开放以来我国经济高速增长。
研究表明，国内生产总值与诸多因素有关.下表是1995-2005 年各季度的国内生产总值和人力资本的数据。
试确定国内生产值和人力资本的关系并进行统计检验和分析。
年份季度人力资本(X)国内产值(Y)年份季度人力资本(X)国内产值(Y) 1995春324.0 14410.3 2000秋618.9 24463.5 1995夏325.0 14660.5 2000冬640.7 24803.7 1995秋348.6 15021.0 2001春656.2 25621.5 1995冬352.9 15198.5 2001夏683.9 26428.8 1996春373.6 15679.1 2001秋710.5 27011.0 1996夏391.1 16478.7 2001冬764.3 27413.8 1996秋400.6 17185.9 2002春788.7 27887.0 1996冬417.9 17794.2 2002夏815.4 29186.3 1997春423.9 18262.3 2002秋838.7 29885.5 1997夏441.2 18761.3 2002冬872.9 30083.2 1997秋456.4 19188.0 2003春886.6 31605.4 1997冬465.6 19743.3 2003夏921.1 32363.1 1998春473.5 20036.4 2003秋941.4 33121.0 1998夏488.4 20365.5 2003冬962.7 33955.7 1998秋498.9 20610.9 2004春1001.4 34878.4 1998冬508.1 21100.6 2004夏1038.5 36710.4 1999春521.6 21513.5 2004秋1076.7 38121.5
设计要求：
1、输入数据，画出相应图形；
2、建立线性回归模型, 对给定的数据拟合回归方程。
3、给出回归系数的估计与检验，画出散点图和回归直线图、残差图。
4、给出方差分析，做出方差分析表，对结果进行简单分析。
5、进行假设检验分析，数据分布的推测；
解答：
（1）设计题目
国内产值与人力资本的回归分析
（2）设计目的
上表是1995-2005 年各季度的国内生产总值和人力资本的数据。
改革开放以来我国经济高速增长。
研究表明，国内生产总值与诸多因素有关.根据这些数
据确定国内生产值和人力资本的关系并进行统计检验和分析。
在分析计算这些
数据的同时重在为了更好的了解概率论与数理统计的知识，熟练掌握概率论与
数理统计在实际问题上的应用，并将所学的知识结合matlab对数据的处理解决
实际问题。
本设计是利用一元线性回归理论对国内产值与人力资本关系建立数
学模型，并用matlab分析工具库中的回归分析软件进行解算并检验。
（3）设计原理
本题是一道确定生产总值和人力资本的关系问题，首先做出该组数据的散点图，由图分析该数据属于线性回归问题，可以利用matlab
“最小二乘法”做直线拟合，用matlab求F临界值，得出结论，然
后根据matlab的输出结果进行线性回归方差分析和相关系数的显著
性检验。
x=[324,325,348.6,352.9,373.6,391.1,400.6,417.9,423.9,441.2,456.4,465 .6,473.5,488.4,498.9,508.1,521.6,541,554.2,571.8,585.4,598.8,618.9,64 0.7,656.2,683.9,710.5,764.3,788.7,815.4,838.7,872.9,886.6,921.1,941.4 ,962.7,1001.4,1038.5,1076.7,1116.5,1153.9,1195.7,1233.6,1290.4];
y=[14410.3,
14660.5 ,15021.0,15198.5,15679.1,16478.7,17185.9,17794.2,18262.3,1876 1.3,19188.0 ,19743.3,20036.4,20365.5,
20610.9,21100.6,21513.5,21885.5,22064.0,22419.3,23118.8,23854.0,24463 .5 ,24803.7,25621.5,26428.8,27011.0,27413.8,27887.0,29186.3,29885.5,3 0083.2,31605.4,32363.1,33121.0,33955.7,34878.4
,36710.4 ,38121.5,39969.6,41353.3,43564.0,44882.8,45580.2];
[p,s]=polyfit(x,y,1)
p =
1.0e+03 *
0.0313 4.5698 s =
R: [2x2 double] df: 42
normr: 4.8289e+03
1、得到回归分析结果： a=4569.8，b=31.3
x 31.34569.8+=Λ
y
（4） 实现过程 （包括Matlab 代码）
为了研究这些数据中所蕴含的规律，将生产总值看做因变量，人力资本看做自变量，画出它们的散点图，可见这些点分布在一条直线附近，所以两组分的百分比之间可能符合一元线性模型。
（5）设计总结、心得体会。