轧钢机matlab仿真

河南工业大学

控制系统仿真

姓名：张天赐

班级：自动化103

学号：201123910415

成绩：

2013年10月30 日

中国于 1871 年在福州船政局所属拉铁厂开始用轧机﹔轧制厚 15mm 以下的铁板﹐ 6 ～ 120mm 的方﹑圆钢。 1890 年汉冶萍公司汉阳铁厂装有蒸汽机拖动的横列双机架 2450mm 二辊中板轧机和蒸汽机拖动的三机架横列二辊式轨梁轧机以及 350/300mm 小型轧机。随着冶金工业的发展﹐现已有多种类型轧机。

1.1轧钢机的原理图及结构

轧钢机原理图

轧钢机的主要设备有工作机座和传动装置作机座（轧辊﹑轧辊轴承﹑机架﹑轨座﹑轧辊调整装置﹑上轧辊平衡装置和换辊装置等组成。）轧辊轴承，支承轧辊并保持轧辊在机架中的固定位置。轧辊轴承工作负荷重而变化大﹐因此要求轴承摩擦系数小﹐具有足够的强度和刚度﹐而且要便于更换轧辊。不同的轧机选用不同类型的轧辊轴承。滚动轴承的刚性大﹐摩擦系数较小﹐但承压能力较小﹐且外形尺寸较大﹐多用于板带轧机工作辊。滑动轴承有半干摩擦与液体摩擦两种。半干摩擦轧辊轴承主要是胶木﹑铜瓦﹑尼龙瓦轴承﹐比较便宜﹐多用于型材轧机和开坯机。液体摩擦轴承有动压﹑静压和静 - 动压三种。优点是摩擦系数比较小﹐承压能力较大﹐使用工作速度高﹐刚性好﹐缺点是油膜厚度随速度而变化。液体摩擦轴承多

用于板带轧机支承辊和其它高速轧机。

轧机机架，由两片“牌坊”组成以安装轧辊轴承座和轧辊调整装置﹐需有足够的强度和钢度承受轧制力。机架形式主要有闭式和开式两种。闭式机架是一个整体框架﹐具有较高强度和刚度﹐主要用于轧制力较大的初轧机和板带轧机等。开式机架由机架本体和上盖两部分组成﹐便于换辊﹐主要用于横列式型材轧机。传动装置（电动机﹑减速机﹑齿轮座和连接轴等组成。齿轮座将传动力矩分送到两个或几个轧辊上。）

轧钢机系统框图

轧钢机调速系统的组成和原理

采用PI调节的单个转速闭环直流调速系统可以在保证系统稳定的前提下实现转速无静差。但是，如果对系统的动态性能要求较高，单闭环系统就难以满足需要，这主要是因为在单闭环系统中不能控制电流和转矩的动态过程。电流截止负反馈环节是专门用来控制电流的，并不能很理想地控制电流的动态波形。

在起动过程中，始终保持电流（转矩）为允许的最大值，使电力拖动系统以最大的加速度起动，到达稳态转速时，立即让电流降下来，使转矩马上与负载相平衡，从而转入稳态运行。

为了实现在允许条件下的最快起动，关键是要获得一段使电流保持为最I的恒流过程。按照反馈控制规律，采用某个物理量的负反馈就可以大值

保持该量基本不变，那么，采用电流负反馈应该能够得到近似的恒流过程。应该在起动过程中只有电流负反馈，没有转速负反馈，达到稳态转速后，又希望只要转速负反馈，不再让电流负反馈发挥作用。

2.1转速、电流双闭环直流调速系统的组成

系统中设置两个调节器，分别调节转速和电流，如图1-1所示。把转速调节器的输出当作电流调节器的输入，再用电流调节器的输出去控制电力电子变换器UPE 。从闭环结构上看，电流环在里面，称作内环；转速环在外边，称作外环。这就形成了转速、电流双闭环调速系统。

3.2双闭环直流调速系统的动态数学模型

双闭环直流调速系统的动态结构图，如图2-1所示。图中

)

(s W ASR 和

)

(s W ACR 分别表示转速调节器和电流调节器的传递函数。

TG n

ASR

ACR

+ -

U n

U i

U c

U d

I d

UPE

图1-1 转速、电流双闭环直流调速系统框图

ASR —转速调节器 ACR —电流调节器 TG —测速发电机

TA —电流互感器 UPE —电力电子变换器

内外

22222

4.3960174.052621--∑=??=+=

s s T h h K n

N ，因为C e =（U N -I N R a ）/n N =（440-220×0.088）/1000=0.234V ?min/r ，α=10V/ n N =10/1800=0.006

V ?r/min ，于是可得ASR 的比例系数为

36.30174

.012.0006.05201

.0234.003.062)1(=???????=+=

∑n m e n RT h T C h K αβ

6.5 转速调节器的实现

取R 0=40k Ω，则R n =K n R 0=3.36×40=134.48k Ω，取1400k Ω；

C n =τn /R n =0.087/（140×103）≈0.621×10-6F=0.621μF ，取0.7μF ；

C on =4T on /R 0=4×0.01/（40×103）=1×10-6

=1μF ，取1μF 。

对双闭环直流调速系统中的电力拖动系统、直流电动机电枢回路、可控硅电路三部分进行串联校正设计，得此开环传递函数为

s T s T T K s T R

s T R s T K s G i s m s m i s s )1)(1(111)(0++=?+?+=

令12.0=s T ，02.0=i T ，m s T K K =

由以上所取的参数可列写出系统反馈系统传递函数可近似为

s s K

s G )02.01)(12.01()(++=

用matlab 软件画出根轨迹为

当阻尼比为0.707时，K 值应小于正负45度直线与根轨迹交点。经计算可以得到

Matlab通信系统仿真实验报告

Matlab通信原理仿真学号： 2142402 姓名：圣斌

实验一Matlab 基本语法与信号系统分析一、实验目的： 1、掌握MATLAB的基本绘图方法； 2、实现绘制复指数信号的时域波形。二、实验设备与软件环境： 1、实验设备：计算机 2、软件环境：MATLAB R2009a 三、实验内容： 1、MATLAB为用户提供了结果可视化功能，只要在命令行窗口输入相应的命令，结果就会用图形直接表示出来。 MATLAB程序如下： x = -pi::pi; y1 = sin(x); y2 = cos(x); %准备绘图数据 figure(1); %打开图形窗口 subplot(2,1,1); %确定第一幅图绘图窗口 plot(x,y1); %以x，y1绘图 title('plot(x,y1)'); %为第一幅图取名为’plot(x,y1)’ grid on; %为第一幅图绘制网格线 subplot(2,1,2) %确定第二幅图绘图窗口 plot(x,y2); %以x，y2绘图 xlabel('time'),ylabel('y') %第二幅图横坐标为’time’,纵坐标为’y’运行结果如下图： 2、上例中的图形使用的是默认的颜色和线型，MATLAB中提供了多种颜色和线型，并且可以绘制出脉冲图、误差条形图等多种形式图： MATLAB程序如下： x=-pi:.1:pi; y1=sin (x); y2=cos (x); figure (1); %subplot (2,1,1); plot (x,y1); title ('plot (x,y1)'); grid on %subplot (2,1,2); plot (x,y2);

matlab常用计算命令

Matlab常用计算命令（部分） by sunny_疑似天人 1.多项式运算： poly2sym函数，将给定的多项式系数向量转化为符号表达式，以降幂排序。 poly函数，得到矩阵的特征多项式（首项系数为1）的系数向量，然后也可以用poly2sym函数转化为多项式的符号表达式。 roots函数，得到方程的根，调用形式为roots（a），其中a 为多项式的系数；也可以直接调用roots（[1 2 1]）。 compan函数与eig函数，通过compan函数建立多项式的伴随矩阵再通过eig函数求伴随矩阵特征值以得到多项式的所有根。效果与roots函数相同；同时这两个函数也可单独使用： compan函数，建立多项式的伴随矩阵，如：a=[1 2 3 ];compan(a) ans = -2 -3 1 0 eig函数，求矩阵的特征值。 conv函数，求多项式的乘积，如：pd=conv(p,d)，其中p和d均为多项式系数向量，得到的同样也是多项式的系数向量。 deconv函数，求多项式的除法。 polyder函数,求多项式的微分。即求一阶导数，如果要求多项式的高阶微分，可以通过循环实现。

polyfit 函数，对数据拟合得到多项式，这个多项式即可大致代表数据变化规律。例如： x=0:pi/20:pi/2; y=sin(x); p=polyfit(x,y,5) x1=0:pi/30:pi*2; y1=sin(x1); y2=polyval(p,x1); plot(x1,y1,'b-',x1,y2,'r*') legend('原曲线','拟合曲线') axis([0 7 -1.2 4]) p = 0.0057 0.0060 -0.1721 0.0021 0.9997 0.0000 1 2 3 4 5 6 7 -1-0.500.511.522.533.54 原曲线拟合曲线 2.向量及其运算 x=linspace(a,b,n)，生成一个向量x ，其中a ，b 分别是生成矢量的第一个和最后一个元素，n 是采样总点数。当n 缺省时默认生成100维的向量。

时域有限差分法的Matlab仿真

时域有限差分法的Matlab仿真关键词: Matlab 矩形波导时域有限差分法摘要：介绍了时域有限差分法的基本原理，并利用Matlab仿真，对矩形波导谐振腔中的电磁场作了模拟和分析。关键词：时域有限差分法；Matlab；矩形波导；谐振腔目前，电磁场的时域计算方法越来越引人注目。时域有限差分（Finite Difference Time Domain，FDTD）法［1］作为一种主要的电磁场时域计算方法，最早是在1966年由K. S. Yee提出的。这种方法通过将Maxwell旋度方程转化为有限差分式而直接在时域求解，通过建立时间离散的递进序列，在相互交织的网格空间中交替计算电场和磁场。经过三十多年的发展，这种方法已经广泛应用到各种电磁问题的分析之中。 Matlab作为一种工程仿真工具得到了广泛应用［2］。用于时域有限差分法，可以简化编程，使研究者的研究重心放在FDTD法本身上，而不必在编程上花费过多的时间。下面将采用FDTD法，利用Matlab仿真来分析矩形波导谐振腔的电磁场，说明了将二者结合起来的优越性。 1FDTD法基本原理时域有限差分法的主要思想是把Maxwell方程在空间、时间上离散化，用差分方程代替一阶偏微分方程，求解差分方程组，从而得出各网格单元的场值。FDTD 空间网格单元上电场和磁场各分量的分布如图1所示。电场和磁场被交叉放置，电场分量位于网格单元每条棱的中心，磁场分量位于网格单元每个面的中心，每个磁场（电场）分量都有4个电场（磁场）分量环绕。这样不仅保证了介质分界面上切向场分量的连续性条件得到自然满足，而且

还允许旋度方程在空间上进行中心差分运算，同时也满足了法拉第电磁感应定律和安培环路积分定律，也可以很恰当地模拟电磁波的实际传播过程。 1.1Maxwell方程的差分形式旋度方程为：将其标量化，并将问题空间沿3个轴向分成若干网格单元，用Δx，Δy和Δz 分别表示每个网格单元沿3个轴向的长度，用Δt表示时间步长。网格单元顶点的坐标（x，y，z）可记为：其中：i，j，k和n为整数。同时利用二阶精度的中心有限差分式来表示函数对空间和时间的偏导数，即可得到如下FDTD基本差分式：由于方程式里出现了半个网格和半个时间步，为了便于编程，将上面的差分式改写成如下形式：

MATLAB计算概率

一、实验名称已知随机向量（X ，Y ）独立同服从标准正态分布，D={(x,y)|a0&&e<6 if e==1

p=erchong(a,b,c,d) end if e==2 p=wangge(a,b,c,d); end if e==3 p=fenbu(a,b,c,d); end if e==4 p=mente(a,b,c,d); end if e==5 [X,Y]=meshgrid(-3:0.2:3); Z=1/(2*pi)*exp(-1/2*(X.^2+Y.^2)); meshz(X,Y,Z); end e=input('请选择: \n'); end % ===============================用二重积分计算function p=erchong(a,b,c,d) syms x y; f0=1/(2*pi)*exp(-1/2*(x^2+y^2)); f1=int(f0,x,a,b); %对x积分 f1=int(f1,y,c,d); %对y积分 p=vpa(f1,9); % ================================等距网格法function p=wangge(a,b,c,d) syms x y ; n=100; r1=(b-a)/n; %求步长 r2=(d-c)/n; za(1)=a;for i=1:n,za(i+1)=za(i)+r1;end %分块 zc(1)=c;for j=1:n,zc(j+1)=zc(j)+r2;end for i=1:n x(i)=unifrnd(za(i),za(i+1));end %随机取点 for i=1:n y(i)=unifrnd(zc(i),zc(i+1));end s=0; for i=1:n for j=1:n s=1/(2*pi)*exp(-1/2*(x(i)^2+y(j)^2))+s;%求和end end p=s*r1*r2;

各种BP学习算法MATLAB仿真

3.3.2 各种BP学习算法MATLAB仿真根据上面一节对BP神经网络的MATLAB设计，可以得出下面的通用的MATLAB程序段，由于各种BP学习算法采用了不同的学习函数，所以只需要更改学习函数即可。 MATLAB程序段如下： x=-4:0.01:4; y1=sin((1/2)*pi*x)+sin(pi*x); %trainlm函数可以选择替换 net=newff(minmax(x),[1,15,1],{'tansig','tansig','purelin'},'trainlm'); net.trainparam.epochs=2000; net.trainparam.goal=0.00001; net=train(net,x,y1); y2=sim(net,x); err=y2-y1; res=norm(err); %暂停，按任意键继续 Pause %绘图，原图（蓝色光滑线）和仿真效果图（红色+号点线） plot(x,y1); hold on plot(x,y2,'r+'); 注意：由于各种不确定因素，可能对网络训练有不同程度的影响，产生不同的效果。如图3-8。标准BP算法（traingd）

图3-8 标准BP算法的训练过程以及结果（原图蓝色线，仿真图+号线）增加动量法(traingdm) 如图3-9。图3-9 增加动量法的训练过程以及结果（原图蓝色线，仿真图+号线）弹性BP算法（trainrp）如图3-10 图3-10 弹性BP算法的训练过程以及结果（原图蓝色线，仿真图+号线）

动量及自适应学习速率法（traingdx）如图3-11。图3-11 动量及自适应学习速率法的训练过程以及结果（原图蓝色线，仿真图+号线）共轭梯度法(traincgf)如图3-12。

matlab与通信仿真实验指导书(上)

《Matlab与通信仿真》实验指导书(上) 刘毓杨辉徐健和煦黄庆东吉利萍编著通信与信息工程学院 2011-1

目录第一章 MALTAB基础知识 (1) 1.1MATLAB基础知识 (1) 1.2MATLAB基本运算 (2) 1.3MATLAB程序设计 (7) 第二章 MATLAB计算结果可视化和确知信号分析 (13) 2.1计算结果可视化 (13) 2.2确知信号分析 (17) 第三章随机信号与数字基带仿真 (23) 3.1基本原理 (23) 3.2蒙特卡罗算法 (30) 第四章模拟调制MATLAB实现 (34) 4.1模拟调制 (34) 4.2信道加性高斯白噪声 (35) 4.3AM调制解调的MATLAB实现 (36) 第五章模拟信号的数字传输 (45) 5.1脉冲编码调制 (45) 5.2低通抽样定理 (45) 5.3均匀量化原理 (46) 5.4非均匀量化 (48) 第六章数字频带传输系统 (52) 6.1数字频带传输原理 (52) 6.2数字频带传输系统的MATLAB实现 (53) 第七章通信系统仿真综合实验 (67) 7.1基本原理 (67) 7.2实验内容 (67)

第一章 MALTAB基础知识本章目标 ●了解MATLAB 程序设计语言的基本特点，熟悉MATLAB软件运行环境 ●掌握创建、保存、打开m文件及函数的方法 ●掌握变量等有关概念，具备初步的将一般数学问题转化为对应的计算机模型并进行处理的能力 1.1 MATLAB基础知识 1.1.1 MATLAB程序设计语言简介 MATLAB，Matrix Laboratory的缩写，是由MathWorks公司开发的一套用于科学工程计算的可视化高性能语言，具有强大的矩阵运算能力。与大家常用的Fortran和C等高级语言相比，MATLAB的语法规则更简单，更贴近人的思维方方式，被称为“草稿纸式的语言”。MATLAB软件主要由主包、仿真系统（simulink）和工具箱（toolbox）三大部分组成。 1.1.2 MATLAB界面及帮助 MATLAB基本界面如图1-1所示，命令窗口包含标题栏、菜单栏、工具栏、命令行区、状态栏、垂直和水平波动条等区域。图1-1 MATLAB基本界面（1）菜单栏

MATLAB仿真教程

一、设计目的通过运用MATLAB对函数进行Z域分析和单边带信号的调制与解调，使我们进一步加深对MATLAB的认识和运用，以实现以下目的： 1.本次试验进一步熟悉了MATLAB软件的使用方法及相关的操作。 2.对Z变换及其反变换函数在MATLAB中的调用有了掌握。 3.理论与实际的仿真相结合，更直观的看到结果。 4.观察了单边带信号调制与解调后的图像，加深认识。二、设计原理 MATLAB是The MathWorks公司在1984年推出的一种商品化软件，它提供了大量丰富的应用函数，并且具有扩充的开放性结构。目前，该软件包涵盖了控制系统应用、数字信号处理、数字图像处理、通讯、神经网络、小波理论分析、优化与统计、偏微分方程、动态系统实时仿真等多学科专业领域。其中单边带调制信号是将双边带信号中的一个边带滤掉而形成的。根据方法的不同，产生单边带调制信号的方法有：滤波和相移法。由于滤波法在技术上比较难实现所以在此我们将用相移法对单边带调制与解调系统进行讨论与设计。三、设计内容和MATLAB图像

1、数字系统的响应源代码如下： b=[0 1 2 1 0]; a=[1 -0.5 0 0.3 -0.005]; subplot(421);zplane(b,a); title('系统的零极点图'); subplot(422);impz(b,a,21); title('单位脉冲响应'); subplot(423);stepz(b,a,21); title('单位阶跃响应');

N=21;n=0:N-1; x=exp(-n); x0=zeros(1,N); y0=[1,-1]; xi=filtic(b,a,y0); y1=filter(b,a,x0,xi); xi0=filtic(b,a,0); y2=filter(b,a,x,xi0); y3=filter(b,a,x,xi); [h w]=freqz(b,a,21); subplot(424);stem(n,y1); title('零输入响应');grid on; subplot(425);stem(n,y2); title('零状态响应');grid on; subplot(426);stem(n,y3); title('系统的全响应');grid on; subplot(427);plot(w,abs(h)); title('幅频特性曲线');grid on; subplot(428);plot(w,angle(h)); title('相频特性曲线');grid on;

利用MATLAB实现循环卷积.doc

一、实验目的 1.利用MATLAB 实现循环卷积。 2.比较循环卷积与线性卷积的区别。二、实验条件 PC 机，MATLAB7.0 三、实验内容 1）循环卷积的定义：两个序列的N 点循环卷积定义为： )0()()()]()([1 0N n m n x m h n x n h N k N N <≤-=?∑-= 利用MATLAB 实现两个序列的循环卷积可以分三个步骤完成： (1)初始化：确定循环点数N ，测量输入2个序列的长度。 (2)循环右移函数：将序列x(n)循环右移，一共移N 次（N 为循环卷积的循环次数），最后将每次循环成的新序列组成一个矩阵V 。 (3)相乘：将x(n)移位后组成的矩阵V 与第二个序列h(n)对应相乘，即得循环卷积结果。程序如下：程序一： clear;close all ; N=10; x1=[6 15 -6 3 5 7 0 1]; x2=[7 1 2 9 4 3 20 6]; xn1=length(x1); xxn1=0:xn1-1; xn2=length(x2); xxn2=0:xn2-1; subplot(3,1,1); stem(xxn1,x1); subplot(3,1,2); stem(xxn2,x2); x11=fft(x1,N);

x12=fft(x2,N); y11=x11.*x12; y1=ifft(y11,N); subplot(3,1,3); n=0:length(y1)-1; stem(n,y1,'.'); title('循环卷积的结果'); xlabel('n');ylabel('y1(n)'); 运行后所得图形如下：观察所得的循环卷积结果发现并没有呈现周期性的序列，因此将程序做下列改变。程序二： clear;close all; N=40; x1=[6 15 -6 3 5 7 0 1]; x2=[7 1 2 9 4 3 20 6]; x2=[x2,x2,x2,x2]; xn1=length(x1); xxn1=0:xn1-1; xn2=length(x2); xxn2=0:xn2-1; subplot(3,1,1);

MATLAB实现通信系统仿真实例

补充内容：模拟调制系统的MATLAB 仿真 1.抽样定理为了用实验的手段对连续信号分析，需要先对信号进行抽样（时间上的离散化），把连续数据转变为离散数据分析。抽样（时间离散化）是模拟信号数字化的第一步。 Nyquist 抽样定律：要无失真地恢复出抽样前的信号，要求抽样频率要大于等于两倍基带信号带宽。抽样定理建立了模拟信号和离散信号之间的关系，在Matlab 中对模拟信号的实验仿真都是通过先抽样，转变成离散信号，然后用该离散信号近似替代原来的模拟信号进行分析的。【例1】用图形表示DSB 调制波形)4cos()2cos(t t y ππ= 及其包络线。 clf %%计算抽样时间间隔 fh=1;%%调制信号带宽(Hz) fs=100*fh;%%一般选取的抽样频率要远大于基带信号频率，即抽样时间间隔要尽可能短。 ts=1/fs; %%根据抽样时间间隔进行抽样，并计算出信号和包络 t=(0:ts:pi/2)';%抽样时间间隔要足够小，要满足抽样定理。 envelop=cos(2*pi*t);%%DSB 信号包络 y=cos(2*pi*t).*cos(4*pi*t);%已调信号 %画出已调信号包络线 plot(t,envelop,'r:','LineWidth',3); hold on plot(t,-envelop,'r:','LineWidth',3); %画出已调信号波形 plot(t,y,'b','LineWidth',3); axis([0,pi/2,-1,1])% hold off% xlabel('t'); %写出图例【例2】用图形表示DSB 调制波形)6cos()2cos(t t y ππ= 及其包络线。 clf %%计算抽样时间间隔 fh=1;%%调制信号带宽(Hz) fs=100*fh;%抽样时间间隔要足够小，要满足抽样定理。 ts=1/fs; %%根据抽样时间间隔进行抽样

matlab仿真模拟

中国地质大学长城学院本科课程设计题目:双手协调机器人系别信息工程系学生姓名专业电气工程及其自动化学号指导教师王密香职称研究生 2015年12 月22日

双手协调机器人摘要多机器人的协同作业是制造业发展的必然要求，双臂机器人就是适应这一要求而开发出的一种新型机器人，相对于单臂机器人它可以大大增强机器人对复杂装配任务的适应性，同时可以提高工作空间的利用效率。当前大多数工业机器人的应用是为单臂机器人独自工作的能力准备的。一般地，单臂机器人只适合于刚性工件的操作，并受制于环境，随着现代工业的发展和科学技术的进步，对于许多任务而言单臂操作是不够的。为了适应任务的复杂性、智能性的不断提高以及系统柔顺性的要求而扩展为双手协调控制。即由两个单臂机器人相互协调、相互配合的去完成某种作业，但由于组成双手协调控制系统的是两个机器人它们不可能是两个单手机器人的简单组合，除了它们各自共同目标的控制实现外，它们相互间的协调控制以及对环境的适应性就成为组合的关键，这样双手协调控制机器人系统的进一步应用就受到了限制。而双臂机器人能完成对于人来说易于实现的功能，它比双手协调机器人更具有实用价值，它的高自律性以及学习性，能够适应许多环境，使其在工业生产、危险处理、国防、航天航空等方面运用广泛采用了Matlab/Simulink 仿真软件，分别用模块法和程序法对双手协调机器人系统进行了PID 控制器的校正仿真设计。关键词：双手协调机器人；Matlab/Simulink；PID 校正；仿真 Using the Matlab/Simulink simulation software, using method of module and the procedural law on hands coordinate robot system has carried on the correction for the simulation of PID controller design simulation results show that these two kinds of design method is not only convenient and quick, and the correction effect is satisfactory to people Keywords：Hands coordinate robot; Matlab/Simulink; PID correction; The simulation

matlab循环语句

matlab 基本语句 1.循环语句for for i=s1:s3：s2 循环语句组 end 解释：首先给i赋值s1；然后，判断i是否介于s1与s2之间；如果是，则执行循环语句组,i=i+s3（否则，退出循环.）；执行完毕后，继续下一次循环。例：求1到100的和，可以编程如下： sum=0 for i=1:1:100 sum=sum+i end 这个程序也可以用while语句编程。注：for循环可以通过break语句结束整个for循环. 2.循环语句while 例：sum=0;i=1; while(i<=100) sum=sum+i;i=i+1; end 3.if语句 if(条件) 语句 end if(条件) 语句 else 语句 end if(条件) 语句 elseif 语句 end 4.关系表达式：

=,>,<,>=,<=,==(精确等于) 5.逻辑表达式：|(或),&(且) 6.[n,m]=size(A)(A为矩阵) 这样可以得到矩阵A的行和列数 n=length(A)，可以得到向量A的分量个数；如果是矩阵，则得到矩阵A的行与列数这两个数字中的最大值。 7.！后面接Dos命令可以调用运行一个dos程序。 8.常见函数： poly():为求矩阵的特征多项式的函数，得到的为特征多项式的各个系数。如 a=[1,0,0;0,2,0;0,0,3]，则poly(a)＝1 -6 11 -6。相当于poly(a)＝1入^3+(-6)入^2+11入+(-6)。 compan():可以求矩阵的伴随矩阵. sin()等三角函数。 MATLAB在数学建模中的应用(3) 一、程序设计概述 MATLAB所提供的程序设计语言是一种被称为第四代编程语言的高级程序设计语言，其程序简洁，可读性很强，容易调试。同时，MATLAB的编程效率比C/C++语言要高得多。 MATLAB编程环境有很多。常用的有： 1.命令窗口 2.word窗口 3.M-文件编辑器，这是最好的编程环境。 M-文件的扩展名为“.m”。M-文件的格式分为两种： ①λ M-脚本文件，也可称为“命令文件”。 ② M-函数文件。这是matlab程序设计的主流。λ 保存后的文件可以随时调用。二、MATLAB程序结构按照现代程序设计的观点，任何算法功能都可以通过三种基本程序结构来实现，这三种结构是：顺序结构、选择结构和循环结构。其中顺序结构是最基本的结构，它依照语句的自然顺序逐条地执行程序的各条语句。如果要根据输入数据的实际情况进行逻辑判断，对不同的结果进行不同的处理，可以使用选择结构。如果需要反复执行某些程序段落，可以使用循环结构。 1 顺序结构顺序结构是由两个程序模块串接构成。一个程序模块是完成一项独立功能的逻辑单元，它可以是一段程序、一个函数，或者是一条语句。看图可知，在顺序结构中，这两个程序模块是顺序执行的，即先执行<程序

内点法matlab仿真doc资料

编程方式实现： 1.惩罚函数 function f=fun(x,r) f=x(1,1)^2+x(2,1)^2-r*log(x(1,1)-1); 2.步长的函数 function f=fh(x0,h,s,r) %h为步长 %s为方向 %r为惩罚因子 x1=x0+h*s; f=fun(x1,r); 3. 步长寻优函数 function h=fsearchh(x0,r,s) %利用进退法确定高低高区间，利用黄金分割法进行求解h1=0;%步长的初始点 st=0.001; %步长的步长 h2=h1+st; f1=fh(x0,h1,s,r); f2=fh(x0,h2,s,r); if f1>f2 h3=h2+st; f3=fh(x0,h3,s,r); while f2>f3 h1=h2; h2=h3; h3=h3+st; f2=f3; f3=fh(x0,h3,s,r); end else st=-st; v=h1; h1=h2; h2=v; v=f1; f1=f2; f2=v; h3=h2+st; f3=fh(x0,h3,s,r); while f2>f3 h1=h2; h2=h3; h3=h3+st; f2=f3;

f3=fh(x0,h3,s,r); end end %得到高低高的区间 a=min(h1,h3); b=max(h1,h3); %利用黄金分割点法进行求解 h1=1+0.382*(b-a); h2=1+0.618*(b-a); f1=fh(x0,h1,s,r); f2=fh(x0,h2,s,r); while abs(a-b)>0.0001 if f1>f2 a=h1; h1=h2; f1=f2; h2=a+0.618*(b-a); f2=fh(x0,h2,s,r); else b=h2; h2=h1; f2=f1; h1=a+0.382*(b-a); f1=fh(x0,h1,s,r); end end h=0.5*(a+b); 4. 迭代点的寻优函数 function f=fsearchx(x0,r,epson) x00=x0; m=length(x0); s=zeros(m,1); for i=1:m s(i)=1; h=fsearchh(x0,r,s); x1=x0+h*s; s(i)=0; x0=x1; end while norm(x1-x00)>epson x00=x1; for i=1:m s(i)=1; h=fsearchh(x0,r,s);

Matlab与通信仿真课程设计报告材料

《MATLAB与通信仿真》课程设计指导老师: 水英、汪泓班级:07通信(1)班学号:E07680104 :林哲妮

目录目的和要求 (1) 实验环境 (1) 具体容及要求 (1) 实验容题目一 (4) 题目容流程图程序代码仿真框图各个参数设置结果运行结果分析题目二 (8) 题目容流程图程序代码仿真框图各个参数设置结果运行结果分析题目三 (17) 题目容流程图程序代码仿真框图各个参数设置结果运行结果分析题目四 (33) 题目容流程图程序代码仿真框图各个参数设置结果运行结果分析心得与体会 (52)

目的和要求通过课程设计，巩固本学期相关课程MATLAB与通信仿真所学知识的理解，增强动手能力和通信系统仿真的技能。在强调基本原理的同时，更突出设计过程的锻炼。强化学生的实践创新能力和独立进行科研工作的能力。要求学生在熟练掌握MATLAB和simulink仿真使用的基础上，学会通信仿真系统的基本设计与调试。并结合通信原理的知识，对通信仿真系统进行性能分析。实验环境 PC机、Matlab/Simulink 具体容及要求基于MATLAB编程语言和SIMULINK通信模块库，研究如下问题： (1)研究BFSK在加性高斯白噪声信道下（无突发干扰）的误码率性能与信噪比之间的关系； (2)研究BFSK在加性高斯白噪声信道下（有突发干扰）的误码率性能与信噪比之间的关系；分析突发干扰的持续时间对误码率性能的影响。 (3)研究BFSK+信道编码（取BCH码和汉明码）在加性高斯白噪声信道下（无突发干扰）的误码率性能与信噪比之间的关系；分析不同码率对误码率性能的影响。比较不同信道编码方式的编码增益性能。 (4)研究BFSK+信道编码（取BCH码和汉明码）在加性高斯白噪声信道下（有突发干扰）的误码率性能与信噪比之间的关系；分析突发干扰的持续时间对误码率性能的影响。分析不同码率对误码率性能的影响。比较不同信道编码方式的编码增益性能。

MATLAB仿真实验全部

实验一 MATLAB 及仿真实验（控制系统的时域分析）一、实验目的学习利用MATLAB 进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性；二、预习要点 1、系统的典型响应有哪些？ 2、如何判断系统稳定性？ 3、系统的动态性能指标有哪些？三、实验方法（一）四种典型响应 1、阶跃响应：阶跃响应常用格式： 1、)(sys step ；其中sys 可以为连续系统，也可为离散系统。 2、),(Tn sys step ；表示时间围0---Tn 。 3、),(T sys step ；表示时间围向量T 指定。 4、),(T sys step Y =；可详细了解某段时间的输入、输出情况。 2、脉冲响应：脉冲函数在数学上的精确定义：0 ,0)(1)(0 ?==?∞ t x f dx x f 其拉氏变换为：)()()()(1 )(s G s f s G s Y s f === 所以脉冲响应即为传函的反拉氏变换。脉冲响应函数常用格式： ① )(sys impulse ； ② ); ,();,(T sys impulse Tn sys impulse ③ ),(T sys impulse Y = （二）分析系统稳定性有以下三种方法： 1、利用pzmap 绘制连续系统的零极点图； 2、利用tf2zp 求出系统零极点； 3、利用roots 求分母多项式的根来确定系统的极点（三）系统的动态特性分析 Matlab 提供了求取连续系统的单位阶跃响应函数step 、单位脉冲响应函数impulse 、零输入响应函数initial 以及任意输入下的仿真函数lsim.

PID控制算法的matlab仿真

PID 控制算法的matlab 仿真 PID 控制算法就是实际工业控制中应用最为广泛的控制算法,它具有控制器设计简单,控制效果好等优点。PID 控制器参数的设置就是否合适对其控制效果具有很大的影响,在本课程设计中一具有较大惯性时间常数与纯滞后的一阶惯性环节作为被控对象的模型对PID 控制算法进行研究。被控对象的传递函数如下: ()1d s f Ke G s T s τ-= + 其中各参数分别为30,630,60f d K T τ===。MATLAB 仿真框图如图1所示。图1 2 具体内容及实现功能 2、1 PID 参数整定 PID 控制器的控制参数对其控制效果起着决定性的作用,合理设置控制参数就是取得较好的控制效果的先决条件。常用的PID 参数整定方法有理论整定法与实验整定法两类,其中常用的实验整定法由扩充临界比例度法、试凑法等。在此处选用扩充临界比例度法对PID 进行整定,其过程如下: 1) 选择采样周期由于被控对象中含有纯滞后,且其滞后时间常数为 60d τ=,故可选择采样周期1s T =。 2) 令积分时间常数i T =∞,微分时间常数0d T =,从小到大调节比例系数K , 使得系统发生等幅震荡,记下此时的比例系数k K 与振荡周期k T 。 3) 选择控制度为 1.05Q =,按下面公式计算各参数:

0.630.490.140.014p k i k d k s k K K T T T T T T ==== 通过仿真可得在1s T =时,0.567,233k k K T ==,故可得: 0.357,114.17,32.62, 3.262p i d s K T T T ==== 0.0053.57 p s i i p d d s K T K T K T K T === = 按此组控制参数得到的系统阶跃响应曲线如图2所示。 01002003004005006007008009001000 0.20.40.60.811.21.41.6 1.8 图2 由响应曲线可知,此时系统虽然稳定,但就是暂态性能较差,超调量过大,且响应曲线不平滑。根据以下原则对控制器参数进行调整以改善系统的暂态过程: 1) 通过减小采样周期,使响应曲线平滑。 2) 减小采样周期后,通过增大积分时间常数来保证系统稳定。 3) 减小比例系数与微分时间常数,以减小系统的超调。改变控制器参数后得到系统的阶跃响应曲线如图3所示,系统的暂态性能得到明显改善、

Matlab与通信仿真课程设计报告

《MATLAB与通信仿真》课程设计指导老师: 张水英、汪泓班级:07通信(1)班学号:E07680104 姓名:林哲妮

目录目的和要求 (1) 实验环境 (1) 具体内容及要求 (1) 实验内容题目一 (4) 题目内容流程图程序代码仿真框图各个参数设置结果运行结果分析题目二 (8) 题目内容流程图程序代码仿真框图各个参数设置结果运行结果分析题目三 (17) 题目内容流程图程序代码仿真框图各个参数设置结果运行结果分析题目四 (33) 题目内容流程图程序代码仿真框图各个参数设置结果运行结果分析心得与体会 (52)

目的和要求通过课程设计，巩固本学期相关课程MATLAB与通信仿真所学知识的理解，增强动手能力和通信系统仿真的技能。在强调基本原理的同时，更突出设计过程的锻炼。强化学生的实践创新能力和独立进行科研工作的能力。要求学生在熟练掌握MATLAB和simulink仿真使用的基础上，学会通信仿真系统的基本设计与调试。并结合通信原理的知识，对通信仿真系统进行性能分析。实验环境 PC机、Matlab/Simulink 具体内容及要求基于MATLAB编程语言和SIMULINK通信模块库，研究如下问题： (1)研究BFSK在加性高斯白噪声信道下（无突发干扰）的误码率性能与信噪比之间的关系； (2)研究BFSK在加性高斯白噪声信道下（有突发干扰）的误码率性能与信噪比之间的关系；分析突发干扰的持续时间对误码率性能的影响。 (3)研究BFSK+信道编码（取BCH码和汉明码）在加性高斯白噪声信道下（无突发干扰）的误码率性能与信噪比之间的关系；分析不同码率对误码率性能的影响。比较不同信道编码方式的编码增益性能。 (4)研究BFSK+信道编码（取BCH码和汉明码）在加性高斯白噪声信道下（有突发干扰）的误码率性能与信噪比之间的关系；分析突发干扰的持续时间对误码率性能的影响。分析不同码率对误码率性能的影响。比较不同信道编码方式的编码增益性能。

MATLAB仿真程序

窗型选择仿真程序： clear,clc bw=3e6;% 信号带宽 T=1e-4;%信号脉冲宽度 A=2;%信号幅度 fs=4*bw; lfft=round(T*fs);%采样点数 lfft=2^nextpow2(lfft); dt=1/fs;%采样间隔 f0=1e6; t=(0:lfft-1)*dt;%时域采样点 q=(0:lfft-1)*2*pi/lfft; s=A*exp(j*2*pi*f0*t+j*pi*bw*t.*t/T);%产生线性调频信号 S=(fft(s));%线性调频信号的傅立叶变换fft H=conj(S);%匹配滤波器的频率响应 Y=S.*H;%线性调频信号的频域匹配滤波输出 y=fftshift(ifft(Y));%线性调频信号的时域匹配滤波输出 %对chirp信号进行时域加权 h1=(triang(lfft))';%三角窗函数 s1=s.*h1;S1=fft(s1);H1=conj(S1); Y1=S1.*H1; y1=fftshift(ifft(Y1));%加三角窗后的线性调频信号的时域匹配滤波输出h2=(hanning(lfft))';%汉宁窗函数 s2=s.*h2;S2=fft(s2);H2=conj(S2); Y2=S2.*H2; y2=fftshift(ifft(Y2));%加汉宁窗后的线性调频信号的时域匹配滤波输出h3=(hamming(lfft))';%海明窗函数 s3=s.*h3;S3=fft(s3);H3=conj(S3); Y3=S3.*H3; y3=fftshift(ifft(Y3));%加海明窗后的线性调频信号的时域匹配滤波输出 figure; subplot(3,1,1), plot(t,real(s)),title('chirp signal'); subplot(3,1,2), plot(q,abs(S)),title('线性调频信号幅度谱'); subplot(3,1,3), plot(q,angle(S)),title('线性调频信号相位谱');

实验一典型环节的MATLAB仿真汇总

实验一典型环节的MATLAB 仿真一、实验目的 1．熟悉MATLAB 桌面和命令窗口，初步了解SIMULINK 功能模块的使用方法。 2．通过观察典型环节在单位阶跃信号作用下的动态特性，加深对各典型环节响应曲线的理解。 3．定性了解各参数变化对典型环节动态特性的影响。二、SIMULINK 的使用 MATLAB 中SIMULINK 是一个用来对动态系统进行建模、仿真和分析的软件包。利用SIMULINK 功能模块可以快速的建立控制系统的模型，进行仿真和调试。 1．运行MATLAB 软件，在命令窗口栏“>>”提示符下键入simulink 命令，按Enter 键或在工具栏单击按钮，即可进入如图1-1所示的SIMULINK 仿真环境下。 2．选择File 菜单下New 下的Model 命令，新建一个simulink 仿真环境常规模板。 3．在simulink 仿真环境下，创建所需要的系统三、实验内容按下列各典型环节的传递函数，建立相应的SIMULINK 仿真模型，观察并记录其单位阶跃响应波形。 ① 比例环节1)(1=s G 和2)(1=s G 实验处理：1)(1=s G SIMULINK 仿真模型

波形图为：实验处理：2)(1=s G SIMULINK 仿真模型波形图为：实验结果分析：增加比例函数环节以后，系统的输出型号将输入信号成倍数放大. ② 惯性环节11)(1+= s s G 和15.01)(2+=s s G 实验处理：1 1 )(1+=s s G SIMULINK 仿真模型

波形图为：实验处理：1 5.01 )(2+= s s G SIMULINK 仿真模型波形图为：实验结果分析：当1 1 )(1+= s s G 时，系统达到稳定需要时间接近5s,当

Matlab与通信仿真课程设计

实验一单边带调幅系统的建模仿真班级：姓名：学号：一、实验目的 1.了解单边带调幅系统的工作原理 2.掌握单边带调幅系统的Matlab和Simulink建模过程二、实验内容 1、Matlab设计一个单边带发信机、带通信道和相应的接收机，参数要求如下。（1）输入话音信号为一个话音信号，采样率8000Hz。话音输入后首先进行预滤波，预滤波器是一个频率范围在[300,3400]Hz的带通滤波器。其目的是将话音频谱限制在3400Hz以下。单边带调制的载波频率设计为10KHz，调制输出上边带。要求观测单边带调制前后的信号功率谱。（2）信道是一个带限高斯噪声信道，其通带频率范围是[10000, 13500]Hz。要求能够根据信噪比SNR要求加入高斯噪声。（3）接收机采用相干解调方式。为了模拟载波频率误差对解调话音音质的影响，设本地载波频率为9.8KHz，与发信机载波频率相差200Hz。解调滤波器设计为300Hz到3400Hz的带通滤波器。程序框图：

设计思想：程序分为三部分：一：SSB调制模块首先从计算机中读入音频信号，作为原信号，读入完成后，对源信号进行参数采集和与滤波处理。进行与滤波之后，对信号进行希尔伯特变换，将原来的信号和载波相乘，将希尔伯特变换后的信号和载波进行希尔伯特变换后的信号相乘之后两者想减，得到SSB调制后的信号。二：信道加噪声模块通过信道，通过设置信道的信噪比来加入相应的噪声三：解调模块：将SSB调制后的信号通过信道加入噪声以后得到新的信号，并将信号和本地载波相乘进行想干解调，得到输出信号，并通过语句输出到相应的目录下。 2、用Simulink方式设计一个单边带传输系统并通过声卡输出接收机解调的结果声音。系统参数参照实例5.9，系统仿真参数设置为50KH 系统设计：

MATLAB仿真技术

MATLAB仿真技术作业合集

第1章习题 5.利用直接输入法和矩阵编辑器创建矩阵A=? ? ? ? ??642531。解：⑴利用直接输入法输入程序 A=[1 3 5;2 4 6] 按Enter 键后，屏幕显示 A = 1 3 5 2 4 6 ⑵用矩阵编辑器创建矩阵，如图1.1所示。图1.1 MATLAB 编辑器 7.用矩阵编辑器创建矩阵a,使a 具有如下矩阵形式。 a=??????642531?a=??????????654321?a=??????????987654321?a=???? ??????098706540321?a=????? ???????00 00 09870654 0321 解：用矩阵编辑器创建矩阵a 的过程如图1.2、1.3、1.4、1.5、1.6所示。图1.2 图1.3 图1.4 图1.5

图1.6 9.已知矩阵B=????? ?? ?????????922518113211912102201304161475231501017，试：①提取矩阵B 的第一行和第二行的第2、4、5个元素组成新矩阵1B ;②提取矩阵B 的第三行和第一行的全部元素组成新矩阵2B ；③使矩阵B 的第一行和第三行的第2；4个元素为0；④标出矩阵B 的第一行中小于5的元素。解：①如上题，用矩阵编辑器生成矩阵B ，再输入程序 B1=B([1,2],[2,4,5]) 按Enter 键后，屏幕显示 B1 = 0 0 15 5 14 16 ②输入程序 B2=B([1,3],:) 按Enter 键后，屏幕显示 B2 = 17 0 1 0 15 4 0 13 0 22 ③第一行和第三行的第2；4个元素原本就为0。 ④输入程序如下 C=B(1,:)<5; %将B 矩阵第一行中小于5 的值标记为1 D=B(1,C) %去B 矩阵第一行中标为1的元素按Enter 键后，屏幕显示 D= 0 1 0 11.已知矩阵a 为4阶魔方阵，令a+3赋值给b ，a+b 赋值给c ，求b 和c 。解：程序如下。 >> a=magic(4) %建立4阶魔方矩阵 a = 16 2 3 13 5 11 10 8 9 7 6 12 4 14 15 1 >> b=a+3 %将a 中各元素加3 b = 19 5 6 16 8 14 13 11 12 10 9 15 7 17 18 4