MATLAB结课作业(2016.12)

计算应用软件

MATLAB程序设计及绘图

班级：

学号：

姓名：

2016.12

题目1：分别利用For循环和While循环求下值：

63

263 0

21222

i

K==++++

∑ 。

For循环代码：

k=0;

for i=0:1:63

k=k+2^i;

end

k

While循环代码:

k=0;

i=0;

while (i<64)

k=k+2^i;

i=i+1;

end

k

结果：k = 1.8447e+019

题目2：已知方程组：

347124

57423

859

652108

x y z w

x y z w

x z w

x y z w

+--=

?

?-++=-

?

?

+-=

?

?-+-+=-

?

，求方程组的解。

代码：

A=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10]; B=[4;-3;9;-8];

X=A\B

结果：X =

-1.4841

-0.6816

0.5337

-1.2429

题目3：某试验试件在受火后，其弹性模量发生变化，数据如下：受火温度T=[20, 100, 200, 300, 400, 500, 600, 700, 800, 900];

受火后弹性模E=[21.95,1.85,1.76,1.70,1.55,1.50,1.42,1.35,1.3,1.32]*1011;

用合适的多项式拟合弹性模量随受火温度变化的规律。

代码：

T=[20,100,200,300,400,500,600,700,800,900];

E=[21.95,1.85,1.76,1.70,1.55,1.50,1.42,1.35,1.3,1.28]*10^11;

plot(T,E,'LineWidth',4)

结果：

题目4：（1）编写一个函数式文件，使其产生如下函数段：

22()50.52646

x

x f x x

x x x <=??

=-<<=??->?

（2）编写命令式文件，绘制x 在[0,10]范围内f (x ).f (x +2)的曲线。

（1）代码： function f=Untitled(x) if x<=2 f=2*x; elseif x>6 f=x-4; else

f=5-0.5*x; end

（2）代码： i=0;

for x=0:0.1:10 i=i+1;

y1(1,i)=Untitled(x); end i=0

for x=0:0.1:10 i=i+1;

y2(1,i)=Untitled(x+2); end y3=y1.*y2; x=0:0.1:10;

plot(x,y3,'LineWidth',5)

题目5：设f(x)=x5-4x4+3x2-2x+6；

（1）取x=[-2,8]之间函数的值（取100个点），画出曲线，看它有几个零点。（提示：用polyval函数）

（2）用roots函数求此多项式的根。

（1）代码：p=[1,-4,3,-2,6];

x=-2:0.1:8;

y=polyval(p,x);

plot(x,y,'LineWidth',5)

（2）代码：

r=roots(p)

n=length(r);

for i=1:n

if r(i)>=-2 && r(i)<=8 if (isreal(r(i)))

r(i)

end

end

end

结果：

ans =

3

ans =

1.6956

题目6：在x =[-10,10]，y =[-10,10]范围内画出函数z =的三维图形。

代码：

x=-10:0.01:10; y=-10:0.01:10;

[X,Y]=meshgrid(x,y);

Z=sin(sqrt(X.^2+Y .^2))./sqrt(X.^2+Y .^2); mesh(X,Y ,Z) 结果：

题目7：由100*rand(30)作为全班成绩S (i )，用饼图分别表示60i S <、6070i S ≤<、

7080i S ≤<、8090i S ≤<、90i S ≤所占比例。

代码： S=100*rand(30);

k1=0;k2=0;k3=0;k4=0;k5=0; for i=1:900 if S(i)<60 k1=k1+1;

elseif S(i)>=60&S(i)<70

k2=k2+1;

elseif S(i)>=70&S(i)<80

k3=k3+1;

elseif S(i)>=80&S(i)<90

k4=k4+1;

else

k5=k5+1;

end

end

x=[k1 k2 k3 k4 k5];

y=[1 1 1 1 1];

pie(x,y,{'Si<60','60<=Si<70','70<=Si<80','80<=Si<90','90<=Si'}); 结果：

题目8：列举“MATLAB程序设计课程”在本专业的一个相关应用实例。

求解任意长度为L的简支梁在任意大小的集中F的作用下的弯矩图和剪力图，最大剪力值:F1max和最大弯矩值Mmax：。

代码：

clc

clear;

disp('输入集中力到A点的距离a,梁的长度l,集中力F')

a=input('a （/m）值：');

L=input('L （/m）值：');

F=input('F （/KN）值：');

M=[];

n=0;

for x=0:0.1:L;

if x>=0 && x<=a

n=n+1;

F1(n)=(L-a)/L*F;

M(n)=(L-a)/F*x;

else

n=n+1;

F1(n)=-a/L*F;

M(n)=a/L*F*(L-x);

end

end

x=0:0.1:L;

y2=0*x;

plot(x,M,'linewidth',3)

title('弯矩图')

xlabel('x /m')

ylabel('M /KN*m') figure(2)

plot(x,F1,'linewidth',3) hold on

plot(x,y2,'linewidth',8) legend('剪力F1','梁y=0') title('剪力图')

xlabel('x /m')

ylabel('F1 /KN')

F1max=F*a/L

Mmax=F*a*(L-a)/L

结果：

F1max =

5 Mmax =

250

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.人在识别事物时是否可以避免错识 2.如果错识不可避免，那么你是否怀疑你所看到的、听到的、嗅 到的到底是真是的，还是虚假的 3.如果不是，那么你依靠的是什么呢用学术语言该如何表示。 4.我们是以统计学为基础分析模式识别问题，采用的是错误概率 评价分类器性能。如果不采用统计学，你是否能想到还有什么合理地分类 器性能评价指标来替代错误率 1.知觉的特性为选择性、整体性、理解性、恒常性。错觉是错误的知觉，是在特定条件下产生的对客观事物歪曲的知觉。认知是一个过程，需要大脑的参与.人的认知并不神秘，也符合一定的规律，也会产生错误 2.不是 3.辨别事物的最基本方法是计算 . 从不同事物所具有的不同属性为出发点认识事物. 一种是对事物的属性进行度量，属于定量的表示方法(向量表示法 )。另一种则是对事务所包含的成分进行分析，称为定性的描述(结构性描述方法)。 4.风险 第二次课堂作业 作为学生，你需要判断今天的课是否点名。结合该问题(或者其它你熟悉的识别问题，如”天气预报”)，说明: 先验概率、后验概率和类条件概率 按照最小错误率如何决策 按照最小风险如何决策 ωi为老师点名的事件,x为判断老师点名的概率 1.先验概率: 指根据以往经验和分析得到的该老师点名的概率,即为先验概率 P(ωi ) 后验概率: 在收到某个消息之后，接收端所了解到的该消息发送的概率称为后验概率。 在上过课之后,了解到的老师点名的概率为后验概率P(ωi|x) 类条件概率:在老师点名这个事件发生的条件下,学生判断老师点名的概率p(x| ωi ) 2. 如果P(ω1|X)>P(ω2|X)，则X归为ω1类别 如果P(ω1|X)≤P(ω2|X)，则X归为ω2类别 3.1)计算出后验概率 已知P(ωi)和P(X|ωi)，i=1,…，c，获得观测到的特征向量X 根据贝叶斯公式计算 j=1,…，x

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) ,基于仿真设计可以在计算机上建立虚拟样机,对产品的外形、结构、强度、动力

西南大学含弘在线课堂PowerPoint多媒体课程第一次作业答案

一、 单项选择题(本大题共5小题，共25.0分) 1. 传播信息的电缆、电磁波等称为： (5.0分 ) A ：媒介 选择： B B ：媒材 C ：媒体 D ：以上三个均不对 2. 在演示文稿制作时，如果在每一页幻灯片上添加版权信息，则可以在下列哪一种视图模式中进行操作设置： (5.0分 ) A ：讲义母版 选择： B B ：幻灯片母版 C ：备注母版 D ： 备注页视图 3. 在演示文稿中，如果仅需要查看幻灯片内的文字信息，而不需要观看其缩略图内容，则可以选择下列哪一个选项卡： (5.0分 ) A ：大纲选项卡 选择： A B ：幻灯片选项卡

C ：属性选项卡 D ：以上三个均不对 4. 在计算机中用构成一幅图的位置、维数、大小和形状的指令来进行描述的是： (5.0分 ) A ：图形 选择： B B ：图像 C ：图形与图像 D ：以上三个均不对 5. 适用于表现比较细节、层次以及丰富色彩的是 ： (5.0分 ) A ：图形 选择： C B ：图像 C ：图形与图像 D ：以上三个均不对 二、 多项选择题(本大题共2小题，共10.0分) 1. 多媒体技术所涉及到以下哪些形式: (5.0分 ) 选择：A,B,C

A ：媒介 B ：媒材 C ：媒体 D ：以上三个均不对 2. 评价多媒体课件质量的指标有 : (5.0分 ) A ：科学性 选择： A,B,C,D B ：技术性 C ：艺术性 D ：教育性 三、 判断题(本大题共3小题，共15.0分) 1. 多媒体课件在开发时其艺术性必须服务于教育性。 (5.0分) √ × 选择：√ 2. 多媒体课件在开发时要求必须使用多种媒体元素。 (5.0分) √ × 选择：× 3. 多媒体课件脚本要严格按照模板制作，不能有个性化调整。 (5.0分) √ × 选择：× 四、

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

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

春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)。

第一次作业

数理经济学第一次作业题 董志勇 1. You can write your answer in either English or in Chinese or in both. 2. You have to do this paper independently; giving identical answers will be assumed cheating, which will lead to severe punishment. 3. Answer to the point; be precise and avoid unnecessary length. 1、单个垄断企业分别在两个市场上实现价格歧视，求证其曲线从下面与各曲线的水平之和相交。 2、设，为两种投入要素，共同生产一种产品，该产品的价格是，生产函数是，的价格是，的价格是，，是使得利润最大化的要素投入量，求证：和不会同时小于0。 3、生产函数是，求出要素需求曲线，对根据该生产函数进行生产而且追求利润最大化的企业进行比较静态分析，证明对该企业而言，交叉效应项的符号为负。 4、一个垄断者在两个市场上销售同一种产品，市场1的收益为，市场2的收益为，两个市场的总成本为，假设在市场1中，每销售一个单位的产品要征税t。 (a)证明t增加将会使得市场1销售得产品量降低。 (b)有关t增加对市场2中的产品的影响，最大化假设意味着什么？ (c)即使假设二阶条件总是成立，市场1中的税收的增加也可能导致总产出，请说明在什么情况下会是这样的。 (d)假设根据政府法律，市场2的产出固定在前面所述的利润最大化水平上，与不加管制的情况相比，在受到政府管制的情况下，市场1中的产出对税收增加得反应在绝对量上较小，请加以证明，并对此进行直观的经济解释。 5、垄断者在两个市场中销售其产品，收入函数分别为和。总产出的函数。在两个市场中，每销售一个单位的产品，征收同样的税收t。 (a)求解，和。如果有解的话，利润最大化隐含了哪些偏导的符号？ (b)假设产出固定不变，求解，能够确定该式的符号？

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大作业

选 题 说 明 本人选做第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课程论文要求

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/6b17424285.html,/pub/wml.txt/980810-2.html,1998-08-16/1998-10-04. [2]万锦坤.中国大学学报论文文摘（1983-1993）.英文版[DB/CD].北京:中国大百科全书出版社，1996.

人教版小学四年级下册语文课堂作业本答案

语文园地二 1、精湛财富属于朴实器重喉结例外赞赏闷闷不乐狼吞虎咽面呈菜色 2、维持生计得到宽容回忆往事动作熟练神情严肃生活拮据 3、(1)行果 (2)言而有信 (3)勿施于人学而不厌，诲人不倦 4、(1)深信(2)自信(3)守信(4)诚信 (1)寡言少语(2)豪言壮语(3)花言巧语(4)胡言乱语 5、(1)严肃道德难题，把车留下来还是还给库伯 （2）小明这次语文考得很差，试卷发下来后，他战战兢兢地递给了妈妈。妈妈的脸顿时由晴转阴，大声嚷道：“你上课都听了些什么呀？这么简单的题目也会错？” 6、(1)不小心，文中指小男孩不小心掉进了大西洋。肯定，一定文中指船长肯定孩子掉进了大海。(2)肯定/确定用力/尽力/奋力的确/确实 (3)他实在游不动了，准备放弃，加上没人发现他掉进大海里，他觉得自己活不下去了认为船长一定会返航去救他的小男孩相信船长一定会去救他，说明船长被孩子所信任 (4)一个人能被他人相信也是一种幸福。他人在绝望时想起你，相信你会给予拯救更是一种无可替代的幸福，我们要真诚地去对待别人。 九、自然之道 1、沙滩侦察蠢事棒球帽愚不可及无遮无拦气喘吁吁 2、啄逐逐啄彻切彻切 4、(1)这个词突出了向导非常不情愿，生怕惹出麻烦，对同伴们善良的呼喊非常理解，但也感到非常无奈 (2)这个词突出了向导悲伤地叹息，对自

己之前的所作所为感到非常羞愧，表达了“自己”自责、悔恨的心理。5、(1)看着说飞走答道抱起 (2)R 艹 5 ① T 页 7 ③ (3)眼前的一幕是指，当幼龟踌躇不前时，一只嘲鸫突然飞来，它用尖嘴啄幼龟的头，企图把它拉到沙滩上去。 (4) ③ 十一、蝙蝠和雷达 1、清朗蝙蝠飞蛾避开蚊子揭开障碍荧光屏 2、bǔ líng qī zhuàng cāng ruì 3、启发灵敏(灵活) 仿照清晰 4、横七竖八地拉了许多绳子系着许多铃铛 5.(1)√×× (2)①波浪②天上的云/春天的风小羊/妈妈的手飘来飘去/吹在脸上真舒服。 （3）蝙蝠用嘴发射超声波，超声波向四周传播，遇到障碍物就反射回来，反射回来的那部分被蝙蝠用耳朵接收，经过大脑处理，就知道哪里有障碍物了。 6、嘴巴超声波耳朵青蛙－电子蛙眼蜻蜓－飞机蝴蝶－卫星螳螂－镰刀顺风耳－电话

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

利用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》期末大作业

给出程序、图、作业分析，程序需加注释。 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)]] %输出矩阵

人教版小学四年级下册语文课堂作业本答案

人教版小学四年级下册语文课堂作业本答案 一、古诗词三首 1、敬亭山洞庭湖青螺 2.×√××√ 3.(1)完了，没有了悠闲孤独众孤因怀才不遇而感到的孤独、寂寞的心情 (2)月光下的洞庭湖洞庭湖的君山翠白青 (3)胜过仿佛，像蓝草，一种植物，叶蓝绿色，可提取青蓝色染料5、江南春唐杜牧千里莺啼绿映红，水村山郭酒旗风。南朝四百八十寺，多少楼台烟雨中。 二、桂林山水 1、无瑕攀登屏障骆驼危峰兀立波澜壮阔奇峰罗列 2、sì juàn dào m?ng 3、波澜壮阔水平如镜峰峦雄伟红叶似火衣服（图画）山峰（山峦） 4、(1)秋天，落叶仿佛像一个会跳舞的精灵，从树上翩翩飞舞而下。 (2)老房子的围墙真危险呐，手一碰，墙壁上的石灰好像要掉下来似的。 5、(1)M 辶 6 ④ D 斤 7 ② (2)静清绿奇秀险 (3) 绿树红花竹筏小舟云雾山水 6、田野上的野花真多啊，多得像天上的星星数也数不清；田野上的野花真美啊，美得仿佛那是一条五彩的锦缎；田野上的野花真香啊，香得引来群群蜜蜂争相采蜜！（我游览过水平如镜的西湖，玩赏过飞泻而下的瀑布，却从没看到过秦淮河的水。秦淮河的水真凉啊，凉得让你想马上和河水来

个亲密的接触；秦淮河的水真清啊，清得就像是一块彩色的大玻璃；秦淮河的水真亮啊，当河边的万家灯火映在水中，像是水底铺满了熠熠生辉的大钻石，使河水变得无比明亮、美丽；秦淮河的水真野啊，刚才还是块大玻璃，现在就波涛起伏，使得河上的画舫也一个劲儿地摇摆，再加上朱雀桥洞内随风而舞的狂草，更体现了秦淮河的粗犷之美；秦淮河的水真神啊，要不，它怎么会吸引来那么多的小雨点儿？想当年大自然在创造世界时，一定挥舞着最大的画笔，而放弃了那些精细的刻刀，毫无雕琢地在南京勾勒出了这样一条粗犷、神奇、美丽的秦淮河吧？） 三、记金华的双龙洞 1、浙江臀部稍微额角擦伤蜿蜒依据石钟乳 2、一簇臀部稍微蜿蜒民族宫殿树梢大碗 3、(1)②① (2)①② 4、（1）小妹妹总是时而笑，时而哭，时而忧，时而喜，真是捉摸不透。（六 月的天气，时而风，时而雨，时而凉，时而热，真是变化多端。）（2）即使这件事情再困难，我们也要尽自己最大的力量去做好它。 5、(1)根据名称欣赏 (2)，，、、， (3)在洞顶的那两条龙，一条黄龙，一条青龙。 (4)蜿蜒在洞顶的双龙变化多端，颜色各异的石钟乳和石笋。 6、怎样美的笔盒呢？让人一看就爱不释手，是这样美的笔盒。（怎样简陋的小屋呢？裂缝的石灰墙，风化的水泥地，一扇破旧的窗，一张油漆脱落的书桌，根本找不到一样值钱的东西，就是这样一间简陋的小屋。） 语文园地一 1、洞庭玩赏扩散泰山拥挤孔隙仰卧油桐水平如镜拔地而起突兀

MATLAB课程论文

基于MATLAB在自动控制频域中稳定性分析的应用 学院：物信学院 班级：08电信二班 姓名：王军祥 学号：281060217

基于MATLAB在自动控制频域中稳定性分析的应用 摘要:自动控制系统主要利用MATLAB高级语言对其进行计算机分析。 MATLAB是一套高性能的数值计算和可视化软件,它集数值分析、矩阵运算、信号处理和图形绘制集于一体,构成了一个方便的、界面友好的用户环境。本文主要介绍了利用MATLAB在自动控制中对频域中系统稳定性的判定，通过MATLAB建立某一系统的模型,并分析该系统的性能。根据响应曲线判断系统的稳定性，当系统的性能不能满足所要求的性能指标时,通过调整系统参数和增添校正装置来改善系统性能并展示方便灵活的动态仿真结果。 关键词:自动控制系统；稳定性；频域分析；频率响应；稳定裕度 引言 频域分析法是应用频域特性研究线性控制系统的一种经典方法，采用这种方法可以直观的表达出系统的频率特性，利用系统的传递函数绘制系统的bode 图、nyquist曲线和nichols图，然后进行系统稳定性的判定，这样可以很明确的判定系统频域的稳定性。频域分析法是自动控制领域中应用又一种数学工具———频率特性来研究系统控制过程性能，即稳定性、快速性及稳态精度的一种方法。这种方法不必直接求解系统的微分方程，而是间接的运用系统的开环频率特性曲线，分析闭环系统的响应，因此它是一种图解的方法。本文介绍了应用MATLAB在bode图、nyquist曲线和nichols图等控制系统频域分析中的主要方法，通过具体实例叙述了MATLAB在频域分析中的应用过程。频域分析里主要用到三种曲线（或叫图）：Bode图、Nyquist曲线图和 Nichols（尼柯尔斯）曲线图。这三种曲线就是频率分析的三种工具。Bode图可以用于分析相角稳定裕度、 -穿越频率、带宽、扰动抑制及其稳定性幅（或模）值稳定裕度、剪切频率、π 等，所以Bode图在频域分析里占有重要的地位。Nyquist与Nichols曲线图在频域分析里也很有用。本文中最重要的函数命令有bode、nyquist、pade、nichols、margin等。 一.频域分析法的基础 1.有关频率分析的几个概念 （1）频率响应 当正弦函数信号作用于线性系统时，系统稳定后输出的稳态分量仍然是同频率的正弦信号，这种过程叫做系统的频率响应。