一、实训目的 1. 熟悉MATLAB语言的使用 2. 了解MATLAB在电子信息课程中的应用 3. 掌握MATLAB矩阵输入、运算以及MATLAB数值的运算功能 4．掌握各种数据的创建、访问、扩建及缩减 5.了解利用MATLAB计算系统响应的方法；验证信号与系统的基本概念、基本理论，掌握信号与系统的分析方法。 二、实训任务和要求 1、熟练利用MATLAB语言编程对系统S平面分析 2、掌握利用MATLAB分析系统时域的方法，熟悉系统的零输入响应，零状态响应及冲击响应的步骤。 3、MATLAB是目前国际上最流行，应用最广泛的科学与工程计算软件，它由MATLAB语言，MATLAB工作环境，MATLAB图像处理系统，MATLAB数据函数库，MATLAB 应用程序接口五大部分组成的集数值计算，图形处理，程序开发为一体的功能强大的系统．它应用于自动控制，数学计算，信号分析，计算机技术，图像信号处理，财务分析，航天工业，汽车工业，生物医学工程，语音处理和雷达工程等各行业，也是国内高校和研究部门进行许多科学研究的重要工具。 MATLAB的基本数据单位是矩阵，它的指令表达式与数学,工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完相同的事情简捷得多。MATLAB是“矩阵实验室”（MATrix LABoratoy）的缩写，它是以矩阵运算为基础的交互式程序语言，能够满足科学、工程计算和绘图的需求。与其它计算机语言相比，其特点是简洁和智能化，适应科技专业人员的思维方式和书写习惯，使得编程和调试效率大大提高。它用解释方式工作，键入持续后立即得出结果，人机交互性能好，易于调试并被科技人员所乐于接受。特别是它可适应多种平台，并且随着计算机硬软件的更新及时升级，因此MATLAB语言在国外的大学工学院中，特别是频繁进行数值计算的电子信息类学科中，已经成为每个学生都掌握的工具了。它大大提高了课程教学、解题作业、分析研究的效率。 4、系统S平面分析的基本理论 MATLAB在系统S平面分析中，一般是求系统函数的零极点分布图，画单位冲激响应和幅频响应的图形，通常按常规方法这些波形很难画出，但是应用MATLAB就可简便快捷的画出图形，使系统的分析更加便捷．例如，已知系统函数为H，利用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 心得体会 本学期通过对MATLAB的系统环境，数据的各种运算，矩阵的分析和处理，程序设计，绘图，数值计算及符号运算的学习，初步掌握了MATLAB的实用方法。通过理论课的讲解与实验课的操作，使我在短时间内学会使用MATLAB，同时，通过上机实验，对理论知识的复习巩固实践，可以自己根据例题编写设计简单的程序来实现不同的功能，绘制出比较满意的二维三维图形，在实践中找到乐趣。 MATLAB是一个实用性很强，操作相对容易，比较完善的工具软件，使用起来比较方便，通过操作可以很快看到结果，能够清晰的感觉到成功与失败，虽然课程中也会出现一些小问题，但是很喜欢这门课程。在为学习这门课前就听说了他的强大，因为现在的很多模型都是需要这些分析软件的。曾经旁听过学校数学建模的课程，当时老师用的是lingo。对那个只需要U盘携带就可以安装的小东西记忆深刻。等到学习matlab时觉得这才是真正的王道啊。 它不仅有强大的运算功能，还有强大的绘图功能，虽然学习了有一个学习，但是我对他的了解额仅仅是一点点，或许连入门都谈不上。因为我学习时了解到一个现实。就是matlab 的学习依赖有比较好的数学功底，其中我看最经常运用到的就是矩阵。我从网上了解到matlab是一门高等数学和计算机技术结合的东西，学习它必须具有相应的数学和计算机知识。然而很可惜，我的书写不是很好。每次讲到这个部分的时候就觉得听说理解无能了。特别是我今年还是大三。虽然这学期的学习的时间短暂，就算时间足够，老师也不能把所有的都讲解给我们，因为一个软件的功能需要我们自己不断的去摸索，老师也不可能知道所有。老师只是个指路人，最终的学习还是要靠自己。而且在摸索的过程中，我们能够发现和体会学习的快乐。痛并快乐着是种常态了吧。 自我感觉学习matlab与其说是学习一门软件，更不如说是学习一门语言。用一种数理的语言描述现象，揭示表象下的规律。此外，我认为matlab中的作图功能很强大，不仅简单的函数现象可以明确画出，而且一些点状物，甚至立体图也可以画出。大一上微积分的时候，老师曾经多次在课件中加入用matlab画出的图来。不论是一维二维三维等等，都能很好的画出来。只要能编写出函数式，在短短的几秒之内，他就会呈现在你眼前。另外就是图形的直观性，这是由阴影的制作的。而且可以根据需要，坐标图上加标题，坐标轴标记，文本注释级栅格等，也可以指定图线形式，比如是虚线。颜色也可以自己来定。可以在同一张图上画，也可以单个显示。 在学习的过程中，因为以前学过access中的select语言，觉得就编写这方面是有共性的，但是matlab的编程语言似乎更多更复杂一点，这是由于涉及的数学模型，数学公式更多的原因。可是今年的这门课真的是让我感到没学到什么，估计也是因为我抱着看一看的随意态度来的吧，也没有那种遇到不懂的就一定要弄懂它的决心和毅力。说什么都是借口了，无法掩饰我没有学好它的事实。事实上，我觉得今年这门课的重点并不是让我们掌握这种软件的具体用法，而是主要向我们展示如何用它去解决一些金融问题，数学问题。这点让我很郁闷，因为我不懂得原理，听起来这门课倍感吃力啊。可是嘛，年轻没有什么不可以，又有谁可以断言我接下来的生活中不能好好学习这个东西为自己的工作，学习，生活，研究兴趣带来方便呢。 从大学开学的见闻到现在学习MATLAB，感觉这是一个很好的软件，语言简便，实用性强。作为一个做新手，想要学习好这门语言，可以说还是比较难的。在我接触这门语言的这些天，除了会画几个简单的图形，其他的还是有待提高。从另一个方面也对我们大学生提出了两个要求——充实的课外基础和良好的英语基础。在现代，几乎所有好的软件都是来自国外，假如不会外语，想学好是非常难的。
信号与系统MATLAB第一次实验报告 一、实验目的 1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。 2.学会运用MATLAB表示常用连续时间信号的方法 3.观察并熟悉一些信号的波形和特性。 4.学会运用MATLAB进行连续信号时移、反折和尺度变换。 5.学会运用MATLAB进行连续时间微分、积分运算。 6.学会运用MATLAB进行连续信号相加、相乘运算。 7.学会运用MATLAB进行连续信号的奇偶分解。 二、实验任务 将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。 三、实验内容 1.MATLAB软件基本运算入门。 1). MATLAB软件的数值计算： 算数运算 向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn为结束值。 矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开； 矩阵的不同行之间必须用分号”;”或者ENTER分开。2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。
举例：计算一个函数并绘制出在对应区间上对应的值。 2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名” 2.MATLAB软件简单二维图形绘制 1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y) 2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p 表示第p个区域，表达为subplot(mnp)或者subplot(m,n,p) 3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin]) 4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’) 5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’) 6).输出：grid on 举例1：
Matlab程序设计 班级 姓名 学号
《MATLAB程序设计》作业 1、考虑如下x-y 一组实验数据： x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10] y=[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2] 分别绘出plot的原始数据、一次拟合曲线和三次拟合曲线，给出MATLAB代码和运行结果。 代码如下： x=[1,2,3,4,5,6,7,8,9,10]; y=[1.2,3,4,4,5,4.7,5,5.2,6,7.2]; plot(x,y); title('原始数据'); p=polyfit(x,y,1); q=polyval(p,x); figure,plot(x,q); title('一次拟合'); p=polyfit(x,y,2); q=polyval(p,x); figure,plot(x,q); title('二次拟合'); 运行结果如下：
1 2 3 4 5 6 7 8 9 10 12 3 4 5 6 7 8 原始数据 123 456789 102 2.5 3 3.54 4.5 55.56 6.57一次拟合 123456789 101 2 3 4 5 6 7 二次拟合 2、在[0，3π]区间，绘制y=sin(x)曲线（要求消去负半波，即(π，2π)区间内的函数值置零），求出曲线y 的平均值，以及y 的最大值及其最大值的位置。给出执行代码和运行结果。 代码如下： clear clc x=(0:0.01:3*pi); y=sin(x); plot(x,y); y1=(y>=0).*y; figure,plot(x,y1);
中南大学MATLAB程序设计实践学长有爱奉献，下载填上信息即可上交，没有下载券的自行百度。所需m文件照本文档做即可，即新建（FILE）→脚本（NEW-Sscript）→复制本文档代码→运行（会跳出保存界面，文件名默认不要修改，保存）→结果。第一题需要把数据文本文档和m文件放在一起。全部测试无误，放心使用。本文档针对做牛顿法求非线性函数题目的同学，当然第一题都一样，所有人都可以用。←记得删掉这段话 班级： ？ 学号： 姓名：
一、《MATLAB程序设计实践》Matlab基础 表示多晶体材料织构的三维取向分布函数（f＝f（φ1，φ，φ2））是一个非常复杂的函数，难以精确的用解析函数表达，通常采用离散 空间函数值来表示取向分布函数，是三维取向分布函数的一个实例。 由于数据量非常大，不便于分析，需要借助图形来分析。请你编写一 个matlab程序画出如下的几种图形来分析其取向分布特征： （1）用Slice函数给出其整体分布特征； " ~ （2）用pcolor或contour函数分别给出(φ2＝0, 5, 10, 15, 20, 25, 30, 35 … 90)切面上f分布情况(需要用到subplot函数)；
(3) 用plot函数给出沿α取向线(φ1=0~90，φ＝45，φ2＝0)的f分布情况。 (
备注：数据格式说明 解： （1）（ （2）将文件内的数据按照要求读取到矩阵f(phi1,phi,phi2)中，代码如 下： fid=fopen(''); for i=1:18 tline=fgetl(fid); end phi1=1;phi=1;phi2=1;line=0; f=zeros(19,19,19); [ while ~feof(fid) tline=fgetl(fid); data=str2num(tline); line=line+1;数据说明部分，与 作图无关此方向表示f随着 φ1从0,5,10,15, 20 …到90的变化而 变化 此方向表示f随着φ 从0,5,10,15, 20 … 到90的变化而变化 表示以下数据为φ2＝0的数据，即f（φ1，φ，0）
MATLAB程序设计语言 实 验 报 告 专业及班级：电子信息工程 姓名：王伟 学号：1107050322 日期 2013年6月20日
实验一 MATLAB 的基本使用 【一】 实验目的 1.了解MATALB 程序设计语言的基本特点，熟悉MATLAB 软件的运行环境； 2.掌握变量、函数等有关概念，掌握M 文件的创建、保存、打开的方法，初步具备将一般数学问题转化为对应计算机模型处理的能力； 3.掌握二维图形绘制的方法，并能用这些方法实现计算结果的可视化。 【二】 MATLAB 的基础知识 通过本课程的学习，应基本掌握以下的基础知识： 一. MATLAB 简介 二. MATLAB 的启动和退出 三. MATLAB 使用界面简介 四. 帮助信息的获取 五. MATLAB 的数值计算功能 六. 程序流程控制 七. M 文件 八. 函数文件 九. MATLAB 的可视化 【三】上机练习 1. 仔细预习第二部分内容，关于MATLAB 的基础知识。 2. 熟悉MATLAB 环境，将第二部分所有的例子在计算机上练习一遍 3. 已知矩阵???? ??????=??????????=123456789,987654321B A 。求A*B ，A .* B ，比较二者结果是否相同。并利用MATLAB 的内部函数求矩阵A 的大小、元素和、长度以 及最大值。 程序代码： >> A=[1 2 3;4 5 6;7 8 9]; >> B=[9 8 7;6 5 4;3 2 1]; >> A*B ans =
30 24 18 84 69 54 138 114 90 >> A.*B ans = 9 16 21 24 25 24 21 16 9 两者结果不同 >> [m,n]=size(A) m = 3 n = 3 >> b=sum(A) b = 12 15 18 >> a=length(A) a = 3 >>max(A) ans =
1.在同一个图形窗口内画出衰减震荡曲线及其包络线,取值范围是[0,6pi]. 代码： t=(0:pi/100:6*pi); y1=exp(-3*t).*cos(t/2); y2=exp(-3*t); t3=pi*(0:9)/9; y3=exp(-3*t3); plot(t,y1,'r--',t,y2,'b',t3,y3,'bo') 2.画出所表示的三维曲面。X,y 的取值范围是[-9,9]。 x=linspace(-9,1,9); y=linspace(-9,1,9); [x,y]=meshgrid(x,y); z=cos((sqrt(2*x.^2+2*y.^2))/sqrt(x.^2+y.^2)); surf(x,y,z);
3.求和当n =100 时的值。 sum=0; for i=1:100 sum=sum+1/factorial(i); end disp(sum); 1.7183 4. 求1000 个元素的随机数向量A 中大于0.5 的元素个数。 A=rand(1,1000); count=0; for i=1:length(A) if A(i)>0.5 count=count+1; end end count count = 514 5.画正态分布2 的概率密度函数曲线，产生个相应的随机数，N (1,4 m 10000 画出直方图和带正态密度曲线的直方图。将随机数的频率曲线与概率密数曲线画在一起进行比对。x=-20:20; y=normpdf(x,1,4); a=1+4*randn(1,10000); [N,h]=hist(a); y=h(1)-h(2); subplot(2,2,1);plot(x,y); subplot(2,2,2);bar(h,N/(10000*y),1); subplot(2,2,3);plot(x,y); hold on; bar(h,N/(10000*y),1);
海南大学本科生 2010—2011学年度第2学期 课程考查论文 学院（中心、所）：信息科学技术学院专业：电子信息工程研究方向：班级： 学生姓名：学生证号： 课程名称：Matlab应用基础 论文题目：学习Matlab的总结与感想 任课老师： （以上由学生填写） 教师评阅： 阅卷教师（签名）：年月日
摘要 本文从计算机语言、数学建模、网络控制系统仿真与结构化思维等方面阐述了半年来学习Matlab的心得体会与感想。由于个人知识有限，在部分细节问题的理解上可能存有偏差，还请杜老师批评指正，不吝赐教。 关键词：Matlab语言数学建模软件网络控制系统仿真结构化思维
学习Matlab 快半个学期了，虽然还有很多问题不是很清楚，但通过实践学习，我对于Matlab 总算有个整体的理解，而且每次上机操作，都会有一定的收获和感想，下面，就谈谈我个人对于Matlab 的一些看法。 （Matlab 语言） Matlab 和其它语言不一样，我这个学期学习的是C 语言，另外，对于Action Script 、HTML 、php 语言也接触过一些。C 语言主要是面向过程的，它的灵活性比较强，可根据自己的意图编辑程序，但所耗费的时间和精力比较大。例如定义变量，就分为int 、float 、char 等类型，十分麻烦，而Action Script 与php 就显得比较随意，不必纠结于哪一种类型的变量，比如，定义Var number=3，Var play=true 即可。相对于前两者而言，Matlab 则显得更为灵活与快捷，它是一门解释性语言，能自动将高级语言翻译成机器语言。比如，求t f 2=，当t=0,1,2,3,4,5时)(t f 的值。如果使用C 语言则需要定义变量，调用math 函数，还要应用for 循环、输出函数，而Matlab 则不然，只需输入t=0:5;f=2.^t ，然后回车即可。 另外，Matlab 还配有许多常用公式，操作起来十分方便，例如，想求出)(2)()(2)(3)(''''t f t f t y t y t y +=++在1)0(=y ，1)0('=y 时的零输入响应，应用dsolve 函数，只需输入x=dsolve('D2y+3*Dy+2*y=0','y(0)=1,Dy(0)=1') 回车，即得结果：x=3*exp(-t)-2*exp(-2*t)。或许，也正是Matlab 语言简洁、优化的特点，才使得它在学术界被广泛应用吧。 （数学建模） 对于数学建模而言，Matlab 是一款相当不错的建模辅助工具，因为 Matlab 中有统计函数，线性分析函数，插值函数，非线性分析函数等等这些数模必备的函数，而且，Matlab 强大的绘图功能可使很多数学演算过程变得可视化。这些对于分析问题都很有帮助。虽然我们学习的Matlab 是电子信息工程方向的，但在下个学期,，班里的大部分同学都要参加数模竞赛，所以掌握好Matlab 的各种函数模式就显得尤为重要了。
MATLAB 数学实验报告 指导老师： 班级： 小组成员： 时间：201_/_/_
Matlab 第二次实验报告 小组成员： 1 题目：实验四，MATLAB 选择结构与应用实验 目的：掌握if 选择结构与程序流程控制，重点掌握break，return , pause语句的应用。 问题：问题1：验证“哥德巴赫猜想” ，即：任何一个正偶数（n>=6）均可表示为两个质数的和。要求编制一个函数程序，输入一个正偶数，返回两个质数的和。 问题分析：由用户输入一个大于6 的偶数，由input 语句实现。由if 判断语句判断是否输入的数据符合条件。再引用质数判断函数来找出两个质数，再向屏幕输出两个质数即可。 编程：function [z1,z2]=gede(n); n=input('please input n')
if n<6 disp('data error'); return end if mod(n,2)==0 for i=2:n/2 k=0; for j=2:sqrt(i) if mod(i,j)==0 k=k+1; end end for j=2:sqrt(n-i) if mod(n-i,j)==0 k=k+1; end end if k==0 fprintf('two numbers are') fprintf('%.0f,%.0f',i,n- i) break end
end end 结果分析 如上图，用户输入了大于6的偶数返回两个质数5和31,通过 不断试验，即可验证哥德巴赫猜想。 纪录：if判断语句与for循环语句联合嵌套使用可使程序结构更加明晰，更快的解决问题。 2题目：实验四，MATLAB选择结构与应用实验 目的：用matlab联系生活实际，解决一些生活中常见的实际问 题。
目标函数 function f=fmin1(x)%x(1)是D外径/mm；x(2)是d内径/mm f=pi*(x(1).^2-x(2).^2)/4; 非线性约束函数 function[f,feq]=noncon(x) f=[6400-30*sqrt(x(1).^2+x(2).^2), 480000*x(1)-pi*(x(1).^4-x(2).^4)]; feq=[] 主程序 A=[-11]; b=[-3.5]; x0=[150140]; [x,fval,exitflag,output,lambda,grad,hessian]= fmincon('fmin1',x0,A,b,[],[],[],[],'noncon') 运算结果 x=152.5893149.0893 fval=829.2823 exitflag=1 output=iterations:3 funcCount:9 lssteplength:1 stepsize:2.5331e-010 algorithm:'medium-scale:SQP,Quasi-Newton,line-search' firstorderopt:9.8475e-010 constrviolation:0 message:[1x788char] lambda= lower:[2x1double] upper:[2x1double] eqlin:[0x1double] eqnonlin:[0x1double] ineqlin:236.9059 ineqnonlin:[2x1double] grad= 239.6867 -234.1889 hessian= 1.0e+006* 1.8964-1.8964 -1.8964 1.8964
实验一 Matlab 程序设计与M 文件 一、实验目的 1. 掌握Matlab 程序设计常用命令，如，循环、选择、暂停、显示输出、输入变量值等。 2. 掌握Matlab 的M 文件的创建于使用，包括脚本文件和函数文件。 3. 掌握常用的编程技巧。 二、实验内容 1. 输入一个百分制成绩，要求输出成绩等级A 、B 、C 、D 、E 。其中90分~100分为A ，80分~89分为B ，79分~79分为C ，60分~69分为D ，60分以下为E 。 要求： (1) 分别用if 语句和switch 语句实现。 (2) 输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。 2. 硅谷公司员工的工资计算方法如下： (1) 工作时数超过120小时者，超过部分加发15%。 (2) 工作时数低于60小时者，扣发700元。 (3) 其余按每小时84元计发。 试编程按输入的工号和该号员工的工时数，计算应发工资。 3. 根据2 222211116123n π=++++，求π的近似值。当n 分别取100、1000、10000时，结果是多少？ 要求：分别用循环结构和向量运算（使用sum 函数）来实现。 4. 考虑以下迭代公式： 1n n a x b x +=+ 其中a 、b 为正的学数。 (1) 编写程序求迭代的结果，迭代的终止条件为|x n+1-x n |≤10-5，迭代初值x 0=1.0，迭代次数不超过500次。 (2) 如果迭代过程收敛于r ，那么r 的准确值是，当(a,b)的值 取(1,1)、(8,3)、(10,0.1)时，分别对迭代结果和准确值进行比较。
5. 已知 12312311021 323 n n n n f n f n f n f f f f n ---==??==??==??=-+>? 求f 1~f 100中： (1) 最大值、最小值、各数之和。 (2) 正数、零、负数的个数。 6. 若两个连续自然数的乘积减1是素数，则称这两个边疆自然数是亲密数对，该素数是亲密素数。例如，2×3-1=5，由于5是素数，所以2和3是亲密数，5是亲密素数。求[2,50]区间内： (1) 亲密数对的对数。 (2) 与上述亲密数对对应的所有亲密素数之和。 7. 设2411()(2)0.1(3)0.01 f x x x =+-+-+，编写一个MATLAB 函数文件fx.m ，使得调用f(x)时，x 可用矩阵代入，得出的f(x)为同阶矩阵。 8. 一物理系统可用下列方程组来表示： 111211 12 220cos sin 0sin 0cos 000sin 000cos 1a m m a m g m N m N m g θθθθθθ--??????????????????=??????-??????-?????? 从键盘输入m 1、m 2和θ的值，求a 1、a 2、N 1和N 2的值。其中g 取9.8，输入θ时以角度为单位。 要求：定义一个求解线性方程组AX=B 的函数文件，然后在命令文件中调用该函数文件。
Matlab程序设计 上交作业要求： 1）电子文档：设计分析报告一份（包括系统建模、系统分析、系统设计思路、程序及其执行结果）。 2）Matlab程序：可执行程序一份（运行程序可显示、输出执行结果） 按班级统一上交。 题目一： 考虑如图所示的倒立摆系统。图中，倒立摆安装在一个小车上。这里仅考虑倒立摆在图面内运动的二维问题。倒立摆系统的参数包括：摆杆的质量（摆杆的质量在摆杆中心）、摆杆的长度、小车的质量、摆杆惯量等。 图倒立摆系统 设计一个控制系统，使得当给定任意初始条件(由干扰引起)时，最大超调量 %≤10%，调节时间ts ≤4s ，使摆返回至垂直位置，并使小车返回至参考位置(x=0)。 要求：1、建立倒立摆系统的数学模型 2、分析系统的性能指标——能控性、能观性、稳定性 3、设计状态反馈阵，使闭环极点能够达到期望的极点，这里所说的期望的极点确定 是把系统设计成具有两个主导极点，两个非主导极点，这样就可以用二阶系统的 分析方法进行参数的确定 4、用MATLAB 进行程序设计，得到设计后系统的脉冲响应、阶跃响应，绘出相应状 态变量的时间响应图。
题目二： 根据自身的课题情况，任意选择一个被控对象，按照上题所示步骤进行分析和设计，并给出仿真程序及其执行结果。 选择的被控对象在第四周上报，通过与教师协商确认，以避免选题不当与选题重复。解： 一、 第一部分单阶倒立摆系统建模 由于此问题为”单一刚性铰链、两自由度动力学问题”，因此，依据经典力学的牛顿定律即可满足要求。 如图1.1所示，设小车的质量为M，倒立摆均匀杆的质量为m，摆长为2l， O为摆摆的偏角为 ，小车的位移为x，作用在小车上的水平方向上的力为F， 1 杆的质心。 根据刚体绕定轴转动的动力学微分方程，转动惯量与角加速度乘积等于作用于刚体主动力对该轴力矩的代数和，则 1）摆杆绕其重心的转动方程为 Jθ=F y lsinθ?F x lcosθ（1-1）2）摆杆重心的水平运动可描述为 (x+lsinθ)（1-2） F x=m d2 dt 3）摆杆重心在垂直方向上的运动可描述为 (lcosθ)（1-3） F y?mg=m d2 dt 4）小车水平方向运动可描述为 F?F x=M d2x （1-4） dt 由式（1-2）和式（1-4）得 (M+m)x?+ml(cosθ?θ?sinθ?θ2)=F（1-5）由式（1-1）、式（1-2）和式（1-3）得 (J+ml2)+mlcosθ?x?=mlgsinθ（1-6）整理式（1-5）和式（1-6），得
Matlab程序设计作业 姓名: 学号: 专业:
? MATLAB 程序设计》作业 1、考虑如下x-y 一组实验数据: x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10] y 二[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2] 分别绘出plot 的原始数据、一次拟合曲线和三次拟合曲线，给出 原始曲线 MATLAB 代码和运行结果。 7 6 5 4 3 2 2 3 4 5 6 7 8 9 10
7 6.5 6 5.5 5 4.5 4 3.5 3 2.5 10 一次拟合 三次拟合
x=[1, 2, 3, 4, 5, 6, 7, 8, 9, 10]; y=[1.2, 3, 4, 4, 5, 4.7, 5, 5.2, 6, 7.2]; figure; plot(x,y) p1=polyfit(x,y,1); y1=polyval(p1,x); figure; plot(x,y1) p2=polyfit(x,y,3); y2=polyval(p2,x); figure; plot(x,y2) 2、在［0, 3n区间，绘制y二Sin(x)曲线(要求消去负半波，即(n 2n)区间内的函数值置零)，求出曲线y 的平均值，以及y 的最大值及其最大值的位置。给出执行代码和运行结果。 x=0:pi/1000:3*pi; y=Sin(x); y1=(y>=0).*y; %消去负半波figure(1); plot(x,y1, 'b' ); a=mean(y1) %求出y1 的平均值 b=max(y1) %求出y1 的最大值b, 以及最大值在矩阵中的位置; d=x(find(y1==b)) >> ex1 a = 0.4243 b = 1 d = 1.5708 7.8540 >>
1.前言 2.matlab的一些特点 3.学习matlab心得体会 4.matlab的一些资源 1.前言 我接触Matlab的时间比较长了，最开始是在大学里面的数学实验课上了解了一些，学了些基础的命令，后来参加过一次数学建模，又自学了点。而后由于所学的专业是生命科学和环境相关的东西，用到matlab的机会不多，主要是一些功能用matlab实现起来不是很方便，而且手边有现成的软件可以做到，例如图像分析，还有DNA序列分析都有现成软件等。本以为不会与其有太多交集。我下决心学习matlab是在经历几件事情之后。当时，在做硕士论文时需要对电泳图片做微生物种群的多样性分析和相似性分析，当时手头的软件只能将电泳图转化为各个泳道的灰度和位置方面的数据，而不能对数据进行分析，而能进行这样分析的软件（Bionumerics）比较贵，只为了这个用几次而买显然很不划算。无奈之下，我查了些文献，了解计算的原理后便用比较熟悉的matlab编程解决这个问题，其实这个程序比较简单--DGGE中条带Shannon多样性指数的计算，在现在看来，根本不值一提，但是在当时自我感觉还是不错的，相当有成就感了。后来在课程（数值分析，微分方程数值解）中matlab经常用到，另外在帮师姐做管理方面的数学模型时用的比较多，便自学了相关方面的知识，主要是看书，自己编程还有上网交流，这时在百度上回答了很多matlab相关的问题，并成为百度matlab技术论坛的副团长，在emuch中蒙前计算模拟区区长cenwanglai 看重，聘为计算模拟版的版主。Matlab涉及的方面非常广，下面我就自己的理解谈下matlab 一些特点和我学习matlab的一点体会，希望能对大家有点帮助，有什么不对的地方，敬请指正！ 2.matlab的一些特点 A．Matlab是一个基于矩阵运算的软件，这恐怕是众所周知的事情了，但是，真正在运用的时候（就是在编程的时候），许多人（特别是初学者）往往没有注意到这个问题，因此，for 循环（包括while循环）嵌套了十几层，这不仅是暴殄天物（没有发挥matlab所长），还浪费了你宝贵的时间，就只见左下角一直busy。 B．友好的界面，易于操作，虽然matlab一打开总看到命令行窗口，其实matlab有很多
附件2 《matlab程序设计与应用》实验指导书 山东建筑大学信息与电气工程学院
前言 一、实验目的 本课程是电气工程及其自动化、自动化、电力工程与管理专业本科生的 学科基础选修课，它在线性代数、信号分析和处理、控制系统设计和仿真等 方面有着广泛的应用。主要是学习MATLAB的语法规则、基本命令和使用环境，使学生掌握MATLAB的基本命令和基本程序设计方法，提高使用该语言 的应用能力，具有使用MATLAB语言编程和调试的能力，以便为后续多门课 程使用该语言奠定必要的基础。 通过上机实验，使学生掌握MATLAB在线帮助功能的使用、熟悉MATLAB运行环境和MATLAB语言的主要特点，掌握MA TLAB语言的基本 语法规则及基本操作命令的使用，学会M文件的建立和使用方法以及应用MATLAB实现二维和三维图形的绘制方法，具有使用MATLAB语言编程和 调试的能力。 二、实验前预习 每次实验前，学生须仔细阅读本实验指导书的相关内容，明确实验目的、要求；明确实验步骤、测试数据及需观察的现象；复习与实验内容有关的理论知识；预习仪器设备的使用方法、操作规程及注意事项；做好预习要求中提出的其它事项。 三、实验注意事项 1．实验开始前，应先检查本组电脑设备是否工作正常，matlab软件工作是否正常。 2．实验时每位同学应独立完成实验任务，避免抄袭。 3．实验后应及时将实验数据进行记录与存盘，避免因电脑故障或其它原因造成实验数据的丢失。 4．实验中严格遵循电脑操作规程，如电脑发生异常现象，应立即切断电源，
报告指导教师检查处理。 5．测量数据或观察现象要认真细致，实事求是。 6．未经许可，不得对电脑进行软件及硬件的更改操作。 7．实验结束后，实验记录交指导教师查看并认为无误后，方可关掉电脑离开。 8．爱护公物，发生仪器设备等损坏事故时，应及时报告指导教师，按有关实验管理规定处理。 9．自觉遵守学校和实验室管理的其它有关规定。 四、实验总结 每次实验后，应对实验进行总结，即实验数据进行整理，绘制波形和图表，分析实验现象，撰写实验报告。实验报告除写明实验名称、日期、实验者姓名、同组实验者姓名外，还包括： 1．实验目的； 2．实验内容 3．实验流程图； 4．程序命令清单； 5．运行结果； 6．实验的收获与体会； 7．回答每项实验的有关问答题。
2.Matlab 编程思想——向量化编程 C 语言的基本元素是单数值（比如单变量或数组里的元素），再加上 其结构化的特点，决定了通常 C 语言程序大都充斥着大量罗嗦的单变量循环和判断语句（注 1）。 而 Matlab 是以向量、矩阵为基本元素的，所以要编写真正的 Matlab 程序必须抛弃【 C 语言那种“单数值、元素化”考虑问题】的思路，转以向量、矩阵为最小单位来考虑问题。也就是说， Matlab 的编程思想是——向量化编程，即面向向量或矩阵。 这样做的好处，至少有两个： （1）代码大大简化，易编程、清晰可读性强； 这样的代码才叫Matlab 代码，否则只能是不伦不类的代码：C 不 C，Matlab 不 Matlab. （2）执行效率也更高； 这是次要的，随着 Matlab 对循环机制的优化，速度差异已不再 那么明显，关键是（ 1）
下面针对 Matlab 中常见的向量化处理问题方法，举例加以说明：（一）整体操作“大块数据” Matlab 为同型的数据块（向量或矩阵）的整体做运算提供了【点 运算】，这里我借用《线性代数》里的说法，矩阵行数、列数相同 称为同型矩阵， Matlab 里矩阵可能不止 2 维。 比如， A.*C 和 A./C 表示 A 与 C的对应位置的各元素做* 和 /运算得到与它们同型的一个新矩阵。 例1.物理实验利用测得的电压电流具体数据，验证欧姆定律R=U/I. 代码 1（C语言风格） U = [0.89, 1.20, 3.09, 4.27, 3.62, 7.71, 8.99, 7.92, 9.70, 10.41]; I = [0.028,0.040,0.100,0.145,0.118,0.258,0.299,0.257, 0.308, 0.345]; L = length(U); S=0; for k = 1:L R(k) = U(k)/I(k); S = S + R(k); end R=S/L 运行结果： R = 30.5247
有限元大作业程序设计 学校：天津大学 院系：建筑工程与力学学院 专业：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学习心得体会（精选3篇） matlab学习心得体会一:matlab学习心得matlab中有丰富的图形处理能力,提供了绘制各种图形、图像数据的函数。他提供了一组绘制二维和三维曲线的函数,他们还可以对图形进行旋转、缩放等操作。matlab内部还包含丰富的数学函数和数据类型,使用方便且功能非常强大。 本学期通过对matlab的系统环境,数据的各种运算,矩阵的分析和处理,程序设计,绘图,数值计算及符号运算的学习,初步掌握了matlab的实用方法。通过理论课的讲解与实验课的操作,使我在短时间内学会使用matlab,同时,通过上机实验,对理论知识的复习巩固实践,可以自己根据例题编写设计简单的程序来实现不同的功能,绘制出比较满意的二维三维图形,在实践中找到乐趣。 matlab是一个实用性很强,操作相对容易,比较完善的工具软件,使用起来比较方便,通过操作可以很快看到结果,能够清晰的感觉到成功与失败,虽然课程中也会出现一些小问题,但是很喜欢这门课程。 matlab学习心得体会二:matlab学习心得(463字) 学习matlab是听说它是一个功能强大的数学软件,但是正被微积分的计算缠身,听说有一个高级的计算器当然高兴,以后可以偷懒了,当然现在不能偷懒。听说关于自动化的计算特别复杂,如果有一种软件能帮忙解题,那是一种极大的解脱,有益于缩短研究时间。目前我只知道有三种数学软件,都是国外的,没有国内的,差距挺大的。matlab学起来挺顺手的,比c语言简单。但是深入学习的时候却困难重重,因为很多知识都没有学习,就算知道那些函数,也没有什么用处。老师布置的作业难度大,写一篇实验,大一什么都不会,写一篇这种论文谈何容易。最多也就会一些数值计算、符号计算、简单绘图,根本不会什么实验。 学习matlab体会最多的是这个软件的功能强大,好多数学题都被轻易的解出。但是有一点遗憾,不知是我不会用,还是它没个功能,已知空间的电荷分布,求空间的电场分布。其中电场分布是无法用函数表达式表示。我知道计算机肯定可以实现,但是这个软件能不能实现就不知道了,我看过许多资料,但是在这方面没有提到相关信息。 总之,这个软件功能强大,不知什么时候国内才有类似的软件。 matlab学习心得体会三:学习matlab的心得(817字) 这是我在学习的过程中的一些技巧,或许对你有帮助,可能字数不你能满足你的要求,但是绝对是精华。