经济统计学上机1：Matlab的基本操作(共6页)

1、Matlab的基本(jbn)操作(cozu)界面(jimin)与帮助界面命令窗；工作空间和工作目录；矩阵编辑器；m文件编辑器。语句后加分号表示计算结果不显示在命令窗中。%表示注释。帮助可使用Matlab的帮助菜单，或在命令窗口中使用help命令。假如准确知道所要求助的主题词或指令名称，那么用help 是获得在线帮助的最简单有效的途径。例如在命令窗口中键入：help length数据输入、保存和导入使用命令窗或矩阵编辑器。Matlab中的数据一般以向量或矩阵的形式存储。例：矩阵的输入步骤。（1）在键盘上输入下列内容(一行中的数据用逗号或空格隔开；行与行用分号隔开)A = 1,2,3; 4,5,6

2、; 7,8,9 （2）按【Enter】键，指令被执行。（3）在指令执行后，MATLAB指令窗中将显示以下结果：A = 1 2 3 4 5 6 7 8 9 例：矩阵的分行输入A=1,2,34,5,67,8,9 （以下是显示结果）A = 1 2 3 4 5 6 7 8 9 练习（将以下语句依次拷到Maltab的命令窗口并回车运行）：x=1 2 3x=1;2;3 x=1:5 x=1:2:9x=-3:0.2:3x=5:-0.3:3A=eye(3)%生成单位阵B=zeros(3,4)C=ones(3,7)D=A B;C%将几组数据(shj)并为一组数据使用(shyng)workspace窗口中的工具按钮

3、将所有(suyu)变量数据保存为文件或从文件中打开。也可选中一个变量后点右键用Save as命令将这个变量单独保存。算术运算符和关系运算符+ - *（乘） / （幂） .* (元素乘)./ (元素除) . (元素幂)A.2A中每个元素均平方； A*B 矩阵A与B相乘； A.*B矩阵A与B的对应元素相乘；A 矩阵A的转置对数据的检索已知20012010年我国的国内生产总值数据如下表：年 份国内生产总值（亿元）2001109655.22002120332.72003135822.82004159878.32005184937.42006216314.42007265810.32008314045.

4、42009340902.82010401202.0用矩阵A表示以上数据：A=2001109655.22002120332.72003135822.82004159878.32005184937.42006216314.42007265810.32008314045.42009340902.82010401202.0用x表示GDP数据：x=109655.2 120332.7 135822.8 159878.3 184937.4 216314.4 265810.3 314045.4 340902.8 401202.0练习以下语句：x(3)%x中第3个元素； A(3,2)%A中第3行第2列的元素；x

5、(1 2 5)%x中第1、个元素； A(1 2 5,2)%A中位于第行第行第5行及第2列的元素；x(3:6)%x中第3个到第个元素；A(3:6,1:2)%A中第3行到第行、第1列到第2列的元素； x(3:end)%x中第个到最后一个元素； A(3:end,2)%A中第行到最后一行、第2列的元素；A(3:6,:)%A中第行到第行、所有列的元素； A(:,2)%A中所有行、第2列的元素；x(1:2:9)%x中第、个元素； A(3:2:9,2)%A中第、行第2列的元素；x(9:-2:1) %x中第9、个元素(yun s)； flipud(A)%将A中数据(shj)上下倒排显示或取出；fliplr(A

6、) %将A中数据左右(zuyu)倒排显示或取出；i,j = find(A150000)%找出A中大于150000的元素，其行号放入i，列号放入j;i,j = find(A=159878.3)%找出A中等于159878.3的元素的行号和列号（即元素的下标）；i= find(A=159878.3)%找出A中等于159878.3的元素的编号（Matlab中矩阵元素除了行号和列号外，还有一个编号，也称为线性索引）；find(A)%显示A中不为0的元素的编号x(find(x=150000&x200000|x150000)%显示出x中大于200000或者小于150000的元素；max(x) %返回x中最大

7、的元素值；max (A) %返回A中各列的最大的元素；min(x) %返回向量x中最小的元素；min(A) %返回A中各列的最小的元素；length(x) %返回向量x的长度，即x中元素的个数；numel(A) %返回矩阵A中元素的个数；sort(x) %将向量x的元素按递增排序；若x是矩阵，对其每列进行递增排序。sort(x, descend) %将向量x的元素按递减排序；若x是矩阵，对其每列进行递减排序。x(3) = 0%把x中的第3个元素重新赋值为0；x(1 4)=1 1%把x的第1、4个元素都赋值为1。常用基本命令和函数基本命令clc 清空命令窗口clear x 清除变量xclear

8、清除所有变量format long 计算结果显示15位数字format long g计算结果显示15位数字且不使用科学计数法format long e计算结果显示15位数字且使用科学计数法format short 计算结果显示5位数字format short g计算结果显示5位数字且不使用科学计数法format short e计算结果显示5位数字且使用科学计数法基本函数令矩阵A和向量x为：A=2001109655.22002120332.72003135822.82004159878.32005184937.42006216314.42007265810.32008314045.4200934

9、0902.82010401202.0 x=109655.2 120332.7 135822.8 159878.3 184937.4 216314.4 265810.3 314045.4 340902.8 401202.0练习(linx)以下语句：sum(x) %求x中所有(suyu)元素的和prod(x) %求x中所有(suyu)元素的积sqrt(x) %求x中各元素的平方根，亦即x.0.5diff(x) %求x中相邻元素之差(可理解为逐期增长量)，此处即各年GDP增长数量cumsum(x) %对x中所有元素累积求和cumprod(x) %对x中所有元素累积求积m,n = size(A) %返

10、回矩阵的行数和列数round(x) %对x中所有元素取整（方法为四舍五入）fix(x) %对x中所有元素取整（方法为去掉小数部分）ceil(x) %对x中所有元素取整（方法为向上取整）floor(x) %对x中所有元素取整（方法为向下取整）r=x(2:end)./x(1:end-1)-1%计算各年GDP增长率（记为r）find(r0.1)%找出增长率小于10%的增长率的序号r(find(r0.1,1)%找出增长率大于10%的年份RA=sortrows(R, 2)%按矩阵R的第2列（即增长率）递增排序RD=sortrows(R, -2)%按矩阵R的第2列（即增长率）递减排序RD(1:3,1)%找

11、出增长最快的3个年份RA(1:3,1)%找出增长最慢的3个年份命令的编写（本小节选学）打开一个空白m文件，键入若干语句；将m文件保存。其文件名即成为一条新的命令，可以在命令窗口中运行。练习：编写一个名为mycomd的命令，该命令生成100个均匀分布的随机数，并求出其中的最小值和最大值。编写后运行它。步骤：在一空白m文件中键入：x=rand(100,1); a=min(x)b=max(x)将文件保存，命名为mycomd。在命令窗口中键入mycomd，回车。随机数的生成A=rand(10,10)%生成10行10列的0,1上均匀分布的随机数，放入A中； B=randn(10,10)%生成10行10列

12、的标准正态分布的随机数，放入B中；画图(hu t)用ezplot、ezmesh、ezsurf画函数(hnsh)图的例子：ezplot(x2 - 2*x + 1,0 5) %x取值范围(fnwi)为从0到5。ezplot(normpdf(x,1,2),-6,8)%画期望值为1，标准差为2的标准正态分布概率密度曲线ezsurf(sqrt(x2 + y2),-6,6) %画二元函数图形 ezmesh(sqrt(x2 + y2),-6,6)%画二元函数图形（网格图）散点图y=109655.2 120332.7 135822.8 159878.3 184937.4 216314.4 265810.3 3

13、14045.4 340902.8 401202.0; plot(y,.)%对y中的数据画散点图。.可用*、d、o、s、v等代替。plot(2001:2010, y,*)%plot(x,y)，以x为横坐标，y为纵坐标，画散点图。xlim(2000 2011)%设定x轴显示范围为2000到2011ylim(90000 500000)%设定y轴显示范围为90000到500000grid on%给图形加网格title(2001年到2010年的中国GDP); %加上标题xlabel(年份); ylabel(GDP);简单的折线图plot(2001:2010, y)%plot(x,y)，以x为横坐标，y为

14、纵坐标，绘制曲线。x=-3:0.1:3; %表示x中的数据从-3开始到3为止，间隔为0.1。plot(x, x.2)%以x为横坐标，画平方函数图 plot(x, normpdf(x)%画标准正态分布的概率密度函数图彩色散点图（本小节选学）彩色散点图函数：scatter(x,y,s,c)，其中x、y为两个向量，x为数据点的横坐标，y为数据点的纵坐标，s为绘图点的大小，c为绘图所使用的色彩，s和c均可以以向量或表达式形式给出。在 scatter函数的前4个参数之后还可以增加第五个参数filled，表示填充绘图点。scatter与plot 的最大差别在于Scatter可以绘制变尺寸、变颜色的点图。例

15、：y=109655.2 120332.7 135822.8 159878.3 184937.4 216314.4 265810.3 314045.4 340902.8 401202.0;scatter(2001:2010,y,y)例：给定数据t=0:pi/10:2*pi, y=sin(t)，观察在不同输入参数时函数的绘图结果。t=0:pi/10:2*pi; y=sin(t)subplot(3,2,1); scatter(t,y)%subplot(3,2,1)表示准备画3行2列共6个图，这是第一个subplot(3,2,2); scatter(t,y,v)subplot(3,2,3); scat

16、ter(t,y,(abs(y)+2).4,filled)subplot(3,2,4); scatter(t,y,30,0:2: 40,v,filled)subplot(3,2,5); scatter(t,y,(t+1).3,y,filled)作业（将相应Matlab语句和要求记录的数据抄在纸上上交）1.写出完成如下操作的Matlab语句：产生(chnshng)10行20列的0,1上均匀分布的随机(su j)数矩阵A和B，将A中各元素(yun s)平方后与B中各元素相乘，得到矩阵C。找出C中大于0.3并且小于或等于0.5的元素值。2、下表反映了新中国成立以来人口的变化。新中国成立以来人口的变化

17、单位：100万人年份人口年份人口年份人口1949542197083019911158195055219718521992117219515631972872199311851952575197389219941199 1953588197490919951211 1954603197592419961224 1955615197693719971236 1956628197795019981248 1957647197896319991258 1958660197997520001267 1959672198098720011276 19606621981100120021285 196165

18、91982101720031292 19626731983103020041300 19636921984104420051308 19647051985105920061314 19657251986107520071321 19667451987109320081328 19677641988111020091335 19687851989112720101341 19698071990114320111347以上数据在Matlab中可写为：x=1949:2011y=542 552 563 575 588 603 615 628 647 660 672 662 659 673 692 705 725 745 764 785 807 830 852 872 892 909 924 937 950 963 975 987 1001 1017 1030 1044 1059 1075 1093 1110 1127 1143 1158 1172 1185 1199 1211 1224 1236 1248 1258 1267 1276 1285 1292 1300 1308 1314 1321 1328 1335 1341 1347 (1)计