经济统计学上机1：Matlab的基本操作

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

2、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%将几组数据并为一组数据使用workspace窗口中的工具按钮将所有变量数据保存为文件或从文件中打

3、开。也可选中一个变量后点右键用Save as命令将这个变量单独保存。3 算术运算符和关系运算符1) + - *（乘） / （幂） .* (元素乘)./ (元素除) . (元素幂)2) A.2A中每个元素均平方； 3) A*B 矩阵A与B相乘； 4) A.*B矩阵A与B的对应元素相乘；5) A' 矩阵A的转置4 对数据的检索已知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列

5、的元素；x(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、个元素； flipud(A)%将A中数据上下倒排显示或取出；fliplr(A) %将A中数

6、据左右倒排显示或取出；i,j = find(A>150000)%找出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&x<200000)%显示出x中大于或等于150000并且小于200000的元素；x(find(x>2000

7、00|x<150000)%显示出x中大于200000或者小于150000的元素；max(x) %返回x中最大的元素值；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(

8、1 4)=1 1%把x的第1、4个元素都赋值为1。5 常用基本命令和函数5.1 基本命令clc 清空命令窗口clear x 清除变量xclear 清除所有变量format long 计算结果显示15位数字format long g计算结果显示15位数字且不使用科学计数法format long e计算结果显示15位数字且使用科学计数法format short 计算结果显示5位数字format short g计算结果显示5位数字且不使用科学计数法format short e计算结果显示5位数字且使用科学计数法5.2 基本函数令矩阵A和向量x为：A=2001109655.22002120332.72

9、003135822.82004159878.32005184937.42006216314.42007265810.32008314045.42009340902.82010401202.0x=109655.2 120332.7 135822.8 159878.3 184937.4 216314.4 265810.3 314045.4 340902.8 401202.0'练习以下语句：sum(x) %求x中所有元素的和prod(x) %求x中所有元素的积sqrt(x) %求x中各元素的平方根，亦即x.0.5diff(x) %求x中相邻元素之差(可理解为逐期增长量)，此处即各年GDP增长

10、数量cumsum(x) %对x中所有元素累积求和cumprod(x) %对x中所有元素累积求积m,n = size(A) %返回矩阵的行数和列数round(x) %对x中所有元素取整（方法为四舍五入）fix(x) %对x中所有元素取整（方法为去掉小数部分）ceil(x) %对x中所有元素取整（方法为向上取整）floor(x) %对x中所有元素取整（方法为向下取整）r=x(2:end)./x(1:end-1)-1%计算各年GDP增长率（记为r）find(r<0.1)%找出增长率小于10%的增长率的序号r(find(r<0.1)%找出增长率小于10%的值sort(r) %将增长率按递增

11、排序sort(r, 'descend') %将增长率按递减排序R=2002:2010' r%将增长率和相应年份放在一起R(R(:,2)>0.1,1)%找出增长率大于10%的年份RA=sortrows(R, 2)%按矩阵R的第2列（即增长率）递增排序RD=sortrows(R, -2)%按矩阵R的第2列（即增长率）递减排序RD(1:3,1)%找出增长最快的3个年份RA(1:3,1)%找出增长最慢的3个年份6 命令的编写（本小节选学）打开一个空白m文件，键入若干语句；将m文件保存。其文件名即成为一条新的命令，可以在命令窗口中运行。练习：编写一个名为mycomd的命令，

12、该命令生成100个均匀分布的随机数，并求出其中的最小值和最大值。编写后运行它。步骤：在一空白m文件中键入：x=rand(100,1); a=min(x)b=max(x)将文件保存，命名为mycomd。在命令窗口中键入mycomd，回车。7 随机数的生成A=rand(10,10)%生成10行10列的0,1上均匀分布的随机数，放入A中； B=randn(10,10)%生成10行10列的标准正态分布的随机数，放入B中；8 画图8.1 用ezplot、ezmesh、ezsurf画函数图的例子：ezplot('x2 - 2*x + 1',0 5) %x取值范围为从0到5。ezplot(&

13、#39;normpdf(x,1,2)',-6,8)%画期望值为1，标准差为2的标准正态分布概率密度曲线ezsurf('sqrt(x2 + y2)',-6,6) %画二元函数图形 ezmesh('sqrt(x2 + y2)',-6,6)%画二元函数图形（网格图）8.2 散点图y=109655.2 120332.7 135822.8 159878.3 184937.4 216314.4 265810.3 314045.4 340902.8 401202.0; plot(y,'.')%对y中的数据画散点图。'.'可用'*

14、'、'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');8.3 简单的折线图p

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

16、表示填充绘图点。scatter与plot 的最大差别在于Scatter可以绘制变尺寸、变颜色的点图。例：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个图，这是第一

17、个subplot(3,2,2); scatter(t,y,'v')subplot(3,2,3); scatter(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')9 作业（将相应Matlab语句和要求记录的数据抄在纸上上交）1.写出完成如下操作的Matlab语句：产生10行20列的0,1上均匀分布的随机数矩阵A和B，

18、将A中各元素平方后与B中各元素相乘，得到矩阵C。找出C中大于0.3并且小于或等于0.5的元素值。2、下表反映了新中国成立以来人口的变化。新中国成立以来人口的变化 单位：100万人年份人口年份人口年份人口1949542197083019911158195055219718521992117219515631972872199311851952575197389219941199 1953588197490919951211 1954603197592419961224 1955615197693719971236 1956628197795019981248 1957647197896319991258 1958660197997520001267 1959672198098720011276 19606621981100120021285 19616591982101720031292 19626731983103020041300 19636921984104420051308 196470519851059200613