广州大学matlab大作业

广州大学机电学院电气101

MATLAB大作业

MATLAB是由美国公司发布主要面对科学计算、可视化以的计算环境。它可以将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个视窗环境中，为科学研究众多科学领域提供了一种全面的解决方案，代表了当今国际科学计算软件的先进水平。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。

MATLAB在以下的领域里解决各种问题是一个十分有效的工具：

? 工业研究与开发。

? 数学教学，特别是线性代数。所有基本概念都能涉及。

? 在数值分析和科学计算方面的教学与研究。能够详细地研究和比较各种算法。? 在诸如电子学、控制理论和物理学等工程和科学学科方面的教学与研究。? 在诸如经济学、化学和生物学等有计算问题的所有其他领域中的教学与研究。

这学期我们做了诸多matlab实验，从符号计算及程序设计到一维、二维数组实验，还有图形显示等实验，我们初步掌握了matlab操作方法。我会在后文中用三个例子在三个应用方面着重汇报我的matlab使用心得。

本报告将以如下顺序进行叙述：

一、MATLAB在线性代数方面的应用

1.简单的矩阵的生成

2.常用矩阵的生成

3. 线性方程求解

二、MATLAB在经济学中的应用

价格平衡模型分析

三、MATLAB在三维图形绘制中的应用

1.函数PLOT3命令

2.如何改变视角

四、心得体会

一、MATLAB在线性代数方面的应用

1980年，MATLAB的首创者Cleve Moler博士在New Mexico大学讲授线性代数课程时，看到了用高级语言编程解决工程计算问题的诸多不便，因而构思开发了用Fortran语言编写而成，集命令翻译、工程计算功能于一身的MATLAB软件。在数学上，矩阵是指纵横排列的二维数据表格，最早来自于方程组的系数及常数所构成的方阵。矩阵概念在生产实践中也有许多应用，比如矩阵图法以及保护个人帐号的矩阵卡系统等等。

在第一次上机实验课中我们就做了简单的矩阵实验，下面稍作探讨：

1.简单的矩阵的生成

在MATLAB中，可以采用多种不同的方式生成矩阵。

(1)直接输入矩阵元素

对于较小的简单的矩阵，从键盘上直接输入矩阵是最常用、最方便和最好的数值矩阵创建方法。直接从键盘输入一系列元素生成矩阵，只要遵循下面几个基本原则：

1、矩阵每一行的元素必须用空格或逗号分开；

2、在矩阵中，采用分号或回车表明每一行的结束；

3、整个输入矩阵必须包含在方括号“[]”中。

输入：A=[4,5,7,8;6,1,2,5;3,5,4,6;4,2,4,8]

显示：A = 4 5 7 8

6 1 2 5

3 5

4 6

4 2 4 8

(2)生成大矩阵

在MATLAB中，可以将小矩阵连接起来生成一个较大的矩阵。事实上，前面直接输入法生成矩阵就是将单个元素连接起来生成矩阵。方括号“[]”就是连接算子。

输入：B=[A,A+1;A+2,A+3]

显示：B =

4 5 7 8 5 6 8 9

6 1 2 5

7 2 3 6

3 5

4 6 4 6

5 7

4 2 4 8

5 3 5 9

6 7 9 10 7 8 10 11

8 3 4 7 9 4 5 8

5 7

6 8 6 8

7 9

6 4 6 10

7 5 7 11

2.常用矩阵的生成

下面介绍一些常用矩阵的生成命令：zeros，ones，eye，rand和randn。(1)zeros生成全0阵

调用格式为:B=zeros(n)：生成n n的全0矩阵；如果n不是标量将给出出错信息。

B=zeros(m,n)：生成m n的全0矩阵。

B=zeros(size(A))：生成与矩阵A大小相同的全0阵。

(2)ones生成全1阵

调用格式为：

Y=ones(n)：生成n n的全1矩阵；如果n不是标量将给出出错信息。

Y=ones(m,n)：生成m n的全1阵。

Y=ones(size(A))：生成与矩阵A大小相同的全1阵。

(3)eye生成单位阵

调用格式为：

Y=eye(n)：生成n n的单位阵。

Y=eye(m,n)：生成m n的矩阵，其对角线元素为1，其它元素为0。

Y=eye(size(A))：生成一个与矩阵A大小相同的单位阵。

3. 线性方程求解

当然，在线性代数中不止矩阵的运算，还有线性方程组运算等。在工程计算

中，一个很重要的问题是线性方程组的求解。在矩阵表示方法中，上述问题可以表述为：给定两个矩阵A和B，求X的唯一解使得：

AX=B或XA=B

MATLAB求解这种问题时并不计算矩阵的逆。尽管在标准数学中没有矩阵除法概念，在此MATLAB求解线性方程组时，采用前面介绍的除法运算“\”和“/”求解。

X=A\B：表示求矩阵方程AX=B的解。

X=B/A：表示求矩阵方程XA=B的解。

对于X=A\B，要求矩阵A和B有相同的行数，X和B有相同的列数，它的行数等于矩阵A的列数。对于X=B/A，则行和列的角色相互交换。

实际上，线性方程形式AX=B比形式XA=B出现的频率更高一些，相应地，左除“\”比右除“/”用得更多一些。并且由于(B/A)T=(AT\BT)，因此下面仅讨论“\”（左除）。

矩阵A并不要求是方阵，如果矩阵A的维数是m n，则有三种情况：

1 m=n：适定方程组，寻求精确解；

2 m>n：超定方程组，寻求最小二乘解；

3 m

针对不同的情况，左除算子采用不同的算法求解。

一般形式的线性方程组为：

Ax=b或AX=B

其中A是个方阵，b是一个列向量，B是个方阵，X是一个与A,B同样大小的方阵。它们的求解命令为：

x=A\b或X=A\B

如果A是奇异的，则AX=B的解或者不存在，或者存在但不唯一。

当A接近奇异时，A\B将给出警告信息，如果发现A是奇异的，一方面给出警告信息，另一方面给出结果为inf。

如果A是正定的，可进行Cholesky分解A=RT*R，从而可对线性方程组A*x=b 进行如下替换；

RT*Rx=b

由于左除算子可以处理三角矩阵，因此可以快速的解出：x=R\(RT\x)

如果A可进行LU分解A=L*U，则线性方程组A*x=b，可如下进行快速计算：x=U\((L\b)

正式凭借MATLAB的这些突出的优势，它现在已成为世界上应用最广泛的工程计算软件。在大学里MATLAB是一种必须掌握的基本工具，而在研究设计单位，更是研究和解决计算问题的一种标准软件。

二、MATLAB在经济学中的应用

经济学在现代社会中起着重要的作用，在分析经济学的历程中，数学上都是用线性代数对其进行建模。但是随着经济发展，人们的经济行为越来越丰富，单靠线形分析有其局限性。而众多经济分析软件的问世，为经济分析提供了很好的工具，因此通过MATLAB来对经济续模型进行研究，对现代经济学分析有着重要意义。

价格平衡模型分析

在 Leontiff 成为诺贝尔奖金获得者的历史中，线性代数曾起过重要的作

用，如今来看他的基本思路；假定一个国家或区域经济可以分解为n 个部门，这些部门都有生产产品或服务的独立功能。设单列n 元向量vx 是这个n 个部门的产出，组成在Rn 空间的产出向量。

先假定该社会是自给自足的经济，这个是一个最简单的情况。因此各经济的本名生产出的产品，完全被自己部门和其他部门所消费。Leontiff 提出的问题是，各生产部门的实际产出的价格P 应该是多少，才能使各部门的收入和消耗相等，以维持持续的生产。

举一个最简单的例子，假如一个自给自足的经济体由三个部门组成，它们是煤炭业、电力业和钢铁业。它们的单位消耗列向量和销售收入列向量P 如下表：

这就是说，电力业产出了个点位的产品，有40 个单位会被煤炭业消耗，

10 单位被自己消耗，而被钢铁业消耗的是个单位，各行业付出的费用为：：

这就是内部消耗的计算方法，把几个部门都算上，可以写出消耗成本

其中

总的价格平衡方程可以写成为：

此等式右端常数项为零，使一个齐次方程。它有非零解的条件是系数行列式等于零，可以用行阶阶梯简化来求解。

用MATLAB 语句写出其解得表达式：

V=[0,0.4,0.6;0.6,0.1,0.2;0.4,0.5,0.2], H=rref([eye(3)-V,zeros(3,1)]) 其中eye 为生成单位矩阵，计算结果如下：

这个结果是合理的，简化行阶梯形式只有两行，说明[I-V]的秩是2，所以它的行列式必定为零。由于现在有三个变量可以作为自由变量。记住H 矩阵中各列的意义，它们分别是原方程中pe,pc,ps, 的系数，所以简化行阶梯矩阵H 表示的是下列方程：

这里取钢铁业价格PS 为自由变量，所以煤炭业和电力业的价格应该为钢铁业价格的0.94 和0.85 倍。如果钢铁业产品价格总计为100 万元，则煤炭业的产品价格总计为94 万，电力业的价格总计为85 万。

采用vm 文件，m 文件名为price pilo vexjuan.m,程序代码为：

V=[0,0.4 ,0.6;0.6,0.1,0.2;0.4,0.5,0.2]; H=rref([eye(3)-V,zeros(3,1)]); Ps=100; Pc=H(1,3)*Ps; Pe=H(2,3)*ps

三、MATLAB在三维图形绘制中的应用

为了显示三维图形，MATLAB 提供了各种各样的函数。有一些函数可在三维空

间中画线，而另一些可以画曲面与线格框架。另外，颜色可以用来代表第四维。

1.函数PLOT3命令

plot3 命令将绘制二维图形的函数plot 的特性扩展到三维空间。函数格式除了包括第三维的信息（比如Z 方向）之外，与二维函数plot 相同。plot3 一般语法调用格式是plot3(x1,y1,z1,S1,x2,y2,z2,S2,…)，这里xn,yn 和zn 是向量或矩阵，Sn 是可选的字符串，用来指定颜色、标记符号和/或线形。

总的来说，plot3 可用来画一个单变量的三维函数。如下为一个三维螺旋线例子：

? t=0:pi/50:10*pi;

? plot3(sin(t),cos(t),t)

?title( ‘Helix ‘),xlabel( ‘sint(t) ‘),ylabel( ‘cos(t) ‘ ),zlabel( ‘ t ‘ )

? text(0,0,0, ‘ Origin ‘ )

? grid

? v = axis

v =

-1 1 -1 1 0 40

输出见图

从上例可明显看出，二维图形的所有基本特性在三维中仍都存在。axis 命令扩展到三维只是返回Z 轴界限（0 和40），在数轴向量中增加两个元素。函数zlabel 用来指定z 轴的数据名称，函数grid 在图底绘制三维网格。函数test(x,y,z, ‘ string ‘ )在由三维坐标x,y,z 所指定的位置放一个字符串。另外，子图和多图形窗口可以直接应用到三维图形中。

2.如何改变视角

注意两个图形，一个是以30 度视角向下看z=0 平面,一个是以37.5 度视角向上看x=0 平面。这是对所有三维图形的缺省视角。与z=0 平面所成的方向角叫仰角，与x=0 平面的夹角叫做方位角。这样，缺省的三维视角方向仰角为30 度，方位角为-37.5 度。而缺省的二维视角仰角为90 度，方位角为0 度。仰角和方位角的概念在图中形象地画出。

在MATLAB 中，函数view 改变所有类型的二维和三维图形的图形视角。view(az,el)和view([az,el)]将视角改变到所指定的方位角az 和仰角el。

经过查阅书籍，得到View函数的用法如下图：

四、心得体会

我在学习的过程中也做了诸如心得体会的记录，希望能记录到最后作为总结。

1.如果不是计算机转业，只是为了方便自己的工作或学习，那么没有必要把matlab教程全部学会，只需要学需要的那部分即可。根据你个人的需要而定，但是基本命令、数据类型、基本的程序结构、文件的IO是必须看的，因为任何一个程序都需要这几个基本的块。

2．需要找一本matlab的函数工具词典，就像汉语词典一样，你要尽量多的熟悉matlab自带的函数，及其作用，因为matlab的自带函数特别多，基本上能够满足一般的数据和矩阵的计算，所以基本上不用你自己编函数（如我们大一学的C语言还有大三的vb中，大部分的函数都需要自己编）。这一点对程序非常有帮助，可以使程序简单，运行效率高。

3.把基本的知识看过之后，就需要找一个实际的程序来动手编一下，不要等所有的知识都学好之后再去编程，要在编程的过程中学习，程序需要什么知识再去补充，最好做一些随手笔记什么的。

作为和Mathematica、Maple并列的三大数学软件。其强项就是其强大的矩阵计算以及仿真能力。要知道Matlab的由来就是Matrix + Laboratory = Matlab，所以这个软件在国内也被称作矩阵实验室。Matlab提供了自己的编译器：全面兼容C++以及 Fortran两大语言。所以Matlab是工程师，科研工作者手上最好的语言，最好的工具和环境。Matlab 已经成为广大科研人员的最值得信赖的助手和朋友！我们更有理由去学习掌握这门技术！

参考资料：

《MATLAB教程R2011a》北京航空航天大学出版社

《经典教程——matlab三维图形的绘制》

《MATLAB 基础与应用简明教程》张平

《线形代数实践及MATLAB入门》陈怀琛,龚杰民

《MATLAB在经济评价中的应用》袁以美,叶合欣

@KungfuLeon

matlab 作业

实验一 1、熟悉MATLAB的窗口结构（命令窗口、历史命令窗口、工作区窗口、当前目录 窗口） 2、掌握命令窗口中基本命令的使用 3、在命令窗口中，给定圆的半径r，求得圆的周长c和面积s，并查看工作区窗 口的变化 4、将r,c,s变量保存到磁盘文件abc.mat中，并删除内存变量r,c,s，查看工作 区窗口的变化 5、将abc.mat文件中变量装入内存，查看工作区窗口的变化 6、将历史命令窗口中的命令再装入命令窗口中使用 7、改变当前目录，查看当前目录窗口的变化 8、掌握命令窗口中 cd,quit,help,date,dir,ls,what,who,clocl,fix(clock),format,save,loa d,clc,clear等命令的使用 9、注意各种MATLAB版本的差别 实验二 1、在命令窗口中，输入长方形的长和宽，求长方形的周长和面积 2、输入三角形的三条边（要满足构成三角形的条件），求三角形的周长和面积 3、掌握MATLAB中各标准函数的使用（sin,cos,sind,fix,mod,…） 4、用fprintf输出各种类型的数据（如fprintf('a=%d\n',123) a=123 >> fprintf('b=%f\n',123.456) b=123.456000 >> fprintf('c=%c\n','A') c=A……) 实验三 1、在编辑窗口中：输入学生成绩，输出该成绩的等级。等级规定如下：[90， 100]为A等，[80，90)为B等，[70，80)为C等，[60，70)为D等，[0，60)为E等。要求用if和 switch两种方法实现。 2、商场购物，100件以下，不优惠，100~199件95折，200~399件90折，400~799 件85折，800~1499件80折，1500件以上，75折。输入所购货物的单价、件数，求实际付款数目。要求用if和 switch两种方法实现（在编辑窗口中实现）。 实验四 1、求两个正整数的最大公约数和最小公倍数（在编辑窗口中实现，命令窗口中 调用）。 2、求100~300内所有素数（在编辑窗口中实现，命令窗口中调用）。

matlab期末大作业

电气学科大类 Modern Control Systems Analysis and Design Using Matlab and Simulink Title: Automobile Velocity Control Name: 巫宇智 Student ID: U200811997 Class：电气0811

电气0811 巫宇智 Catalogue Preface (3) The Design Introduction (4) Relative Knowledge (5) Design and Analyze (6) Compare and Conclusion (19) After design (20) Appendix (22) Reference (22)

Automobile Velocity Control 1.Preface: With the high pace of human civilization development, the car has been a common tools for people. However, some problems also arise in such tendency. Among many problems, the velocity control seems to a significant challenge. In a automated highway system, using the velocity control system to maintain the speed of the car can effectively reduce the potential danger of driving a car and also will bring much convenience to drivers. This article aims at the discussion about velocity control system and the compensator to ameliorate the preference of the plant, thus meets the complicated demands from people. The discussion is based on the simulation of MATLAB. Key word: PI controller, root locus

系统仿真结课作业

系统仿真导结课作业 一、概述 建模与仿真技术已成功地应用于航空航天、生产制造、交通运输、信息、生物、医学、材料、能源、教育、军事、社会、经济等众多领域;并成功地应用于产品研制的全生命周期,包括需求分析、方案论证、概念设计、初步设计、详细设计、生产制造、试验试飞、运行、维护、训练等各个阶段。仿真科学与技术正是从其广泛的应用中获得了日益强大的生命力，而仿真技术的发展反过来使得其得到愈来愈广泛的应用。 广义而言,仿真是采用建模的方法和物理的方法对真实环境客观事物进行 抽象、映射、描述和复现。基于系统原理、理论、定律、系统数据等应用计算机技术、软件技术和信息技术建立仿真环境(虚拟环境) ,在仿真环境中对客观事物进行研究。客观事物包括真实环境中的实体/系统、自然环境(地形、大气、海洋、空间)、和人的行为(操作、决策、推理)。仿真环境包括模型、数据、软件、物理效应设备、计算机等。 计算机仿真的三要素是系统、模型、计算机，三个基本活动是模型设计、模型执行、模型分析。计算机仿真的三要素和三个基本活动的关系关系如图1 所示。 系统 模型设计模型分析 模型执行 模型计算机 图1 从模型设计到模型分析经历的过程，即对实物进行仿真可概括为以下几个方面： 1) 问题的描述; 2) 建立概念模型; 3) 建立仿真模型; 4) 收集数据; 5) 编写程序; 6) 在计算机进行模型试验; 7) 模型和数据的验证; 8) 仿真结果显示; 9) 仿真结果分析和评估。 仿真是建立模型在计算机上运行,但这属于数学仿真,随着技术的发展,许多应用 领域建立仿真系统时除了模型外还要求将实物和人员包含在仿真回路中。 由此可见，基于仿真设计与传统设计的方法和流程两者有很大区别(图2) ,基于仿真设计可以在计算机上建立虚拟样机,对产品的外形、结构、强度、动力

数学应用软件作业2Matlab作图

注意：上机作业文件夹以自己的姓名学号命名，文件夹中包括如下上机报告和Matlab程序。

5、用surf ，mesh 绘制曲面22 2z x y =+，]3,3[],3,3[-∈-∈y x 。 6、用polar 绘制阿基米德螺线r a θ=和三叶玫瑰线cos3r a θ=。（a=100） 7、在同一平面的两个窗口中分别画出心形线和马鞍面。 三. 上机方法与步骤 给出相应的问题分析及求解方法，并写出Matlab 程序。并有上机程序显示。 第1题：要在同一坐标系中分别画出四个函数的图形，可以利用plot 来将这四条曲线画在一起。 Matlab 程序： x=linspace(15,200,50); y1=(1+1./x).^x; y2=(1+1./x).^(x+1); y3=[1+1./(x+1)].^x; y4=2.7183; plot(x,y1,'r',x,y2,'g',x,y3,'b',x,y4,'co') 第2题：直接用ezplot 命令绘制函数的图形。 Matlab 程序： ezplot('（exp(x*y)-sin(x+y)）',[-3,3])

第3题：取a=1，直接用ezplot命令绘出这两条曲线。 Matlab程序： ezplot('[t-sin(t)]','[1-cos(t)]',[0,2*pi]) 第4题：利用subplot命令将这四个函数的图像画在两行两列的同一个图形上。Matlab程序： x=[0:0.1:5]; subplot(2,2,1); plot(x,sin(5*x)); grid on; title('plot-y=sin(5*x)'); subplot(2,2,2); plot(x,cos(3*x));

春MATLAB仿真期末大作业

MATLAB仿真 期末大作业 姓名：班级：学号：指导教师：

2012春期末大作业 题目：设单位负反馈控制系统前向通道传递函数由)()(21s G s G 和串联，其中： ) 1(1)()(21++==s A s G s K s G A 表示自己学号最后一位数（可以是零），K 为开环增益。要求： （1）设K=1时，建立控制系统模型，并绘制阶跃响应曲线（用红色虚线，并标注坐标和标题）；求取时域性能指标，包括上升时间、超调量、调节时间、峰值时间； （2）在第（1）问中，如果是在命令窗口绘制阶跃响应曲线，用in1或者from workspace 模块将命令窗口的阶跃响应数据导入Simulink 模型窗口，用示波器显示阶跃响应曲线；如果是在Simulink 模型窗口绘制阶跃响应曲线，用out1或者to workspace 模块将Simulink 模型窗口的阶跃响应数据导入命令窗口并绘制阶跃响应曲线。 （3）用编程法或者rltool 法设计串联超前校正网络，要求系统在单位斜坡输入信号作用时，速度误差系数小于等于0.1rad ，开环系统截止频率s rad c /4.4''≥ω，相角裕度大于等于45度，幅值裕度大于等于10dB 。

仿真结果及分析： （1）、（2）、将Simulink模型窗口的阶跃响应数据导入命令窗口并绘制阶跃响应曲线 通过在Matlab中输入命令： >> plot(tout,yout,'r*-') >> title('阶跃响应曲线') 即可得出系统阶跃响应曲线，如下： 求取该控制系统的常用性能指标：超调量、上升时间、调节时间、峰值时间的程序如下： G=zpk([],[0,-1],5)。 S=feedback(G,1)。

科技写作结课作业(时域有限差分法的Matlab仿真开题报告)

开题报告 论文题目：基于matlab的时域有限差分法的电磁仿真研究（10分） 学院：电气工程及其自动化学院学号：1103000105姓名：__杨志刚___ 一、论文选题的目的和意义（300字以内；15分） 时域有限差分法，因具有多种优点被运用到电磁场理论研究的各个方面，而且其使用成效和应用领域还在迅速扩大和提高，在现代电磁场理论研究中具有很大的重要性和很强的可操作性。但是同时这种方法也存在一定的缺陷，主要表现在对无边界问题需要吸收边界条件处理，有色散误差，消耗内存大等方面。本课题在利用时域有限差分法对一些实际的算例进行实验仿真和验证，同时对这种方法在解决实际问题的缺陷进行一定程度的研究和分析。 Matlab作为一种工程仿真工具得到了广泛应用。用于时域有限差分法，可以简化编程，使研究者的研究重心放在FDTD法本身上，而不必在编程上花费过多的时间。 二、国内外关于该论题的研究现状和发展趋势（500字以内；15分） 时域有限差分方法作为一种典型的全波时域分析方法，因其原理直观、编程简便、实用性强在目前的计算电磁学领域内被人们广泛深入地研究，并取得巨大应用成功的方法。时域数值技术的一个突出优点是可以给出关于问题空间的丰富的时域信息，而且经过简单的时频变换，即可得到宽带范围的频域信息，相对频域方法显著地节约了计算量。最近几十年，是电磁场数值计算时域技术蓬勃发展的时期，各具优势和特色的新颖时域算法层出不穷。 但是到目前为止国内关于时域有限差分法中的PML 算法文献较少，其中绝大多数文献集中在综述和应用方面。而在国际的学报和杂志上对于这方面的文献非常多。时域有限差分法经过了三十年多年的高速发展之后，仍然还是计算电磁学制高点的研究热潮，而且其应用的范围和成效还在迅速的扩大和提高。本课题正是利用时域有限差分法的基础理论，利用matlab对一些实际的电磁场问题进行仿真研究。 三、论文的主攻方向、主要内容、研究方法及技术路线（1000字左右；40分） 通过对时域有限差分法理解基础之上，利用matlab仿真软件按照这种方法编程，实现对三种情况下的电磁场情况的仿真研究。

matlab教程课后作业

【例1.3-5】图示复数i z i z 21,3421+=+=的和。 z1=4+3*i;z2=1+2i; z12=z1+z2 clf,hold on plot([0,z1,z12],'-b','LineWidth',3) plot([0,z12],'-r','LineWidth',3) plot([z1,z12],'ob','MarkerSize',8) hold off,grid on axis equal axis ([0,6,0,6]) text(3.5,2.3,'z1') text(5,4.5,'z2') text(2.5,3.5,'z12') xlabel('real') ylabel('image') shg z12 = a=-8; r_a=a^(1/3) p=[1,0,0,-a]; R=roots(p) MR=abs(R(1)); t=0:pi/20:2*pi; x=MR*sin(t); y=MR*cos(t); plot(x,y,'b:'),grid on hold on plot(R(2),'.','MarkerSize',30,'Color','r') plot(R([1,3]),'o','MarkerSize',15,'Color','b') axis([-3,3,-3,3]),axis square

hold off r_a = 1.0000 + 1.7321i R = -2.0000 1.0000 + 1.7321i 【例1.3-10】画出衰减振荡曲线t e y t 3sin 3-=，t 的取值范围是]4,0[π。 t=0:pi/50:4*pi; y=exp(-t/3).*sin(3*t); plot(t,y,'r','LineWidth',2) axis([0,4*pi,-1,1]) xlabel('t'),ylabel('y')

MATLAB大作业

选 题 说 明 本人选做第2、4、5、9、11、12、13、14、16、19、24 题。 作业内容题目2：问题描述：在[0 ， 2π]范围内绘制二维曲线图y=cos(5x)*sin(x) （1）问题分析 这是一个二维绘图问题，先写出x的取值范围，再用plot函数画出y的图像。 （2）软件说明及源代码 >> x = 0:pi/100:2.*pi; y=cos(5*x).*sin(2*x); >> plot(x,y) （3）实验结果 题目4：问题描述：创建符号函数并求解，要求写出步骤和运行结果 (1)创建符号函数f=ax2+bx+c

(2)求f=0的解 （1）问题分析 这是符号计算问题，首先要确定符号变量，然后创建符号函数，最后利用subs函数求解特值。 （2）软件说明及源代码 >> syms a b c x f; f=a*x^2+b*x+c; subs(f,0) （3）实验结果 ans = c 题目5：问题描述：求积分 （1）问题分析 这是符号计算的积分求解问题，首先需要确定符号变量，然后利用int函数计算积分。 （2）软件说明及源代码 >> syms x y; y=sqrt(1-2*sin(2*x)); >> int(y,x,0,pi/2) （3）实验结果 ans = ellipticE(-pi/4, 4)*1i - ellipticE(pi/4, 4)*1i - ellipticE(-pi/6, 4)*2i + ellipticE(pi/6, 4)*2i 题目9：问题描述：按水平和竖直方向分别合并下述两个矩阵：

（1）问题分析 这是考查矩阵的基本操作，首先定义矩阵，然后合并矩阵。 （2）软件说明及源代码 >> A=[1,0,0;1,1,0;0,0,1]; B=[2,3,4;5,6,7;8,9,10]; >> a=[A,B],b=[A;B] （3）实验结果 a = 1 0 0 2 3 4 1 1 0 5 6 7 0 0 1 8 9 10 b = 1 0 0 1 1 0 0 0 1 2 3 4 5 6 7 8 9 10 题目11：问题描述：计算z＝yx2+3y2x+2y3的和： （1）问题分析 这是符号计算问题，首先确定符号变量，然后构造函数，最后利用diff函数进行求导。 （2）软件说明及源代码 >> syms x y z; >> z=y*x^2+3*y^2*x+2*y^3; >> diff(z,y,1),diff(diff(z,y,1),x,1) （3）实验结果 ans = x^2 + 6*x*y + 6*y^2

MATLAB结课作业

4.10 上机操作步骤 1在MatLab 的命令窗口输入如下命令序列: clf subplot(1,2,1) hold on grid on n=1:1000; m=1./n.*cos(n*pi/2); plot(n,m,'k.') 观察数列的散点图22,当n 趋于无穷大时,数列趋于 0 subplot(1,2,2) hold on grid on n=500:10000; m=1./n.*cos(n*pi/2); plot(n,m,'k.') fplot('0.001',[500,10000]) fplot('-0.001',[500,10000]) axis([500,10000,-0.005,0.005]) 观察图23,当001.0=ε时,可以取N= 1000 ,当n>N 时有επε<< -2 co s n 1n . 图22 图23 2 在MatLab 的命令窗口输入如下命令序列: clf subplot(1,2,1) hold on grid on fplot('x.*x',[1,3])

观察函数图24, 当2x →时,2x y =的极限是 4 subplot(1,2,2) hold on grid on fplot('x.*x',[1.9,2.1]) fplot('4.001',[ 1.9,2.1]) fplot('3.999',[ 1.9,2.1]) axis([1.9997,2.0005,3.9989,4.0011]) % 调整显示图形的范围是该实验的重点 观察图25,当001.0=ε时, δ取 0.003 δ<-<2 0x 时,001.04<-y ? 图24 图25 3 在MatLab 的命令窗口输入: syms x limit((2.^x-log(2.^x)-1)./(1-cos(x)),x,0) 运行结果为 ans = log(2)^2 理论上用洛必达法则计算该极限: x x x cos 112ln 2lim 0x ---→= 1 4 在MatLab 的命令窗口输入如下命令序列: （1）syms x y=sqrt(x+2)*(3-x)^4/(x+1)^5 diff(y,x) ％求一阶导数 运行结果 ＝y'1/2/(x+2)^(1/2)*(3-x)^4/(x+1)^5-4*(x+2)^(1/2)*(3-x)^3/(x+1)^5-5*(x+2)^(1/2)*(3-x )^4/(x+1)^6 x=1; eval(y) ％求导数在x ＝1处的值 运行结果 1'=x y = 0.8660

期末大作业报告

期末大作业报告 课程名称：数字图像处理 设计题目：车牌识别 学院：信息工程与自动化学院 专业：计算机科学与技术 年级：xxxxx 学生姓名：xxxxxxx（学号xxxxxxxxxxxxx） 指导教师：xxxx 日期：20XX.6.10 教务处制 车牌识别 摘要:数字图像处理技术是20世纪60年代发展起来的一门新兴学科，随着图像处理理论和方法的进一步完善，使得数字图像处理技术在各个领域得到了广泛应用，并显示出广阔的应用前景。MATLAB既是一种直观、高效的计算机语言,同时又是一个科学计算平台。它为数据分析和数据可视化、算法和应用程序开发提供了最核心的数学和高级图形工具。根据它提供的500多个数学和工程函数，工程技术人员和科学工作者可以在它的集成环境中交互或编程以完成各自的计算。MATLAB中集成了功能强大的图像处理工具箱。由于MA TLAB语言的语法特征与C语言极为相似，而且更加简单，更加符合科技人员对数学表达式的书写格式，而且这种语言可移植性好、可扩展性强，再加上其中有丰富的图像处理函数，所以MA TLAB在图像处理的应用中具有很大的优势。车牌识别技术是智能交通系统的重要组成部分，在近年来得到了很大的发展。本文从预处理、边缘检测、车牌定位、字符分割、字符识别五个方面，具体介绍了车牌自动识别的原理。并用MATLAB软件编程来实现每一个部分,最后识别出汽车牌照。 关键词:车牌识别、数字图像处理、MATLAB

一、设计原理 车辆牌照识别系统的基本工作原理为：将摄像头拍摄到的包含车辆牌照的图像通过视频卡输入到计算机中进行预处理，再由检索模块对牌照进行搜索、检测、定位，并分割出包含牌照字符的矩形区域，然后对牌照字符进行二值化并将其分割为单个字符，然后输入JPEG或BMP 格式的数字，输出则为车牌号码的数字。牌照自动识别是一项利用车辆的动态视频或静态图像进行牌照号码、牌照颜色自动识别的模式识别技术。其硬件基础一般包括触发设备、摄像设备、照明设备、图像采集设备、识别车牌号码的处理机等，其软件核心包括车牌定位算法、车牌字符分割算法和光学字符识别算法等。某些牌照识别系统还具有通过视频图像判断车辆驶入视野的功能称之为视频车辆检测。一个完整的牌照识别系统应包括车辆检测、图像采集、牌照识别等几部分。当车辆检测部分检测到车辆到达时触发图像采集单元，采集当前的视频图像。牌照识别单元对图像进行处理，定位出牌照位置，再将牌照中的字符分割出来进行识别，然后组成牌照号码输出。 二、设计步骤 1. 提出总体设计方案: （1）车牌图像预处理方法 因为车牌图像都是在室外拍摄的，所以不可避免地会受到光照、气候等因素的影响，而且拍摄者的手部抖动与车辆的移动会造成图像的模糊。要去除这些干扰就得先对车牌图像进行预处理。由于当前数码相机的像素较高，原始图像的数据一般比较大，输入的彩色图像包含大量颜色信息，会占用较多的存储空间，且处理时也会降低系统的执行速度。因此对图像进行识别等处理时，常将彩色图像转换为灰度图像，以加快处理速度。对图像进行灰度化处理后常用的方法是图像二值化、去除背景图像、增强处理、边缘检测、滤波等处理等。

有限元大作业matlab课程设计例子

有 限 元 大 作 业 程 序 设 计 学校：天津大学 院系：建筑工程与力学学院 专业：01级工程力学 姓名：刘秀 学号：\\\\\\\\\\\ 指导老师： 连续体平面问题的有限元程序分析 [题目]： 如图所示的正方形薄板四周受均匀载荷的作用，该结构在边界 上受正向分布压力， m kN p 1=，同时在沿对角线y 轴上受一对集中压 力，载荷为2KN ，若取板厚1=t ，泊松比0=v 。 [分析过程]： 由于连续平板的对称性， 只需要取其在第一象限的四分之一部分

参加分析，然后人为作出一些辅助线将平板“分割”成若干部分，再为每个部分选择分析单元。采用将此模型化分为4个全等的直角三角型单元。利用其对称性，四分之一部分的边界约束，载荷可等效如图所示。 [ 用和单元信息文件DATA.OUT。 位移模式：用用线性位移模式 载荷类型：节点载荷，非节点载荷应先换算为等效节点载荷 材料性质：弹性体由单一的均匀材料组成 约束方式：为“0”位移固定约束，为保证无刚体位移，弹性体至少应有对三个自由度的独立约束 方程求解：针对半带宽刚度方程的Gauss消元法 输入文件：由手工生成节点信息文件NODE.IN，和单元信息文件ELEMENT.IN 结果文件：输出一般的结果文件DATA.OUT 程序的原理如框图：

（1） ID ： ID=2时为平面应变问题 （平面问题） ，LJK_ELE(I,1)，LJK_ELE(I,2)， X(I),Y(I)分别存放节点I 的x ，y 表示第I 个作用有节点载荷的节点x,y 方向的节点载荷数值 存放节点载荷向量，解方程后该矩 （2 READ_IN ： 读入数据 BAND_K ： 形成半带宽的整体刚度矩阵 FORM_KE ： 计算单元刚度矩阵 FORM_P ： 计算节点载荷 CAL_AREA ：计算单元面积 DO_BC ： 处理边界条件 CLA_DD ： 计算单元弹性矩阵 SOLVE ： 计算节点位移 CLA_BB ： 计算单元位移……应变关系矩阵 CAL_STS ：计算单元和节点应力 （3）文件管理： 源程序文件： chengxu.for 程序需读入的数据文件：

《科学计算与MATLAB》期末大作业

杭州电子科技大学信息工程学院《科学计算与MATLAB》期末大作业

给出程序、图、作业分析，程序需加注释。 1. 试编写名为fun.m 的MATLAB 函数，用以计算下述的值： ?? ? ??-<->=t t n t t t n t f 的)4/sin()(si 对所有)4/sin(其他情况)sin(的)4/sin()(si 对所有)4/sin()(ππππ 绘制t 关于函数f(t)的图形，其中t 的取值范围为ππ66≤≤-t ，间距为10/π。 function y=fun()%定义函数 % t=-6*pi:pi/10:6*pi; %定义变量范围 y = (sin(pi/4)).*(sin(t)>sin(pi/4))+(sin(-pi/4)).*(sin(t)=sin(-pi/4)));%函数表示 plot(t,y); %画图 end

2.解以下线性方程组 ??? ??=+=++=--3 530 42231 321321x x x x x x x x A=[2 -1 -1;1 1 4;3 0 5];%输入矩阵 B=[2;0;3]; %输入矩阵 X = A\B %计算结果 3.已知矩阵? ? ??? ???? ???=44434241 3433323124232221 14131211A 求： (1)A(2:3,2:3) (2)A(:,1:2) (3)A(2:3,[1,3]) (4)[A,[ones(2,2)；eye(2)]]

A=[11 12 13 14;21 22 23 24;31 32 33 34;41 42 43 44];%输入矩阵A(2:3,2:3) %输出矩阵 A(:,1:2) %输出矩阵 A(2:3,[1,3]) %输出矩阵 [A,[ones(2,2);eye(2)]] %输出矩阵

利用Matlab实现Romberg数值积分算法----系统建模与仿真结课作业

利用Matlab 实现Romberg 数值积分算法 一、内容摘要 针对于某些多项式积分，利用Newton —Leibniz 积分公式求解时有困难，可以采用数值积分的方法，求解指定精度的近似解，本文利用Matlab 中的.m 文件编写了复化梯形公式与Romberg 的数值积分算法的程序，求解多项式的数值积分，比较两者的收敛速度。 二、数值积分公式 1.复化梯形公式求解数值积分的基础是将区间一等分时的Newton —Cotes 求积公式: I =(x)[f(a)f(b)]2 b a b a f dx -≈ +? 其几何意义是，利用区间端点的函数值、与端点构成的梯形面积来近似(x)f 在区间[a,b]上的积分值，截断误差为： 3" (b a)()12 f η-- (a,b)η∈ 具有一次的代数精度，很明显，这样的近似求解精度很难满足计算的要求，因而，可以采用将积分区间不停地对分，当区间足够小的时候，利用梯形公式求解每一个小区间的积分近似值，然后将所有的区间加起来，作为被求函数的积分，可以根据计算精度的要求，划分对分的区间个数，得到复化梯形公式： I =1 1 (b a)(b a) (x)dx [f(a)f(b)2(a )]2n b a k k f f n n -=--≈+++∑? 其截断误差为：

2" (b a)h ()12 R f η--= (a,b)η∈ 2.Romberg 数值积分算法 使用复化的梯形公式计算的数值积分，其收敛速度比减慢，为此，采用Romberg 数值积分。其思想主要是，根据I 的近似值2n T 加上I 与2n T 的近似误差，作为新的I 的近视，反复迭代，求出满足计算精度的近似解。 用2n T 近似I 所产生的误差可用下式进行估算： 12221 ()3 n n n I T T T -?=-=- 新的I 的近似值： 122 n n j T T -=?+ j =（0 1 2 ….） Romberg 数值积分算法计算顺序 i=0 (1) 002T i=1 (2) 102T (3) 012T i=2 (4) 202T (5) 112T (6) 022T i=3 (7) 302T (8) 212T (9) 122T (10) 032T i=4 (11) 402T (12) 312T (13) 222T (14) 132T … … … … 其中，第一列是二阶收敛的，第二列是四阶收敛的，第三列是六阶收敛的，第四列是八阶收敛的，即Romberg 序列。

MATLAB期末大作业模板

MATLAB应用技术 期末大作业 专业： 姓名： 学号： 分数

一、在一个图形窗口中以子图形式同时绘制正弦、余弦、正切、余切曲线。请写下完整代码，展示图形结果。（请标注题图和坐标轴，用不同颜色和不同线型分别绘制以上曲线）。（15分） 二、某公司员工的工资计算方法如下。 （1）工作时数超过120小时者，超过部分加发15%。 （2）工作时数低于60小时者，扣发700元。 （3）其余按每小时84元发。 根据员工的工时数，计算应发工资。请写下完整的程序代码，并任意输入一工时数（使用input 函数），将结果展示（使用disp 函数）利用该代码进行计算工资，请写下计算结果。（15分） 三、编写一个函数文件，使其能够产生如下的分段函数： ?? ? ??≥<<≤-=66225.0,25.05.15.0)(x x x x x x f 请编写完整的函数文件（保存函数文件名为hanshu.m ），并编写脚本文件代码，任意输入x 值（使用input 函数），在脚本文件中调用函数文件求)(x f ，展示结果（使用disp 函数），请写下计算结果。（15分） 四、将5个学生的6门功课的成绩存入矩阵P 中，进行如下处理： （1）分别求每门课的最高分、最低分及相应学生的序号。 （2）分别求每门课的平均分和标准差。 （3）5门课总分的最高分、最低分及相应学生序号。 （4）将5门课总分按从大到小顺序存入score 中，相应学生序号存入num 。 请将各小题的运行代码完整写下来，并写下运行结果。（20分） 五、请利用所学的MATLAB 知识，自主设计一个图形用户界面，请完整记录它的设计过程，需提供文字、代码和图片，以充分说明设计的图形用户界面可实现

matlab结课论文

精选文库 山西大同大学matlab课程结课作业MATLAB程序应用 姓名： 课程序号： 2 班级： 学号： 2013年12月

1.实验内容：已知!123n n =?????L ，编写一个程序求满足100!10n ≤的 最大的n 值以及此时!n 的值。 function n n=2;m=1; while m<=10^100 m=m.*n;n=n+1; end m=m/(n-1)；n=n-2; m n m = 1.7112e+098 n =69 2.设)15113111191715131 1(22Λ+--++--+=π，试根据公式编出计算pi 的Mat lab 主程序文件，pi 的精度为0.00001。 程序： k=0;n=1;b=0;a=0; while abs((pi-a))>0.00001 a=2*sqrt(2)*k; k=( bcos( *pi/2)+sin(b*pi/2))/n+k; n=n+2; b=b+1; end a 输出a=3.141602572083633 ; a-pi= 9.918493839577991e-006 3.有两个矩阵A 和B 如下：????????????---=771175420132861-1A ，????????????------=0162310013125673B ， 将A 中所有等于-1的元素改为-2，将B 中所有小于0的元素改为1，然后将B 中等于0的元素的值改为A 的相应位置元素的值。请用Matlab 函数文件实现上述运算。 clear;

clc; A=[1 -1 6 8;2 3 -1 0;-2 4 5 7;1 -1 7 7]; B=[-3 -7 6 -5;-2 1 3 -1;0 0 1 3;2 6 -1 0]; C=A;A(A==-1)=-2;U=A; D=B;B(B<0)=1;V=B; A=C;B=D;[i,j]=find(B==0);A(i,j)=0;W=A; A=C;B=D; A,B,W,U,V %用函数文件实现矩阵中元素的变换。 %A、B为输入变量。 %U、V、W分别存放A、B中间变换结果。 ; 4.用matlab主程序文件产生动画：呈现一小圆（半径为1）在一大圆（半径为3）的圆周外部滚动的动画，要求连续滚动20周。 clea close;clc;r; axis([-6 6 -6 6],'equal','manual');hold on; ezplot('x^2+y^2-9'); h=ezplot('x^2+y^2-1'); x=get(h,'xdata'); y=get(h,'ydata'); for t=1:7200 set(h,'xdata',x+4*cosd(t),'ydata',y+4*sind(t)); drawnow; end

MATLAB大作业

M A T L A B大作业作业要求： （1）编写程序并上机实现，提交作业文档，包括打印稿（不含源程序）和电 子稿（包含源程序），以班为单位交，作业提交截止时间6月24日。 （2）作业文档内容：问题描述、问题求解算法（方案）、MATLAB程序、结果 分析、本课程学习体会、列出主要的参考文献。打印稿不要求MATLAB程序，但电 子稿要包含MATLAB程序。 （3）作业文档字数不限，但要求写实，写出自己的理解、收获和体会，有话 则长，无话则短。不要抄袭复制，可以参考网上、文献资料的内容，但要理解，要变成自己的语言，按自己的思路组织内容。 （4）从给出的问题中至少选择一题(多做不限，但必须独立完成，严禁抄袭)。 （5）大作业占过程考核的20%，从完成情况、工作量、作业文档方面评分。 第一类：绘制图形。（B级） 问题一：斐波那契（Fibonacci）螺旋线，也称黄金螺旋线（Golden spiral），是根据斐波那契数列画出来的螺旋曲线，自然界中存在许多斐波那契螺旋线的图案，是自然界最完美的经典黄金比例。斐波那契螺旋线，以斐波那契数为边的正方形拼成的长方形，然后在正方形里面画一个90度的扇形，连起来的弧线就是斐波那契 螺旋线，如图所示。 问题二：绘制谢尔宾斯基三角形（Sierpinskitriangle）是一种分形，由波兰数学家谢尔宾斯基在1915年提出,它是一种典型的自相似集。其生成过程为：取一个实心的三角形（通常使用等边三角形），沿三边中点的连线，将它分成四个小三

角形，然后去掉中间的那一个小三角形。接下来对其余三个小三角形重复上述操作，如图所示。 问题三：其他分形曲线或图形。分形曲线还有很多，教材介绍了科赫曲线，其他还有皮亚诺曲线、分形树、康托（G. Cantor）三分集、Julia集、曼德布罗集合（Mandelbrot set），等等。这方面的资料很多（如），请分析构图原理并用MATLAB 实现。 问题四：模拟掷骰子游戏：掷1000次骰子，统计骰子各个点出现的次数，将结果以下表的形式显示，并绘制出直方图。 点数 1 2 3 4 5 6 出现次数166 150 164 162 184 174 问题五：利用MATLAB软件绘制一朵鲜花，实现一定的仿真效果。 提示：二维/三维绘图，对花瓣、花蕊、叶片、花杆等的形状和颜色进行详细设置。 第二类：插值与拟合。（B级） 问题一：有人对汽车进行了一次实验，具体过程是，在行驶过程中先加速，然后再保持匀速行驶一段时间，接着再加速，然后再保持匀速，如此交替。注意，整个实验过程中从未减速。在一组时间点上测得汽车的速度如表所示。 （1）分别使用最近点插值、线性插值、三次埃尔米特插值和三次样条插值进行计算[0,110]时间段50个时间点的速度。 （2）绘制插值图形并标注样本点。

(完整版)有限元大作业matlab---课程设计例子

有限元大作业程序设计 学校：天津大学 院系：建筑工程与力学学院 专业：01级工程力学 姓名：刘秀 学号：\\\\\\\\\\\ 指导老师：

连续体平面问题的有限元程序分析 [题目]： 如图所示的正方形薄板四周受均匀载荷的作用，该结构在边界 上受正向分布压力， m kN p 1=，同时在沿对角线y 轴上受一对集中压 力，载荷为2KN ，若取板厚1=t ，泊松比0=v 。 [分析过程]： 由于连续平板的对称性，只需要取其在第一象限的四分之一部分参加分析，然后人为作出一些辅助线将平板“分割”成若干部分，再为每个部分选择分析单元。采用将此模型化分为4个全等的直角三角型单元。利用其对称性，四分之一部分的边界约束，载荷可等效如图所示。

[程序原理及实现]： 用FORTRAN程序的实现。由节点信息文件NODE.IN和单元信息文件ELEMENT.IN，经过计算分析后输出一个一般性的文件DATA.OUT。模型基本信息由文件为BASIC.IN生成。 该程序的特点如下： 问题类型：可用于计算弹性力学平面问题和平面应变问题 单元类型：采用常应变三角形单元 位移模式：用用线性位移模式 载荷类型：节点载荷，非节点载荷应先换算为等效节点载荷 材料性质：弹性体由单一的均匀材料组成 约束方式：为“0”位移固定约束，为保证无刚体位移，弹性体至少应有对三个自由度的独立约束 方程求解：针对半带宽刚度方程的Gauss消元法

输入文件：由手工生成节点信息文件NODE.IN，和单元信息文件ELEMENT.IN 结果文件：输出一般的结果文件DATA.OUT 程序的原理如框图：

matlab综合大作业(附详细答案)

m a t l a b综合大作业(附详细 答案) 标准化文件发布号：（9312-EUATWW-MWUB-WUNN-INNUL-DQQTY-

《MATLAB语言及应用》期末大作业报告1．数组的创建和访问（20分，每小题2分）： 1)利用randn函数生成均值为1，方差为4的5*5矩阵A；实验程序：A=1+sqrt(4)*randn(5) 实验结果： A = 2)将矩阵A按列拉长得到矩阵B； 实验程序：B=A(:) 实验结果： B =

3)提取矩阵A的第2行、第3行、第2列和第4列元素组成2*2的 矩阵C； 实验程序：C=[A(2,2),A(2,4);A(3,2),A(3,4)] 实验结果： C = 4)寻找矩阵A中大于0的元素；] 实验程序：G=A(find(A>0)) 实验结果： G = 5)求矩阵A的转置矩阵D； 实验程序：D=A' 实验结果： D = 6)对矩阵A进行上下对称交换后进行左右对称交换得到矩阵E； 实验程序：E=flipud(fliplr(A)) 实验结果： E =

7)删除矩阵A的第2列和第4列得到矩阵F； 实验程序：F=A; F(:,[2,4])=[] 实验结果： F = 8)求矩阵A的特征值和特征向量； 实验程序：[Av,Ad]=eig(A) 实验结果： 特征向量Av = + - + - - + + - 特征值Ad = 0 0 0 0 0 + 0 0 0 0 0 - 0 0 0 0 0 0 0 0 0 0 9)求矩阵A的每一列的和值； 实验程序：lieSUM=sum(A) 实验结果： lieSUM = 10)求矩阵A的每一列的平均值； 实验程序：average=mean(A) 实验结果： average = 2．符号计算（10分，每小题5分）： 1)求方程组20,0 uy vz w y z w ++=++=关于,y z的解； 实验程序：S = solve('u*y^2 + v*z+w=0', 'y+z+w=0','y,z'); y= S. y, z=S. z

matlab结课论文

山西大同大学matlab课程结课作业MATLAB程序应用 姓名： 课程序号： 2 班级： 学号： 2013年12月

1.实验内容：已知!123n n =????? ，编写一个程序求满足100!10n ≤的 最大的n 值以及此时!n 的值。 function n n=2;m=1; while m<=10^100 m=m.*n;n=n+1; end m=m/(n-1)；n=n-2; m n m = 1.7112e+098 n =69 2.设)15113111191715131 1(22 +--++--+=π，试根据公式编出计算pi 的Mat lab 主程序文件，pi 的精度为0.00001。 程序： k=0;n=1;b=0;a=0; while abs((pi-a))>0.00001 a=2*sqrt(2)*k; k=( bcos( *pi/2)+sin(b*pi/2))/n+k; n=n+2; b=b+1; end a 输出a=3.141602572083633 ; a-pi= 9.918493839577991e-006 3.有两个矩阵A 和B 如下：????????????---=771175420132861-1A ，????????????------=0162310013125673B ， 将A 中所有等于-1的元素改为-2，将B 中所有小于0的元素改为1，然后将B 中等于0的元素的值改为A 的相应位置元素的值。请用Matlab 函数文件实现上述运算。

clear; clc; A=[1 -1 6 8;2 3 -1 0;-2 4 5 7;1 -1 7 7]; B=[-3 -7 6 -5;-2 1 3 -1;0 0 1 3;2 6 -1 0]; C=A;A(A==-1)=-2;U=A; D=B;B(B<0)=1;V=B; A=C;B=D;[i,j]=find(B==0);A(i,j)=0;W=A; A=C;B=D; A,B,W,U,V %用函数文件实现矩阵中元素的变换。 %A、B为输入变量。 %U、V、W分别存放A、B中间变换结果。 ; 4.用matlab主程序文件产生动画：呈现一小圆（半径为1）在一大圆（半径为3）的圆周外部滚动的动画，要求连续滚动20周。 clea close;clc;r; axis([-6 6 -6 6],'equal','manual');hold on; ezplot('x^2+y^2-9'); h=ezplot('x^2+y^2-1'); x=get(h,'xdata'); y=get(h,'ydata'); for t=1:7200 set(h,'xdata',x+4*cosd(t),'ydata',y+4*sind(t)); drawnow; end