选修课MATLAB作业

Matlab程序设计

作业

学院：信息科学与技术学院

专业：通信工程

姓名：刘潇

学号：201313070404

班级：通信4班

教师：徐静

日期：2014年11月9日

文件1

%学号：201313070404

%姓名：刘潇

clear all %清除工作空间的所有变量

clf %清除图形窗口的内容

mfilename('fullpath') %返回当前正在运行的函数所在文件的文件名

echo on

%本m文件运行后会输出6个图形

subplot(2,3,1) %subplot（m,n,p）是将多个图画到一个平面上的工具。其中，m表示是图排成m行，n表示图排成n列，p表示图所在的位置。

t = 0:0.1:10;

z = impulse(1, [1 1 1], t);

stairs(t(1:5:end),z(1:5:end)) %stairs函数用于绘制阶梯状图，在图中为蓝色的那根线hold on %hold on 表示保留原来的图

plot(t,z,'r') %绘图函数，在图中为红色的那根连续的线

plot([0 t(end)], [0 0], 'k:') %绘图函数，在图中为水平虚线，左边对准了0刻度线

title('Impulse Response - (STAIRS)') %在图形上方添上标题

subplot(2,3,2) %此函数表示接下来的图形绘制在第二个位置

theta = 2*pi*(0:74)/75;

x = cos(theta);

y = sin(theta);

z = abs(fft(ones(10,1), 75))';

stem3(x, y, z) %stem3函数绘制离散序列数据的三维杆图

title('Polar FFT - (STEM3)') %在图形上方添上标题

subplot(2,3,3) %此函数表示接下来的图形绘制在第三个位置

[X,Y,Z] = peaks(-2:0.25:2);

[U,V] = gradient(Z, 0.25);

contour(X,Y,Z,10); %该函数用于绘制矩阵的等高线

hold on %hold on 表示保留原来的图

quiver(X,Y,U,V); %使用箭头绘制矢量图

title('Surface Gradient - (CONTOUR & QUIVER)') %在图形上方添上标题

theta = 0:0.1:4*pi;

[x,y] = pol2cart(theta(1:5:end), theta(1:5:end));

subplot(2,3,4) %此函数表示接下来的图形绘制在第四个位置

polar(theta,theta)

axis([-13 13 -12.5 14.5])

title('Spiral Plot - (POLAR)') %在图形上方添上标题

subplot(2,3,5) %此函数表示接下来的图形绘制在第五个位置

compass(x,y) %在一个圆形带网格的区域内，使用带箭头的线条绘制从原点出发的若干向量。绘制出来的效果图就像一个罗盘

axis([-13 13 -12.5 14.5])

title('Direction Vectors - (COMPASS)') %在图形上方添上标题

subplot(2,3,6) %此函数表示接下来的图形绘制在第六个位置

feather(x(1:19),y(1:19)) %feather复平面图形，把复数矩阵中的元素的相角和幅值显示成沿横轴等间隔辐射的箭头

axis([1 21 -5 10])

title('Direction Vectors - (FEATHER)') %在图形上方添上标题

set(gcf,'Position', [64 111 887 564])

echo off

shg

文件1运行后输出的图形为

文件2

%学号：201313070404

%姓名：刘潇

clear all %清除工作空间的所有变量

clf %清除图形窗口的内容

echo on

%本m文件运行后会输出6个图形

data = [10 2 3 5; 5 8 10 3; 9 7 6 1; 3 5 7 2; 4 7 5 3];

subplot(2,3,1) %subplot（m,n,p）是将多个图画到一个平面上的工具。其中，m表示是图排成m行，n表示图排成n列，p表示图所在的位置。

bar(data, 'stacked'); %绘制直方图

title('Bar Graph - (BAR, ''stacked'')'); %在图形上方添上标题

subplot(2,3,2) %此函数表示接下来的图形绘制在第二个位置

bar3h(data); %绘制三维水平直方图

title('Horizontal Bar Graph - (BAR3H, ''grouped'')'); %在图形上方添上标题

subplot(2,3,3) %此函数表示接下来的图形绘制在第三个位置

hist(randn(1000,3)); %绘制直方图函数

title('Histogram - (HIST)'); %在图形上方添上标题

subplot(2,3,4) %此函数表示接下来的图形绘制在第四个位置

area(data); %area是填充画图区域的命令

title('Area Plot - (AREA)'); %在图形上方添上标题

subplot(2,3,5) %此函数表示接下来的图形绘制在第五个位置

pie3(sum(data), [0 0 1 0]); %画3维饼图的函数

title(['3-D Pie Chart';' (PIE3) ']); %在图形上方添上标题

subplot(2,3,6) %此函数表示接下来的图形绘制在第六个位置

rose(5/3*randn(1000,1), 18); %使用rose函数可以绘制玫瑰图

title('Polar Histogram - (ROSE)'); %在图形上方添上标题

set(gcf,'Position',[184 248 740 424]) %设置对象的属性

echo off

shg

文件2运行后输出的图形为

文件3

%学号：201313070404

%姓名：刘潇

echo on

%本m文件运行后会输出6个图形

subplot(2,3,1) %subplot（m,n,p）是将多个图画到一个平面上的工具。其中，m表示是图排成m行，n表示图排成n列，p表示图所在的位置。

x = -3:0.3:3; y = x;

[X,Y]=meshgrid(x,y); % meshgrid（）生成绘制3-D图形所需的网格数据

[theat,R] = cart2pol(X,Y); %在MATLAB中，该函数用于把笛卡尔坐标转换为极坐标（柱坐标）

Z = sinc(R);

contourf(peaks(30), 10) %contourf（）填充二维等高线图。即先画出不同等高线，然后相邻的等高线之间用同一颜色进行填充。填充用的颜色决定于当前的色图颜色colorbar

grid on %打开网格线

title('Peaks Function - (CONTOURF & COLORBAR)') %在图形上方添上标题

subplot(2,3,2) %此函数表示接下来的图形绘制在第二个位置

plot3(X,Y,Z) %plot3(x,y,z)用来绘制3维曲线图

grid on %打开网格线

axis([-3 3 -3 3 -1 1])

title('Sinc Function - (PLOT3)') %在图形上方添上标题

subplot(2,3,3) %此函数表示接下来的图形绘制在第三个位置

waterfall(membrane(1)); %waterfall函数，绘制三维表面网格图，只是它的网格线在x轴方向出现，具有瀑布效果

title('L-shaped Membrane - (WATERFALL)') %在图形上方添上标题

subplot(2,3,4) %此函数表示接下来的图形绘制在第四个位置

contour3(peaks(30), 25); %三维空间等高线图，该命令生成一个定义在矩形格栅上曲面的三维等高线图

title('Peaks Function - (CONTOUR3)') %在图形上方添上标题

subplot(2,3,5) %此函数表示接下来的图形绘制在第五个位置

mesh(X,Y,Z) %生成由X，Y和Z指定的网线面

axis([-3 3 -3 3 -1 1])

title('Sinc Function - (MESH)') %在图形上方添上标题

subplot(2,3,6) %此函数表示接下来的图形绘制在第六个位置

surf(membrane(1)) %画三维曲面图

title('L-shaped Membrane - (SURF)') %在图形上方添上标题

set(gcf,'Position',[211 248 713 413]) %设置对象的属性

echo off

shg

文件3运行后输出的图形为

部分函数的具体用法

subplot函数

使用方法：subplot（m,n,p）或者subplot（m n p）。

subplot是将多个图画到一个平面上的工具。其中，m表示是图排成m行，n表示图排成n列，也就是整个figure中有n个图是排成一行的，一共m行，如果m=2就是表示2行图。p表示图所在的位置，p=1表示从左到右从上到下的第一个位置。

plot函数

plot（x）：x为长度为n的数值向量，坐标的纵坐标是向量x，横坐标为matlab系统根据x向量的元素序号自动生成的从1开始的向量。

plot（x，y）：纵坐标是y，横坐标是x。x，y向量的数目必须相等。

plot（x）：x为m*n的矩阵，plot(x)命令为矩阵的每一列画出一条线，共n条，曲线的横坐标为向量1：m，m是矩阵的行数。

plot(z): 可以复制复向量曲线它与（plot(real(z),imag(z))等效，以实部为横坐标，虚部为纵坐标。

plot(x1,y1,x2,y2...) ：绘制多条曲线。

stem3函数

函数绘制离散序列数据的三维杆图，常用格式为： stem3(z) ，stem3(x,y,z)。

axis函数

axis(v):v是一维向量，形式可以是[xmin xmax ymin ymax],定义二维图形x轴和y轴坐标轴的范围，其中v向量中的元素必须满足：xmin

polar函数

polar()进行极坐标系的绘图

polar（theta，rho）：以弧度为theta和半径为rho在极坐标中进行绘图，其中两个输入参数满足极坐标下rho=f （theta），其作用相当于直坐标中的plot（x,y）。

hist函数

hist（）绘制直方图

hist（y,n）将y的取值范围分成等差的n段，然后将所有元素分别分类到该n段中，根据每段元素个数绘制直方图的高度，n可以省略，默认值为10。

Hist（y,x）将y的取值分成等差的length（x）份，将y 中的元素放到各个分段中，然后以x中元素指定的位置为中心的直方图。

Area（x,y）函数

area（x,y）：若x和y是向量，则以x中的元素为横坐标，y中的元素为纵坐标，并且填充线条和x轴之间的空间，如

果y是矩阵，则绘制y中每一列的和。

bar(x,y)或barh(x,y)函数

在指定的横坐标x上画出y，参数x为单调增加的。如果y为矩阵，则将每个行向量画出。

contour()绘制等高线图，contourf（）绘制经过填充的等高线图

contour（z）：绘制矩阵z的等高线。

Contour（x,y,z）:在坐标（x,y）的值为矩阵z的等高线

Contour（z,n）：绘制n条等高线。

Contour（x,y,z,[v v]）:绘制高度为v的等高线。

compass函数：罗盘图

用法说明：

compass(x,y)函数绘制一个由原点出发由(x,y)组成的向量箭头图形；

compass(z)等价于compass（real(z),imag(z))；

compass(...,LineSpec)函数用参量LineSpec指定箭头的线型、标记符号、颜色等属性；

h = compass(...)函数返回line对象的句柄给h。

feather函数

feather复平面图形，把复数矩阵中的元素的相角和幅值显示成沿横轴等间隔辐射的箭头，格式feather(z),feather(x,y)等价于feather(x+y*i),feather(z,str),str是确定的线形绘制箭头。

rose函数

使用rose函数可以绘制玫瑰图，rose函数使用的基本方法是rose（t），其中t是随机数，绘制出来的图形为线性。

plot3函数

plot3（x1，y1，z1，选项1，x2，y2，z2，选项2，…）其中每一组x，y，z组成一组曲线的坐标参数，选项的定义和plot的选项一样。当x，y，z是同维向量时，则x，y，z对应元素构成一条三维曲线。当x，y，z是同维矩阵时，则以x，y，z对应列元素绘制三维曲线，曲线条数等于矩阵的列数。

mesh（x，y，z，c）函数

surf（x，y，z，c）函数

一般情况下，x，y，z是维数相同的矩阵，x，y是网格坐标矩阵，z是网格点上的高度矩阵，c用于指定在不同高度下的颜色范围。c省略时，Matlab认为c=z，也即颜色的设定是正比于图形的高度的。这样就可以得到层次分明的三维图形。当x，y省略时，把z矩阵的列下标当作x轴的坐标，把z矩阵的行下标当作y轴的坐标，然后绘制三维图形。当x，y是向量时，要求x的长度必须等于z矩阵的列，y的长度必须等于必须等于z的行，x，y向量元素的组合构成网格点的x，y坐标，z坐标则取自z矩阵，然后绘制三维曲线。

pie3函数

绘制三维饼图，常用格式为：

pie3（x）

x为向量，用x中的数据绘制一个三维饼图。

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

MATLAB结课论文设计.

MATLAB程序设计(论文) 基于MATLAB实现语音信号的去噪 院（系）名称电子与信息工程学院 专业班级通信工程 学号 学生姓名 任课教师

论文任务

摘要 滤波器设计在数字信号处理中占有极其重要的地位，FIR数字滤波器和IIR 滤波器是滤波器设计的重要组成部分。利用MATLAB信号处理工具箱可以快速有效地设计各种数字滤波器。课题基于MATLAB有噪音语音信号处理的设计与实现，综合运用数字信号处理的理论知识对加噪声语音信号进行时域、频域分析和滤波。通过理论推导得出相应结论，再利用MATLAB 作为编程工具进行计算机实现。在设计实现的过程中，使用窗函数法来设计FIR数字滤波器，用巴特沃斯、切比雪夫和双线性变法设计IIR数字滤波器，并利用MATLAB 作为辅助工具完成设计中的计算与图形的绘制。通过对对所设计滤波器的仿真和频率特性分析，可知利用MATLAB信号处理工具箱可以有效快捷地设计FIR和IIR数字滤波器，过程简单方便，结果的各项性能指标均达到指定要求。 关键词数字滤波器 MATLAB 窗函数法巴特沃斯切比雪夫双线性变换

目录 第1章绪论 (1) 1.1数字信号处理的意义 (1) 1.2语音去噪设计要求 (2) 第2章语音去噪方案设计 (3) 2.1语音去噪的应用意义 (3) 2.2 语音去噪设计框图 (3) 2.3设计原理 (4) 第3章程序分析 (5) 3.1 语音去噪采样过程 (5) 3.2 语音去噪方案 (6) 第 4 章总结 (8) 参考文献 (9) 附录 (10)

第1章绪论 1.1数字信号处理的意义 数字信号处理是利用计算机或专用处理设备，以数值计算的方法对信号进行采集、抽样、变换、综合、估值与识别等加工处理，借以达到提取信息和便于应用的目的。它在语音、雷达、图像、系统控制、通信、航空航天、生物医学等众多领域都获得了极其广泛的应用。具有灵活、精确、抗干扰强、设备尺寸小、造价低、速度快等优点。 数字滤波器, 是数字信号处理中及其重要的一部分。随着信息时代和数字技术的发展，受到人们越来越多的重视。数字滤波器可以通过数值运算实现滤波，所以数字滤波器处理精度高、稳定、体积小、重量轻、灵活不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊功能。数字滤波器种类很多，根据其实现的网络结构或者其冲激响应函数的时域特性，可分为两种，即有限冲激响应( FIR，Finite Impulse Response)滤波器和无限冲激响应( IIR，Infinite Impulse Response)滤波器。 FIR滤波器结构上主要是非递归结构，没有输出到输入的反馈，系统函数H (z)在处收敛，极点全部在z = 0处（因果系统），因而只能用较高的阶数达到高的选择性。FIR数字滤波器的幅频特性精度较之于IIR数字滤波器低，但是线性相位，就是不同频率分量的信号经过fir滤波器后他们的时间差不变，这是很好的性质。FIR 数字滤波器是有限的单位响应也有利于对数字信号的处理，便于编程，用于计算的时延也小，这对实时的信号处理很重要。FIR滤波器因具有系统稳定，易实现相位控制，允许设计多通带（或多阻带）滤波器等优点收到人们的青睐。 IIR滤波器采用递归型结构，即结构上带有反馈环路。IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成，可以组合成直接型、正准型、级联型、并联型四种结构形式，都具有反馈回路。同时，IIR数字滤波器在设计上可以借助成熟的模拟滤波器的成果，如巴特沃斯、契比雪夫和椭圆滤波器等，有现成的设计数据或图表可查，在设计一个IIR数字滤波器时，我们根据指标先写出模拟滤波器的公式，然后通过一定的变换，将模拟滤波器的公式转换成数字滤波器的公式。 滤波器的设计可以通过软件或设计专用的硬件两种方式来实现。随着MATLAB软

系统仿真结课作业

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

春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语言与应用课程作业MATLAB Simulink在电路暂态分析中的应用学生姓名陈志豪所在专业轮机工程（陆上）所在班级陆上1102 指导教师徐国保（博士） MATLAB Simulink在电路暂态分析中的应用（广东海洋大学轮机工程（陆上）1102 陈志豪）摘要本文通过引入举了实际的例子，简要介绍了Matlab语言在电工学电路暂态分析中的应用；并先使用普通方法分析暂态电路，然后再用Matlab Simulink来仿真暂态电路；通过Matlab Simulink 仿真技术，可以使得暂态分析可视化。关键词：MATLAB；Simulink仿真；电工学；暂态分析1，引言MATLAB是Matrix Laboratory的缩写,事实上MATLAB最初就是纯粹的矩阵计算软件。如今MATLAB既表示一种交互式的数值计算软件,又表示一门高级科学计算语言,是一套功能十分强大的工程计算及数据分析软件,其应用范围涵盖了数学、工业技术、电子科学、医疗卫生、建筑、金融、数字图像处理等各个领域。它把计算、图示 和编程集成到一个易用的交互式环境中,用大家熟悉的数学表达式来描述问题和求解方法,从而使许多用C 或FORTRAN实现起来十分复杂和费时的问题用MATLAB可以轻松地解决。许多工程师和研究人员发现,MATIAB能迅速测试其构思,综合评测系统性能,并能借此快速设计出更多的解决方案,达到更高的技术要求。[1]MATLAB因为提供了非常方便的绘图功能和强大的图形图像处理能力，以及强大的仿真技术,所以收到了广泛的欢迎。2，问题背景《电工学》是一门非电专业的技术基础课，通过本课程的学习，学生掌握电工技术的基本理论、基本定律、基本概念及基本分析方法和理论的实际应用。它的内容广泛，理 论性和系统性也很强。采用传统的教学模式，往往只能在理论上进行论述、推导、验证和证明，并借助 一些公式来阐述问题，很难给学生较直观的印象[2]，教学效果不理想。若把MATLAB应用到学习中，利用其强大的数值计算功能、绘图功能、可视化的仿真功能，可以很好的弥补传统教学的不足，使一些不 容易理解的抽象、复杂的变化过程，通过MATLAB仿真比较直观的的显示出来，便于学生理解和应用。 同时，可以随机修改电路和参数，即时观察输出结果，从而加深学生对电路本质的理解，全面掌握教学 内容[3]。下面通过实例探讨MATLAB SIMULINK在电工学暂态分析中的应用。图1所示电路是一个一阶电路。已知R=20Ω，U=6V，U=10V，C=O.02F。假s0 设在t=O时开关S从闭合在a端换路闭合到b端，求t>O时，电容电压u和电 c 容电流i。c图1，一阶电路的电路图 3，理论推导根据一阶电路暂态分析的三要素法有：（1）确定初始值由换路前的电路求得u(0)=U=10V C0再由换路后的电路求得 (??)??????????????????i(0)===?0.2A C??????(2)确定稳态值有电路图易知：i（∞）=0A Cu(∞)=6V C (3)确定时间常数τ=RC =20×0.02=0.4s (4),求出待求响应????????i= i（∞）+[ i(0)? i （∞）]e =?0.2 e ????.??CCCC????????u= u(∞)+[ u(0) ?uC(∞)] e =6+4e ????.??CCC运用Matlab 编程画出ic和u波形图； C 其代码如下所示：subplot(1,2,1); fplot(‘6+4*exp(-x/0.4)’,[0,6]); subplot(1,2,2); fplot(‘-0.2*exp(-x/0.4)’,[0,6]); 其运行结果图2所示 图2，电容电压u和电容电流ic波形图C4，应用MATLAB进行仿真图1电路对应的仿

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

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

matlab课程论文要求

matlab课程论文要求 一、时间安排 （一）2016年X月X日之前必须提交纸质版（时间待定，另行通知，尽早完成，以免影响其他科目的复习考试）。 （二）电子版统一写清楚学号（学号在前）+姓名+专业发送给学委。打包文件夹发送给我，不接受单独发给我的。 二、选题 （一）选题要紧密结合本学科专业的教学科研和MATLAB，符合专业培养目标的要求。 （二）论文一般为一人一题，严格控制与往年的重复率。 三、成绩评定 平时成绩（0.3）＋课程论文（0.7）=最终成绩。 四、论文写作规范要求 （一）封面：封面要使用统一格式。 （二）目录：“目录”两字黑体小二号、居中，“目录”两字间空四格、与正文空一行。各部分名为宋体小四号字，各小部分名间有缩进。 （三）题目：题目要对论文的内容有高度的概括性，简明、易读，字数应在20个字以内，论文题目用黑体三号字。 （四）署名：论文署名的顺序为：专业学号学生姓名指导老师姓名，用宋体小四号字。可用以下表示： 专业：XXXXX 学号：XXXXX 学生姓名：XXXXX 指导老师姓名：XXXX （五）内容摘要：中文内容摘应简要说明所研究的内容、目的、实验方法、主要成果和特色，一般为200－300字，用宋体小四号字，其中“内容摘要”四个字加粗。 （六）关键词：一般为3－6个，用分号隔开，用宋体小四号字，其中“关键词”三个字加粗。 （七）正文：正文要符合一般学术论文的写作规范，统一用宋体小四号字，行距为1.5倍。字数一般要求为不得少于5000字。

内容要理论联系实际，涉及到他人的观点、统计数据或计算公式的要注明出处（引注），涉及计算内容的数据要求准确。标题序号从大到小的顺序为：“1”“1.1”“1.1.1”……。 （八）注释：论文中所引用文献按学术论文规范注明出处，注序要与文中提及的序号一致。注释方法参见参考文献顺序。 （九）参考文献：论文后要标注参考文献和附录，参考文献按照以下格式排列： 1．专著、论文集、学位论文、报告 [序号]主要责任者．文献题名[文献类型标识]．出版地：出版者，出版年．起止页码。 [1]刘国钧，陈绍业，王凤．图书馆目录[M]．北京：高等教育出版社，1957．10-12. [2]辛希孟．信息技术与信息服务国际研讨会论文集：A集[C]．北京：中国社会科学出版社，1994．12-13. [3] 查正军.《基于机器学习方法的视觉信息标注研究》.[D].北京.中国科技大学.2010年.32-35 2．期刊文章 [序号]主要责任者．文献题名[J]．刊名，年卷（期）：起止页码． [1]何龄修．读顾城《南明史》[J]．中国史研究，1998（3）：12-13. [2]金显贸，王昌长，王忠东等．一种用于在线检测局部放电的数字滤波技术 [J]．清华大学学报（自然科学版），1993（4）:12-13. 3．电子文献 [序号]主要责任者．电子文献题名[电子文献及载体类型标识] ．电子文献的出处或可获得地址，发表或更新日期／引用日期（任选）． [1]王明亮.关于中国学术期刊标准化数据库系统工程的进展[EB/OL]. https://www.wendangku.net/doc/8e4287572.html,/pub/wml.txt/980810-2.html,1998-08-16/1998-10-04. [2]万锦坤.中国大学学报论文文摘（1983-1993）.英文版[DB/CD].北京:中国大百科全书出版社，1996.

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

期末大作业报告

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

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

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

《科学计算与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结课论文

MATLAB结课论文 题目：基于matlab的双音频电话机的图形界面 装 订 线 学院信息工程学院 专业通信工程 学号 姓名 任课教师赵晓燕 2013年 5 月28 日

摘要 MATLAB是美国MathWorks公司出品的商业数学软件，本学期通过对matlab的学习，我们了解了matlab的基本使用方法，并能很好的利用matlab进行信号与系统等课程的分析和学习，对我们今后的学习和工作有很大的帮助，本文是matlab的结课论文，题目要求是创建双音频电话机的图形用户界面（phone）、创建演示抽样定理的图形用户界面或者用matlab设计电子音乐。我选择了设计双音频电话机的图形用户界面。本文叙述的是制作双音频电话机图形界面的主要过程。 关键字：matlab 双音频图形用户界面

一、matlab简介 MATLAB（矩阵实验室）是MATrix LABoratory的缩写，是一款由美国The MathWorks 公司出品的商业数学软件。MATLAB是一种用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。除了矩阵运算、绘制函数/数据图像等常用功能外，MATLAB还可以用来创建用户界面及与调用其它语言（包括C，C++和FORTRAN）编写的程序。 尽管MATLAB主要用于数值运算，但利用为数众多的附加工具箱（Toolbox）它也适合不同领域的应用，例如控制系统设计与分析、图像处理、信号处理与通讯、金融建模和分析等。另外还有一个配套软件包Simulink，提供了一个可视化开发环境，常用于系统模拟、动态/嵌入式系统开发等方面。1970年代末到80年代初，时任美国新墨西哥大学教授的克里夫·莫勒尔为了让学生更方便地使用LINPACK及EISPACK（需要通过FORTRAN 编程来实现，但当时学生们并无相关知识），独立编写了第一个版本的MATLAB。这个版本的MATLAB只能进行简单的矩阵运算，例如矩阵转置、计算行列式和本征值，此版本软件分发出大约两三百份。 1984年，杰克·李特、克里夫·莫勒尔和斯蒂夫·班格尔特合作成立了MathWorks 公司，正式把MATLAB推向市场。MATLAB最初是由莫勒尔用FORTRAN编写的，李特和班格尔特花了约一年半的时间用C重新编写了MATLAB并增加了一些新功能，同时，李特还开发了第一个系统控制工具箱，其中一些代码到现在仍然在使用。C语言版的面向MS-DOS 系统的MATLAB 1.0在拉斯维加斯举行的IEEE决策与控制会议（IEEE Conference on Decision and Control）正式推出，它的第一份订单只售出了10份拷贝，而到了现在，根据MathWorks自己的数据，目前世界上100多个国家的超过一百万工程师和科学家在使用MATLAB和Simulink。 1992年，学生版MATLAB推出；1993年，Microsoft Windows版MATLAB面世；1995年，推出Linux版。 MATLAB的主要提供以下功能： ①可用于技术计算的高级语言 ②可对代码、文件和数据进行管理的开发环境 ③可以按迭代的方式探查、设计及求解问题的交互式工具 ④可用于线性代数、统计、傅立叶分析、筛选、优化以及数值积分等的数学函数 ⑤可用于可视化数据的二维和三维图形函数

利用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软件的各菜单、工具栏及常用命令的使用。掌握MatLab有关矩阵的创建方法、矩阵的基本运算符、矩阵的计算函数。掌握MatLab的符号运算。熟练掌握二维、三维图形的绘制；掌握简单动画的制作；了解分形几何学，绘制Koch 雪花曲线和Minkowski“香肠”曲线。熟炼掌握MatLab程序设计的顺序、分支和循环结构;熟炼掌握脚本M文件和自定义函数的设计和使用；复习高等数学中有关函数极限、导数、不定积分、定积分、二重积分、级数、方程近似求解、常微分方程求解的相关知识.通过作图和计算加深对数学概念：极限、导数、积分的理解.学会用MatLab 软件进行有关函数极限、导数、不定积分、级数、常微分方程求解的符号运算；了解数值积分理论,学会用MatLab软件进行数值积分;会用级数进行近似计算.复习线性代数中有关行列式、矩阵、矩阵初等变换、向量的线性相关性、线性方程组的求解、相似矩阵及二次型的相关知识.学会用MatLab软件进行行列式的计算、矩阵的基本运算、矩阵初等变换、向量的线性相关性的判别、线性方程组的求解、二次型化标准形的运算. 二作图应用 1描点作图 2显函数作图（fplot） 3隐函数作图（ezplot）

4参数方程作图(ezplot) 5极坐标作图 6 空间曲面作图 三高等数学应用 1 极限问题 2 求导数问题 3 求不定积分问题 4 求定积分问题 5 求偏导数问题 6 求二重积分问题 7 求级数和的问题 8求函数的泰勒展开式问题 9 求常微分方程的通解和特解问题 四线性代数应用 1 行列式问题 2 矩阵运算问题（包括加，减，乘，转置，求逆，求秩） 3 一般的线性方程组求解问题(jfch) 4 向量组的线性相关性问题 5 矩阵的特征值和特征向量问题 6 二次型化标准型问题

MATLAB语言课程论文

《MATLAB语言》课程论文 基于MATLAB在高等数学基础部分中 的应用

基于MATLAB的控制系统分析 摘要：运用MATLAB软件处理和仿真,分析所建立的控制系统模型的可行性。利用MATLAB平台分析系统传递函数的稳定性，利用MATLAB软件得出系统的单位阶跃响应和脉冲响应，分析说明MATLAB在自动控制系统方面的应用。 关键词：控制系统MA TLAB 稳定性Simulink 0引言 自动控制技术已经广泛应用于工业、农业、交通运输业、航空及航天业等众多产业部门，极大地提高了社会劳动生产率，改善了人们的劳动条件，丰富与提高了人们的生活水平。在当今的社会生活中，自动化装置无所不在，为人类文明进步做出了重要贡献。随着科学技术的发展，控制系统变得越来越复杂，控制理论和系统的分析如果仅靠人工计算已经远远不能满足学习和研究的要求了。MATLAB是高性能的数值计算和可视化软件，他集数值分析、矩阵运算、信号处理和图形显示于一体，构成了一个方便、界面友好的用户环境。MATLAB的出现给控制系统的分析提供了极大的方便。 MATLAB(Matix Laboratory,即“矩阵实验室”)是美国MathWorks公司开发的用于概念设计、算法开发、建模仿真、实时实现的理想的集成环境。从1984年MATLAB诞生到现在，它已经集成了许多工具箱，例如控制系统工具箱(Control System Toolbox)、Simulink工具箱等。为此，MATLAB软件在控制工程领域已获得了广泛的应用，使得MATLAB成为控制系统仿真分析的必要工具 1系统的稳定性分析 稳定是控制系统的重要性能，也是系统能够正常运行的首要条件。在分析控制系统时，首先遇到的问题就是系统的稳定性。对线性系统来说，如果一个系统的所有几点都位于左半s平面，则该系统是稳定的。对于离散系统来说，如果一个系统的全部极点都在单位圆内，则该系统可以被认为是稳定的。由此可见，线性系统的稳定性完全取决于系统的极点在根平面上的位置。 判断一个线性系统稳定性的一种最有效的方法是直接求出系统所用的极