MATLAB编程与Geoda应用(1)

研究生学位课程试卷

院（系、所）商学院专业数量经济学

考试科目Matlab编程与Geoda应用第二学期研究生姓名余帮正学号 162200627

说明

一、凡学位课程考试试题、试卷必须与本封面一起装订。阅卷导师务必用红笔批卷，并在本封面规定位置打分、写完评语后在二周（论文考试一个月）内交院（系、所）办公室教务员，教务员及时做好成绩登记，在学期结束前或第二学期初将成绩单交研究生处统一整理归档。试题、试卷由院（系、所）办公室保管。

二、学位课程考试用纸除计算机专用打字纸、16开小方格稿子纸外，一律使用研究生处统一印制的“学位课程考试纸”。

三、该封面请用A4纸双面打印，将此说明打印于封面背面。

硕士研究生课程作业

作业题目: 房价对我国居民储蓄率的分析——基于VAR模型实证研究

作业类型：学术论文

课程名称: MATLAB编程与Geada应用

授课老师: 崔百胜

专业班级: 数量经济学

研究生姓名: 余帮正

研究生学号: 162200627

完成时间2017 年8 月

2017年8月

摘要

购房已作为我国的刚性需求，造成了我国房价不断上升，居民大量储蓄买房而减少消费。在此社会背景下，高房价是否造就了居民高储蓄率。本文针对这一问题，从理论和实证出发来探讨我国房价如何影响我国居民储蓄率，并试图通过研究能够合理解释我国高房价、高储蓄、低消费的现象。通过研究发现，我国房价对居民储蓄率的影响显著，但是房价与我国居民储蓄率之间还存在着其他内生性以及非对称的互动影响，比如居民收入水平，收入水平同样对于居民储蓄率有着很大影响。产业结构的升级会改变社会经济中的收入分布，那么对于居民的储蓄也同样产生影响。对于居民储蓄的影响因素较多，我们不能单独来分析房价对于居民储蓄率的影响。因此，有必要在考虑到相关影响因素下再进行研究。

本文在介绍VAR 模型的基础上，使用公共因子的提取以及脉冲响应函数来选取相关特定几个变量来进行分析，很好的解决了VAR模型容纳变量过少的问题又能够保证研究影响因素的相关性。使用我国2000 年至2015年的宏观经济数据，运用 VAR 模型对我国房地产市场价格影响居民储蓄率的有效性进行了分析，并得到了一些结论。

关键词：房价消费居民储蓄率互动影响 VAR模型

Abstract

Purchase has been as China's rigid demand, resulting in rising house prices in China, a large number of residents to buy a house and reduce consumption. In this social context, whether high housing prices have created a high savings rate of residents. In this paper, from the theoretical and empirical study to explore how China's housing prices affect the savings rate of our residents, and trying to study through a reasonable explanation of China's high prices, high savings, low consumption phenomenon. It is found that the impact of house price on household savings rate is significant, but there are other endogenous and asymmetric interaction effects between house price and household savings rate, such as the income level of residents and the income level. Big influence. The upgrading of the industrial structure will change the income distribution in the socio-economy, so the savings for the residents also have an impact. There are many factors that affect the savings of residents, and we can not analyze the impact of house prices on household savings rates alone. Therefore, it is necessary to study the relevant factors in consideration of the relevant factors.

Based on the introduction of the V AR model, this paper uses the extraction of the common factor and the impulse response function to select the relevant variables to analyze the variables, which solves the problem that the V AR model contains too few variables and ensures the correlation of the influencing factors Sex. Using the macroeconomic data of China from 2000 to 2015, the VAR model is used to analyze the effectiveness of the real estate market price in China, and some conclusions are obtained.

Key words:price consumption, household savings rate, interaction, V AR model

目录

摘要 ................................................................................................................................................... I Abstract ............................................................................................................................................ II 目录 ................................................................................................................................................ III 第一章导论. (1)

1.1研究背景与意义 (1)

1.2 研究方法和思路 (1)

1.3 研究技术路线图 (2)

第二章文献综述与理论基础 (3)

2.1 文献综述 (3)

2.2 理论基础 (4)

2.2.1 居民储蓄行为理论回顾 (4)

2.3 相关概念界定 (6)

2.3.1 居民储蓄 (6)

2.3.2 商品房价格 (6)

第三章我国房价与居民储蓄率的分析 (7)

3.1 我国高房价情况分析 (7)

3.2 我国居民高储蓄分析 (8)

第四章房价波动影响居民储蓄的实证分析 (8)

4.1数据说明及分析方法 (8)

4.2 房价影响居民储蓄的实证研究 (10)

第五章结论与建议 (18)

5.1 研究结论 (18)

5.2 研究建议 (19)

参考文献 (20)

第一章导论

1.1研究背景与意义

投资、消费和进出口是我国经济增长的三驾马车。而我国房地产是投资、消费和进出口三大增长发动机的主要推动力。房地产业是经济发展的重要组成部分，是各国居民消费和投资的重要组成部分，也是个人与机构持有的主要实物资产。我国城市家庭财产调查显示房产在家庭总资产中所占比重已接近一半。这表明，房产己逐渐成为我国普通居民家庭价值量最大的财产。房产具有消费和投资的双重属性，一方面作为耐用消费来说，其价格的上涨通过影响居民的消费意愿和消费数量来影响居民的储蓄行为，这主要取决于住房的需求弹性。另一方面作为投资品来说，其价格的上涨，将引起居民预期收益和总财富的变化，从而影响居民的消费一储蓄行为。而房价变动又可以通过储蓄一投资渠道对实体经济产生影响。因此，研究房价波动对储蓄的影响及其影响的渠道，对我们清楚的认识房地产价格波动对国民经济发展的影响具有重要意义。

我国房地产市场持续繁荣多年，房价持续上涨，政府刺激住宅消费，扩大内需的政策不能达到预期效果，与此同时居民消费率不断下降，储蓄意愿强烈，储蓄率长期居高不下，高房价与高储蓄并存的现状成为政府和学者们关注的焦点。深入研究房地产市场价格波动影响储蓄的机理，有助于全面认识房地产市场与居民储蓄的内在联系，认识房地产价格过度波动。通过储蓄一投资渠道对实体经济所产生的影响对居民购房行为的动机分析以及房价上涨，对不同微观主体储蓄行为的影响分析可为政府管理部门监控房地产市场、利用房市调控居民消费一储蓄行为，促进社会和谐发展提供理论指导与实证支持。对它国经验的研究还为我国解决居民住房问题，稳定房价提供经验借鉴。在当前房地产市场价格持续攀升和居民储蓄总量连年上涨并存的背景下，本文的研究对我国房地产市场管理和宏观经济稳定具有现实意义。

1.2 研究方法和思路

本文拟采用的研究方法主要有三个：一是理论分析。本论文从生命周期假说

和预防性储蓄理论入手，考虑我国经济政策，制度变迁，消费习惯等特点来深入分析住房价格对储蓄的影响，建立房价与居民储蓄之间的理论联系。二是模型分析。本文将运用住房消费一投资需求模型和使用成本模型，结合我国现阶段的实际情况分析我国居民购房动机，并构建了作为耐用消费品，居民购房一储蓄过程模型,来分析消费型购房者的储蓄行为，另外运用居民储蓄一购房选择模型分析房价对不同微观主体储蓄行为的影响。三是实证分析。运用图表分析和基于VAR 模型的计量方法，利用中国的时间序列数据对模型所反映的统计规律进行计量检验，探寻通过房地产渠道对居民储蓄的影响，为我国通过调控房地产市场引导居民储蓄行为，促进消费、扩大内需、提供实证支持。

1.3 研究技术路线图

第二章文献综述与理论基础

2.1 文献综述

本文文献综述主要从两个方面来分析对于储蓄率的影响，一是房价对于储蓄率的影响；二是人口结构对于储蓄率的影响。

通过文献研究发现，房价对于储蓄率的影响研究主要在于三个方面：一是房价与储蓄率呈现正相关关系。相关研究有陈崇和葛杨基（2011）通过我国 1997～2008年省际面板数据模型，分析表明从我国整体和中西部地区的角度上来看，房价与储蓄的关系显著为正。陈斌开和杨汝岱（2013）利用 2002～2007年国家统计局发布的中国城镇住户调查（UHS）的微观数据分析发现，住房价格上升的过高会抑制年轻人和老年人的消费需求，促进储蓄率的升高。陈彦斌和邱哲圣（2011）通过构建含有房价因素的Bewley模型，通过内生化住房需求的手段以及根据相关生命周期的特征，结果表明，年轻家庭会因为房价的上涨而储蓄。

二是房价与储蓄负相关。骆祖春、李勇刚和李祥（2014）先把住房消费加入到世代交叠模型，然后通过构建面板联立方程，结果理论和实证均表明房价与储蓄的关系显著负相关。三是房价与储蓄率关系不显著。田志华和田艳芳（2014）使用2006～2011年省际面板数据，运用系统GMM方法检验房价与储蓄率的关系，结果房价增长对社会居民总体储蓄率影响微弱且不显著。

有关人口结构对储蓄率影响的文献国内外都有不同程度的研究，其中最典型的莫过于生命周期假说（life cycle hypothesis，简称 LCH）。该理论认为，为了使一生效用最大化，理性人会根据预期收入平滑人生三个阶段（少年，成年和老年）的消费，少年和老年是低储蓄，成年是高储蓄阶段。Leff（1969）使用74个国家（包括27个发达和47个欠发达国家）的跨国数据，发现抚养比（少儿抚养比和老年抚养比）与国民储蓄率呈反方向变化。汪伟（2009）根据生命周

期理论，并利用中国面板数据进行实证，研究表明抚养比与储蓄率呈反方向变化。董丽霞和赵文哲（2011）利用中国面板数据，通过面板自回归（pavar）模型把人口结构内生化，实证表明抚养比与储蓄率呈反方向变化。Higgins和Williamson（1997）以及袁志刚和宋铮（2000）的研究均支持LCH的假说。当然也有部分研究不支持生命周期理论假设，王领和张余（2015）通过构造全国时间序列和省际面板数据，结果时间序列支持生命周期理论，而省际面板数据实证表明少儿抚养比与储蓄率负相关，而老年抚养比与储蓄率成正相关。刘铠豪和刘渝琳（2015）利用系统广义矩估计和门槛回归证明少儿抚养比与储蓄率负相关，而老年抚养比与储蓄率正相关。杨继军和张二震（2013）在面板模型中考虑养老保险制度，结果表明由于家庭支出的“远期化”导致少儿抚养比与储蓄正相关，老年抚养比与储蓄率负相关。李祥和李勇刚（2013）通过面板数据联立方程研究人口抚养比和居民消费的关系，实证表明只有少儿抚养比与居民消费显著负相关，而老年抚养比却与居民消费关系不显著，从反面证实了不支持生命周期理论。钟水映和李魁（2009）、刘世杰和张士斌（2009）以及 Goldberger（1973）的研究也均不支持生命周期假说。

2.2 理论基础

2.2.1 居民储蓄行为理论回顾

储蓄是社会进步与经济增长的促进力量，国外学者们一直致力于储蓄理论的研究。本文主要介绍以下四种储蓄理论：

（1）生命周期一持久收入模型（LC/PIH）

生命周期一持久收入模型（LC/PIH）一直是研究储蓄行为（消费行为）的主要理论框架。它是美国经济学家莫迪里安尼在费雪的效用最大化基础上提出的。该理论认为消费者的消费行为受一生中收入的影响，消费者为了实现整个生命周期内的效用最大化，就要通过在年轻时进行储蓄而在退休以后进行消费的方法来平滑一生的消费。

（2）预防性储蓄理论（Precautionary Saving Hypothesis）

预防性储蓄理论对生命周期一持久收入储蓄理论进行了修正和补充。预防性储蓄假说吸收了生命周期理论中的理性消费者和效用最大化的假说，引入不确定性问题及消费者跨时选择分析，认为消费者储蓄不仅仅只是将财富平均分配于整

个生命周期，还在于防范不确定事件的发生，如未来收人的波动等。首次对预防动机的储蓄模型进行分析的是利兰德，他定义预防性储蓄为由未来不确定性收入而引起的额外的储蓄。

（3）流动性约束理论（Liquidity Constrain Theory）

流动性约束理论（Liquidity Constrain Theory）认为，生命周期一永久收入储蓄理论假设个人能够在相同的利率水平上借贷和储蓄与现实不符，实际中存在着流动性约束。流动性约束从两个途径增加储蓄，首先，不论流动性约束何时发生，它都会使一个人的消费比他想要的消费少。其次，流动性约束不在现期发生，但它将在未来发生的事实同样会降低现期消费。例如，假设在下期预期收入降低，若没有任何流动性约束，并且下期实际收入降低了，那么个人可以通过借款来避免消费锐减，但如果存在流动性约束在个人没有储蓄情况下，那么收入下降就会引起消费下降。因此，流动性约束的存在会导致个人将增加储蓄作为保险以减少收入在未来下降产生的影响。

（4）缓冲存货储蓄理论（Butter-Stock Saving Theory）

缓冲存货储蓄模型是由盯等提出的。如果消费者存在预防性储蓄动机，并且如果预期未来有更高收入，同时又是不耐心的，他们将倾向于选择大于当一前收入的消费。在这种情况下消费者就会表现出到所称的缓冲存货储蓄行为。这种理论认为，储蓄相当于一种缓冲存货，以便在境况艰难时维持消费而在境况如意时增加消费。缓冲存货储蓄者一般有一个财富对持久收入的目标比率，如果财富低于目标，预防性储蓄动机将战胜不耐心从而加大储蓄。在相反的情况下，不耐心就会占上风从而使消费者选择负储蓄。而认为，每个人都有一系列拇指法则用来指导他们的消费行为。这些拇指法则的许多例子表明花费一个人的现期收入是合理的，但资产急剧下降只能发生在意外的环境中。这样的拇指法则可以导致消费者通过储蓄和借贷来拉平短期收入波动，从而使消费在短期里较合理地遵循永久收入假说的预言，但它们还可以使消费在长期中与收入轨迹相当紧密。我们从以上储蓄理论可以看出，虽然储蓄理论得到了不断的发展对居民储蓄行为的解释力不断加强，但是各种储蓄理沦更多的是重视收入主要是指工资性收入的不确定性对储蓄行为的影响，而没有涉及资产价格变动对居民储蓄行为的影响。不过，这些理论为我们深入研究房地产价格变动对居民储蓄行为的影响提供了许多有用

的观察视角和思路。

2.3 相关概念界定

为了更明确地开展我们后面的研究，这里首先对有关的概念进行界定，主要包括以下几个概念。

2.3.1 居民储蓄

居民储蓄可分为两个层次：广义的和狭义的居民储蓄。广义居民储蓄是指居民可支配收入减去消费与转移支付后的差额，包括实物性储蓄和金融性储蓄，其实物性储蓄包括城乡个人固定资产投资、购置的耐用消费品和居民购买的商品房等金融性储蓄包括银行存款、手持现金、有价证券、保险费和外币储蓄等。而狭义的居民储蓄是指居民银行储蓄存款。由于银行存款有较好的流动性是居民持有的主要金融资产，孙凤（2014）研究认为我国的居民银行存款占总储蓄的比值。狭义的居民储蓄按照储蓄对象的不同又可以分为城镇居民储蓄存款和农村居民储蓄存款。本文研究对象为我国的住房价格和居民储蓄，由于我国城镇和农村实行不同的住房制度，房价波动影响的主要是城镇居民储蓄，因此本论文中居民储蓄指的是狭义的储蓄中的城镇居民储蓄存款，包括城镇居民活期存款和定期存款。

2.3.2 商品房价格

商品房在中国兴起于80年代，它是指在市场经济条件下，具有经营资格的房地产开发公司（包括外商投资企业）通过出让方式取得土地使用权后经营的住宅，均按市场价出售。商品房价格由成本、税金、利润、代收费用以及地段、层次、朝向、质量、材料差价等组成。另外，从法律角度来分析，商品房是指按法律、法规及有关规定可在市场上自由交易，不受政府政策限制的各类商品房屋，包括新建商品房、二手房（存量房）等。商品房根据其销售对象的不同，可以分为外销商品房和内销商品房两种。本文的研究对象是指狭义的房地产市场，即商品房市场，而较少涉及商铺和写字楼等房地产市场。如无特殊说明，下文中所涉及的房地产均指商品房住宅资产。

第三章我国房价与居民储蓄率的分析

3.1 我国高房价情况分析

高房价带来的经济影响是中国内需缺乏，其主要原因是民众的购买力低下。而导致购买力不够的一个重要原因就是民众在住房上投入了太多的资金，这一点从房价收入比就可以看出。2000 年全国住宅均价为2112 元，2015 年全国住宅均价为6793元，上涨了221.6%。我国目前的房价，对于多数人来说都是难以坦然和承受的。在许多城市，购买一套房相当于吞噬掉一个中等收入家庭数十年的全部收入。这对于一个家庭其他消费的影响是显而易见的。住房是和民生息息相关的行业，这一定会给住房民生带来危害，造成社会的不稳定，给中国经济的发展带来了巨大的伤害。高房价、高收入比极大地降低了大多数人的真实收入，可以说是房地产商对普通居民和终端消费工业的剥夺。

有调查显示，2005年底北京、上海、天津等10个城市平均月供收入比为35%。有31.75%的人房租、月供占到了其收入的50%以上，而这一比例到2009年又进一步提高，严重超过了国际警戒线1/3的危险水平。高房价造成家庭住房消费过度透支，成为房奴的人不得不挤压其它消费，缩减各项开支，对其它产品的需求被迫降低，这自然就影响了更多工业的发展。高房价对产业结构的影响高房价对产业结构的扭曲破坏是深远和严重的。一方面，高房价刺激的庞大房地产业带动了高能耗、高污染的粗钢生产的一路发展，延缓慢了部分基础工业的结构调整。有数据显示，2008年在GDP构成中，钢筋水泥占了57%，到了2009年占了67%，房地产业的拉动，让我国成了全球最大的钢铁消费市场，并且给部分产能落后、规模小、效率低、污染重、对资源和生态环境构成严重威胁的中小企业创造了很大的生存空间，使我国付出了巨大的环境代价。另一方面，高房价、高收入比严重挤压了其他工业的生存空间，使得原本就过剩的消费工业生存状况更加严峻；再加上超高利润的诱惑，工业资金大量进入房地产业。从2007年底开始，格力、海尔、长虹等多家老牌家电企业，高调介入地产领域，其它企业更是纷纷效仿。大量工业资金的流失，不但抑制了工业的发展和技术改造，延缓了产业结构的调整速度，破坏了工业长期的发展潜力和国际竞争力，而且进一步推高了房价和投机的欲望，导致我国整体经济结构严重扭曲。

3.2 我国居民高储蓄分析

高储蓄对我国经济会带来很多的负面影响，负面影响主要有以下几个方面：一是储蓄的高增长加大了银行的压力，增加了银行的经营成本。而居民的高债权极易形成和加剧银行的系统性风险。由于目前我国金融改革还没有到位，金融市场规模偏小，投资渠道狭窄，银行除了贷款给企业，很少有其他渠道来消化存款。从而加大了企业高负债对银行经营的威胁，一旦出现储蓄挤兑，银行会出现支付风险。二是储蓄转化投资的效率不高，就有可能导致银行业的风险累积。储蓄的高增长使我国形成了以银行贷款为主的直接融资格局，而债券和股票市场等间接融资却发展总体滞后，导致直接融资比例过大。结果企业融资高度依赖于银行体系，银行承担了一些本应由金融市场承担的风险，金融风险向银行业集中。三是我国经济中存在投资率较高、出口依赖性较强、消费率偏低等问题，归根结底是因为我国的储蓄率过高。有关数据显示，2012年前11个月，我国投资增加27.8%，依然保持高位。而过去几年，我国的储蓄率一直占到GDP的45%左右，这就意味着有很大的资金闲置,造成资金的使用率低下和一定程度上的浪费。四是储蓄与消费反差过大，会严重影响生产和消费。居民储蓄的高增长，意味着消费必将相应的减少，消费品在市场上出现供过于求，生产将出现相应的缩减，如果持续如此将出现市场疲软，经济萎靡。因此，根据中国的经济发展水平和国情，我们应该减小储蓄比例、增加内需，扩大投资，保持我国经济在未来长期稳定的增长。

第四章房价波动影响居民储蓄的实证分析在本文中，前部分主要以理论分析了房价与居民储蓄的问题，房价上涨与居民储蓄之间的关系存在不确定性。本章中，主要以实证分析我国房价对居民储蓄的影响，并建立各变量之间的模型，基于VAR模型进行深入分析。首先说明要使用的方法以及变量数据选择及来源，然后是模型的构建与设定，最后对实证结构进行解释说明。考虑到数据的可得性，选取比较有代表性的解释变量进行分析。最后利用操作软件R、Eviews9.0 和MATLAB等来建模进行说明。

4.1数据说明及分析方法

（1）指标选取与数据来源

本文研究的是我国房价对于居民储蓄的影响，考虑到数据的可得性和连续

性，本文中以我国2000-2015年的相关数据进行实证分析。数据来源于《中国统计年鉴》和《地区统计年鉴》。

鉴于对于居民储蓄率的影响因素较多，通过相关参考文献研究和实际理论分析，本文先通过选取一些有代表性的变量，再进一步分析选取有针对性变量来实证分析。本文中被解释变量为居民储蓄率，解释变量有：居民收入房价、城镇化率、教育经费、产业结构和人口年龄结构抚养比。相关变量解释说明和衡量方法如下表4-1所示：

表4-1 变量说明与描述统计

数据来源：中国统计年鉴计算整理得到

（2）分析方法

选取的数据为时间序列数据，由于数据的非平稳性，所以在进行构造向量自回归模型（VAR）和普通最小二乘法（OLS）之前有必要先对数据变量的平稳性进行检验。由于选取的数据解释变量较多，先基于主成分分析，得出特定几个解释变量来考虑，其余为控制变量来考虑。实证部分先是观察时序图变化，再进行单位根检验，之后做脉冲响应分析和建立相关模型来分析。

4.2 房价影响居民储蓄的实证研究

本文选取的数据为我国2000-2015年的16年的相关时序数据，在进行建模之前，首先利用操作R软件先进行主成分分析，根据分析结果选取变量来考虑。之后利用Eviews软件进行单位根检验，协整检验，脉冲响应分析等，最后利用MATLAB建立VAR模型。

（1）主成分分析

通过软件操作可得主成分分析KMO值0.87，适合做主成分分析。通过图4-1 的碎石土可以得出一个因子能够解释大部分信息。由于选取一个因子分析不适合建立模型分析，为了下一步建模分析选取三个因子来说明。通过表2的旋转前因子载荷矩阵和旋转之后的载荷矩阵分析，可以得到三个因子可以有效解释说明。在表4-2 中，旋转前的因子对于房价的解释力度,是0.99，对于城镇化率的解释力度也是0.99，对教育经费支出解释力度是0.97，对产业结构的解释是0.98，对人口年龄结构总抚养比-0.92，对人均国内生产总值解释力度是0.98。说明因子能够解释大部分解释变量。这些因子存在很大的内生相关性，下一步的建模分析中，容易产生共线性问题。

图4-1 碎石图分析

表 2 因子载荷矩阵

通过软件R操作实现得到。

（2）单位根检验

对于序列数据的平稳性检验，利用 ADF 检验法。其中对变量房价、教育经费支出和居民收入水平取对数之后，再对序列数据sr、lnhp、ub、lnes、lnes、ps、ns和lnpgdp进行单位根检验，可根据赤池（AIC ）信息准则和施瓦茨（ SC ）信息准则来确定滞后阶数，利用 Eviews9. 0 得到检验结果如下：得到变量数据几乎是平稳的，通过一阶差分或者两阶差分达到平稳。

通过表4-3单位根检验，居民储蓄率原序列是平稳的，其他解释变量原序列数据都是不平稳的。变量房价在取对数之后，一阶差分平稳；城镇化率原序列数据不平稳，一阶差分平稳；教育经费支出取对数之后在二阶差分达到平稳；产业结构在一阶差分滞后期为1的情况

下达到平稳；人口年龄结构总抚养比在二阶差分滞后期为2的情况下达到平稳；居民收入水平取对数在二阶差分滞后期为1的情况下达到平稳。

表4-3 各变量的单位根检验

（3）脉冲响应分析

根据VAR模型单位根全部落在单位圆内，说明VAR模型稳定，可以对VAR 模型进行一个标准差的脉冲响应函数分析。下图4-2、图4-3和图4-4分布显示VAR模型的单位根全部落在单位圆内，说明选择的变量建立的模型是稳定的。

图4-2是关于房价，教育经费和产业结构对于居民储蓄的影响。在除考虑房价对于居民储蓄的影响外，教育经费和产业结构对于居民储蓄率的影响重大。政府财政支出中的教育支出规模越大，则居民的教育投资的储蓄比值下降。产业结构的升级那个改变社会经济中的收入分布，进而对居民储蓄产生较大影响。针对这几个变量分析，可以反映在非城镇地区。

图4-3是关于房价，人口年龄抚养比和居民收入水平对于居民储蓄率的影响研究。在出考虑房价对于居民储蓄率的影响外，人口结构和居民收入水平对于储蓄率的影响重大。本文中人口结构是以少儿抚养比和老年抚养比来衡量的，非劳动年龄人口比重越高，居民储蓄率越低，家庭支出会更加大。居民收入水平是影响房价和居民储蓄的重要因素，显而易见的是城镇人口人均收入水平较高而非城镇人口收入水平低。选取的这几个变量能够反映城镇地区情况。

图4-4是关于房价、城镇化率和居民收入水平对于储蓄率的影响。本文中的城镇是以城市常住人口占比总人口来衡量的。我国不同城镇化率地区房价也明显

存在差异，一二线城市房价越来越高，进几年表现出持续上升且无下降的趋势。而四五线城市，房价虽然不高，但是买房的欲望仍然很低，而且有着下降的趋势。造成这种房价极端的情况，还是与所在城市的不同。所以，针对这一情况，为了说明房价对于储蓄率的影响，有必要对我国不同地区加以考虑。

-1.5

-1.0

-0.5

0.0

0.5

1.0

1.5

Inverse Roots of AR Characteristic Polynomial

图4-2关于sr 、lnhp 、lnes 和ps 的脉冲响应分析。

-1.5

-1.0

-0.5

0.0

0.5

1.0

1.5

Inverse Roots of AR Characteristic Polynomial

图4-3 关于sr 、ns 、lnhp 和pgdp 的脉冲响应分析

《MATLAB程序设计与应用(刘卫国)》(第二版) 答案

《MATLAB程序设计与应用（刘卫国）》（第二版）实验一MATLAB运算基础 1.（1） z1=2*sin(pi*85/180)/(1+exp(2)) （2）x=[2,1+2i;-0.45,5]; z2=log(x+sqrt(1+x.^2))/2 （3）a=-3.0:0.1:3.0; z3=0.5*(exp(0.3*a)-exp(-0.3*a)).*sin(a+0.3)+log((0.3+a)/2) （4）t=0:0.5:2.5; z4=(t>=0&t<1).*t.^2+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3).*(t.^2-2*t+1) 2. A=[12,34,-4;34,7,87;3,65,7]; B=[1,3,-1;2,0,3;3,-2,7]; （1）a=A+6*B b=A-B+eye(size(A)) %I=eye(size(A)) （2）c=A*B d=A.*B （3）e=A^3 f=A.^3 （4）g=A/B h=B\A （5）m=[A,B] n=[A([1,3],:);B^2] 3. 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,11]; （1）C=A*B （2）D=C(3:end,2:end) 4.（1） a=100:999; b=rem(a,21)==0; c=find(b); d=length(c) （2）ch='Just as Bianhaiman said,Xiehong is ...'; e=find(ch>='A'&ch<='Z'); ch(e)=[] 实验二 MATLAB矩阵分析与处理 1.E=eye(3);

Matlab编程与应用习题和一些参考答案

Matlab 上机实验一、二 3.求下列联立方程的解???????=+-+-=-+=++-=--+4 1025695842475412743w z y x w z x w z y x w z y x >> a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10]; >> b=[4;4;9;4]; >> c=a\b 4．设???? ??????------=81272956313841A ，??????????-----=793183262345B ，求C1=A*B’;C2=A’*B;C3=A.*B,并求上述所有方阵的逆阵。 >> A=[1 4 8 13;-3 6 -5 -9;2 -7 -12 -8]; >> B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7]; >> C1=A*B' >> C2=A'*B >> C3=A.*B >> inv(C1) >> inv(C2) >> inv(C3) 5．设 ?? ????++=)1(sin 35.0cos 2x x x y ，把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线。 >> x=linspace(0,2*pi,101); >> y=cos(x)*(0.5+(1+x.^2)\3*sin(x)); >> plot(x,y,'r') 6．产生8×6阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。并求该矩阵全体数的平均值和均方差。 （mean var ） a=randn(8,6) mean(a) var(a) k=mean(a) k1=mean(k) i=ones(8,6) i1=i*k1 i2=a-i1 i3=i2.*i2 g=mean(i3) g2=mean(g)

刘卫国版MATLAB程序设计与应用课后实验六八九

实验六 高层绘图操作 %第一题： 程序代码如下： x=linspace(0,2*pi,101); y=(0.5+3*sin(x)./(1+x.^2)).*cos(x); plot(x,y) 01234567 -1 -0.5 0.5 1 1.5 %第二题： %（1） 程序代码如下： x=linspace(-2*pi,2*pi,100); y1=x.^2; y2=cos(2*x); y3=y1.*y2; plot(x,y1,'b-',x,y2,'r:',x,y3,'y--'); text(4,16,'\leftarrow y1=x^2'); text(6*pi/4,-1,'\downarrow y2=cos(2*x)'); text(-1.5*pi,-2.25*pi*pi,'\uparrow y3=y1*y2');

-8 -6 -4 -2 2 4 6 8 -30-20 -10 10 20 30 40 %（2） 程序代码如下： x=linspace(-2*pi,2*pi,100); y1=x.^2; y2=cos(2*x); y3=y1.*y2; subplot(1,3,1);%分区 plot(x,y1); title('y1=x^2');%设置标题 subplot(1,3,2); plot(x,y2); title('y2=cos(2*x)'); subplot(1,3,3); plot(x,y3); title('y3=x^2*cos(2*x)');

-10 10 0510 15202530 35 40y1=x 2 -10 10 -1-0.8 -0.6 -0.4-0.200.20.4 0.6 0.8 1y2=cos(2*x) -10 10 -30-20 -10 10 20 30 40 y3=x 2*cos(2*x) %（3） 程序代码如下： x=linspace(-2*pi,2*pi,20); y1=x.^2; subplot(2,2,1);%分区 bar(x,y1); title('y1=x^2的条形图');%设置标题 subplot(2,2,2); stairs(x,y1); title('y1=x^2的阶梯图'); subplot(2,2,3); stem(x,y1); title('y1=x^2的杆图'); subplot(2,2,4); fill(x,y1,'r');%如果少了'r'则会出错 title('y1=x^2的填充图'); %其他的函数照样做。

matlab程序设计与应用第二版习题答案

matlab程序设计与应用第二版习题答案【篇一：matlab程序设计与应用(第二版)实验答案】 %实验一 matlab运算基础 %第1题 %（1） z1=2*sin(85*pi/180)/(1+exp(2)) %（2） x=[2,1+2i;-0.45,5]; z2=0.5*log(x+sqrt(1+x.^2)) %(3) a=-3.0:0.1:3.0; z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2) %(4) t=0:0.5:2.5; z4=t.^2.*(t=0t1)+(t.^2-1).*(t=1t2)+(t.^2-2*t+1).*(t=2t3) %第2题 a=[12 34 -4;34 7 87;3 65 7]; b=[1 3 -1;2 0 3;3 -2 7]; a+6*b a-b+eye(size(a)) a*b a.*b a^3 a.^3 a/b b\a [a,b] [a([1,3],:);b^2] %第3题 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 11] c=a*b f=size(c) d=c(f(1)-2:f(1),f(2)-1:f(2)) whos %第4题 %(1):

a=100:999; b=rem(a,21); c=length(find(b==0)) %(2): a=lsdhksdlkklsdkl; k=find(a=aa=z); a(k)=[] %实验二 matlab矩阵分析与处理 %第1题 e=eye(3); r=rand(3,2); o=zeros(2,3); s=diag([2,3]); a=[e,r;o,s]; a^2 b=[e,(r+r*s);o,s^2] %第2题 h=hilb(5) p=pascal(5) hh=det(h) hp=det(p) th=cond(h) tp=cond(p) %第3题 a=fix(10*rand(5)) h=det(a) trace=trace(a) rank=rank(a) norm=norm(a) %第4题 a=[-29,6,18;20,5,12;-8,8,5] [v,d]=eig(a) %数学意义略 %第5题方法一 %(1): a=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; b=[0.95,0.67,0.52]; x=inv(a)*b %(2):

MATLAB程序设计与应用(第二版)

阅阅读读时时：：请请选选择择““视视图图||文文档档结结构构图图””，，弹弹出出文文档档中中的的标标题题链链接接。。数学软件 MATLAB 程序设计与应用

第1章MATLAB系统环境 1.1 MATLAB概貌 1.2 MATLAB环境的准备 1.3 MATLAB操作界面 1.4 MATLAB帮助系统 自上世纪80年代以来，出现了科学计算语言，亦称数学软件。 MATLAB Mathematica Mathcad Maple LINDO LINGO 1.1 MATLAB概貌p3 MATLAB 是MATrix LABoratory（矩阵实验室）的缩写。 1984年由Math Works 公司推出，现已成为国际公认的优秀的工程应用开发环境，是影响最大，流行最广的科学计算语言。 1.1.1 MATLAB的发展

1. 从MATLAB 4.2c开始，每个版本增加了一个建造编号； 2. 例如MATLAB7.6的建造编号是R2008a。说明MATLAB7.6与MATLAB2008a是等同的； 3. 对于建造编号，正规化以后，每年出两个版本。一般来说。a是测试版，b是正式版。a是前半 年出，b是后半年出。 教材采用MATLAB7.0（R14，2004） 实验室采用MATLAB 7.8（R2009a，2009.3，汉化） 1.1.2 MATLAB的主要功能p4 ◆数值计算和符号计算功能 ◆绘图功能 ◆语言体系 ◆MATLAB工具箱 (1) 数值计算和符号计算功能 MATLAB以矩阵作为数据操作的基本单位，还提供了十分丰富的数值计算函数。 MATLAB先后和著名的符号计算语言Maple与MuPAD（从MATLAB 2008b开始使用MuPAD）相结合，使得MATLAB具有符号计算功能。 (2) 绘图功能 可以绘制二维和三维图形。 MATLAB提供了两个层次的绘图操作： ●对图形句柄进行的低层绘图操作； ●建立在低层绘图操作之上的高层绘图操作。 (3) 语言体系 MATLAB具有程序结构控制、函数调用、数据结构、输入输出、面向对象等程序语言特征，而且简单易学、编程效率高。 MATLAB是解释性语言，不能脱离MATLAB环境而独立运行。 (4) MATLAB工具箱 MATLAB包含两部分内容：基本部分和各种可选的工具箱。 MATLAB工具箱分为两大类：功能性工具箱和学科性工具箱。 基本部分 构成MATLAB的核心内容，也是使用和构造工具箱的基础。 功能性工具箱 主要用来扩充其符号计算功能、可视建模仿真功能及文字处理功能等。 学科性工具箱 Control System Toolbox 控制系统工具箱 Signal Processing Toolbox 信号处理工具箱 Neural Network Toolbox 神经网络工具箱 Optimization Toolbox 最优化工具箱 Financial Toolbox 金融工具箱 Statistics Toolbox 统计学工具箱 开始→工具箱 MATLAB具备很强的开放性 除内部函数外，所有MATLAB基本文件和各工具箱文件都是可读、可改的源文件，用户可通过对源文件的修改或加入自己编写的文件去构成新的专用工具箱。

MATLAB程序设计及应用(第二版)课后实验答案

。 Matlab课后实验题答案 《 实验一 MATLAB运算基础 1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。 (1) 0 12 2sin85 1 z e = + (2) 21 ln( 2 z x =，其中 212 0.455 i x + ??=?? -?? (3) 0.30.3 3 0.3 sin(0.3)ln, 3.0, 2.9,,2.9,3.0 22 a a e e a z a a - -+ =++=--

(4) 22 42011 122123t t z t t t t t ?≤

) 3. 设有矩阵A 和B 1234 53 166789101769,11 121314150 23416171819209 7021222324254 1311A B ???? ????-??? ?????==-??? ? ???????????? (1) 求它们的乘积C 。 (2) 将矩阵C 的右下角3×2子矩阵赋给D 。 (3) 查看MATLAB 工作空间的使用情况。 4. 完成下列操作： (1) 求[100,999]之间能被21整除的数的个数。 (2) 建立一个字符串向量，删除其中的大写字母。 解：(1) 结果：

MATLAB程序设计与应用课后习题答案

西安科技大学MATLAB程序设计 专业：信息与计算科学 班级：1001班 学号：1008060129 姓名：刘仲能 2012年6月27日

实验一 2.已知： ??? ? ? ?????-= 765 3 8773443412A ，???? ? ?????--=72 3 302131 B 求下列表达式的值： （1）A+6*B 和A-B+I （其中I 为单位矩阵） （2）A*B 和A.*B （3）A^3和A.^3 （4）A/B 及B\A （5）[A,B]和 [A([1,3],:);B^2]

3.设有矩阵A 和B ????? ?? ? ????????= 2524 23 22 21 2019181716151413121110987654321A ，??????? ? ????????--=1113 4 079423096171603 B （1） 求它们的乘积 C 。 （2） 将矩阵C 的右下角3×2子矩阵赋给D 。 （3） 查看MATLAB 工作空间的使用情况 (1) (2) （3）

4.完成下列操作 （1）求[100,999]之间能被21整除的数的个数。（2）建立一个字符串向量，删除其中的大写字母。 (1)(2)

实验二 3.建立一个5×5矩阵，求它的行列式值、迹、秩和范数。 运行截图： A 矩阵的行列式值、迹、秩分别如下： 范数如下： 4.已知 ???? ? ?????--= 58 8 1252018629A 求A 的特征值及特征向量，并分析其数学意义。 运行截图：

5.下面是一个线性方程组：???? ??????= ???? ? ????????? ? ?????52.067.095.06/15 /14 /15 /14/13 /14/13/12 /1321x x x （1）求方程的解； （2）将方程右边向量元素改为0.53，在求解，并比较的变化和解的相对 变化； （3）计算系数矩阵A 的条件数并分析结论。 （2） 变大，其解中，相对未变化前的的解：x1变大，x2变小，x3变大。 （3） 由于A 矩阵的条件数很大，故当线性方程组中的b 变大时，x 也将发生很大的变

Matlab程序设计与应用第二版刘卫国

实验 T1： %%第一小题 z1=2*sin(85*pi/180)/(1+exp(2)) %%第二小题x=[2,1+2i;-0.45,5]; z2 %%第三小题 a=-3.0:0.1:3.0; 0.3)+log((0.3+a)/2) %%第四题 t=0:0.5:2.5 z4=(t>=0&t<1).*(t.A2)+(t>=1 &t<2).*(t 八2-1)+(t>=2&t<3).*(t.A2-2*t+1) T2： A=[12,34,-4;34,7,87;3,65,7] B=[1,3,-1;2,0,3;3,-2,7] disp ('A+6*B='); disp(A+6*B); disp('A-B+I=');disp(A-B+eye(3)); disp('A*B=');

disp(A*B); disp('A.*B='); disp(A.*B); disp('A A3='); disp(AT); disp('A.A3二')； disp(A.A3); disp('A/B='); disp(A/B); disp('B\A='); disp(B\A); disp('[A,B]='); disp([A,B]); disp('[A([1,3],:);BA2]='); disp([A([1,3],:);BA2]); T3: z=1:25; A=reshape(z,5,5)'; B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11];C=A*BD=C(3:5,2:3) T4-1: a=100:999; b=find(rem(a,21)==0); c=length(b)

MATLAB程序设计与应用(第二版)实验参考答案

MATLAB程序设计与应用（第二版）实验参考答案 %实验一MATLAB运算基础 %第一题 %（1） z1=2*sin(85*pi/180)/(1+exp(2)) %（2） x=[2,1+2i;-0.45,5]; z2=0.5*log(x+sqrt(1+x.^2)) %(3) a=-3.0:0.1:3.0; z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2) %(4) t=0:0.5:2.5; z4=t.^2.*(t>=0&t<1)+(t.^2-1).*(t>=1&t<2)+(t.^2-2*t+1).*(t>=2&t<3) %第二题 A=[12 34 -4;34 7 87;3 65 7]; B=[1 3 -1;2 0 3;3 -2 7]; A+6*B A-B+eye(size(A)) A*B A.*B A^3 A.^3 A/B B\A [A,B] [A([1,3],:);B^2] %第三题 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 11] C=A*B F=size(C) D=C(F(1)-2:F(1),F(2)-1:F(2)) whos %第四题 %(1): A=100:999; B=rem(A,21); C=length(find(B==0)) %(2): A='lsdhKSDLKklsdkl';

k=find(A>='A'&A<='Z'); A(k)=[] %实验二MATLAB矩阵分析与处理 %第一题 E=eye(3); R=rand(3,2); O=zeros(2,3); S=diag([2,3]); A=[E,R;O,S]; A^2 B=[E,(R+R*S);O,S^2] %第二题 H=hilb(5) P=pascal(5) Hh=det(H) Hp=det(P) Th=cond(H) Tp=cond(P) %第三题: A=fix(10*rand(5)) H=det(A) Trace=trace(A) Rank=rank(A) Norm=norm(A) %第四题: A=[-29,6,18;20,5,12;-8,8,5] [V,D]=eig(A) %数学意义略 %第五题方法一: %(1): A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/6]; b=[0.95,0.67,0.52]'; x=inv(A)*b %(2): B=[0.95,0.67,0.53]'; x=inv(A)*B %(3): cond(A) %第五题方法二: A=hilb(4) A(:,1)=[] A(4,:)=[] B=[0.95,0.67,0.52]';

(完整word版)MATLAB程序设计及应用(第二版)课后实验答案

Matlab课后实验题答案 实验一 MATLAB运算基础 1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。 (1) 0 12 2sin85 1 z e = + (2) 21 ln( 2 z x =+，其中 212 0.455 i x + ??=?? -?? (3) 0.30.3 3 0.3 sin(0.3)ln, 3.0, 2.9,,2.9,3.0 22 a a e e a z a a - -+ =++=--L (4) 2 2 4 2 01 112 2123 t t z t t t t t ?≤< ? =-≤< ? ?-+≤< ? ，其中t=0:0.5:2.5 2. 已知：

1234413134787,2033657327A B --???? ????==???? ????-???? 求下列表达式的值： (1) A+6*B 和A-B+I （其中I 为单位矩阵） (2) A*B 和A.*B (3) A^3和A.^3 (4) A/B 及B\A (5) [A,B]和[A([1,3],:);B^2] 解： 3. 设有矩阵A 和B 1234 53 166789101769,11 121314150 23416171819209 7021222324254 1311A B ???? ????-??? ?????==-??? ? ???????????? (1) 求它们的乘积C 。 (2) 将矩阵C 的右下角3×2子矩阵赋给D 。 (3) 查看MATLAB 工作空间的使用情况。

4. 完成下列操作： (1) 求[100,999]之间能被21整除的数的个数。 (2) 建立一个字符串向量，删除其中的大写字母。 解：(1) 结果： (2). 建立一个字符串向量 例如： ch='ABC123d4e56Fg9';则要求结果是： 实验二 MATLAB 矩阵分析与处理 1. 设有分块矩阵33 322322E R A O S ?????? =? ??? ，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩

电子书_MATLAB程序设计与应用(第二版)

阅阅读读时时：：请请选选择择““视视图图 || 文文档档结结构构图图””，，弹弹出出文文档档中中的的标标题题链链接接。。 数学软件 MATLAB 程序设计与应用

第1章MATLAB系统环境 1.1 MATLAB概貌 1.2 MATLAB环境的准备 1.3 MATLAB操作界面 1.4 MATLAB帮助系统 自上世纪80年代以来，出现了科学计算语言，亦称数学软件。 MATLAB Mathematica Mathcad Maple LINDO LINGO 1.1 MATLAB概貌p3 MATLAB 是MATrix LABoratory（矩阵实验室）的缩写。 1984年由Math Works 公司推出，现已成为国际公认的优秀的工程应用开发环境，是影响最大，流行最广的科学计算语言。 1.1.1 MATLAB的发展 2. 例如MATLAB7.6的建造编号是R2008a。说明MATLAB7.6与MATLAB2008a是等同的；

3. 对于建造编号，正规化以后，每年出两个版本。一般来说。a是测试版，b是正式版。a是前半 年出，b是后半年出。 教材采用MATLAB7.0（R14，2004） 实验室采用MATLAB 7.8（R2009a，2009.3，汉化） 1.1.2 MATLAB的主要功能p4 ◆数值计算和符号计算功能 ◆绘图功能 ◆语言体系 ◆MATLAB工具箱 (1) 数值计算和符号计算功能 MATLAB以矩阵作为数据操作的基本单位，还提供了十分丰富的数值计算函数。 MATLAB先后和著名的符号计算语言Maple与MuPAD（从MATLAB 2008b开始使用MuPAD）相结合，使得MATLAB具有符号计算功能。 (2) 绘图功能 可以绘制二维和三维图形。 MATLAB提供了两个层次的绘图操作： ●对图形句柄进行的低层绘图操作； ●建立在低层绘图操作之上的高层绘图操作。 (3) 语言体系 MATLAB具有程序结构控制、函数调用、数据结构、输入输出、面向对象等程序语言特征，而且简单易学、编程效率高。 MATLAB是解释性语言，不能脱离MATLAB环境而独立运行。 (4) MATLAB工具箱 MATLAB包含两部分内容：基本部分和各种可选的工具箱。 MATLAB工具箱分为两大类：功能性工具箱和学科性工具箱。 基本部分 构成MATLAB的核心内容，也是使用和构造工具箱的基础。 功能性工具箱 主要用来扩充其符号计算功能、可视建模仿真功能及文字处理功能等。 学科性工具箱 Control System Toolbox 控制系统工具箱 Signal Processing Toolbox 信号处理工具箱 Neural Network Toolbox 神经网络工具箱 Optimization Toolbox 最优化工具箱 Financial Toolbox 金融工具箱 Statistics Toolbox 统计学工具箱 开始→工具箱 MATLAB具备很强的开放性 除内部函数外，所有MATLAB基本文件和各工具箱文件都是可读、可改的源文件，用户可通过对源文件的修改或加入自己编写的文件去构成新的专用工具箱。

Matlab程序设计与应用第二版刘卫国

实验一：T1： %%第一小题 z1=2*sin(85*pi/180)/(1+exp(2)) %%第二小题 x=[2,1+2i;-0.45,5]; z2 %%第三小题 a=-3.0:0.1:3.0; 0.3)+log((0.3+a)/2) %%第四题 t=0:0.5:2.5 z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t .^2-1)+(t>=2&t<3).*(t.^2-2*t+1) T2： A=[12,34,-4;34,7,87;3,65,7] B=[1,3,-1;2,0,3;3,-2,7] disp ('A+6*B='); disp(A+6*B); disp('A-B+I=');disp(A-B+eye(3)); disp('A*B=');

disp(A*B); disp('A.*B='); disp(A.*B); disp('A^3='); disp(A^3); disp('A.^3='); disp(A.^3); disp('A/B='); disp(A/B); disp('B\A='); disp(B\A); disp('[A,B]='); disp([A,B]); disp('[A([1,3],:);B^2]='); disp([A([1,3],:);B^2]); T3: z=1:25; A=reshape(z,5,5)'; B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11];C=A*BD=C(3:5,2:3) T4-1: a=100:999;

b=find(rem(a,21)==0); c=length(b) T4-2: a=input('请输入一个字符串：','s'); b=find(a>='A'&a<='Z'); a(b)=[]; disp(a); 实验二： T1： E=eye(3),R=rand(3,2),O=zeros(2,3),S=d iag([1,2]); A=[E,R;O,S] disp('A^2='); disp(A^2); disp('[E,R+RS;O,S^2]'); B=[E,R+R*S;O,S^2] T2： H=hilb(5) P=pascal(5)Hh=det(H) Hp=det(P) Th=cond(H)

最新MATLAB程序设计与应用(第二版)课后实验答案

Matlab 课后实验题答案 实验一 MATLAB 运算基础 1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。 (1) 0 12 2sin851z e =+ (2) 221 ln(1)2z x x = ++，其中2120.45 5i x +??=??-?? (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+= ++=-- (4) 22 42011 122123t t z t t t t t ?≤=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3) .*(t.^2-2*t+1) 2. 已知：

1234413134787,2033657327A B --???? ????==???? ????-???? 求下列表达式的值： (1) A+6*B 和A-B+I （其中I 为单位矩阵） (2) A*B 和A.*B (3) A^3和A.^3 (4) A/B 及B\A (5) [A,B]和[A([1,3],:);B^2] 解： M 文件： A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7]; A+6.*B A-B+eye(3) A*B A.*B A^3 A.^3 A/B B\A [A,B] [A([1,3],:);B^2] 3. 设有矩阵A 和B 1234 53 166789101769,11 121314150 23416171819209 7021222324254 1311A B ???? ????-??? ?????==-??? ? ???????????? (1) 求它们的乘积C 。 (2) 将矩阵C 的右下角3×2子矩阵赋给D 。 (3) 查看MATLAB 工作空间的使用情况。 解：. 运算结果： E=(reshape(1:1:25,5,5))';F=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11]; C= E*F H=C(3:5,2:3) C = 93 150 77

MATLAB程序设计与应用

实验一 MATLAB 运算基础 第二题：已知 A=??????????7653877344-3412 B=??? ? ??????72-33021-31 求下列问题： 1) A+6*B 和A-B+I 2) A*B 和A.*B 3) A^3 A.^3 4) A/B B\A 5) [A,B] [A(1,3),:B.^2] 解： >> A=[12 34 -4;34 7 87;3 65 7]; B=[1 3 -1;2 0 3;3 -2 7]; I=[1 0 0;0 1 0;0 0 1]; （1） >> A+6*B ans = 18 52 -10 46 7 105 21 53 49 >> A-B+I ans = 12 31 -3 32 8 84 0 67 1 （2） >> A*B ans = 68 44 62 309 -72 596 154 -5 241 >> A.*B

ans = 12 102 4 68 0 261 9 -130 49 （3） >> A^3 ans = 37226 233824 48604 247370 149188 600766 78688 454142 118820 >> A.^3 ans = 1728 39304 -64 39304 343 658503 27 274625 343 （4） >> A/B ans = 16.4000 -13.6000 7.6000 35.8000 -76.2000 50.2000 67.0000 -134.0000 68.0000 >> B\A ans = 109.4000 -131.2000 322.8000 -53.0000 85.0000 -171.0000 -61.6000 89.8000 -186.2000 （5） >> [A,B] ans =

matlab程序设计与应用(第二版)第三章部分课后答案解析

第三章1. (1A=eye(3 (2A=100+100*rand(5,6 (3A=1+sqrt(0.2*randn(10,50 (4B=ones(size(A (5A+30*eye(size(A (6B=diag(diag(A 2. B=rot90(A C=rot90(A,-1 3. B=inv(A ;A的逆矩阵 C=det(A ;A的行列式的值 D=A*B E=B*A D=E 因此A与A-1是互逆的。 4. A=[4 2 -1;3 -1 2;12 3 0]; b=[2;10;8];

x=inv(A*b x = -6.0000 26.6667 27.3333 5. (1 diag(A ;主对角线元素ans = 1 1 5 9 triu(A ;上三角阵 ans = 1 -1 2 3 0 1 -4 2 0 0 5 2

0 0 0 9 tril(A ;下三角阵ans = 1 0 0 0 5 1 0 0 3 0 5 0 11 15 0 9 rank(A ;秩ans = 4 norm(A ;范数 ans = 21.3005 cond(A ;条件数ans = 11.1739 trace(A ;迹 ans = 16 (2略

6. A=[1 1 0.5;1 1 0.25;0.5 0.25 2] A = 1.0000 1.0000 0.5000 1.0000 1.0000 0.2500 0.5000 0.2500 2.0000 [V,D]=eig(A V = 0.7212 0.4443 0.5315 -0.6863 0.5621 0.4615 -0.0937 -0.6976 0.7103 D = -0.0166 0 0 0 1.4801 0 0 0 2.5365

MATLAB程序设计和应用课后习题答案解析

专业技术资料分享 西安科技大学MATLAB程序设计 专业：信息与计算科学 班级： 1001班 学号：1008060129 姓名：刘仲能 2012年6月27日

实验一 2.已知： ?? ??? ?????-=76538773443412A ，???? ? ?????--=723302131B 求下列表达式的值： （1）A+6*B 和A-B+I （其中I 为单位矩阵） （2）A*B 和A.*B （3）A^3和A.^3 （4）A/B 及B\A （5）[A,B]和 [A([1,3],:);B^2]

3.设有矩阵A 和B ????? ?? ?????????=25242322212019181716151413121110987654321A ，??????? ?????????--= 11134079423096171603B （1） 求它们的乘积C 。 （2） 将矩阵C 的右下角3×2子矩阵赋给D 。 （3） 查看MATLAB 工作空间的使用情况 (1) (2) （3）

4.完成下列操作 （1）求[100,999]之间能被21整除的数的个数。 （2）建立一个字符串向量，删除其中的大写字母。 (1) (2) 实验二 3.建立一个5×5矩阵，求它的行列式值、迹、秩和范数。

运行截图： A 矩阵的行列式值、迹、秩分别如下： 范数如下： 4.已知 ?? ?? ? ?????--=5881252018629A 求A 的特征值及特征向量，并分析其数学意义。 运行截图：

5.下面是一个线性方程组：?? ??? ?????=???? ????????????????52.067.095.06/15/14/15/14/13/14/13/12/1321x x x （1） 求方程的解； （2） 将方程右边向量元素 改为0.53 ，在求解，并比较 的变化和解的相 对变化； （3） 计算系数矩阵A 的条件数并分析结论。 （2） 变大，其解中，相对未变化前的 的解：x1变大，x2变小，x3变大。 （3） 由于A 矩阵的条件数很大，故当线性方程组中的b 变大时，x 也将发生很大的变化，即数值稳定性较差。 实验三 3.硅谷公司员工的工资计算方法如下： （1）工作时数超过120小时者，超过部分加发15%； （2）工作时数低于60小时者，扣发700元； （3）其余按每小时84元计发。

Matlab程序设计及应用实验讲义

Matlab程序设计及应用实验讲义 自编 电子科学与工程系 2013.02

实验一 MATLAB 环境与命令窗口 1、实验目的 1）熟悉MATLAB 的操作环境及基本操作方法； 2）掌握MATLAB 的搜索路径及其设置方法； 3）熟悉MATLAB 帮助信息的查阅方法； 2、实验主要仪器设备和材料 计算机PC 一台 2010a 软件 3、实验内容和原理 1、先建立自己的工作目录，再将自己的工作目录设置到MATLAB 的搜索路径下，再试验用help 命令能否查询到自己的工作目录。 2、在MATLAB 环境下验证下面几个例子，并总结MATLAB 的主要优点。 1）绘制正弦曲线和余弦曲线 2）求方程432379230x x x ++-=的全部根 3）求积分()1 0ln 1x x dx +? 4）求解线性方程组234832245917x y z x y z x y z -+=?? ++=??+-=? 3、利用MATLAB 的帮助功能分别查询inv 、plot 、 max 、round 等函数的功能与用法 4、完成下列操作： 1）在MATLAB 命令窗口输入以下命令： x=0:pi/10:2*pi;

y=sin(x); 2）在工作空间窗口选择变量y，再在工作空间窗口选择绘图菜单命令或在工具栏中单击绘图命令按钮，绘制变量的图形，并分析图形的含义。 5、访问MathsWorks公司的主页，查询有关MATLAB的产品信息。 思考与练习 1、如何启动和退出MATLAB的集成环境？ 2、简述MATLAB的主要功能。 3、如果一个MATLAB命令包含的字符很多，需要分成多行输入，该如何处理？ 4、help命令和look for命令有何区别？ 5、在MATLAB环境下，建立一个变量fac,同时又在当前目录下建立了一个M 文件fac.m，如果需要运行fac.m文件，该如何处理？

MATLAB程序设计与应用(刘卫国编)课后实验答案

实验一 MATLAB 运算基础 1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。 (1) 012 2sin851z e =+ (2) 21ln(2 z x =+，其中2 120.45 5i x +??=? ?-?? (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+= ++=-- (4) 2242011 122123t t z t t t t t ?≤

- 2. 已知： 1234413134787,2033657327A B --???? ????==???? ????-???? 求下列表达式的值： (1) A+6*B 和A-B+I （其中I 为单位矩阵） (2) A*B 和A.*B (3) A^3和A.^3 (4) A/B 及B\A

(5) [A,B]和[A([1,3],:);B^2] @ 解： 运算结果：

3. 设有矩阵A和B

1234 53 166789101769,11 121314150 23416171819209 7021222324254 1311A B ???? ????-??? ?????==-??? ? ???????????? (1) 求它们的乘积C 。 (2) 将矩阵C 的右下角3×2子矩阵赋给D 。 ( (3) 查看MATLAB 工作空间的使用情况。 解：. 运算结果： ! 4. 完成下列操作： (1) 求[100,999]之间能被21整除的数的个数。 (2) 建立一个字符串向量，删除其中的大写字母。 解：(1) 结果：