专业基础实践Matlab运算与二阶动态电路分析

1、武汉理工大学专业基础实践课程设计说明书课程设计任务书学生姓名： 专业班级：电信1306班指导教师： 工作单位：信息工程学院课程设计名称：专业基础实践课程设计题目：matlab运算与二阶动态电路分析初始条件：1. matlab7.1以上版本软件；2 .专业基础实践辅导资料：“matlab语言基础及使用入门” 、“ matlab及在电子信息课程中的应用”、“线性 代数”及“信号处理类”相关书籍等；3 .先修课程：高等数学、线性代数、电路、 matlab应用实践及信号处理类基础课程等。要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）1 .实践内容：根据指导老师给定的7

2、套题目，按规定选择其中 1套独立完成；2 .本专业基础实践统一技术要求：研读辅导资料对应章节，对选定的设计题目进行理论分析，完成针对 具体设计部分的原理分析、建模、必要的推导和可行性分析，画出程序设计框图，编写程序代码（含 注释），上机调试运行程序，记录实验结果（含计算结果和图表等），并对实验结果进行分析和总结。具体设计要求包括：初步了解matlab、熟悉matlab界面、进行简单操作等；matlab的数值计算：创建矩阵、矩阵运算、多项式运算、线性方程组、数值统计等；基本绘图函数：了解 plot, plot3, mesh, surf等，要求掌握以上绘图函数的用法、简单图形标注、简 单颜色设定等

3、； 使用文本编辑器编辑 m文件，会函数调用等； 能完成简单电路的 matlab编程分析；按要求参加专业基础实践的实验演示和答辩等。3 .课程设计说明书按学校“课程设计工作规范”中的“统一书写格式”撰写，具体包括：目录； 与设计题目相关的理论分析、归纳和总结；与设计内容相关的原理分析、建模、推导、可行性分析；程序设计框图、程序代码（含注释）、程序运行结果和图表、实验结果分析和总结；课程设计的心得体会（至少 500字）；参考文献（不少于5篇）；其它必要内容等。时间安科卜：1周指导教师签名： 系主任（或责任教师）签名:武汉理工大学专业基础实践课程设计说明书inimatlab运算与二阶动态电路分析 3

4、1 . matlab 概述 31.1 matlab 简介 31.2 matlab 的功能 31.3 matlab 的应用 42 .设计内容 52.1 题 一 52.2 题二 62.3 题三 72.4 题四 82.5 题五 142.6 题六 152.7 题七 172.8 题八 192.9 题九 202.10 题十 222.11 题 h一 243 .心得体会 294 .参考文献 30武汉理工大学专业基础实践课程设计说明书matlab运算与二阶动态电路分析1. matlab 概述1.1 matlab 简介matlab是美国mathworks公司出品的商业数学软件，用于算法开发、数据可视化、数 据分析

5、以及数值计算的高级技术计算语言和交互式环境，主要包括 matlab和simulink两 大部分。matlab将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等 诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效 数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式 程序设计语言(如c、fortran)的编辑模式，代表了当今国际科学计算软件的先进水平。1.2 matlab的功能(1)强大的科学计算机数据处理能力matlab是一个包含大量计算算法的集合。其拥有600多个工程中要用到的数学运算函 数，可以方便的实现用户所需的

6、各种计算功能。matlab勺这些函数集包括从最简单最基本的函数到诸如矩阵，特征向量、快速傅立叶变换的复杂函数。函数所能解决的问题其大致包括 矩阵运算和线性方程组的求解、微分方程及偏微分方程的组的求解、符号运算、傅立叶变换 和数据的统计分析、工程中的优化问题、稀疏矩阵运算、复数的各种运算、三角函数和其他 初等数学运算、多维数组操作以及建模动态仿真等。(2)出色的图形处理功能图形处理功能matlab自产生之日起就具有方便的数据可视化功能， 以将向量和矩阵用 图形表现出来，并且可以对图形进行标注和打印。高层次的作图包括二维和三维的可视化、 图象处理、动画和表达式作图。可用于科学计算和工程绘图。应用广

7、泛的模块集合工具箱matlab对许多专门的领域都开发了功能强大的模块集和工具箱。它们都是由特定领域的专家开发的，用户可以直接使用工具箱学习、应用和评估不同的方法而不需要自己编写代码。目前，matlab已经把工具箱延伸到了科学研究和工程应用的诸多领域，诸如数据采集数据库接口、概率统计、样条拟合、优化算法、偏微分方程求解、神经网络、小波分析、信 号处理、图像处理、系统辨识、控制系统设计、lmi控制、鲁棒控制、模型预测、模糊逻辑、金融分析、地图工具、非线性控制设计、实时快速原型及半物理仿真、嵌入式系统开发、定 点仿真、dsp与通讯、电力系统仿真等，都在工具箱家族中有了自己的一席之地。(4)实用的程序

8、接口和发布平台新版本的matlab可以利用matlab编译器和c/c+数学库以及图形库，将自己的 matlab程序自动转换为独立于 matlab运行的c和c+的代码。允许用户编写可以和 matlab进行交互的c或c+语言程序。matlab的一个重要特色就是具有一套程序扩展 系统和一组称之为工具箱的特殊应用子程序。(5)应用软件开发在开发环境中，使用户更方便地控制多个文件和图形窗口；在编程方面支持了函数嵌套,有条件中断等；在图形化方面，有了更强大的图形标注和处理功能，包括对性对起连接注释 等；在输入输出方面，可以直接向 excel和hdf5进行连接。1.3 matlab的应用matlab产品族可

9、以用来进行以下各种工作： 数值分析 数值和符号计算 工程与科学绘图控制系统的设计与仿真数字图像处理技术数字信号处理技术 通讯系统设计与仿真 财务与金融工程 管理与调度优化计算(运筹学)matlab的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、 财务建模和分析以及计算生物学等众多应用领域。 附加的工具箱扩展matlab环境，以解决 这些应用领域内特定类型的问题。2.设计内容2.1题一：求100,200之间第一个能被 21整除的整数2.1.1分析利用for循环语句控制，并结合rem函数求得符合条件的数2.1.2程序编码 c:u写e5adminkt同口。.2000语句最终达到

10、题目要求，求得y=20702.4.12.4 题四：分析下面每条指令的功能并运行，观察执行结果x=0:0.25:3;x ), ylabel( y ,x：ie();y=x*exp(y=x.*exp(-x);plot(x,y),xlabel(武汉理工大学专业基础实践课程设计说明书2.4.1.1 分析x=0:0.25:3;得到以0.25为步长，取得0到3内的所有实数，确定自变量 x的范围。利 用plot函数画出函数y的图像，并设置图像的一些基本性质。2.4.1.2结果显示图2.4.1题4.1的结果显示2.4.2题4.2程序编码c;usersadministratardeskto4雪u ntitled4

11、2,m*ledit text go cell tools debug desktop window help0乂,触，也 al 目惹11a=zeros(2f 5)a :i：5l=abs(a)l islogical d x=a(l)先生成一个2行5列的全。矩阵务将且重新哦值常将a中绝对值大于4的数置1，其余为。常判断l是否为逻辑书机为逻辑数，将a中相同位置的元素期出组成一个新的矩驻scriptbi 5 col 430v- i图2.4.2题4.2的程序编码2.4.2.1 分析zeros数建立一个2*5全0矩阵赋给a,再将-4到5内的整数赋给a,形成新的矩阵。l=abs(a) 4,获得新的矩阵l,通

12、过逻辑函数islogical得到新的矩阵赋给x武汉理工大学专业基础实践课程设计说明书2.4.2.2结果显示图2.4.3题4.2的结果显示2.4.3 题4.3程序编码 c:usersadministratordesktop专业塞朴实任程序4huntitled43,m* ；*图2.4.4题4.3的程序编码2.4.3.1 分析先生成1至5和6至10的两行五列矩阵赋给 a,再利用pow2函数的使用获得新的矩阵 pow2(a)的使用方法为a=2aa2.4.3.2结果显示武汉理工大学专业基础实践课程设计说明书图2.4.5题4.3的结果显示2.4.4 题4.4程序编码图2.4.6题4.4的程序编码2.4.4

13、.1 分析zeros数建立一个3*2全0矩阵赋给a,再将1到6内的整数赋给a,形成新的矩阵。将a中的每个元素以a=a*(1+i)进行运算，获得新的矩阵a,通过转置运算和共腕转置运算得到新的矩阵赋给a1和b1。2.4.4.2结果显示图2.4.7题4.4的结果显示(1)武汉理工大学专业基础实践课程设计说明书command windowa =1426362.4.5.1 分析ones3数建立一个2*3全1矩阵赋给a,建立一个2*2全1矩阵赋给b,再建立一个4阶单位矩alblloooo2.00003. 00001. 0000工 0000l 00004. 0000+ 1,00001+ 2,00001+ 3

14、, ooooi+ 1,ooooi+ l ooooil ooooi-4. ooooi工 00005, 00006. 0000l 00005. 0000工 00005. 0000+ 工 ooooi+ 5. ooooi4- 6, ooooi+ 2.ooooi+ 5. ooooiz ooooi5.000013.0000 + 3.ooooi6. 0000 + 6. ooooi3. 000。- 3. ooooi5. 0000 - 6. ooooia图2.4.8题4.4的结果显示（2）2.4.5 题4.5程序编码ca,ur$ad mi 口 i strat口 rde&kt叩t专业实注民用第4毁u ntitle

15、cms.m生成一个2*3的全为1的矩阵%生成一个2*2的全为1的矩阵需生成一个4阶单位矩阵%即匚矩陈中的对角元素作为列向里e=rei：nztc. r ?/将u矩阵复制为1*3块的矩胜口 |回file edrt text go cell tools debug desktop window help*目目目i 目lox噱遨qa=ones 3】b=ones(2)c=5ye(4)d=diag coscriptcol 33ovr图2.4.9题4.5的程序编码武汉理工大学专业基础实践课程设计说明书阵赋给c；禾ij用diag函数对c取对角元素作列向量，形成新的矩阵 d。利用repmata数将c矩阵复制为1

16、*3块的矩阵，获得新的矩阵e。2.4.5.2结果显示图2.4.10题4.5的结果显示(1)command windowc =1000010000100001d =111110001000100001000100010q001000100010a000100010001图2.4.11题4.5的结果显示(2)2.4.6 结果分析matlab中有多种生成矩阵的方法，以及矩阵的一些操作函数等。写代码的时候要注意分号的使用，若一个输出表达式后有分号，则只计算而不在command窗口中显示出来。2.5 题五：利用 whileend1环结构编程求由1+3+5+101的值2.5.3 分析while语句的一般格

17、式为：while（条件）循环体语句end。其执行过程为：若条件成立,则执行循环体语句，执行后再判断条件是否成立，如果不成立则跳出循环2.5.4 程序编码显不c:ljsersadministrator口e5ktop%专业基牡实践程序5sjllntitled5.m* l 1 | 回图2.5.1题五的程序编码2.5.5 程序流程图图2.5.2题五的程序流程2.5.6 结果显不corrirriand window sum -42601一图2.5.3题五的结果显示2.5.7 结果分析该题目从i=1开始，每循环一次其值加2,由sum=sun+i起到累加器的效果，一直到i101 结束循环，最终结果为sum=

18、2601。2.6 题六：随机生成一个 10000刈的矩阵，并绘制其柱状图。2.6.3 分析利用rand函数生成本题矩阵。该题生成一个 10000*3的矩阵赋给a,再利用bar函数画 一个a的二维条形图，即本题要求的柱状图。2.6.4 程序编码显示c:usersadministratordesktop-专业基社实笈程序6cuntitledsl.mfile edit text go cell tools debug desktop window help学：8 ci | - lo | + | + ii * | 靖遨 q1 a=rand(10000, 3)上2 -bar (a)回 |il -titl

19、ed电信j306班陈雪燕)|二 i3 u&/ga，of 启 found scriptln m20 0vr .二 1图2.6.1题六的程序编码由于10000过大，在有限的区域内显示10000*3的柱状图不清晰、不明显，为了方便观 察柱状图的分布，将以上程序稍加修改，利用axis函数显示横轴由0到10,纵轴由0到1的区域的柱状图。修改程序显示如图2.6.2所示：武汉理工大学专业基础实践课程设计说明书图2.6.4题六改进的结果显示图2.6.2题六改进的程序编码2.6.3柱状图显示图2.6.3题六结果显示武汉理工大学专业基础实践课程设计说明书2.6.4结果分析随机生成矩阵可用rand (m,n)来表示

20、，其中m代表行数，n代表列数。绘制柱状图可用 bar函数表示。由图2.6.3可观察得10000*3的柱状图全局，但因其容量过大无法清晰观察， 为此我截取一部分数据进行柱状图展示，如图 2.6.4方便观察。2.7题七：利用switch cas编写以下程序某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price来表示)：price200无折扣200price5005% 折扣500price100010% 折扣price 之 100015% 折扣输入所售商品的价格，求其实际销售价格。2.7.1 分析switch语句根据表达式的取值不同，分别执行不同的语句，其语句格式为： switch表达

21、式case表达式1语句组1case表达式m语句组motherwise语句组nend当表达式的值等于表达式1的值时，执行语句组1,当表达式的值等于表达式2的值时, 执行语句组2,，当表达式的值等于表达式 m的值时，执行语句组 m,当表达式的值不等 于case所列的表达式的值时，执行语句组 no当任意一个分支的语句执行完成后，直接执行 switch语句的下一句。这个题需要建立mfile文件，来将进行求解。武汉理工大学专业基础实践课程设计说明书2.7.2 m-file文件程序廿 c:usersadministratarde与ktop,、专业基叱实法理序作7题untitled7.mfile edit

22、text go cell tools debug desktop window help学目白。|总值含t* ，用|见母人用 唱唱| - l0 | + | + l1 * |嘴蟆 q1 -pric=input。请输入商品价格)；口2 - switch fix(price/100)3 -case【口，1)4 -rate=0 ：5 -case % 3, i6 -racte=5/100;7 -case nuji2cell (5： 9)8 -rate=i0/100;9 -otherwise10 -tate=15/100：11 -end12 -price=prxce*(1-rate)scriptln 12

23、 col 21ovr .：图2.7.1题七的程序编码2.7.3结果显示+ 口 x清输入商品伯格100price 100请输入商品价格35price =327. 7500请输入商品价格78。price -702清输久商品价格150。price =1275代请输入商品伯格2.7.4结果分析当输入 100,巾x(price/100)为 1,符合 case 0,1, rate 等于 0, price=100*(1-0)=100。当输入 345jix(price/100)为 3,符合 case 2,3,4, rate等于 5%,price=345*(1-5%)=327.7500。当输入 780 , fi

24、x(price/100)为 7 ,符合 case num2cell(5:9) , rate 等于 10% , price=780*(1-10%)=702。当输入 1500,fix(price/100)为 15,符合 otherwise,rate等于 15%,price=1500*(1-15%)=1275。2.8题八：利用函数调用编写函数文件，求半径 r分别为0.5与5时圆 的面积与周长。2.8.1分析在matlab中新建function-m-file，建立一个函数并保存到 matlab文件夹中，输入程序求得结果。2.8.2函数文件c:usersadministratordesktop专业基壮实

25、法理序策8sl.m*, , ,山图2.8.1题八的程序定义图2.8.2题八的程序编码武汉理工大学专业基础实践课程设计说明书2.8.3结果显不command windowidans0.78543.14160.785478.539831.4159ans78,53981=一图2.8.3题八的结果显示2.8.4结果分析本题先进行函数定义，利用 function定义求面积和求周长公式，再来计算相关值时的结 果。由此可见，半径r为0.5时，面积为0.7854,周长为3.1416;半径r为5时，面积为78.5398, 周长为31.41592.9题九:已知y=e-1cos(t),其中t的范围是0 10,计算y

26、的微分和积分 并给由曲线图。武汉理工大学专业基础实践课程设计说明书2.9.1分析分别利用diff,以及int函数求的函数的微分和积分，利用ezplot画图函数得到函数的图像。注意此处ea(-1),用exp(-1)表示。2.9.2程序编码图2.9.1题九的程序编码2.9.3结果显不command windowans 二e283903570ss4e7*sm(10)/22517s9sl36s524s828390957088497 sin(lo)2251799s136s524sans 二-(328390s57083437*costt)/2251799813635248828390857083487 c

27、os(t)2251799813685248a:图2.9.2题九的结果显示武汉理工大学专业基础实践课程设计说明书2.9.4曲线图展示图2.9.3题九的图形显示2.9.5结果分析本题利用求导函数diff求得已知函数的微分，其结果如图2.9.2的第一个结果所示，由此 可见该操作是正确的；再利用int函数求得已知函数的积分，其结果如图 2.9.2的第二个结果 所示，由此可见该操作也是正确的。2.10 题十：编写程序，用plot(x,y)绘制下图所示曲线a)b)2.10.1 分析观察第一个函数图像可以得到函数解析式为 y=3*x+1 , x的范围为：0=x=5。画图时要标明x, y坐标轴，x轴坐标最小单

28、位是0.5, y轴坐标最小单位是2观察第二个函数图像，由两条直线组成：y1=-x+1, y2=2*x。纵坐标最小单位为：0.2,且范围为1,3,横坐标最小单位为0.1,范围为0,1。利用hold on语句在同一个坐标下画出 两个函数的图像。2.10.2 程序编码2.10.2.1 第一个函数图像的程序编码c:usersadministratordesktcip,1.专业基此实建逞序10asurtitledlo.file edit text go cell tools debug desktop window help学图2.10.1题10.1的程序编码2.10.2.2第二个函数图像的程序编程图2

29、.10.2题10.2的程序编码2.10.3 结果显不图2.10.4题10.2的结果显示图2.10.3题10.1的结果显示2.10.4 结果分析本题我先求得两图的表达式，再根据表达式编写程序，给定 x的范围，再编辑表达式， 最后利用plot函数绘制二维图形得到图2.10.3和图2.10.4与题目所给图形一致，因此该题编 程正确。2.11 题h一：在如图所示的二阶电路中， r=1kq , c=2nf , l=2.5h ,电容 原先已充电。且uc(0j=1cv。在t=0时开关s闭合。试求uc(t), i(t), ul。)， 绘出各曲线图；开关s闭合后的濡为多大？c 丰 ucljut(t)2.11.1

30、 分析由电路图及已知条件展开理论分析，当t0时，电路接通。此时电路的微分方程为:初始条件为:设uc对时间的一阶导数为lcdu d71dudu2 rcdudt(2.(1)i(0)(2.(2)uc=10vc(2.(3)p,则可得微分方程的特征方程为:2lcp2 rcp 1 -0由求根公式求得其根为:(2.5)(2.6)r _ (r)2_12l 2l lc又由已知条件可知，当r0sin6 ,如下图所示。 ot根据欧拉公式:可得:所以有:根据上式可得:图2.11.1 口,日,露露的关系图j c cos 二 j sin 1e = cos jsinp = fejp2 -ojuc u；(p2eplt -p2

31、，)f2-puoj： / -： j )，ji( -：j)=- oe eoe e -2juoo _：trej(tle(t：-e u 00aeqsin(6t 十日)ulu o-?ti = e sin( t)lu o- 0 tt二e sin( t -)oimax(2.(9)(2.(10)(2.(11)(2.(12)(2.(13)(2.(14)(2.(15)(2.(16)2.11.2程序编码图2.11.2题h一的程序编码2.11.3结果显不图2.11.3题十一的结果显示2.11.4曲线图展示图2.11.4题十一的图形显示2.11.5结果分析在2.11.1节中，已对电路进行分析计算，得到了u/t), i(t), 5(t)等的表达式，继而将表达式带入程序，通过编程与运行来计算电流的最大值，最后得出开关s闭合后的imax=0.0051。 再利用plot函数绘制二维图形如图2.11.4所示。3 .心得体会经历了一个星期的查阅资料和理论分析，终于完成了专业基础实践的编程和报告。经历 了这次课程设计，大大的提高了我的操作能力以及分析问题的能力，从中也学到了很多书面 上所没有搞清楚的问题，也熟悉了应用 matlab这个软件来进行程序编程。通过