Matlab课程大作业要求
MATLAB大作业
1.设计一个程序，要求用户用键盘输入10次，输入内容可以是数字、英文、汉字、或者
其他文字或者符号。
如果输入的是数字，则把存入一个数字数组当中；如果输入的是英文的字符，则存入一个英文字符数组当中；如果是其他字符，则存入一个字符串数组中。
当10次内容都输入完毕后，先输出英文数组中的内容，然后输出字符串数组中的内容。
如果数字数组中元素数量大于3个，则由数字数组中的数字组成一个a0 X n+a1X n-1+………an=f（X），其中a0…an是数组中第一到最后一个元素，n为元素的个数，并且画出这个函数在区间[0,20]之间的曲线。
如果数字数组中的元素数量小于等于3个，则由数字数组中的数字组成一个a0 X n+a1Y n-1+anZ=0的方程，并且画出这个函数在[0,10]之间的曲线。
具体要求：
（1）画出程序的结构图
（2）画出程序的流程图
（3）每一行或者每一小段都需要标注
（4）附上完整的程序
2.按如下条件，设计一个完整的程序。
（1）有一片草原，草原上有人、狼、羊、鸡和蝗虫生活；
（2）被吃过的草，三年之后可以修复；
（3）一个人一年吃一只羊，10个人一年杀一只狼，一个人一年吃一只鸡；
（4）一头狼一年吃一只羊，一头狼三年成活一个幼崽；
（5）一只羊一年吃掉100平方米的草原；
（6）一只鸡一年吃掉1000只蝗虫；
（7）一只蝗虫一年吃掉一平方米的草原。
按如上条件，模拟出草原及其上述的生物从现在起500年间的状态。
具体要求：
（1）写出具体的解题思路
（2）列出详细的变量表，写出每一个变量的含义
（3）写出完整的程序，最好有图。
matlab大作业实验报告,《Matlab程序设计》期末实验报告-大作业2015.doc
matlab⼤作业实验报告,《Matlab程序设计》期末实验报告-⼤作业2015.doc《MATLAB程序设计》实验报告学院： 学号： 姓名：⼀、题⽬：1、(10分)已知矩阵，⽤Matlab代码实现以下要求：(1)将矩阵赋给变量A，并在屏幕上显⽰A；(2)将A按列进列逆序重排，重排后的矩阵赋给变量B，并在屏幕上显⽰B；(3)⽤reshape命令将A重排为⼀个2⾏6列矩阵并赋给变量C；(4)将A重排为⼀个列向量，将其赋给变量D，求D的平均值，在屏幕上显⽰D和它的平均值；(5)⽤命令查看变量A的维数，并显⽰运⾏结果。
2、(10分)写代码实现以下要求：构造菜单项‘Plot’，菜单项Plot有两个⼦菜单项Plot sin(选择此项后执⾏画出曲线，线型为虚线，线条颜⾊为红⾊)和Plot cos(选择此项后执⾏画出曲线 ，线型为实线，线条宽度为2)。
3、(20分)已知，实现下列操作：(1)在同⼀个图形窗⼝，同⼀坐标系下⽤不同的颜⾊和线型绘制三条曲线，并添加图例来区分三条曲线(5分)。
(2)⽤subplot命令，以⼦图的⽅式绘制三条曲线，图形排列⽅式为三⾏⼀列(5分)。
(3) 分别⽤直⽅图(bar)、棒状图(stem)和填充图(fill)绘制三条曲线，以⼦图⽅式绘制，排列⽅式为3⾏3列，共9幅⼦图(10分)。
4、(10分)⽤surf命令绘制曲⾯图形，⽤shading interp命令进⾏插值着⾊处理并添加垂直颜⾊棒。
5、(15分)⾃2011年9⽉1⽇起，我国实⾏新的个⼈所得税征收办法，起征点为3500元，请⽤If-else if-else-end结构实现⼈⼯输⼊⽉收⼊后能计算出个⼈所得税的缴纳额并显⽰⽉收⼊10000元时应缴纳的税款。
级数应纳税所得额x(元)税率备注1x<=15003%x指⽉收⼊扣除起征点3500元之后的余额；215008000045%同上6. (10分)⽤while-end循环结构计算级数和的值，输⼊n值，能计算出f的值，并显⽰结果。
2021年华科电气大四matlab大作业w
《MATLAB课程设计》班级学号姓名时间12月30日目录一．设计目 (2)二．设计要求 (2)三．设计分析 (2)四．基于根轨迹法相位滞后校正设计 (5)五．相位超前校正设计 (8)六．相位滞后校正设计 (12)七．总结 (16)参考文件 ....................................................................................... 错误!未定义书签。
反馈控制系统设计—铣床控制一、 设计目铣床是用铣刀对工件进行铣削加工机床。
铣床除能铣削平面、 沟槽、 轮齿、 螺纹和花键轴外, 还能加工比较复杂型面, 效率较刨床高, 在机械制造和修理部门得到广泛应用。
很多情况下, 铣床被用来进行高精度加工, 因为对其含有较高精度要求, 那么一个可靠反馈控制系统也不可缺乏。
此次课程设计, 我将经过利用MATLAB 软件设计铣床反馈控制系统, 使其含有良好反馈校正, 确保一个较高精度。
二、 设计要求1、 单位斜坡输入21()R s s下速度误差小于18。
2、 阶跃输入响应下超调量小于20%。
三、 设计分析传输函数简化表示下铣床闭环控制方框图以下图所表示:图1. 铣床闭环控制系统其中: Gc(s)为开环控制系统函数, G(s)为铣床本身传输函数, D(s)为外部扰动函数, N(s)为测量噪声干扰。
铣床传输函数为:2()(1)(5)G s s s s =++系统开环稳态误差()1()()()()()()1()1()G s E s R s Y s R s R s R s G s G s =-=-=++代入2()(1)(5)G s s s s =++, 21()R s s=可求得系统单位斜坡响应稳态误差:220011(1)(5)151lim ()lim ()lim lim 2(1)(5)2281(1)(5)＞→∞→→→++=====++++++ss t s s s s s s e e t sE s ss s s s ss s s s 此时速度误差不满足要求, 需要加入反馈控制系统。
控制工程基础matlab大作业要求及模板2019秋季-机设专业
控制工程基础大作业MATLAB软件应用2019年秋季学期专业名称：机械设计制造及其自动化专业班级：姓名：学号：成绩：授课教师：一、教学目的：使学生能够掌握现代工程工具MATLAB软件使用的基本方法，能够应用MATLAB软件对控制系统进行建模及性能分析。
二、内容要求：1.控制系统建模（1）确定所研究的闭环反馈控制系统，清晰表述系统的具体工作原理及参数条件；（同学们可以通过查阅相关的文献资料、生活或者工程实践中的实际案例确定自己所研究的闭环反馈控制系统）（2）绘制闭环反馈控制系统的职能方框图、函数方框图，并建立系统的传递函数形式的数学模型。
2．应用MATLAB软件进行控制系统性能分析针对所选定的闭环控制系统，应用MATLAB软件完成以下工作：（1）控制系统频域特性分析分别使用nyquist函数和bode函数绘制系统的开环奈奎斯特图和开环波德图，并附程序清单。
（2）控制系统稳定性分析判定控制系统的稳定性，并进行控制系统相对稳定性分析，计算稳定性裕量，并附程序清单。
（3）控制系统时域特性分析使用step函数绘制控制系统的单位阶跃响应曲线，分析控制系统响应的快速性指标，分析比较结构参数变化对系统性能的影响，并附程序清单。
三、作业书写注意事项：1．封皮格式按照此模板内容，不必更改，完整填写相应的个人信息；2．正文按照第二部分内容要求的顺序分项书写，给出运行结果并附上完整的编写程序清单（同时提交电子版程序）；3．本模板及要求保留，另起一页书写正文的内容成果，A4纸双面打印，左侧装订；4．杜绝抄袭，如果雷同，按照零分计；5．抽查答辩，采用百分制记分，占期末总评成绩的10%。
《控制工程基础》MATLAB大作业评分标准控制工程基础大作业MATLAB软件应用2019年秋季学期专业名称：机械设计制造及其自动化专业班级：机设177 姓名：于子雄学号： 170456 成绩：授课教师：一、教学目的：使学生能够掌握现代工程工具MATLAB软件使用的基本方法，能够应用MATLAB软件对控制系统进行建模及性能分析。
matlb课程设计作业
matlb课程设计作业一、教学目标本课程的教学目标是使学生掌握MATLAB基本语法、编程技巧以及应用方法，培养学生解决实际问题的能力。
具体目标如下：1.知识目标：（1）理解MATLAB的基本概念，如变量、数据类型、运算符等。
（2）掌握MATLAB编程的基本语法，如矩阵操作、函数定义与调用、循环结构、条件语句等。
（3）熟悉MATLAB与其他软件（如Mathematica、Python等）的接口转换。
（4）了解MATLAB在工程领域中的应用，如信号处理、控制系统、图像处理等。
2.技能目标：（1）能够运用MATLAB进行简单的数学计算、数据分析及图形绘制。
（2）具备编写MATLAB脚本文件和函数文件的能力。
（3）学会使用MATLAB解决实际问题，如编写程序实现线性方程组求解、最优化问题求解等。
（4）掌握MATLAB在实验数据处理、仿真实验等方面的应用。
3.情感态度价值观目标：（1）培养学生对科学探究的兴趣，提高其创新意识。
（2）培养学生团队协作、沟通交流的能力。
（3）培养学生具备良好的编程习惯和职业道德。
二、教学内容本课程的教学内容主要包括以下几个部分：1.MATLAB基本概念：变量、数据类型、运算符等。
2.MATLAB编程语法：矩阵操作、函数定义与调用、循环结构、条件语句等。
3.MATLAB高级应用：数组运算、图像处理、控制系统、信号处理等。
4.MATLAB与其他软件的接口转换。
5.实践项目：利用MATLAB解决实际问题，如线性方程组求解、最优化问题求解等。
三、教学方法本课程采用讲授法、案例分析法、实验法等多种教学方法相结合，以提高学生的学习兴趣和主动性。
1.讲授法：用于讲解MATLAB基本概念、语法和应用。
2.案例分析法：通过分析实际案例，使学生掌握MATLAB在各个领域的应用。
3.实验法：让学生亲自动手实践，培养其运用MATLAB解决实际问题的能力。
四、教学资源1.教材：选用《MATLAB教程》作为主要教材，辅助以相关参考书籍。
MATLAB大作业——三种函数绘制项目说明书
MATLAB大作业------三种函数的绘制项目说明书201111202013-04MATLAB项目说明书一、设计方案1.设计思路（1）、界面设计：界面空间的布局、空间大小等进行设计。
（2）、利用GUIDE的外观编辑功能在界面绘制需要的控件，并设定合适的位置和大小。
（3）、设置控件的属性：包括按钮，标签,文本和弹出菜单等。
（4）、根据要求功能对回调函数以及其他函数进行m语言编程。
（5）、图形处理：使用axes控件进行图像的生成，并给每个图像注上相应的标签。
（6）、程序的调整测试，对相关部分进行优化改进。
二、实践过程1、界面设计根据题目要求，完成对题目所涉及的信息在脑海中进行构思，首先在脑海中形成一个大概的模型，可以在草稿纸上简单地勾画，初步完成设计思路；2、利用GUIDE进行界面布局（1）、打开一个全新的空白GUI；（2）、调节画布尺寸；（3）、添加控件对象此时形成如图2-2-3：图2-2-3（4）、添加相同控件对象的时候可以使用右键 -> Duplicate 选项生成大小一致的控件，并对齐控件；3、设置控件属性（1）、Name属性设置界面名称（2）、设置按钮1的String属性为Surf。
同样设置按钮2，按钮3，静态文本框的的String属性分别为： Mesh、Contour和Select Data如图2-3-2：图2-3-2（3）、设置弹出菜单的属性，包含弹出菜单项：peaks,membrane,sinc（4）、设置弹出菜单的属性，包含弹出菜单项：peaks,membrane, sinc设置完上述内容,界面如图3-4所示.图2-3-4（5）设置Tag属性：分别对弹出菜单，Surf，Mesh和Contour按钮设置Tag属性为：plot_popup,surf,mesh,contour.这样这些控件的回调函数(Callback)分别以Tag名开头，便于识别。
4、对回调函数及其他函数编程（1）、Opening Function 界面打开时候就执行的程序。
MATLAB大作业
MATLAB大作业作业要求：（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次骰子，统计骰子各个点出现的次数，将结果以下表的形式显示，并绘制出直方图。
matlab课程设计书写规范
matlab课程设计书写规范一、教学目标本课程的目标是使学生掌握MATLAB课程设计的基本知识和技能，能够熟练地运用MATLAB进行数学建模和数据分析。
通过本课程的学习，学生应该能够：1.理解MATLAB的基本语法和操作，包括矩阵运算、函数定义、脚本编写等。
2.掌握MATLAB在数学建模中的应用，包括线性方程组求解、微分方程求解、优化问题求解等。
3.熟悉MATLAB在数据分析中的应用，包括数据导入、数据可视化、统计分析等。
4.培养学生的创新能力和团队合作精神，能够独立完成课程设计任务。
二、教学内容本课程的教学内容主要包括MATLAB的基本语法和操作、数学建模方法、数据分析方法和课程设计实践。
具体安排如下：1.MATLAB基本语法和操作：包括矩阵运算、函数定义、脚本编写等，共4课时。
2.数学建模方法：包括线性方程组求解、微分方程求解、优化问题求解等，共6课时。
3.数据分析方法：包括数据导入、数据可视化、统计分析等，共4课时。
4.课程设计实践：包括选题、方案设计、编程实现和报告撰写，共8课时。
三、教学方法为了激发学生的学习兴趣和主动性，本课程将采用多种教学方法，包括讲授法、案例分析法、实验法和讨论法。
1.讲授法：用于讲解MATLAB的基本语法和操作、数学建模方法、数据分析方法等知识点。
2.案例分析法：通过分析典型的数学建模和数据分析案例，使学生更好地理解和掌握相关方法。
3.实验法：安排上机实验，让学生亲自动手操作MATLAB，巩固所学知识。
4.讨论法：学生进行小组讨论，共同完成课程设计任务。
四、教学资源为了支持教学内容和教学方法的实施，本课程将准备以下教学资源：1.教材：选用《MATLAB课程设计》作为主要教材，辅助以相关参考书。
2.多媒体资料：制作PPT和教学视频，用于讲解和演示知识点。
3.实验设备：提供计算机实验室，让学生进行上机实验。
4.网络资源：利用校园网提供相关学习资料和在线交流平台，方便学生学习和交流。
matlab课程设计内容要求
matlab课程设计内容要求一、教学目标本课程的教学目标是让学员掌握 MATLAB 的基本语法、操作和编程技巧，能够熟练使用 MATLAB 进行数学计算、数据分析、图形绘制等操作。
通过本课程的学习，学员将能够：1.知识目标：理解 MATLAB 的基本概念、语法和编程思想，掌握MATLAB 的数据类型、矩阵运算、函数调用、脚本编写等基本操作。
2.技能目标：能够独立使用MATLAB 进行简单的数学计算、数据分析、图形绘制等任务，具备基本的 MATLAB 编程能力。
3.情感态度价值观目标：培养学员对 MATLAB 编程的兴趣和热情，提高学员的自主学习能力、团队协作能力和创新思维能力。
二、教学内容根据教学目标，本课程的教学内容主要包括以下几个部分：1.MATLAB 概述：介绍 MATLAB 的发展历程、特点和应用领域，让学员对 MATLAB 有一个整体的认识。
2.MATLAB 基本语法：讲解 MATLAB 的数据类型、变量、运算符、函数等基本语法元素。
3.MATLAB 矩阵运算：介绍矩阵的基本操作、矩阵运算符、矩阵函数等，让学员掌握 MATLAB 矩阵运算的基本技巧。
4.MATLAB 编程技巧：讲解 MATLAB 编程的常用技巧，如循环、条件语句、函数编写等。
5.MATLAB 图形绘制：介绍MATLAB 图形绘制的基本方法，包括线图、散点图、柱状图等。
6.MATLAB 应用案例：通过案例分析，让学员学会将 MATLAB 应用于实际问题中，提高解决问题的能力。
三、教学方法本课程采用讲授法、案例分析法、实验法等多种教学方法，以激发学员的学习兴趣和主动性。
1.讲授法：通过讲解 MATLAB 的基本概念、语法和操作，让学员掌握MATLAB 的基本知识。
2.案例分析法：通过分析实际案例，让学员学会将 MATLAB 应用于实际问题中，提高解决问题的能力。
3.实验法：通过上机实验，让学员亲手操作MATLAB，巩固所学知识，提高实际编程能力。
MATLAB大作业
MATLAB大作业一、作业说明MATLAB，即matrix laboratory，作为强大的数学处理软件，在化学化工领域中也有着极为广泛的应用。
而作为化工专业的学生，我认为本课学习MATLAB最直接和有效的意义在于，学生能够利用MATLAB处理实验数据和图形。
而对于MATLAB本身数学逻辑与语言的理解和掌握相信在之前所上交的“高斯消元法”和“二分法改进”两次小作业中得到体现。
因此，本次大作业我选择利用MATLAB处理化工原理实验中已经做过的离心泵实验来进行数据处理和图形制作。
其中实验原始数据皆由本人进行化工原理实验时所记录，水的密度和粘度表由文献查得。
二、编程思路整个M文件的目标是，利用现有的实验做得数据，进行实验数据的处理以及图形绘制，主要分为三个步骤完成：第一部分利用插值和拟合求实验温度下水的密度和粘度值；第二部分利用公式和矩阵变换计算离心泵的扬程He和轴功率Pzhou；第三部分作图：离心泵的扬程和轴功率随流量变化的曲线（双y轴）。
具体思路已在主程序文件中进行注释。
三、MATLAB主程序clearclose allformat compact,format shortecho off% 处理化工原理实验数据及图形制作% 第一部分求水的密度和粘度%已知条件：（1）实验温度为23.1℃% （2）网上查得15℃-30℃下水的密度和粘度表，记录为excel表% （3）此部分数据存放在原始数据表的sheet1和sheet2中%%主要方法：（1）MATLAB数据读入% （2）矩阵初等变换% （3）函数拟合% （4）函数插值%[num1]=xlsread('原始数据表');%读取excel表格中sheet1的数据x0=[num1(:,1)];%x0x表示温度值，sheet1的第1列y0=[num1(:,2)];%y0表示密度值，sheet1的第2列x1=x0';%进行转置，将列向量转换为相应的行向量y=y0';%进行转置，将列向量转换为相应的行向量%下面开始拟合密度随温度的变化函数n = length(y);ss = '';for k = 1:n;yy = num2str(y(k));ss = strcat(ss,['+''(' yy ')''*''x''^' num2str(n+1-k)]);%将数值转换为字符串格式，在MATLAB内将密度与温度的关系拟合为n次多项式endpause,ss%输出密度随温度变化函数式%%进行图像表示%x轴为温度，y轴为密度，表示在figure窗口的上半部分pause,subplot(2,1,1),line(x1,y,'linewidth',2),title('T-density'),xlabel('温度T/℃'),ylabel('密度/(kg?m^-3)')axis([15,30,995,1000])pause, x2=23.1,y1=interp1(x1,y,x2,'spline');%读取excel表格中sheet2的数据，接下来步骤与密度处理方法相同[num2]=xlsread('原始数据表',2);t0=[num2(:,1)];z0=[num2(:,2)];t1=t0';z=z0';n = length(z);rr= '';for k = 1:n;zz = num2str(z(k));rr = strcat(rr,['+''(' zz ')''*''t''^' num2str(n+1-k)]);endpause,rr%x轴为温度，y轴为粘度，表示在figure窗口的下半部分pause,subplot(2,1,2),line(t1,z,'linewidth',2),title('T-viscosity'), xlabel('温度T/℃'),ylabel('粘度/(Pa?s)')axis([15,30,0.0008 0.0012])pause,t2=23.1,z1=interp1(t1,z,t2,'spline');y1,z1%显示实验温度下流体的密度与粘度值，y1为密度，z1为粘度pause,close all%关闭图像窗口pause,clc% 第二部分计算离心泵的扬程He和轴功率Pzhou%%已知条件：（1）所需数据存放在原始数据表的sheet3中% （2）从左至右分别为“入口压强、出口压强、入口速度、出口速度、流量、电机功率”% （3）计算公式：He=pout-pin+0.2+(uout.^2-uin.^2)./(9.81*2);% Pzhou=P.*0.9;%%主要方法：（1）MATLAB数据读入% （2）矩阵初等运算与变换% （3）利用MATLAB进行实验数据处理%读取excel表格中sheet3的数据%“pin、pout、uin、uout、qv、P、He、Pzhou”分别表示为%“入口压强、出口压强、入口速度、出口速度、流量、电机功率、扬程和轴功率”[num3]=xlsread('原始数据表',3);pout0=[num3(:,2)];pin0=[num3(:,1)];uout0=[num3(:,4)];uin0=[num3(:,3)];qv0=[num3(:,5)];pout=pout0';pin=pin0';uout=uout0';uin=uin0';qv=qv0';He=pout-pin+0.2+(uout.^2-uin.^2)./(9.81*2);P0=[num3(:,6)];P=P0';Pzhou=P.*0.9;pause,clc%输出计算结果pause,HePzhouclose allpause,clcpause,%第三部分作图：离心泵的扬程和轴功率随流量变化的曲线（双y轴）%%pause,s1=He;s2=Pzhou;x5=qv;yWidth = 60;%设置两个y轴的间隔（像素）axesPosition = [110 40 200 200];%设置图像生成位置x5Limit = [min(x5) max(x5)];%设定自变量范围x5Offset=-yWidth*diff(x5Limit)/axesPosition(3);%建立图像，设置双变量与双坐标轴h1 = axes('Units','pixels','Position',axesPosition,...'Color','w','XColor','k','YColor','r',...'XLim',x5Limit,'YLim',[10 22],'NextPlot','add');h2 = axes('Units','pixels','Position',axesPosition+yWidth.*[-1 0 1 0],...'Color','none','XColor','k','YColor','m',...'XLim',x5Limit+[x5Offset 0],'YLim',[0.3 0.7],...'XTick',[],'XTickLabel',[],'NextPlot','add');xlabel(h1,'qv');ylabel(h2,'Pzhou');%给坐标轴命名pause,plot(h2,x5,s2,'*m');pause,plot(h1,x5,s1,'^r');四、工作日志density_viscosity（密度_粘度）clearclose allformat compact,format shortecho offss =+(999.126)*x^31+(999.05)*x^30+(998.97)*x^29+(998.888)*x^28+(998.802)* x^27+(998.714)*x^26+(998.623)*x^25+(998.53)*x^24+(998.433)*x^23+(998. 334)*x^22+(998.232)*x^21+(998.128)*x^20+(998.021)*x^19+(997.911)*x^18 +(997.799)*x^17+(997.685)*x^16+(997.567)*x^15+(997.448)*x^14+(997.327 )*x^13+(997.201)*x^12+(997.074)*x^11+(996.944)*x^10+(996.813)*x^9+(99 6.679)*x^8+(996.542)*x^7+(996.403)*x^6+(996.262)*x^5+(996.119)*x^4+(9 95.974)*x^3+(995.826)*x^2+(995.676)*x^1x2 =23.1000rr =+(0.0011404)*t^16+(0.0011111)*t^15+(0.0010828)*t^14+(0.0010559)*t^13+ (0.0010299)*t^12+(0.001005)*t^11+(0.000981)*t^10+(0.0009579)*t^9+(0.0009358)*t^8+(0.0009142)*t^7+(0.0008937)*t^6+(0.0008737)*t^5+(0.000854 5)*t^4+(0.000836)*t^3+(0.000818)*t^2+(0.0008007)*t^1t2 =23.1000y1 =997.5432z1 =9.3361e-04He =Columns 1 through 1021.7000 20.9069 19.9348 19.1688 18.7950 18.3412 17.9656 17.1311 16.2075 15.6581Columns 11 through 1415.1075 14.2847 13.4611 12.6543Pzhou =Columns 1 through 100.3240 0.3600 0.4050 0.4500 0.4770 0.5040 0.5310 0.5670 0.5940 0.6120Columns 11 through 140.6300 0.6390 0.6570 0.6750diary off五、附录附录1. 温度在15-30℃时，水的密度表附录2. 温度在15-30℃时，水的粘度表附录3. 离心泵实验原始数据记录表附录4. 生成图像。
matlab课程设计大作业
matlab课程设计大作业一、教学目标本课程的教学目标是使学生掌握MATLAB基本语法、编程技巧以及MATLAB 在工程计算和数据分析中的应用。
通过本课程的学习，学生将能够熟练使用MATLAB进行简单数学计算、线性方程组求解、函数图像绘制等。
1.掌握MATLAB基本语法和编程结构。
2.了解MATLAB在工程计算和数据分析中的应用。
3.熟悉MATLAB的函数库和工具箱。
4.能够使用MATLAB进行简单数学计算。
5.能够使用MATLAB求解线性方程组。
6.能够使用MATLAB绘制函数图像。
7.能够利用MATLAB进行数据分析和处理。
情感态度价值观目标：1.培养学生对计算机辅助设计的兴趣和认识。
2.培养学生团队合作和自主学习的能力。
二、教学内容本课程的教学内容主要包括MATLAB基本语法、编程技巧以及MATLAB在工程计算和数据分析中的应用。
1.MATLAB基本语法：介绍MATLAB的工作环境、基本数据类型、运算符、编程结构等。
2.MATLAB编程技巧：讲解MATLAB的函数调用、脚本编写、函数文件编写等编程技巧。
3.MATLAB在工程计算中的应用：介绍MATLAB在数值计算、线性方程组求解、图像处理等方面的应用。
4.MATLAB在数据分析中的应用：讲解MATLAB在数据采集、数据分析、数据可视化等方面的应用。
三、教学方法本课程采用讲授法、案例分析法、实验法等多种教学方法相结合的方式进行教学。
1.讲授法：通过讲解MATLAB的基本语法、编程技巧以及应用案例，使学生掌握MATLAB的基本知识和技能。
2.案例分析法：通过分析实际工程案例，使学生了解MATLAB在工程计算和数据分析中的应用。
3.实验法：安排上机实验，使学生在实际操作中巩固所学知识，提高实际编程能力。
四、教学资源本课程的教学资源包括教材、实验设备、多媒体资料等。
1.教材：选用《MATLAB教程》作为主要教材，辅助以相关参考书籍。
2.实验设备：为学生提供计算机实验室，配备有MATLAB软件的计算机。
Matlab语言及应用大作业
MATLAB语言及应用大作业姓名：xxx学号：xxxxxxxxxxxx班级：xxxx1.利用plot函数在区间[0,2π]同时绘制x=sin(t)和y=cos(t),z=sin(t)+cos(t)的图形。
要求：1)对图像x轴和y轴分别标注“时间”和“函数值”>> t=[0:0.01:2*pi];x=sin(t);y=cos(t);z=sin(t)+cos(t);plot(t,x,t,y,t,z);xlabel('时间','fontsize',20,'fontweight','bold');ylabel('函数值','fontsize',20,'fontweight','bold');2)对曲线加图例标注,图例位置自动定位（12分）。
t=[0:0.01:2*pi];x=sin(t);y=cos(t);z=sin(t)+cos(t);plot(t,x,t,y,t,z);xlabel('时间','fontsize',20,'fontweight','bold'); ylabel('函数值','fontsize',20,'fontweight','bold'); gtext('x=sin(t)');gtext('y=cos(t)');gtext('z=sin(t)+cos(t)');2.求2!+4!+6!+8!+10! （12分）s=sum(factorial(2:2:10))s =36698663．编写m 函数，计算函数值（12分）221,10002,,0≥<≤<≤<⎪⎪⎩⎪⎪⎨⎧-=x x x x x x y 首先创建.m 文件function y=funn(x) if x<0y=0; elseif x>=0&&x<1 y=x; elseif 1<=x&&x<2 y=2-x;elsey=0;end在MA TLAB命令窗口中输入>> funn(50)得：ans =4．计算下面矩阵的特征值、特征向量、迹和秩。
MATLAB大作业
M A T L A B大作业作业要求：（1）编写程序并上机实现，提交作业文档，包括打印稿（不含源程序）和电子稿（包含源程序），以班为单位交，作业提交截止时间6月24日。
（2）作业文档内容：问题描述、问题求解算法（方案）、MATLAB程序、结果分析、本课程学习体会、列出主要的参考文献。
打印稿不要求MATLAB程序，但电子稿要包含MATLAB程序。
（3）作业文档字数不限，但要求写实，写出自己的理解、收获和体会，有话则长，无话则短。
90问题五：利用MATLAB软件绘制一朵鲜花，实现一定的仿真效果。
提示：二维/三维绘图，对花瓣、花蕊、叶片、花杆等的形状和颜色进行详细设置。
第二类：插值与拟合。
（B级）问题一：有人对汽车进行了一次实验，具体过程是，在行驶过程中先加速，然后再保持匀速行驶一段时间，接着再加速，然后再保持匀速，如此交替。
注意，整个实验过程中从未减速。
在一组时间段50个时间点的速度。
（2）绘制插值图形并标注样本点。
问题二：估算矩形平板各个位置的温度。
已知平板长为5m，宽为3m，平板上3×5栅格点上的温度值为44，25，20，24，30；42，21，20，23，38；25，23，19，27，40。
（1）分别使用最近点插值、线性插值和三次样条插值进行计算。
（2）用杆图标注样本点。
（3）绘制平板温度分布图。
对a,b,c,d的值。
提示：曲线拟合并绘图分析第三类：定积分问题。
（B级）问题一：地球密度随着离中心（r=0）距离的变化而变化，不同半径处的密度如表所示，试估问题二：河道平均流量Q（m3/s）可使用速度和深度的乘积的积分来计算（河道横截面不规则），公式如下。
其中V(x)是离岸x（m）距离处的水速（m/s）,H(x)是离岸x距离处的水深（m）。
根据收集到过5（1（2（3（Q，单位是m（1（2（1（2（3）将节点1的力改为方向向上，计算这种改变对H2和V2的影响。
（4）将节点1的力撤销，而在节点1和2处施加1500N的水平外力，求节点3处垂直反作用力（V3）。
MATLAB 大作业
MATLAB 大作业请各位同学，自己完成matlab 大作业的内容。
禁止相互抄袭，如有雷同，零分计算。
大作业的格式按照实验报告的格式书写，务必标明题号，作业完成后，将生成的报告打印出来提交。
正文的字号以宋体五号字，1.5倍行距的格式打印。
请与18周五前将大作业报告由班级负责人统一收齐交给我，过期不收。
1、 试编写名为test01.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/π。
（提示：注意要产生一系列的点，这里可考虑t 的输入是向量形式，可以利用find 函数找出所需限定值的元素的位置，对其按需要赋值后，再进行绘图；其次，另外一种思路，也可考虑使用循环的形式来实现）2、 编写函数，在同一窗口的4个子图中利用plot 等语句绘制y=at 2图像，其中a=[1 2 5 10],t错误！未找到引用源。
3、 求函数32)(3-+=x x x f 在区间[-5,5]上的最大值和最小值。
4、 求解函数⎰12dx e x 的数值积分和符号积分，并比较结果。
5、 求解微分方程3|;1|2)1(002='='=''+==x x y y y x y x 的精确解和解析解，并绘制图形。
假设求解区间为[0,10] 。
6、 说说你对MATLAB 及应用这门课程学习后的体会，另外请说明在所学章节中哪一章的内容你最感兴趣，为什么？哪一章的内容你认为是没有必要学习的，为什么？如果可以选择MATLAB 的学习的内容的话，谈谈你所期望学到的知识类别的前三种。
1.实验名称 MATLAB 大作业2.实验目的(1)理解Matlab数据对象的特点；(2)掌握基本Matlab运算规则；(3)掌握Matlab帮助的使用方法；3.实验设备及条件计算机一台（带有MATLAB7.0以上的软件环境）4.实验内容及要求请各位同学，自己完成matlab大作业的内容。
matlab大作业-平衡杆小球位置控制系统-华中科技大学精选全文
可编辑修改精选全文完整版MATLAB语言、控制系统分析与设计大作业平衡杆小球位置控制系统设计与仿真专业：电气工程及其自动化班级：设计者：学号：华中科技大学电气与电子工程学院2008年1月平衡杆小球位置控制系统设计与仿真一、问题描述与实验要求A ball is placed on a beam, see figure below, where it is allowed to roll with 1 degree of freedom along the length of the beam. A lever arm is attached to the beam at one end and a servo gear at the other. As the servo gear turns by an angle theta, the lever changes the angle of the beam by alpha. When the angle is changed from the vertical position, gravity causes the ball to roll along the beam. A controller will be designed for this system so that the ball's position can be manipulated.For this problem, we will assume that the ball rolls without slipping and friction between the beam and ball is negligible. The constants and variables for this example are defined as follows:M mass of the ball 0.11 kgR radius of the ball 0.015 mD lever arm offset 0.03 mG gravitational acceleration 9.8 m/s^2L length of the beam 1.0 mJ ball's moment of inertia 9.99e-6 kgm^2R ball position coordinate Alpha beam angle coordinate Theta servo gear angleSystem EquationsThe Lagrangian equation of motion for the ball is given by the following:Linearization of this equation about the beam angle, alpha = 0, gives us the following linear approximation of the system:The equation which relates the beam angle to the angle of the gear can be approximated as linear by the equation below:Substituting this into the previous equation, we get:Design requirementsThe design criteria for this problem are: •Settling time less than 3 seconds • Overshoot less than 5%1．数学模型的建立将上面推导的简化式做拉普拉斯变换，得到： 22()()()J d m R s s mg s R Lθ+=- 变化后得到：22()()()dmgR s lJ s m s Rθ=-+ 化简为：222()()()R s mgdR s J mR ls θ=-+2．设计目标希望能精确小球的位置，即要求小球的稳态误差为零，同时希望因扰动引起的稳态误差也能为零。
传热学MATLAB温度分布大作业完整版
传热学大作业（第四章）姓名：张宝琪学号：03110608一、题目及要求1.各节点的离散化的代数方程2.源程序3.不同初值时的收敛快慢4.上下边界的热流量（λ=1W/(m℃））5.计算结果的等温线图6.计算小结题目：已知条件如下图所示：二、方程及程序(1)各温度节点的代数方程ta=(300+b+e)/4 ; tb=(200+a+c+f)/4; tc=(200+b+d+g)/4; td=(2*c+200+h)/4 te=(100+a+f+i)/4; tf=(b+e+g+j)/4; tg=(c+f+h+k)/4 ; th=(2*g+d+l)/4ti=(100+e+m+j)/4; tj=(f+i+k+n)/4; tk=(g+j+l+o)/4; tl=(2*k+h+q)/4tm=(2*i+300+n)/24; tn=(2*j+m+p+200)/24; to=(2*k+p+n+200)/24; tp=(l+o+100)/12 (2)源程序【G-S迭代程序】【方法一】函数文件为：function [y,n]=gauseidel(A,b,x0,eps)D=diag(diag(A));L=-tril(A,-1);U=-triu(A,1);G=(D-L)\U;f=(D-L)\b;y=G*x0+f;n=1;while norm(y-x0)>=epsx0=y;y=G*x0+f;n=n+1;end命令文件为：A=[4,-1,0,0,-1,0,0,0,0,0,0,0,0,0,0,0;-1,4,-1,0,0,-1,0,0,0,0,0,0,0,0,0,0;0,-1,4,-1,0,0,-1,0,0,0,0,0,0,0,0,0;0,0,-2,4,0,0,0,-1,0,0,0,0,0,0,0,0;-1,0,0,0,4,-1,0,0,-1,0,0,0,0,0,0,0;0,-1,0,0,-1,4,-1,0,0,-1,0,0,0,0,0,0;0,0,-1,0,0,-1,4,-1,0,0,-1,0,0,0,0,0;0,0,0,-1,0,0,-2,4,0,0,0,-1,0,0,0,0;0,0,0,0,-1,0,-1,0,4,0,0,0,-1,0,0,0;0,0,0,0,0,-1,0,0,-1,4,-1,0,0,-1,0,0;0,0,0,0,0,0,-1,0,0,-1,4,-1,0,0,-1,0;0,0,0,0,0,0,0,-1,0,0,-2,4,0,0,0,-1;0,0,0,0,0,0,0,0,-2,0,0,0,24,-1,0,0;0,0,0,0,0,0,0,0,0,-2,0,0,-1,24,-1,0;0,0,0,0,0,0,0,0,0,0,-2,0,0,-1,24,-1;0,0,0,0,0,0,0,0,0,0,0,-1,0,0,-1,12];b=[300,200,200,200,100,0,0,0,100,0,0,0,300,200,200,100]';[x,n]=gauseidel(A,b,[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]',1.0e-6) xx=1:1:4;yy=xx;[X,Y]=meshgrid(xx,yy);Z=reshape(x,4,4);Z=Z'contour(X,Y,Z,30)Z =139.6088 150.3312 153.0517 153.5639108.1040 108.6641 108.3119 108.1523 84.1429 67.9096 63.3793 62.4214 20.1557 15.4521 14.8744 14.7746 【方法2】>> t=zeros(5,5);t(1,1)=100;t(1,2)=100;t(1,3)=100;t(1,4)=100;t(1,5)=100;t(2,1)=200;t(3,1)=200;t(4,1)=200;t(5,1)=200;for i=1:10t(2,2)=(300+t(3,2)+t(2,3))/4 ;t(3,2)=(200+t(2,2)+t(4,2)+t(3,3))/4;t(4,2)=(200+t(3,2)+t(5,2)+t(4,3))/4;t(5,2)=(2*t(4,2)+200+t(5,3))/4;t(2,3)=(100+t(2,2)+t(3,3)+t(2,4))/4;t(3,3)=(t(3,2)+t(2,3)+t(4,3)+t(3,4))/4; t(4,3)=(t(4,2)+t(3,3)+t(5,3)+t(4,4))/4; t(5,3)=(2*t(4,3)+t(5,2)+t(5,4))/4;t(2,4)=(100+t(2,3)+t(2,5)+t(3,4))/4;t(3,4)=(t(3,3)+t(2,4)+t(4,4)+t(3,5))/4;t(4,4)=(t(4,3)+t(4,5)+t(3,4)+t(5,4))/4;t(5,4)=(2*t(4,4)+t(5,3)+t(5,5))/4;t(2,5)=(2*t(2,4)+300+t(3,5))/24;t(3,5)=(2*t(3,4)+t(2,5)+t(4,5)+200)/24;t(4,5)=(2*t(4,4)+t(3,5)+t(5,5)+200)/24;t(5,5)=(t(5,4)+t(4,5)+100)/12;t'endcontour(t',50);ans =100.0000 200.0000 200.0000 200.0000 200.0000 100.0000 136.8905 146.9674 149.8587 150.7444 100.0000 102.3012 103.2880 103.8632 104.3496 100.0000 70.6264 61.9465 59.8018 59.6008 100.0000 19.0033 14.8903 14.5393 14.5117【Jacobi迭代程序】函数文件为：function [y,n]=jacobi(A,b,x0,eps)D=diag(diag(A));L=-tril(A,-1);U=-triu(A,1);B=D\(L+U);f=D\b;y=B*x0+f;n=1;while norm(y-x0)>=epsx0=y;y=B*x0+f;n=n+1;end命令文件为：A=[4,-1,0,0,-1,0,0,0,0,0,0,0,0,0,0,0;-1,4,-1,0,0,-1,0,0,0,0,0,0,0,0,0,0; 0,-1,4,-1,0,0,-1,0,0,0,0,0,0,0,0,0; 0,0,-2,4,0,0,0,-1,0,0,0,0,0,0,0,0;-1,0,0,0,4,-1,0,0,-1,0,0,0,0,0,0,0; 0,-1,0,0,-1,4,-1,0,0,-1,0,0,0,0,0,0; 0,0,-1,0,0,-1,4,-1,0,0,-1,0,0,0,0,0;0,0,0,-1,0,0,-2,4,0,0,0,-1,0,0,0,0;0,0,0,0,-1,0,-1,0,4,0,0,0,-1,0,0,0;0,0,0,0,0,-1,0,0,-1,4,-1,0,0,-1,0,0;0,0,0,0,0,0,-1,0,0,-1,4,-1,0,0,-1,0;0,0,0,0,0,0,0,-1,0,0,-2,4,0,0,0,-1;0,0,0,0,0,0,0,0,-2,0,0,0,24,-1,0,0;0,0,0,0,0,0,0,0,0,-2,0,0,-1,24,-1,0;0,0,0,0,0,0,0,0,0,0,-2,0,0,-1,24,-1;0,0,0,0,0,0,0,0,0,0,0,-1,0,0,-1,12];b=[300,200,200,200,100,0,0,0,100,0,0,0,300,200,200,100]'; [x,n]=jacobi(A,b,[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]',1.0e-6); xx=1:1:4;yy=xx;[X,Y]=meshgrid(xx,yy);Z=reshape(x,4,4);Z=Z'contour(X,Y,Z,30)n =97Z =139.6088 150.3312 153.0517 153.5639108.1040 108.6641 108.3119 108.152384.1429 67.9096 63.3793 62.421420.1557 15.4521 14.8744 14.7746三、不同初值时的收敛快慢1、[方法1]在Gauss 迭代和Jacobi 迭代中，本程序应用的收敛条件均为norm(y-x0)>=eps ，即使前后所求误差达到e 的-6次方时，跳出循环得出结果。
matlab数字图像处理 课程设计论文 期中大作业要点
《数字图像处理》期中大作业姓名：**学号：***********序号：23南湖学院机电系2013年4月第一部分：IPT函数的应用1、读取并显示一张彩色图像，然后将其灰度化，并将灰度化后的结果存入计算机中，最后再将其二值化；2、解答：（1）程序f=imread('pic.jpg') %读入图像meimshow(f) %显示图像mea=rgb2gray(f) %调用灰度处理函数对头像进行灰度处理figure,imshow(a) %显示灰度处理后的图像b=im2bw(a) %调用二值转换函数将图像转换成二值图像figure,imshow(b) %显示灰度处理后的二值图像imwrite(a,'pic-gai.jpg') %将灰度图像写入计算机imwrite(b,'pic-gai.jpg') %将二值图像写入计算（2）运行结果图1：原图pic.jpg 图2：灰度处理后的pic.jpg图3：灰度处理加二值处理后的pic-gai.jpg（3）结果分析调用灰度系统函数对图像进行灰度化，然后进行二值处理。
灰度处理后去掉了彩色，二值处理让深灰色变成黑色浅灰色变成了白色。
2、使用亮度变换函数完成一张灰度图片的亮度调整：解答：（1）程序：f2=rgb2gray(imread('pic.jpg')) %读入图像并进行灰度处理imshow(f2) %显示灰度处理后的图像F=fft2(f2) %对图像进行DFT计算S=abs(F) %得到频谱图figure,imshow(S,[]) %显示频谱图Fc=fftshift(F) %进行图像的平移figure,imshow(abs(Fc),[]) %显示平移后图像S2=log(1+abs(Fc)) %进行对数变换figure,imshow(S2,[]) %显示对数变换后的图像（3）运行结果图1 灰度图像图2 f4直方图图3 进行规定化之后的效果g4 图4 g4直方图（3）结果分析通过调用函数fft2进行DFT计算，然后利用fftshift将变换的原点平移到频率矩阵中心，可以明显看到图像中心点的变化（图3），然后利用对数变化增强视觉效果（图4），结果变化明显，一目了然。