实验一 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内所有素数（在编辑窗口中实现，命令窗口中调用）。
电气学科大类 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
复数的计算和图示 表3 MATLAB关于复数运算的函数 {范例3_1}复数的加减法 设有两个复数z1 = 1 + 2i和z2 = 4 + 3i，其中i是虚数单位求两个复数的和z1 + z2和差z2–z1。 [解析]复数有三种表示形式 (1)代数式 z = x + i y(3_1_1) (2)三角式 z = r(cosθ + isinθ) (3_1_2) 其中r是复数的模，θ是复角。代数式与三角式的换算关系是 r=arctan y x θ(3_1_3) x = r cosθ，y = r sinθ(3_1_4) (3)指数式 z = r e iθ(3_1_5) 其中利用了欧拉公式 e iθ= cosθ + isinθ(3_1_6) 设有两个复数 z1 = x1 + i y1，z2 = x2 + i y2(3_1_7) 复数加法是 z= z1 + z2 = (x1 + x2) + i(y1 + y2) (3_1_8) 复数减法是 z= z1 - z2 = (x1 - x2) + i(y1 - y2) (3_1_9) [程序]P3_1plus.m如下。 %复数的加减法 clear %清除变量 z1=1+2i; %第1个复数(1) x1=real(z1); %取第1个复数的实部(2) y1=imag(z1); %取第1个复数的虚部(2) x2=4; %第2个复数的实部 y2=3; %第2个复数的虚部 z2=x2+i*y2; %形成第2个复数(3) z=z1+z2; %两复数之和(4) x=real(z); %取复数的实部 y=imag(z); %取复数的虚部 figure %创建图形窗口 quiver(0,0,x1,y1,0) %在复平面画第1个复数(5)
系统仿真导结课作业 一、概述 建模与仿真技术已成功地应用于航空航天、生产制造、交通运输、信息、生物、医学、材料、能源、教育、军事、社会、经济等众多领域;并成功地应用于产品研制的全生命周期,包括需求分析、方案论证、概念设计、初步设计、详细设计、生产制造、试验试飞、运行、维护、训练等各个阶段。仿真科学与技术正是从其广泛的应用中获得了日益强大的生命力，而仿真技术的发展反过来使得其得到愈来愈广泛的应用。 广义而言,仿真是采用建模的方法和物理的方法对真实环境客观事物进行 抽象、映射、描述和复现。基于系统原理、理论、定律、系统数据等应用计算机技术、软件技术和信息技术建立仿真环境(虚拟环境) ,在仿真环境中对客观事物进行研究。客观事物包括真实环境中的实体/系统、自然环境(地形、大气、海洋、空间)、和人的行为(操作、决策、推理)。仿真环境包括模型、数据、软件、物理效应设备、计算机等。 计算机仿真的三要素是系统、模型、计算机，三个基本活动是模型设计、模型执行、模型分析。计算机仿真的三要素和三个基本活动的关系关系如图1 所示。 系统 模型设计模型分析 模型执行 模型计算机 图1 从模型设计到模型分析经历的过程，即对实物进行仿真可概括为以下几个方面： 1) 问题的描述; 2) 建立概念模型; 3) 建立仿真模型; 4) 收集数据; 5) 编写程序; 6) 在计算机进行模型试验; 7) 模型和数据的验证; 8) 仿真结果显示; 9) 仿真结果分析和评估。 仿真是建立模型在计算机上运行,但这属于数学仿真,随着技术的发展,许多应用 领域建立仿真系统时除了模型外还要求将实物和人员包含在仿真回路中。 由此可见，基于仿真设计与传统设计的方法和流程两者有很大区别(图2) ,基于仿真设计可以在计算机上建立虚拟样机,对产品的外形、结构、强度、动力
实验五 MATLAB 计算的可视化（一） 实验目的 1. 熟练掌握MATLAB 二维曲线的绘制 2.掌握图形的修饰 3.掌握三维图的绘制 4.了解各种特殊图形的绘制 内容与步骤 1.在同一幅图形窗口中分别绘制y1=sin(t)和y2=cos(t)二条函数曲线，t 的取值范围为[0,10]。y1用红色虚线表示，y2用蓝色实线表示，横坐标轴名称为“时间t ”，纵坐标轴名称为“正弦、余弦”，整个图形的标题为“正弦和余弦曲线”。在坐标(1.7*pi ，-0.3)处添加文字“sin(t)”， 在坐标(1.6*pi ，0.8)处添加文字“cos(t)”，并在右上角添加图例，其运行界面图如下图所示。之后并尝试修改坐标轴刻度。 2.用subplot 命令在同一个窗口的不同子窗口绘制曲线y=sin(t)，y1=sin(t+0.25) y2=sin(t+0.5)，其中t=[0 10]。 3.绘制三维曲线：?? ? ??=≤≤==)cos()sin()200() cos()sin(t t t z t t y t x π （注意：用plot3命令） 4.三维网线图：绘制z=sin(y)cos(x) 三维网线图。 5. 三维曲面图 绘制22y x z +=的三维曲面图，x 在[-5,5]范围，y 在[-5,5]范围。将曲面图颜色用shading 命令连续变化，并用颜色标尺显示色图（使用函数colorbar 生成）。生成的图形如下图所示。
6.请绘制一个饼形图，数据如下表所示 7. 用semilogx命令绘制传递函数为1//(s+1)(0.5s+1)的对数幅频特性曲线，横坐标为w，纵坐标为Lw，w的范围为10-2-103，按对数分布。
Matlab 特殊图形和高维可视化 2009-10-20 01:06 7.4 特殊图形和高维可视化 7.4.1 特殊图形指令例示 7.4.1.1 面域图area 【* 例7.4.1 .1-1 】面域图指令area 。该指令的特点是：在图上绘制多条曲线时，每条曲线（除第一条外）都是把“前”条曲线作基线，再取值绘制而成。因此，该指令所画的图形，能醒目地反映各因素对最终结果的贡献份额。注意：（ 1 ）area 的第一输入宗量是单调变化的自变量。第二输入宗量是“各因素”的函数值矩阵，且每个“因素”的数据取列向量形式排放。第三输入宗量是绘图的基准线值，只能取标量。当基准值为0 （即以x 轴为基准线）时，第三输入宗量可以缺省。（2 ）本例第<4> 条指令书写格式x' , Y' ，强调沿列方向画各条曲线的事实。 clf;x=-2:2 % 注意：自变量要单调变化 Y=[3,5,2,4,1;3,4,5,2,1;5,4,3,2,5] % 各因素的相对贡献份额 Cum_Sum=cumsum(Y) % 各曲线在图上的绝对坐标 area(x',Y',0) %<4> legend(' 因素A',' 因素B',' 因素C'),grid on,colormap(spring) x = -2 -1 0 1 2 Y = 3 5 2 4 1 3 4 5 2 1 5 4 3 2 5 Cum_Sum = 3 5 2 4 1 6 9 7 6 2 11 13 10 8 7
图 7.4.1 .1-1 面域图表现各分量的贡献 7.4.1.2 各种直方图bar, barh, bar3, bar3h 【 * 例 7.4.1 .2-1 】二维直方图有两种图型：垂直直方图和水平直方图。而每种图型又有两种表现模式：累计式：分组式。本例选其两种加以表现。 x=-2:2; % 注意：自变量要单调变化 Y=[3,5,2,4,1;3,4,5,2,1;5,4,3,2,5]; % 各因素的相对贡献份额 subplot(1,2,1),bar(x',Y','stacked') % “累计式”直方图 xlabel('x'),ylabel('\Sigma y'),colormap(cool)% 控制直方图的用色legend(' 因素 A',' 因素 B',' 因素 C') subplot(1,2,2),barh(x',Y','grouped') % “分组式”水平直方图 xlabel('y'),ylabel('x') 图 7.4.1 .2-1 二维直方图 clf;x=-2:2; % 注意：自变量要单调变化 Y=[3,5,2,4,1;3,4,5,2,1;5,4,3,2,5]; % 各因素的相对贡献份额 subplot(1,2,1),bar3(x',Y',1) % “队列式”直方图 xlabel(' 因素 ABC'),ylabel('x'),zlabel('y') colormap(summer) % 控制直方图的用色 subplot(1,2,2),bar3h(x',Y','grouped') % “分组式”水平直方图 ylabel('y'),zlabel('x')
开题报告 论文题目：基于matlab的时域有限差分法的电磁仿真研究（10分） 学院：电气工程及其自动化学院学号：1103000105姓名：__杨志刚___ 一、论文选题的目的和意义（300字以内；15分） 时域有限差分法，因具有多种优点被运用到电磁场理论研究的各个方面，而且其使用成效和应用领域还在迅速扩大和提高，在现代电磁场理论研究中具有很大的重要性和很强的可操作性。但是同时这种方法也存在一定的缺陷，主要表现在对无边界问题需要吸收边界条件处理，有色散误差，消耗内存大等方面。本课题在利用时域有限差分法对一些实际的算例进行实验仿真和验证，同时对这种方法在解决实际问题的缺陷进行一定程度的研究和分析。 Matlab作为一种工程仿真工具得到了广泛应用。用于时域有限差分法，可以简化编程，使研究者的研究重心放在FDTD法本身上，而不必在编程上花费过多的时间。 二、国内外关于该论题的研究现状和发展趋势（500字以内；15分） 时域有限差分方法作为一种典型的全波时域分析方法，因其原理直观、编程简便、实用性强在目前的计算电磁学领域内被人们广泛深入地研究，并取得巨大应用成功的方法。时域数值技术的一个突出优点是可以给出关于问题空间的丰富的时域信息，而且经过简单的时频变换，即可得到宽带范围的频域信息，相对频域方法显著地节约了计算量。最近几十年，是电磁场数值计算时域技术蓬勃发展的时期，各具优势和特色的新颖时域算法层出不穷。 但是到目前为止国内关于时域有限差分法中的PML 算法文献较少，其中绝大多数文献集中在综述和应用方面。而在国际的学报和杂志上对于这方面的文献非常多。时域有限差分法经过了三十年多年的高速发展之后，仍然还是计算电磁学制高点的研究热潮，而且其应用的范围和成效还在迅速的扩大和提高。本课题正是利用时域有限差分法的基础理论，利用matlab对一些实际的电磁场问题进行仿真研究。 三、论文的主攻方向、主要内容、研究方法及技术路线（1000字左右；40分） 通过对时域有限差分法理解基础之上，利用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)。
[收稿日期]20050228 [作者简介]黄书先(1962),女,1983年大学毕业,硕士,副教授,现主要从事石油勘探开发研究与科研管理工作。 基于MATLAB 平台的可视化图像处理系统 设计方法 黄书先 (长江大学科学技术处,湖北荆州434023) [摘要]以MATLAB 为工作语言和开发环境,开发了一个在M AT LAB 平台下的可视化图像处理系统, 可实现包括对测井图像在内的一般图像的精细处理,并能和用户开发的程序接口。为M AT LAB 的再开发 和可视化系统的设计作了有益的探索。 [关键词]MATLAB;图像处理;可视化GUI 界面 [中图分类号]TP 311111 [文献标识码]A [文章编号]16731409(2005)04015803 MAT LAB 的图像处理工具为自然科学各学科领域的学者、研究人员和工程师提供了一个直观的灵活的环境,用以解决复杂的图像处理问题。用MAT LAB 语言开发的图像处理算法可以在所有支持MAT LAB 的平台上共享。也可以将m 语言算法和现存的C 程序集成在一起或者将MATLAB 开发的m 语言算法和GU Is 编译为C/C++代码,供其他程序调用,或者发布为一个独立的应用程序。下面笔者提出一个基于MATLAB 平台的可视化图像处理系统设计方法,可用于包括测井图像资料在内的一般图像的精细处理。 1 系统总体设计 在MAT LAB 中有个重要的图像处理工具包[1],该工具包是由一系列支持图像处理操作的函数组成的,按功能可以分为以下几类:图像显示;图像文件输入与输出;几何操作;像素值统计;图像增强;图像识别;图像滤波;图像变换;邻域和块操作;二值图像操作;颜色映射和颜色空间转换;图像格式转换等。和其他工具包一样,用户还可以根据需要书写自己的函数,以满足特定的需要;也可以将这个工具包和信号处理工具包或小波工具包等其他工具包联合起来使用。 MAT LAB 提供了交互式的GU I 开发环境[2~4],用户只需要设置各个对象相应的属性,系统自动生成与之对应的界面,大大减少了开发的难度。本设计的MAT LAB 图像处理系统由封面、主界面和各个子功能界面组成,其框图如图1所示。程序总流程图如图2所示。 2 封面界面设计 首先用MAT LAB 编辑封面的脚本文件,生成系统的封面界面,然后再进入处理的主界面。 在设计封面时,要做好封面的总体布局,力求完美。要插入背景,首先要找到所用的函数(这在一般的书中很难找到);其次要注意设置axes 的属性。在显示背景图像时,不能用imshow ()函数,如果用此函数,背景图像只能显示在封面的一部分,不能整屏覆盖;而要用imagesc ()函数。MATLAB 提供修改文本中文字的字体,它支持华文中宋、华文彩云、华文仿宋、楷书、黑体等,功能较强。 在设计封面的过程中,最主要的技术是要解决时间的显示,如果只用MATLAB 中的clock ()函数,则显示的是静态时间,而不会显示和电脑同步的动态时间。为了解决这个问题,需做一个循环判断语句:while find (get (0,'c hildren'))==h0。成立的条件是,只要是当前窗口循环,否则停止,这样可以减少CPU 的负担。同时,要实时提取clock ()函数,可用fix (clock),使提取的时间更美观。#158#长江大学学报(自科版)2005年4月第2卷第4期/理工卷第2卷第2期 Jour nal of Yangtze University (Nat Sci Edit)Apr 12005Vol 12No 14/Sci &Eng V,Vol 12No 12
选 题 说 明 本人选做第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
p345 subplot(2,2,1) t1=0:0.1:2; y1=sin(2*pi*t1); plot(t1,y1); title('y=sin(2\pit)') 练习： subplot(2,2,2) t2=0:0.1:2; y2=[exp(-t2);exp(-2*t2);exp(-3*t2)]; plot(t2,y2) axis([0 2 -0.2 1.2]); title('y=e-t,y=e-2t,y=e-3t') 练习： subplot(2,2,3); t3=[0 1 1 2 2 3 4]; y3=[0 0 2 2 0 0 0]; plot(t3,y3); axis([0 4 -0.5 3]); title('脉冲信号') 练习： subplot(2,2,4); t4=0:0.1:2*pi; plot(sin(t4),cos(t4));
axis([-1.2 1.2 -1.2 1.2]); axis equal; title('圆') 练习： P346 x=0:0.1:20; zeta=0 y1=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt( 1-zeta^2)*x+acos(zeta)); plot(x,y1) zeta=0.3; y2=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt( 1-zeta^2)*x+acos(zeta)); hold on plot(x,y2,'r:') zeta=0.5; y3=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt( 1-zeta^2)*x+acos(zeta)); plot(x,y3,'g*') zeta=0.707; y4=1-1/sqrt(1-zeta^2)*exp(-zeta*x).*sin(sqrt( 1-zeta^2)*x+acos(zeta)); plot(x,y4,'m-') title('二阶系统曲线') legend('\zeta=0','\zeta=0.3','\zeta=0.5','\zeta=0. 707') grid on gtext('\zeta=0') gtext('\zeta=0.3') gtext('\zeta=0.5') gtext('\zeta=0.707') ginput(3) zeta = ans = 2.6037 0.9035 13.1106 2.0029 4.2166 1.0380 P347 h_fig=gcf h_axis=gca h_line1=gco h_title=get(gca,'title') h_text2=findobj(h_fig,'string','\zeta=0.3') h_fig = 1 h_axis = 151.0018 h_line1 = 1 h_title = 152.0018 h_text2 = Empty matrix: 0-by-1 set(h_line1,'linewidth',5)
利用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》期末大作业
给出程序、图、作业分析，程序需加注释。 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)
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 的两种基本绘图功能：二维平面图形和三维立体图形。 5.1 二维平面图形 5.1.1 基本图形函数 plot 是绘制二维图形的最基本函数，它是针对向量或矩阵的列来绘制曲线的。也就是说，使用plot 函数之前，必须首先定义好曲线上每一点的x 及y 坐标，常用格式为：（1）plot(x) 当x 为一向量时，以x 元素的值为纵坐标，x 的序号为横坐标值绘制 曲线。当x 为一实矩阵时，则以其序号为横坐标，按列绘制每列元素值相对于其序号的曲线， 当x 为m× n 矩阵时，就由n 条曲线。 （2）plot(x,y) 以x 元素为横坐标值，y 元素为纵坐标值绘制曲线。 （3）plot(x,y1,x,y2,…) 以公共的x 元素为横坐标值，以y1,y2,… 元素为纵坐标值绘 制多条曲线。 例5.1.1 画出一条正弦曲线和一条余弦曲线。 >> x=0:pi/10:2*pi; >> y1=sin(x); >> y2=cos(x); >> plot(x,y1,x,y2) 图5.1.1 函数plot 绘制的正弦曲线 在绘制曲线图形时，常常采用多种颜色或线型来区分不同的数据组，MATLAB 软件专 门提供了这方面的参数选项（见表5.1.1），我们只要在每个坐标后加上相关字符串，就可实 现它们的功能。 - 2 - 表5.1.1 绘图参数表 色彩字符颜色线型字符线型格式标记符号数据点形式标记符号数据点形式 y 黄- 实线. 点<小于号 m 紫：点线o 圆s 正方形 c 青-. 点划线x 叉号 d 菱形 r 红- - 虚线+ 加号h 六角星 g 绿* 星号p 五角星 b 蓝v 向下三角形 w 白^ 向上三角形 k 黑>大于号 例如，在上例中输入 >> plot(x,y1,'r+-',x,y2,'k*:') 则得图5.1.2 图5.1.2 使用不同标记的plot 函数绘制的正弦曲线 5.1.2 图形修饰 MATLAB 软件为用户提供了一些特殊的图形函数，用于修饰已经绘制好的图形。 表5.1.2 图形修饰函数表
实验六MATLAB数据可视化 一、实验目的 掌握MATLAB 二维、三维图形绘制，掌握图形属性的设置和图形修饰；掌握图像文件的读取和显示。 二、实验内容 (1) 二维图形绘制。 (2) 三维曲线和三维曲面绘制。 三、实验步骤 1．二维图形绘制 (1) 二维图形绘制主要使用函数plot。 >> clear all; >> x=linspace(0,2*pi,100); >> y1=sin(x); >> plot(x,y) >> hold on >> y2=cos(x) >> plot(x,y) >> hold off
注：hold on 用于保持图形窗口中原有的图形，hold off解除保持。 (2) 函数plot 的参数也可以是矩阵。 >> close all >> x=linspace(0,2*pi,100); >> y1=sin(x); >> y2=cos(x); >> A=[y1 ; y2]'; >> B=[x ; x]' >> plot(B,A)
(3) 选用绘图线形和颜色。>> close all >> plot(x,y1,'g+',x,y2, 'r:') >> grid on
(4) 添加文字标注。 >> title('正弦曲线和余弦曲线') >> ylabel('幅度') >> xlabel('时间') >> legend('sin(x)', 'cos(x)') >> gtext('\leftarrowsinx')
(5) 修改坐标轴范围。 >> axis equal >> axis normal >> axis([0 pi 0 1.5]) 程序如下： x=linspace(0,2*pi,100); y1=sin(x); y2=cos(x); A=[y1 ; y2]'; B=[x ; x]' plot(B,A) plot(x,y1,'g+',x,y2, 'r:') axis equal axis normal axis([0 pi 0 1.5])
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
有限元大作业程序设计 学校：天津大学 院系：建筑工程与力学学院 专业：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绘图 摘要：MATLAB 可以表达出数据的二维、三维和四维的图形。通过对图形的线型、立面、色彩、光线、视角等属性的控制，可把数据的内在特征表现得更加细腻完善。二维图形的绘制是MATLAB 语图形处理基础，也是在绝大多数数值计算中广泛应用的图形之一。为了显示三维图形，MATLAB提供了各种各样的函数。有一些函数可在三维空间中画线，而另一些可以画曲面与线格框架。 关键词：二维绘图三维绘图特殊绘图极坐标绘图柱坐标绘图球坐标绘图四维表现图 引言 MATLAB，即“矩阵实验室”，它具有强大的绘图功能。语言丰富的图形表现方法，使得数学计算结果可以方便地、多样性地实现了可视化，这是其它语言所不能比拟的。不仅能绘制几乎所有的标准图形，而且其表现形式也是丰富多样的。matlab语言不仅具有高层绘图能力，而且还具有底层绘图能力——句柄绘图方法。在面向对象的图形设计基础上，使得用户可以用来开发各专业的专用图形。本文从图形出发，详尽地介绍MATLAB的绘图功能及其用法。 一、二维绘图 二维图形的绘制是MATLAB语言图形处理的基础，它包括一般二维图形的绘制和特殊二维图形的绘制。 （一）、一般的二维图形的绘制 1．使用plot作图 【命令】 plot 【调用格式】 plot(x) %向量绘图 plot(x,y) %基本函数绘图 plot(x,y,'cs') %自定义样式基本函数绘图,c 颜色,s线型 plot(x,y,'cs','markersize',n) %自定义样式与大小格式 plot(x1,y1,x2,y2,…) %多曲线绘图格式 plot(x1,y1,'cs1',x2,y2,'cs2',…) %自定义样式多曲线绘图 plot(x1,y1,'cs1',x2,y2,'cs2',…,' markersize',n)%自定义样式与大小多曲线绘图格式 subplot(m,n,p) % 子图分割，m代表行;n代表列;p代表绘图序号此命令绘制不同的线型、点标和颜色的图形，其中cs为字符，可以代表不同的属性，包括线型、点标和颜色。 plot绘图函数属性参数
山西大同大学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