matlab的程序设计基础
matlab优化设计
MATLAB优化设计 学院：机电学院 专业：机械设计制造及其自动化 班级：072&&&-** 学号：20131****** 姓名：大禹 指导老师：祯 2015年10月25日
题目 1 1、求解如下最优化问题 步骤一：对已有的数学模型matlab 编程 1. 编写.m 文件并保存： h=[2 ,-2;-2, 4]; %实对称矩阵 f=[-2;-6]; %列向量 a=[1, 1;-1, 2]; %对应维数矩阵 b=[2;2]; %列向量 lb=zeros(2, 1); [x,value]=quadprog(h, f, a ,b ,[] ,[], lb) 2. 运行.m 文件结果如图1.0所示： subject to 2 21≤+x x 22-21≤+x x 0 21≥x x ，2 2 2121212262)(m in x x x x x x x f +-+--=
图1.0题目一文件运行结果 步骤二：matlab运行结果分析阶段 由图1.0知，当x1=0.8，x2=1.2时，min f (x)= -7.2。 题目 2 2、某农场拟修建一批半球壳顶的圆筒形谷仓，计划每座谷仓容积为300立方米，圆筒半径不得超过3米，高度不得超过10米。半球壳顶的建筑造价为每平方米150元，圆筒仓壁的造价为每平方米120元，地坪造价为每平方米50元，求造价最小的谷仓尺寸为多少？
步骤一：题目分析阶段 设：圆筒的半径为R,圆筒的高度为H 。 谷仓的容积为300立方米，可得： 3003 232=+R H R ππ 圆筒高度不得超过10米，可得： 100≤≤H 圆筒半径不得超过3米，可得： 30≤≤R 当造价最小时： 2225021202150),(m in R H R R H R f πππ+?+?= 步骤二：数学模型建立阶段 2 225021202150),(m in R H R R H R f πππ+?+?=
matlab课程设计题目
课题一： 连续时间信号和系统时域分析及MATLAB实现 课题要求： 深入研究连续时间信号和系统时域分析的理论知识。利用MATLAB强大的图形处理功能、符号运算功能以及数值计算功能，实现连续时间信号和系统时域分析的仿真波形。 课题内容： 一、用MATLAB实现常用连续时间信号的时域波形（通过改变参数，分析其时域特性）。 1、单位阶跃信号， 2、单位冲激信号， 3、正弦信号， 4、实指数信号， 5、虚指数信号， 6、复指数信号。 二、用MATLAB实现信号的时域运算 1、相加， 2、相乘， 3、数乘， 4、微分， 5、积分 三、用MATLAB实现信号的时域变换（参数变化，分析波形变化） 1、反转， 2、使移（超时，延时）， 3、展缩， 4、倒相， 5、综合变化 四、用MATLAB实现信号简单的时域分解 1、信号的交直流分解， 2、信号的奇偶分解 五、用MATLAB实现连续时间系统的卷积积分的仿真波形 给出几个典型例子，对每个例子，要求画出对应波形。 六、用MATLAB实现连续时间系统的冲激响应、阶跃响应的仿真波形。 给出几个典型例子，四种调用格式。 七、利用MATLAB实现连续时间系统对正弦信号、实指数信号的零状态响应的仿真波形。 给出几个典型例子，要求可以改变激励的参数，分析波形的变化。 课题二： 离散时间信号和系统时域分析及MATLAB实现。 课题要求： 深入研究离散时间信号和系统时域分析的理论知识。利用MATLAB强大的图
形处理功能、符号运算功能以及数值计算功能，实现离散时间信号和系统时域分析的仿真波形。 课题内容： 一、用MATLAB绘制常用信号的时域波形（通过改变参数分析其时域特性） 1、单位序列， 2、单位阶跃序列， 3、正弦序列， 4、离散时间实指数序列， 5、离散时间虚指数序列， 6、离散时间复指数序列。 二、用MATLAB实现信号的时域运算 1、相加， 2、相乘， 3、数乘。 三、用MATLAB实现信号的时域变换（参数变化，分析波形的变化） 1、反转， 2、时移（超时，延时）， 3、展缩， 4、倒相。 四、用MATLAB实现离散时间系统卷积和仿真波形 给出几个典型例子，对每个例子要求画出e(k)，h(k)，e(i)，h(i)，h(-i)，Rzs(k)波形。 五、用MATLAB实现离散时间系统的单位响应，阶跃响应的仿真波形 给出几个典型例子，四中调用格式。 六、用MATLAB实现离散时间系统对实指数序列信号的零状态响应的仿真波形 给出几个典型例子，要求可以改变激励的参数，分析波形的变化。 课题三： 连续时间信号傅里叶级数分析及MATLAB实现。 课题要求： 深入研究连续时间信号傅里叶级数分析的理论知识，利用MATLAB强大的图形处理功能，符号运算功能以及数值计算功能，实现连续时间周期信号频域分析的仿真波形。 课题内容： 一、用MATLAB实现周期信号的傅里叶级数分解与综合 以周期矩形波信号为例，绘出包含不同谐波次数的合成波形，观察合成波形与原矩形 波形之间的关系及吉布斯现象。
matlab滤波器设计
长安大学 数字信号处理综合设 计 专业_______电子信息工程_______ 班级__24030602___________ 姓名_______张舒_______ 学号2403060203 指导教师陈玲 日期_______2008-12-27________
一、课程设计目的： 1. 进一步理解数字信号处理的基本概念、基本理论和基本方法； 2.熟悉在Windows环境下语音信号采集的方法； 3.学会用MATLAB软件对信号进行分析和处理； 4.综合运用数字信号处理理论知识，掌握用MATLAB软件设计FIR和IIR数字滤波器的方法； 5. 提高依据所学知识及查阅的课外资料来分析问题解决问题的能力。 二、课程设计内容： 1.语音信号的采集 利用windows下的录音机录制一段自己的话音，时间控制在1秒左右；并对语音信号进行采样，理解采样频率、采样位数等概念。 2.语音信号的频谱分析 利用函数fft对采样后语音信号进行快速傅里叶变换，得到信号的频谱特性。 3.设计数字滤波器 采用窗函数法和双线性变换法设计下列要求的三种滤波器，根据语音信号的特点给出有关滤波器的性能指标： 1）低通滤波器性能指标，fp=1000Hz，fc=1200Hz， As=100dB，Ap=1dB； 2）高通滤波器性能指标，fc=4800Hz，fp=5000Hz ，As=100dB，Ap=1dB； 3）带通滤波器性能指标，fp1=1200Hz，fp2=3000Hz，fc1=1000Hz，fc2=3200Hz，As=100dB，Ap=1dB。 4.对语音信号进行滤波 比较用两种方法设计的各滤波器的性能，然后用性能好的滤波器分别对采集的语音信号进行滤波；并比较滤波前后语音信号的波形及频谱，分析信号的变化。 5.回放语音信号，感觉滤波前后的声音变化。 三、实验原理 （一）基于双线性Z变换法的IIR数字滤波器设计 由于的频率映射关系是根据推导的，所以使jΩ轴每隔2π/Ts便映射到单位圆上一周，利用冲激响应不变法设计数字滤波器时可能会导致上述的频域混叠现象。为了克服这一问题，需要找到由s平面到z平面的另外的映射关系，这种关系应保证： 1) s平面的整个jΩ轴仅映射为z平面单位圆上的一周； 2) 若G(s)是稳定的，由G(s)映射得到的H(z)也应该是稳定的； 3) 这种映射是可逆的，既能由G(s)得到H(z)，也能由H(z)得到G(s)； 4) 如果G(j0)=1，那么。 双线性Z变换满足以上4个条件的映射关系，其变换公式为
MATLAB课设报告
课程设计任务书 学生姓名：董航专业班级：电信1006班 指导教师：阙大顺，李景松工作单位：信息工程学院 课程设计名称：Matlab应用课程设计 课程设计题目：Matlab运算与应用设计5 初始条件： 1.Matlab6.5以上版本软件； 2.课程设计辅导资料：“Matlab语言基础及使用入门”、“Matlab及在电子信息课程中的应 用”、线性代数及相关书籍等； 3.先修课程：高等数学、线性代数、电路、Matlab应用实践及信号处理类相关课程等。 要求完成的主要任务：（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求） 1.课程设计内容：根据指导老师给定的7套题目，按规定选择其中1套完成； 2.本课程设计统一技术要求：研读辅导资料对应章节，对选定的设计题目进行理论分析， 针对具体设计部分的原理分析、建模、必要的推导和可行性分析，画出程序设计框图，编写程序代码（含注释），上机调试运行程序，记录实验结果（含计算结果和图表），并对实验结果进行分析和总结。具体设计要求包括： ①初步了解Matlab、熟悉Matlab界面、进行简单操作； ②MATLAB的数值计算：创建矩阵矩阵运算、多项式运算、线性方程组、数值统计； ③基本绘图函数：plot, plot3, mesh, surf等，要求掌握以上绘图函数的用法、简单图形 标注、简单颜色设定等； ④使用文本编辑器编辑m文件，函数调用； ⑤能进行简单的信号处理Matlab编程； ⑥按要求参加课程设计实验演示和答辩等。 3.课程设计说明书按学校“课程设计工作规范”中的“统一书写格式”撰写，具体包括： ①目录； ②与设计题目相关的理论分析、归纳和总结； ③与设计内容相关的原理分析、建模、推导、可行性分析； ④程序设计框图、程序代码（含注释）、程序运行结果和图表、实验结果分析和总结； ⑤课程设计的心得体会（至少500字）； ⑥参考文献（不少于5篇）； ⑦其它必要内容等。 时间安排：1.5周（分散进行） 参考文献： [1](美)穆尔,高会生,刘童娜,李聪聪．MA TLAB实用教程(第二版) . 电子工业出版社,2010. [2]王正林,刘明．精通MATLAB(升级版) .电子工业出版社,2011. [3]陈杰. MA TLAB宝典(第3版) . 电子工业出版社,2011. [4]刘保柱,苏彦华,张宏林. MATLAB 7.0从入门到精通(修订版) . 人民邮电出版社,2010. 指导教师签名：年月日 系主任（或责任教师）签名：年月日
基于matlab的数字滤波器设计
淮北煤炭师范学院 2009届学士学位论文 基于MA TLAB的数字滤波器设计 学院、专业物理与电子信息学院 电子信息科学与技术 研究方向基于MATLAB的数字滤波器设计 学生姓名耿博 学号200513432024 指导教师姓名邹锋 指导教师职称讲师 2009 年4 月18
基于MATLAB的数字滤波器设计 耿博 （淮北煤炭师范学院物理与电子信息学院235000） 摘要随着信息时代和数字世界的到来，数字信号处理已成为今一门极其重要的学科和技术领域。数字信号处理在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理应用中，数字滤波器十分重要并已获得广泛应用。 数字滤波是数字信号处理的重要内容，数字滤波器可分为IIR和FIR两大类。对于IIR数字滤波器的设计，需要借助模拟原型滤波器，再将模拟滤波器转化为数字滤波器，文中采用的设计方法是脉冲响应不变法、双向性变换法和完全函数设计法；对于FIR数字滤波器的设计，可以根据所给定的频率特性直接设计，文中采用的设计方法是窗函数法。本文根据IIR滤波器和FIR滤波器的特点，在MATLAB坏境下分别用双线性变换法设计IIR和用窗函数设计FIR数字滤波器，并对采集的语音信号进行分析，最后给出了IIR和FIR对语音滤波的效果。 关键词数字滤波器；IIR ；FIR ；MATLAB
The Design of Digital Filter based on MATLAB Geng Bo School of Physics and Electronics Information, Huaibei Coal Industry Teachers? College, 235000 ABSTRACT Along with the information age and the digital world arrival, the digital signal processing has become a now extremely important discipline and the area of technology.The digital signal processing in the correspondence, the multitudinous domains the pronunciation such as the image, the automatic control, the radar, the military, the aerospace, the medical service and the domestic electric appliances and so on have obtained the widespread application.In the digital signal processing application, the digital filter are extremely important and have obtained the widespread application. The digital filter are the digital signal processing important content, the digital filter may divide into IIR and the FIR two main kinds. As for the IIR digital filter design, we need the help of analog prototype filter, and then transform analog filter into digital filter. In the paper we use the design of the pulse response invariable method, the bilinear method and full function design; as for the FIR filter, we can design it directly based on the giving frequency, in the paper it uses the design of the window function.This article according to the IIR filter and the FIR filter characteristic, uses the bilinearity method of transformation under the MATLAB bad boundary to design IIR and to design the FIR numeral filter separately with the window box number, and carries on the analysis to the gathering pronunciation signal, and finally gives IIR and FIR to the pronunciation filter effect. Keywords Digtial Filter;IIR;FIR;MATLAB
基于MATLAB的优化设计
基于MATLAB的曲柄摇杆机构优化设计 1.问题的提出 根据机械的用途和性能要求的不同，对连杆机构设计的要求是多种多样的，但这些设计要求可归纳为以下三种问题：（1）满足预定的运动规律要求；（2）满足预定的连杆位置要求；（3）满足预定的轨迹要求。在在第一个问题 里按照期望函数设计的思想，要求曲柄摇杆机构的曲柄与摇杆转角之间按照φ=f(?)（称为期望函数）的关系实现运动，由于机构的待定参数较少，故一 般不能准确实现该期望函数，设实际的函数为φ=F(?)（称为再现函数），而再 现函数一般是与期望函数不一致的，因此在设计时应使机构再现函数φ=F(?) 尽可能逼近所要求的期望函数φ=f(?)。这时需按机械优化设计方法来设计曲 柄连杆，建立优化数学模型，研究并提出其优化求解算法，并应用于优化模型的求解，求解得到更优的设计参数。 2.曲柄摇杆机构的设计 在图1所示的曲柄摇杆机构中，l1、l2、l3、l4分别是曲柄AB、连杆BC、摇杆CD和机架AD的长度。这里规定?0为摇杆在右极限位置φ0时的曲柄起始 位置角，它们由l1、l2、l3和l4确定。 图1曲柄摇杆机构简图 设计时，可在给定最大和最小传动角的前提下，当曲柄从?0转到?0+90?时，要求摇杆的输出角最优地实现一个给定的运动规律f(?)。这里假设要求： (?-?0)2（1）φE=f(?)=φ0+2 3π
s=30；qb=1；jj=5；fx=0； fa0=acos（（（qb+x（1））^2-x（2）^2+jj^2）/（2*（qb+x（1））*jj））； %曲柄初始角 pu0=acos（（（qb+x（1））^2-x（2）^2-jj^2）/（2*x（2）*jj））；%摇杆初始角for i=1:s fai=fa0+0.5*pi*i/s； pui=pu0+2*（fai-fa0）^2?（3*pi）； ri=sqrt（qb^2+jj^2-2*qb*jj*cos（fai））； alfi=acos（（ri^2+x（2）^2-x（1）^2）/（2*ri*x（2）））； bati=acos（（ri^2+jj^2-qb^2）（/2*ri*jj））； if fai>0&fai<=pi psi=pi-alfi-bati; elseif fai>pi&fai<=2*pi psi=pi-alfi+bati； end fx=fx+（pui-psi）^2； end f=fx； （2）编写非线性约束函数M文件confun.m function［c，ceq］=confun（x）； qb=1；jj=5；m=45*pi/180；n=135*pi/180； c（1）=x（1）^2+x（2）^2-（jj-qb）^2-2*x（1）*x（2）*cos（m）； %最小传动角约束c（2）=-x（1）^2-x（2）^2+（jj+qb）^2+2*x（1）*x（2）*cos（n）； %最大传动角约束ceq=［］； （3）在MATLAB命令窗口调用优化程序 x0=［6；4］； lb=［1；1］； ub=［］； %线性不等式约束 a=［-1-1；1-1；-11］；b=［-6；4；4］；［x，fn］=fmincon（@optimfun， x0，a，b，［］，［］，lb，ub，@confun）； （4）运行结果
Matlab课程设计报告
自控系统仿真软件课程设计报告 MATLAB 设计题目：牛顿摆球 姓名： 学号： 院系： 班级：1203 指导教师： 2014年12月20日
一．课程设计目的 1、熟悉课程设计的基本流程； 2、掌握MATLAB语法结构及调试方法； 3、熟悉MATLAB函数调用，熟练二维画图; 4、掌握MATLAB语言在控制方面的运用； 5、学会用MATLAB进行基本仿真； 6、掌握MATLAB编程技巧，提高编程水平。 二．系统分析 1.题目的描述： （1）牛顿摆球原理描述 五个质量相同的球体由吊绳固定，彼此紧密排列。当摆动最右侧的球并在回摆时碰撞紧密排列的另外四个球，最左边的球将被弹出，并仅有最左边的球被弹出。当然此过程也是可逆的，当摆动最左侧的球撞击其它球时，最右侧的球会被弹出。当最右侧的两个球同时摆动并撞击其他球时，最左侧的两个球会被弹出。同理相反方向同样可行，并适用于更多的球。 为了更接近现实，在这里我将考虑重力及空气阻力的影响，摆球将不会永无止境的运动下去，由于外界因素的影响，摆球运动一段时间后将回归静止状态。（2）通过MATLAB动画程序制作软件，实现下述过程 当运行程序时，把最右边的小球拉到一定的高度放下，让其碰撞其余四个小球，仅让最左边的小球被弹出，当最左边小球回摆碰撞其它球时，最右边小球又被弹出，如此循环。由于是非理想条件下，摆球的摆动幅度会随摆动次数的增加越来越小，直到静止。 时间停顿两秒，把右边两小球一起拉到一定高度放下，让其碰撞其余三个球，同样仅让左边两球被弹出，当球回摆再次碰撞时，最右边两球又被同时弹出，如此循环，因为外界因素的影响，最终五个球都会静止下来。 （3）整个实验看似简单，但要在MATLAB上完成这样一个动画过程，还是需要下点功夫，克服困难的。经过自己的努力，终于实现了整个过程，这也是一种不小的收获。 2.设计要求： （1）能够实现有阻尼摆动，即摆幅随摆动次数增加越来越小，直到静止。（2）能够让摆球弧线摆动。 三．系统设计 1.系统设计过程 （1）通过函数axis建立坐标系 （2）在坐标系范围内通过函数line画各个支架 （3）通过函数title添加标题“动量守恒实验”、函数text添加标注“牛顿摆球” （4）通过函数line画出五个球，并设定其初始位置，颜色,大小,线条的擦拭方式
数字滤波器的MATLAB设计与DSP上的实现
数字滤波器的MAT LAB设计与 DSP上的实现 数字滤波器的MATLAB 设计与DSP上的实现 公文易文秘资源网佚名2007-11-15 11:56:42我要投稿添加到百度搜藏 摘要：以窗函数法设计线性相位FIR数字滤波器为例，介绍用MATLAB工具软件设计数字滤波器的方法和在定点DSP上的实现。实现时，先在CCS5000仿真开发，然后将程序加载到TMS320VC5409评估板上实时运行，结果实现了目标要求。文中还讨论了定标、误差、循环寻址等在DSP上实现的关键问题。关键词 摘要：以窗函数法设计线性相位 FIR数字滤波器为例，介绍用 MATLAB工具软件设计数字滤波器的方法和在定点DSP上的实现。实现时，先在 CCS5000仿真开发，然后将程序加载到 TMS320VC5 409评估板上实时运行，结果实现了目标要求。文中还讨论了定标、误差、循环寻址等在DSP上实 现的关键问题。 关键词：数字滤波器MATLAB DSP 引言 随着信息时代和数字世界的到来，数字信号处理已成为今一门极其重要的学科和技术领域。数字信号处理在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应
用。在数字信号处理应用中，数字滤波器十分重要并已获得广泛应用。 1数字滤波器的设计 1.1数字滤波器设计的基本步骤 数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应（IIR ）滤波器和有限长冲激响应（FIR ）滤波器。IIR滤波器的特征是，具有无限持续时间冲激响应。种滤波器一般需要用递归模型来实现，因而有时也称之为递归滤波器。FIR滤波器的冲激响应只能延续一定时间， 在工程实际中可以采用递归的方式实现，也可以采用非递归的方式实现。数字滤波器的设计方法有多种，如双线性变换法、窗函数设计法、插值逼近法和Chebyshev逼近法等等。随着 MATLAB软件尤 其是MATLAB的信号处理工作箱的不断完善，不仅数字滤波器的计算机辅助设计有了可能，而且还可以使设计达到最优化。 数字滤波器设计的基本步骤如下： (1确定指标 在设计一个滤波器之前，必须首先根据工程实际的需要确定滤波器的技术指标。在很多实际应用中，数字滤波器常常被用来实现选频操作。因此，指标的形式一般在频域中给岀幅度和相位响应。幅度指标主要以两种方式给岀。第一种是绝对指标。它提供对幅度响应函数的要求，一般应用于FI R滤波器的设计。第二种指标是相对指标。它以分贝值的形式给岀要求。在工程实际中，这种指标最受欢迎。对于相位响应指标形式，通常希望系统在通频带中人有线性相位。运用线性相位响应指标进行滤波器设计具有如下优点：①只包含实数算法，不涉及复数运算；②不存在延迟失真，只有固定数量的延迟；③长度为N的滤波器(阶数为N-1)，计算量为N/2数量级。因此，本文中滤波器的设计就以线性相位FIR滤波器的设计为例。 (2)逼近
matlab课程设计拟定题目
第一类：单位转换 1.长度单位换算的设计与实现 2.面积单位换算的设计与实现 3.体积单位换算的设计与实现 4.容积单位换算的设计与实现 5.质量单位换算的设计与实现 6.时间单位换算的设计与实现 7.温度单位换算的设计与实现 7.压强单位换算的设计与实现 8.角度单位换算的设计与实现 8.功率单位换算的设计与实现 第二类：曲线绘制 1.直线的自动绘制和相关计算 2.椭圆的自动绘制和相关计算 3.双曲线的自动绘制和相关计算 4.抛物线的自动绘制和相关计算 5.心脏线的自动绘制和相关计算 6.渐开线的自动绘制和相关计算 7.滚圆线的自动绘制和相关计算 8.三叶玫瑰线的自动绘制和相关计算9.四叶玫瑰线的自动绘制和相关计 10.阿基米德螺线的自动绘制和相关计算第三类：曲面绘制 1.球面的自动绘制和相关计算 2.椭球面的自动绘制和相关计算 3.单叶双曲面的自动绘制和相关计算 4.双叶双曲面的自动绘制和相关计算 5.抛物面的自动绘制和相关计算 6.双曲抛物面的自动绘制和相关计算 7.双曲柱面的自动绘制和相关计算 8.椭圆柱面的自动绘制和相关计算 9.抛物柱面的自动绘制和相关计算 10.圆锥面的自动绘制和相关计算 第四类：线性回归 1.男士身高体重相关计算经验公式 2.女士身高体重相关计算经验公式 3.男士胖瘦等级的确定 4.女士胖瘦等级的确定 5.男士身高脚长相关计算经验公式 6.女士身高脚长相关计算经验公式 7.父子身高相关性研究 8.母子身高相关性研究 9.父女身高相关性研究 10.母女身高相关性研究 第五类：学习成绩 1.期末总评自动计算的设计与实现 2.成绩等级自动评定的设计与实现 3.成绩分段自动统计的设计与实现 4.成绩分布折线自动绘制的设计与实现 5.成绩自动统计分析的设计与实现 6.试卷分布自动分析的设计与实现 7.试卷难度自动分析的设计与实现 8.考试成绩名次自动生成的设计与实现
机械优化设计MATLAB程序文件
机械优化设计作业1.用二次插值法求函数()()()22 ?极小值，精度e=0.01。 t t =t 1- + 在MATLAB的M文件编辑器中编写的M文件，如下： f=inline('(t+1)*(t-2)^2','t') a=0;b=3;epsilon=0.01; t1=a;f1=f(t1); t3=b;f3=f(t3); t2=0.5*(t1+t3);f2=f(t2); c1=(f3-f1)/(t3-t1); c2=((f2-f1)/(t2-t1)-c1)/(t2-t3); t4=0.5*(t1+t3-c1/c2);f4=f(t4); k=0; while(abs(t4-t2)>=epsilon) if t2
运行结果如下： 迭代计算k= 7 极小点坐标t= 2 函数值f=0.0001 2.用黄金分割法求函数()32321+-=t t t ?的极小值，精度e=0.01。 在MATLAB 的M 文件编辑器中编写的M 文件，如下： f=inline('t^(2/3)-(t^2+1)^(1/3)','t'); a=0;b=3;epsilon=0.01; t1=b-0.618*(b-a);f1=f(t1); t2=a+0.618*(b-a);f2=f(t2); k=1; while abs(b-a)>=epsilon if f1 Matalab课后作业 学院：电气信息工程及其自动化 班级： 学号： 姓名： 完成日期： 2012年12月23日 1、 matlab 软件主要功能是什么？电气工程及其自动化专业本科生主要用到哪 些工具箱，各有什么功能？ 答：（1）主要功能：工业研究与开发； 数学教学，特别是线性代数；数值分析和科学计算方面的教学与研究；电子学、控制理论和物理学等工程和科学学科方面的教学与研究； 经济学、化学和生物学等计算问题的所有其他领域中的教学与研究；符号计算功能；优化工具；数据分析和可视化功能；“活”笔记本功能；工具箱；非线性动态系统建模和仿真功能。 （2）常用工具箱： （a ） MATLAB 主工具箱：扩充matlab 的数值计算、符号运算功能、图形建模仿真功能、文字处理功能以及与硬件实时交互功能。 （b ）符号数学工具箱：符号表达式、符号矩阵的创建；符号可变精度求解；因式分解、展开和简化；符号代数方程求解；符号微积分；符号微分方程。 （c ） SIMULINK 仿真工具箱： Simulink 是用于动态系统和嵌入式系统的多领域仿真和基于模型的设计工具。对各种时变系统，包括通讯、控制、信号处理、视频处理和图像处理系统，Simulink 提供了交互式图形化环境和可定制模块库来对其进行设计、仿真、执行和测试。 （d ）信号处理工具箱：数字和模拟滤波器设计、应用及仿真；谱分析和估计；FFT 、DCT 等 变换；参数化模型。 （e ）控制系统工具箱：连续系统设计和离散系统设计；状态空间和传递函数以及模型转换；时域响应（脉冲响应、阶跃响应、斜坡响应）；频域响应（Bode 图、Nyquist 图）；根轨迹、极点配置。 2、设y=23e t 4-sin(43t+3 ),要求以0.01秒为间隔，求出y 的151个点，并求出其导数的值和曲线。 程序如下： clc clear x=0:0.01:1.5; y=sqrt(3)/2*exp(-4*x).*sin(4*sqrt(3)*x+pi/3); y1=diff(y); subplot(2,1,1) plot(x,y) subplot(2,1,2) plot(x(1:150),y1) 曲线如下图所示： %要求设计一butterworth低通数字滤波器，wp=30hz，ws=40hz，rp=0.5，rs=40，fs=100hz。>>wp=30;ws=40;rp=0.5;rs=40;fs=100; >>wp=30*2*pi;ws=40*2*pi; >> [n,wn]=buttord(wp,ws,rp,rs,'s'); >> [z,p,k]=buttap(n); >> [num,den]=zp2tf(z,p,k); >> [num1,den1]=impinvar(num,den); Warning: The output is not correct/robust. Coeffs of B(s)/A(s) are real, but B(z)/A(z) has complex coeffs. Probable cause is rooting of high-order repeated poles in A(s). > In impinvar at 124 >> [num2,den2]=bilinear(num,den,100); >> [h,w]=freqz(num1,den1); >> [h1,w1]=freqz(num2,den2); >>subplot(1,2,1); >>plot(w*fs/(2*pi),abs(h)); >>subplot(1,2,2); >>plot(w1*fs/(2*pi),abs(h1)); >>figure(1); >>subplot(1,2,1); >>zplane(num1,den1); >>subplot(1,2,2); >>zplane(num2,den2); 机械优化设计在matlab中的应用 东南大学机械工程学院** 一优化设计目的： 在生活和工作中，人们对于同一个问题往往会提出多个解决方案，并通过各方面的论证从中提取最佳方案。最优化方法就是专门研究如何从多个方案中科学合理地提取出最佳方案的科学。由于优化问题无所不在，目前最优化方法的应用和研究已经深入到了生产和科研的各个领域，如土木工程、机械工程、化学工程、运输调度、生产控制、经济规划、经济管理等，并取得了显著的经济效益和社会效益。 " 二优化设计步骤： 1.机械优化设计的全过程一般可以分为如下几个步骤： 1)建立优化设计的数学模型； 2)选择适当的优化方法； 3)编写计算机程序； : 4)准备必要的初始数据并伤及计算； 5)对计算机求得的结果进行必要的分析。 其中建立优化设计数学模型是首要的和关键的一步，它是取得正确结果的前提。优化方法的选取取决于数学模型的特点，例如优化问题规模的大小，目标函数和约束函数的性态以及计算精度等。在比较各种可供选用的优化方法时，需要考虑的一个重要因素是计算机执行这些程序所花费的时间和费用，也即计算效率。 2.建立数学模型的基本原则与步骤 ①设计变量的确定； — 设计变量是指在优化设计的过程中，不断进行修改，调整，一直处于变化的参数称为设计变量。设计变量的全体实际上是一组变量，可用一个列向量表示： x=。 ②目标函数的建立； 选择目标函数是整个优化设计过程中最重要的决策之一。当对某以设计性能有特定的要求，而这个要求有很难满足时，则针对这一性能进行优化会得到满意的效果。目标函数是设计变量的函数，是一项设计所追求的指标的数学反映，因此它能够用来评价设计的优劣。 目标函数的一般表达式为： 。 f(x)=，要根据实际的设计要求来设计目标函数。 ③约束条件的确定。 一个可行性设计必须满足某些设计限制条件，这些限制条件称为约束条件，简称约束。 由若干个约束条件构成目标函数的可行域，而可行域内的所有设计点都是满足设计要求的，一般情况下，其设计可行域可表示为 高级数字滤波器设计及Matlab 实现 利用Parks-McClellan 算法设计线性相位FIR 滤波器 一、 算法原理 长度为2n+1的线性相位数字滤波器的传输函数为：20 ()n k k k G z h Z -== ∑。当 Z=exp(j T ω)=exp(2j F π)时，可得到频率响应： ()exp(2)cos 2n k k G F j nF d k F ππ==-∑ exp(2)()j F H F π=- 其中2n k k d h -=，k=0，……，n-1，0n d h = max ()()()W F D F H F - 对于一个理想的低通滤波器上式中的H(F)可以表示为： 误差加权函数： 允许设计者自己给定通带和阻带内的误差范围。令p s B B A = ，设计长度为2n+1的线性相位低通滤波器只要找到k d 并使得m a x ()()()W F D F H F -最小。 设空间A 是[0,1/2]的封闭子空间，为了使0 ()cos 2n k k H F d kF π== ∑是D(F) 在A 上唯一的最佳逼近，加权误差方程()()[()()]E F W F D F H F =-在A 上至少要有n+2个交错点。因此1()()i i E F E F E -=-=±，011,n i F F F F A +<∈ ， max ()E E F =。算法的流程如图1所示。 对于给定的n+2个频率点，需要计算n+2个方程： ()(()())(1)k k k k W F H F D F ρ-=-- 写成矩阵的形式就是： 图 1. 0000001 1 1 1 1011 1 1 1 11 1cos 2cos 4cos 2()()11cos 2cos 4cos 2()()()(1)1cos 2cos 4cos 2()n n n n n n n F F nF W F d D F d F F nF D F W F d D F F F nF W F ππππππρπππ++++++?? ???? ?? ?? ??-????????????=??????????????????-?? ????? ? 通过该方程组可得： 其中： 利用拉格朗日插值公式可得： 这里 利用求得的H(F)求出误差函数E(F)。如果对所有的频率都有()E F ρ≤，说明ρ是纹波极值，交错频率点121,n F F F + 是交错频率点。若存在某些频率使得()E F ρ>，说明初始交错点组中的某些点需要交换。 对于上次确定的121,n F F F + 中每一点，都检查其附近是否存在某一频率 ()E F ρ>，如果有再在该点附近找出局部极值点，并用该，点代替原来的 第三章1. (1)A=eye(3) (2)A=100+100*rand(5,6) (3)A=1+sqrt(0.2)*randn(10,50) (4)B=ones(size(A)) (5)A+30*eye(size(A)) (6)B=diag(diag(A)) 2. B=rot90(A) C=rot90(A,-1) 3. B=inv(A) ;A的逆矩阵 C=det(A) ;A的行列式的值 D=A*B E=B*A D=E 因此A与A-1是互逆的。 4. A=[4 2 -1;3 -1 2;12 3 0]; b=[2;10;8]; x=inv(A)*b x = -6.0000 26.6667 27.3333 5. (1) diag(A) ；主对角线元素 ans = 1 1 5 9 triu(A) ；上三角阵 ans = 1 -1 2 3 0 1 -4 2 0 0 5 2 0 0 0 9 tril(A) ；下三角阵 ans = 1 0 0 0 5 1 0 0 3 0 5 0 11 15 0 9 rank(A) ；秩 ans = 4 norm(A) ；范数 ans = 21.3005 cond(A) ；条件数 ans = 11.1739 trace(A) ；迹 ans = 16 （2）略 6. A=[1 1 0.5;1 1 0.25;0.5 0.25 2] A = 1.0000 1.0000 0.5000 1.0000 1.0000 0.2500 0.5000 0.2500 2.0000 [V,D]=eig(A) V = 0.7212 0.4443 0.5315 -0.6863 0.5621 0.4615 -0.0937 -0.6976 0.7103 D = -0.0166 0 0 0 1.4801 0 0 0 2.5365 课程设计任务书 学生姓名：专业班级：电信 指导教师：工作单位：信息工程学院 题目:MATLAB运算与应用设计2 初始条件： 1.MATLAB6.5以上版本软件； 2.课程设计辅导资料：“MATLAB语言基础及使用入门”、“MATLAB及在电子信息课程中的 应用”等； 3.先修课程：信号与系统、数字信号处理、MATLAB应用实践及信号处理类课程等。 要求完成的主要任务: （包括课程设计工作量及其技术要求，以及说 明书撰写等具体要求） （1）选择一本《MATLAB教程》，学习该教程的全部内容，包括使用方法、数组运算、矩阵运算、数学运算、程序设计、符号计算、图形绘制、GUI设计等内容； （2）对该套综合题的10道题，进行理论分析，针对具体设计部分的原理分析、建模、必要的推导和可行性分析，画出程序设计框图，编写程序代码（含注释），上机调试运行程序，记录实验结果（含计算结果和图表）。 （3）对实验结果进行分析和总结； （4）要求阅读相关参考文献不少于5篇； （5）根据课程设计有关规范，按时、独立完成课程设计说明书。 时间安排: (1) 布置课程设计任务,查阅资料,学习《MATLAB教程》十周； (2) 进行编程设计一周； (3) 完成课程设计报告书一周； 指导教师签名：年月日 系主任（或责任教师）签名：年月日 目录 1 MATLAB概述 (3) 1.1MATLAB简介 (4) 1.2MATLAB的功能 (4) 1.3MATLAB 的典型应用 (6) 2设计题目：MATLAB运算与应用设计套题二 (6) 3设计内容 (8) 3.1 题一 (8) 3.2 题二 (8) 3.3 题三 (9) 3.4 题四 (10) 3.5 题五 (15) 3.6 题六 (15) 3.7 题七 (15) 3.8 题八 (16) 3.9 题九 (17) 3.10题十 (18) 4 课程设计心得 (20) 5参考文献 (21) 6 本科生课程设计成绩评定表 (22) 数字信号处理课程设计 《数字信号处理》 课程设计报告 FIR数字滤波器设计及MATLAB实现 专业：通信工程 班级：通信1101班 组次：第9组 姓名及学号： 姓名及学号： 目录 一、设计目的 (3) 二、设计任务 (3) 三、设计原理 (3) 3.1窗函数法 (3) 3.2频率采样法 (4) 3.3最优化设计 (5) 3.3.1等波纹切比雪夫逼近准则 (5) 3.3.2仿真函数 (6) 四、设计过程 (7) 五、收获与体会 (13) 参考文献 (13) FIR数字滤波器设计及MATLAB实现 一、设计目的 FIR滤波器：有限长单位冲激响应滤波器，是数字信号处理系统中最基 本的元件，它可以在保证任意幅频特性的同时具有严格的线性相频特性， 同时其单位抽样响应是有限长的，因而滤波器是稳定的系统。因此，FIR 滤波器在通信、图像处理、模式识别等领域都有着广泛的应用。滤波器设 计是根据给定滤波器的频率特性，求得满足该特性的传输函数。 二、设计任务 FIR滤波器设计的任务是选择有限长度的() H e满足一定 h n，使传输函数()jw 的幅度特性和线性相位要求。由于FIR滤波器很容易实现严格的线性相位，所以FIR 数字滤波器设计的核心思想是求出有限的脉冲响应来逼近给定的频率响应。 设计过程一般包括以下三个基本问题： （1）根据实际要求确定数字滤波器性能指标； （2）用一个因果稳定的系统函数去逼近这个理想性能指标； （3）用一个有限精度的运算去实现这个传输函数。 三、设计原理 FIR滤波器设计的任务是选择有限长度的() H e满足一定 h n，使传输函数()jw 的幅度特性和线性相位要求。由于FIR滤波器很容易实现严格的线性相位，所以FIR数字滤波器设计的核心思想是求出有限的脉冲响应来逼近给定的频率响应。 设计过程一般包括以下三个基本问题： （1）根据实际要求确定数字滤波器性能指标； （2）用一个因果稳定的系统函数去逼近这个理想性能指标； （3）用一个有限精度的运算去实现这个传输函数。 3.1窗函数法 设计FIR数字滤波器的最简单的方法是窗函数法，通常也称之为傅立叶级数法。FIR数字滤波器的设计首先给出要求的理想滤波器的频率响应()jw H e,设计 dmatlab课程设计题目全
matlab数字滤波器设计程序
matlab(四连杆优化设计)
高级数字滤波器设计及Matlab实现
matlab程序设计与应用(第二版)第三章部分课后答案
MATLAB课程设计任务书
FIR数字滤波器设计及MATLAB使用要点
相关文档
- 【精品】PPT课件 MATLAB程序设计基础
- matlab 程序设计基础课程
- 第三章 Matlab语言程序设计基础 (1)
- matlab的程序设计基础
- 3matlab程序设计基础
- MATLAB程序设计基础介绍
- MATLAB程序设计基础--精品PPT课件
- MATLAB程序设计基础25634
- 第三章 MATLAB 的程序设计基础
- 最新matlab程序设计基础教程
- 2019年matlab程序设计基础教程_图文.ppt.ppt
- MATLAB程序设计基础
- Matlab程序设计教程(第二版)课后参考答案
- MATLAB程序设计与应用
- matlab程序设计与应用(第二版)第三章部分课后答案
- 专题四MATLAB程序设计基础精品
- 第五讲 Matlab程序设计基础PPT课件
- MATLAB程序设计基础
- 《MATLAB程序设计》复习题