MATLAB对降压斩波器的仿真实验
苏洪（09电力牵引一班03）
（华东交通大学电气与电子工程学院，南昌）
摘要：利用matlab的simulink工具建立模拟电路对buck斩波电路进行仿真，并对输出结果进行分析，与理论的波形相仿。
关键词：matlab；simulink；buck；直流-直流变流电路
直流-直流变流电路的功能是将直流电
变为另一固定电压或可调电压的直流
电，包括直接直流变流电路和间接直
流变流电路。直接直流变流电路也称
斩波电路，他的功能是将直流电变为
另一固定电压或可调电压的直流电，
一般是指直接将直流电变成为另一直
流电，这种情况下输入与输出之间不
隔离。间接直流变流电路是在直流变
交流电路中增加了交流环节，在交流
环节中通常采用变压器实现输入输出
间的隔离，因此也称为带隔离的直流-
直流变流电路或直-交-直电路。习惯上
DC-DC变换器包括以上两种情况，且
甚至更多地指后一种情况。
1 降压斩波电路原理
降压斩波电路原理图
该电路使用一个全控型器件IGBT。为在V 关断时给负载中电感电流提供通道，设置了续流二极管VD。斩波电路主要用于电子电路的供电电源，也可拖动直流电动机或带蓄电池负载等，后两种情况下负载中均会出现反电动势，如图中EM。
降压斩波电路的工作波形如下
a）电流连续时的工作波形
b）电流断续时的工作波形
在t=0时刻驱动V导通，电源E向负载供电，
负载电压u0=E，负载电流i0=按指数曲线上
升。
当t=t1时刻，控制V关断，负载电流经二极
管VD续流，负载电压u0近似为零，负载电
t
t
O
O
O
b)
T
E
i
G
t
on
t
off
i
o i1i
2
I
10
I
20
t
1
u
o
O
O
O t
t
t
T
E E
i
G
i
G
t
on
t
off
i
o
t
x
i
1i
2
I
20
t
1
t
2
u
o
E
M
E
V
+
-
M
R
L
VD
a)
i
o
E
M
u
o
i
G
t
O
O
O
E
O
t
t
t
E M
i G
t
t
T
i G
t on
t off
i o i
1i
2
I10I20
t1
u o
a)b)
O
O T
E E
i G
t on
t off
i o t x
i1
i2
I20
t1t2
u o
流呈指数曲线下降。为了使负载电流连续且脉动小，通常使串联的电感L 值较大。 至一个周期T 结束，再驱动V 导通，重复上一周期的过程。当电路工作于稳态时，负载电流在一个周期的初值和终值相等。负载电压的平均值为
U 0=t on E/（t on +t off ）=t on E/T=aE
式中，t on 为V 处于通态的时间；t off 为V 处于断态的时间；T 为开关周期；a 为导通占空比。简称占空比或导通比。
2 利用MATLAB 进行仿真电路
根据电路原理图，选用适当的器件在simulink 中画出电路模拟图。a 为占空比，
b 为周期。
其中IGBT 的信号输出要用Bus selector 进行分解，得到管压和流经的电流。
参数的设置E=100V ，E m =30V ，负载中电阻R=1，L=0.01，C=inf 。Configuration Parameters 中，算法采用软件推荐的ode23tb ，扫描时间设置0到0.4秒，其余为
默认。在示波器设置中勾选save data to workspace ，同时取消掉上面的选项。 对此电路进行仿真，在matlab 命令窗口中输入subplot(4,1,1)
plot(buck1.time,buck1.signals(1,1).values) subplot(4,1,2)
plot(buck1.time,buck1.signals(1,2).values) subplot(4,1,3)
plot(buck1.time,buck1.signals(1,3).values) subplot(4,1,4)
plot(buck1.time,buck1.signals(1,4).values) 此时采用占空比为50%，T 为0.02s 。 得
到
如
下
的
图
形
：
其形式基本符合理论的波形。
将占空比调整为30%，
20%。波形如下
可以看出负载中电压明显减小，即占空比的调节可以调节输出电压平均值大小。以上的负载电流都是断续的，为了得到连续的电流，将L值改为0.1，50%占空比得到如下波形
。
经过本学期的多次实验，对该电路的仿真已经很熟练了，但是真正的去调试出正常的波形，还需要更多的时间。得到的结论虽然与理论的相比有些不整齐，但是基本的形式是没有问题的，对这些数据的处理还需要更多的实验来锻炼自己。
参考文献：
[1]王兆安，黄俊. 电力电子技术[M]. 4版，北京：
机械工业出版社，2000
例6.11 利用卡尔曼滤波估计一个未知常数 题目： 设已知一个未知常数x 的噪声观测集合，已知噪声v(n)的均值为零， 方差为 ，v(n)与x 不相关，试用卡尔曼滤波估计该常数 题目分析： 回忆Kalman 递推估计公式 由于已知x 为一常数，即不随时间n 变化，因此可以得到： 状态方程： x(n)=x(n-1) 观测方程： y(n)=x(n)+v(n) 得到A(n)=1,C(n)=1, , 将A(n)=1,代入迭代公式 得到：P(n|n-1)=P(n-1|n-1) 用P(n-1)来表示P(n|n-1)和P(n-1|n-1),这是卡尔曼增益表达式变为 从而 2v σ1??(|1)(1)(1|1)(|1)(1)(1|1)(1)()()(|1)()[()(|1)()()]???(|)(|1)()[()()(|1)](|)[()()](|1)H w H H v x n n A n x n n P n n A n P n n A n Q n K n P n n C n C n P n n C n Q n x n n x n n K n y n C n x n n P n n I K n C n P n n --=----=----+=--+=-+--=--2()v v Q n σ=()0w Q n =(|1)(1)(1|1)(1)()H w P n n A n P n n A n Q n -=----+21 ()(|1)[(|1)]v K n P n n P n n σ-=--+22(1)()[1()](1)(1)v v P n P n K n P n P n σσ-=--=-+
MATLAB 仿真实验报告 课题名称：MATLAB 仿真——图像处理 学院：机电与信息工程学院 专业：电子信息科学与技术 年级班级：2012级电子二班 一、实验目的 1、掌握MATLAB处理图像的相关操作，熟悉相关的函数以及基本的MATLAB语句。 2、掌握对多维图像处理的相关技能，理解多维图像的相关性质 3、熟悉Help 命令的使用，掌握对相关函数的查找，了解Demos下的MATLAB自带的原函数文件。 4、熟练掌握部分绘图函数的应用，能够处理多维图像。 二、实验条件
MATLAB调试环境以及相关图像处理的基本MATLAB语句，会使用Help命令进行相关函数查找 三、实验内容 1、nddemo.m函数文件的相关介绍 Manipulating Multidimensional Arrays MATLAB supports arrays with more than two dimensions. Multidimensional arrays can be numeric, character, cell, or structure arrays. Multidimensional arrays can be used to represent multivariate data. MATLAB provides a number of functions that directly support multidimensional arrays. Contents ： ●Creating multi-dimensional arrays 创建多维数组 ●Finding the dimensions寻找尺寸 ●Accessing elements 访问元素 ●Manipulating multi-dimensional arrays操纵多维数组 ●Selecting 2D matrices from multi-dimensional arrays从多维数组中选择二维矩 阵 (1)、Creating multi-dimensional arrays Multidimensional arrays in MATLAB are created the same way as two-dimensional arrays. For example, first define the 3 by 3 matrix, and then add a third dimension. The CAT function is a useful tool for building multidimensional arrays. B = cat(DIM,A1,A2,...) builds a multidimensional array by concatenating（联系起来）A1, A2 ... along the dimension DIM. Calls to CAT can be nested（嵌套）. (2)、Finding the dimensions SIZE and NDIMS return the size and number of dimensions of matrices. (3)、Accessing elements To access a single element of a multidimensional array, use integer subscripts（整数下标）. (4)、Manipulating multi-dimensional arrays
Matlab 在电磁场中的 应用 专业: 电气信息与自动化 班级：2012级自动化3班 学号：12012242065 学院：物电学院 指导老师：李虹 完成日期：2013年12月15日
Matlab 在电磁场中的应用 摘要 Matlab是美国Mathworks公司于80年代推出的大型数学软件，通过多年的升级换代，现在已发展成为集数值计算、符号计算、可视化功能以及诸多的工具箱为一体的大型科学计算软件，它已广泛应用于科研院所、工程技术等各个部门，并成为大学生、研究生必备的工具软件。 电磁学是物理学的一个分支，是研究电场和电磁的相互作用现象。电磁学从原来互相独立的两门科学(电学、磁学)发展成为物理学中一个完整的分支学科，主要是基于电流的磁效应和变化的磁场的电效应的发现。这两个实验现象，加上麦克斯韦关于变化电场产生磁场的假设，奠定了电磁学的整个理论体系，发展了对现代文明起重大影响的电工和电子技术。 针对电磁场学习理论性强、概念抽象等特点，利用Matlab强大的数值计算和图形技术，通过具体实例进行仿真，绘制相应的图形，使其形象化，便于对其的理解和掌握。将Matlab引入电磁学中，利用其可视化功能对电磁学实验现象进行计算机模拟，可以提高学习效率于学习积极性，使学习效果明显。 本文通过Matlab软件工具，对点电荷电场、线电荷产生的电位、平面上N 个电荷之间的库仑引力、仿真电荷在变化磁场中的运动等问题分别给出了直观形象的的仿真图，形实现了可视化学习，丰富了学习内容，提高了对电磁场理论知识的兴趣。 关键词：Matlab 电磁学仿真计算机模拟 一、点电荷电场 问题描述： 真空中，两个带正电的点电荷，在电量相同和电量不同情况下的电场分布。根据电学知识，若电荷在空间激发的电势分布为V，则电场强度等于电势梯度的
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);
卫星轨迹 一．问题提出 设卫星在空中运行的运动方程为： 其中是k 重力系数(k=401408km3/s)。卫星轨道采用极坐标表示，通过仿真，研究发射速度对卫星轨道的影响。实验将作出卫星在地球表面（r=6400KM ，θ=0）分别以v=8KM/s,v=10KM/s,v=12KM/s 发射时，卫星绕地球运行的轨迹。 二．问题分析 1．卫星运动方程一个二阶微分方程组，应用Matlab 的常微分方程求解命令ode45求解时，首先需要将二阶微分方程组转换成一阶微分方程组。若设，则有： 2．建立极坐标如上图所示，初值分别为：卫星径向初始位置，即地球半径：y(1,1)=6400；卫星初始角度位置：y(2,1)=0；卫星初始径向线速度：y(3,1)=0；卫星初始周向角速度：y(4,1)=v/6400。 3．将上述一阶微分方程及其初值带入常微分方程求解命令ode45求解，可得到一定时间间隔的卫星的径向坐标值y(1)向量；周向角度坐标值y(2)向量；径向线速度y(3)向量；周向角速度y(4)向量。 4．通过以上步骤所求得的是极坐标下的解，若需要在直角坐标系下绘制卫星的运动轨迹，还需要进行坐标变换，将径向坐标值y(1)向量；周向角度坐标值y(2)向量通过以下方程转换为直角坐标下的横纵坐标值X,Y 。 5．卫星发射速度速度的不同将导致卫星的运动轨迹不同，实验将绘制卫星分别以v=8KM/s ，v=10KM/s ，v=12KM/s 的初速度发射的运动轨迹。 三．Matlab 程序及注释 1．主程序 v=input('请输入卫星发射速度单位Km/s ：\nv=');%卫星发射速度输入。 axis([-264007000-1000042400]);%定制图形输出坐标范围。 %为了直观表达卫星轨迹，以下语句将绘制三维地球。 [x1,y1,z1]=sphere(15);%绘制单位球。 x1=x1*6400;y1=y1*6400;???????-=+-=dt d dt dr r dt d dt d r r k dt r d θ θθ2)(2 22222θ==)2(,)1(y r y ?????????????**-=**+*-===)1(/)4()3(2)4()4()4()1()1()1()3()4()2() 3()1(y y y dt dy y y y y y k dt dy y dt dy y dt dy ???*=*=)] 2(sin[)1(Y )]2(cos[)1(X y y y y
实验一 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.
Matlab 与电磁场模拟 一单电荷的场分布： 单电荷的外部电位计算公式： q φ= 4πε0r 等位线就是连接距离电荷等距离的点，在图上表示就是一圈一圈的圆，而电力线就是由点向 外辐射的线。 MATLAB 程序： theta=[0:.01:2*pi]'; r=0:10; x=sin(theta*r; y=cos(theta*r; plot(x,y,'b' x=linspace(-5,5,100; for theta=[-pi/4 0 pi/4] y=x*tan(theta; hold on ; plot(x,y; end grid on 单电荷的等位线和电力线分布图： 二多个点电荷的电场情况： 模拟一对同号点电荷的静电场 设有两个同号点电荷, 其带电量分别为 +Q1和+Q2(Q1、Q2>0 距离为 2a 则两 电荷在点P(x, y处产生的电势为: 由电场强度可得E = -?U, 在xOy 平面上, 电场强度的公式为: 为了简单起见, 对电势U 做如下变换:
。 Matlab 程序： q=1; xm=2.5; ym=2; x=linspace(-xm,xm; y=linspace(-ym,ym; [X,Y]=meshgrid(x,y; R1=sqrt((X+1.^2+Y.^2; R2=sqrt((X-1.^2+Y.^2; U=1./R1+q./R2; u=1:0.5:4; figure contour(X,Y,U,u grid on legend(num2str(u' hold on
plot([-xm;xm],[0;0] plot([0;0],[-ym;ym] plot(-1,0,'o' , 'MarkerSize' ,12 plot(1,0,'o' , 'MarkerSize' ,12 [DX,DY] = gradient(U; quiver(X,Y,-DX,-DY; surf(X,Y,U; 同号电荷的静电场图像为： 50 40 30 20 10 0-2 2
实验1：上采样与内插 一、实验目的 1、了解上采样与内插的基本原理和方法。 2、掌握上采样与内插的matlab程序的设计方法。 二、实验原理 上采样提高采样频率。上采样使得周期降低M倍，即新采样周期Tu和原有采样周期Ts的关系是T u=T s/M,根据对应的连续信号x(t),上采样过程从原有采样值x(kT s)生成新采样值x(kT u)=x(kT s/M)。操作的结果是在每两个采样值之间放入M-1个零值样点。 更实用的内插器是线性内插器，线性内插器的脉冲响应定义如下： 上采样值x(kT u)=x(kT s/M)通过与线性内插器的脉冲响应的卷积来完成内插。 三、实验内容 仿真正弦波采样和内插，通过基本采样x（k）,用M=6产生上采样x u(k)，由M=6线性内插得到样点序列x i（k）。 四、实验程序 % File: c3_upsampex.m M = 6; % upsample factor h = c3_lininterp(M); % imp response of linear interpolator t = 0:10; % time vector tu = 0:60; % upsampled time vector x = sin(2*pi*t/10); % original samples xu = c3_upsamp(x,M); % upsampled sequence subplot(3,1,1) stem(t,x,'k.') ylabel('x') subplot(3,1,2) stem(tu,xu,'k.') ylabel('xu') xi = conv(h,xu); subplot(3,1,3) stem(xi,'k.') ylabel('xi') % End of script file. % File: c3_upsample.m function out=c3_upsamp(in,M)
《电磁场与电磁波》仿真实验 2016年11月 《电磁场与电磁波》仿真实验介绍 《电磁场与电磁波》课程属于电子信息工程专业基础课之一，仿真实验主要目的在于使学生更加深刻的理解电磁场理论的基本数学分析过程，通过仿真环节将课程中所学习到的理论加以应用。受目前实验室设备条件的限制，目前主要利用 MATLAB 仿真软件进行，通过仿真将理论分析与实际编程仿真相结合，以理论指导实践，提高学生的分析问题、解决问题等能力以及通过有目的的选择完成实验或示教项目，使学生进一步巩固理论基本知识，建立电磁场与电磁波理论完整的概念。 本课程仿真实验包含五个内容： 一、电磁场仿真软件——Matlab的使用入门 二、单电荷的场分布 三、点电荷电场线的图像 四、线电荷产生的电位 五、有限差分法处理电磁场问题 目录 一、电磁场仿真软件——Matlab的使用入门……………............................................... .4 二、单电荷的场分
布 (10) 三、点电荷电场线的图像 (12) 四、线电荷产生的电位 (14) 五、有限差分法处理电磁场问题 (17) 实验一电磁场仿真软件——Matlab的使用入门 一、实验目的 1. 掌握Matlab仿真的基本流程与步骤； 2. 掌握Matlab中帮助命令的使用。 二、实验原理 （一）MATLAB运算 1.算术运算 (1)．基本算术运算 MATLAB的基本算术运算有：＋(加)、－(减)、*(乘)、/(右除)、\(左除)、 ^(乘方)。
注意，运算是在矩阵意义下进行的，单个数据的算术运算只是 一种特例。 (2)．点运算 在MATLAB中，有一种特殊的运算，因为其运算符是在有关算术运算符前面加点，所以叫点运算。点运算符有.*、./、.\和.^。两矩阵进行点运算是指它们的对应元素进行相关运算，要求两矩阵的维参数相同。 例1：用简短命令计算并绘制在0≤x≦6范围内的sin(2x)、sinx2、sin2x。 程序：x=linspace(0,6) y1=sin(2*x),y2=sin(x.^2),y3=(sin(x)).^2; plot(x,y1,x, y2,x, y3) （二）几个绘图命令 1. doc命令：显示在线帮助主题 调用格式：doc 函数名 例如：doc plot，则调用在线帮助，显示plot函数的使用方法。 2. plot函数：用来绘制线形图形 plot(y),当y是实向量时，以该向量元素的下标为横坐标，元素值为纵坐标画出一条连续曲线，这实际上是绘制折线图。 plot(x,y),其中x和y为长度相同的向量，分别用于存储x坐标和y 坐标数据。 plot(x,y,s)
补充内容：模拟调制系统的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仿真实验
增量调制(DM)实验 一、实验目的 （1）进一步掌握MATLAB的应用。 （2）进一步掌握计算机仿真方法。 （3）学会用MATLAB软件进行增量调制(DM)仿真实验。 二、实验原理 增量调制是由PCM发展而来的模拟信号数字化的一种编码方式，它是PCM的一种特例。增量调制编码基本原理是指用一位编码，这一位码不是表示信号抽样值的大小，而是表示抽样幅度的增量特性，即采用一位二进制数码“1”或“0”来表示信号在抽样时刻的值相对于前一个抽样时刻的值是增大还是减小，增大则输出“1”码，减小则输出“0”码。输出的“1”，“0”只是表示信号相对于前一个时刻的增减，不表示信号的绝对值。 增量调制最主要的特点就是它所产生的二进制代码表示模拟信号前后两个抽样值的差别(增加、还是减少)而不是代表抽样值本身的大小，因此把它称为增量调制。在增量调制系统的发端调制后的二进制代码1和0只表示信号这一个抽样时刻相对于前一个抽样时刻是增加(用1码)还是减少(用0码)。收端译码器每收到一个1码，译码器的输出相对于前一个时刻的值上升一个量化阶，而收到一个0码，译码器的输出相对于前一个时刻的值下降一个量化阶。 增量调制(DM)是DPCM的一种简化形式。在增量调制方式下,采用1比特量化器,即用1位二进制码传输样值的增量信息,预测器是
一个单位延迟器,延迟一个采样时间间隔。预测滤波器的分子系数向量是[0,1],分母系数为1。当前样值与预测器输出的前一样值相比较,如果其差值大于零,则发1码,如果小于零则发0码。 三、实验内容 增量调制系统框图如图一所示,其中量化器是一个零值比较器,根据输入的电平极性,输出为 δ,预测器是一个单位延迟器,其输出为前一个采样时刻的解码样值,编码器也是一个零值比较器,若其输入为负值,则编码输出为0,否则输出为1。解码器将输入1,0符号转换为 δ,然后与预测值相加后得出解码样值输出,同时也作为预测器的输入 输入样值 e n e n =δsgn(e n ) 传输 n ) n n-1+δsgn(e n ) x n + - + + 预测输出 + n-1 + 预测输出 解码样值输出 x n-1 预测输入x n =x n-1+δsgn(e n ) 图一 增量调制原理框图 设输入信号为: x(t)=sin2π50t+0.5sin 2π150t 增量调制的采样间隔为1ms,量化阶距δ=0.4,单位延迟器初始值为0。建立仿真模型并求出前20个采样点使客商的编码输出序列以 解码 编码 二电平量化 单位延迟 单位 延迟
《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 仿真实验要求 1. 仿真题目 （1） 线性频谱搬移电路仿真 根据线性频谱搬移原理，仿真普通调幅波。 基本要求：载波频率为8kHz ，调制信号频率为400Hz ，调幅度为0.3；画出调制信号、载波信号、已调信号波形，以及对应的频谱图。 扩展要求1：根据你的学号更改相应参数和代码完成仿真上述仿真；载波频率改为学号的后5位，调制信号改为学号后3位，调幅度设为最后1位/10。（学号中为0的全部替换为1，例如学号2010101014，则载波为11114Hz ，调制信号频率为114，调幅度为0.4）。 扩展要求2：根据扩展要求1的条件，仿真设计相应滤波器，并获取DSB-SC 和SSB 的信号和频谱。 （2） 调频信号仿真 根据调频原理，仿真调频波。 基本要求：载波频率为30KHz ，调制信号为1KHz ，调频灵敏度32310f k π=??，仿真调制信号，瞬时角频率，瞬时相位偏移的波形。 扩展要求：调制信号改为1KHz 的方波，其它条件不变，完成上述仿真。 2. 说明 （1） 仿真的基本要求每位同学都要完成，并且记入实验基本成绩。 （2） 扩展要求可以选择完成。
1.0 >> ma = 0.3; >> omega_c = 2 * pi * 8000; >> omega = 2 * pi * 400; >> t = 0 : 5 / 400 / 1000 : 5 / 400; >> u_cm = 1; >> fc = cos(omega_c * t); >> fa = cos(omega * t); >> u_am = u_cm * (1 + fa).* fc; >> U_c =fft(fc,1024); >> U_o =fft(fa,1024); >> U_am =fft(u_am, 1024); >> figure(1); >> subplot(321);plot(t, fa, 'k');title('调制信号');grid;axis([0 2/400 -1.5 1.5]); >> subplot(323);plot(t, fc, 'k');title('高频载波');grid;axis([0 2/400 -1.5 1.5]); >> subplot(325);plot(t, u_am, 'k');title('已调信号');grid;axis([0 2/400 -3 3]); >> fs = 5000; >> w1 = (0:511)/512*(fs/2)/1000; >> subplot(322);plot(w1, abs([U_am(1:512)']),'k');title('调制信号频谱');grid;axis([0 0.7 0 500]); >> subplot(324);plot(w1, abs([U_c(1:512)']),'k');title('高频载波频谱');grid;axis([0 0.7 0 500]); >> subplot(326);plot(w1, abs([U_am(1:512)']),'k');title('已调信号频谱');grid;axis([0 0.7 0 500]); 1.1 >> ma = 0.8; >> omega_c = 2 * pi * 11138; >> omega = 2 * pi * 138; >> t = 0 : 5 / 400 / 1000 : 5 / 400; >> u_cm = 1; >> fc = cos(omega_c * t);
Matlab 与电磁场模拟 一 单电荷的场分布： 单电荷的外部电位计算公式： 等位线就是连接距离电荷等距离的点，在图上表示就是一圈一圈的圆，而电力线就是由点向外辐射的线。 MATLAB 程序： theta=[0:.01:2*pi]'; r=0:10; x=sin(theta)*r; y=cos(theta)*r; plot(x,y,'b') x=linspace(-5,5,100); for theta=[-pi/4 0 pi/4] y=x*tan(theta); hold on ; plot(x,y); end grid on 单电荷的等位线和电力线分布图： r q 04πεφ=
二多个点电荷的电场情况： 模拟一对同号点电荷的静电场 设有两个同号点电荷,其带电量分别为+Q1和+Q2(Q1、Q2>0 )距离为2a则两电荷在点P(x, y)处产生的电势为: 由电场强度可得E = -?U,在xOy平面上,电场强度的公式为: 为了简单起见,对电势U做如下变换: 。 Matlab程序：
q=1; xm=; ym=2; x=linspace(-xm,xm); y=linspace(-ym,ym); [X,Y]=meshgrid(x,y); R1=sqrt((X+1).^2+Y.^2); R2=sqrt((X-1).^2+Y.^2); U=1./R1+q./R2; u=1::4; figure contour(X,Y,U,u) grid on legend(num2str(u')) hold on plot([-xm;xm],[0;0]) plot([0;0],[-ym;ym]) plot(-1,0,'o','MarkerSize',12) plot(1,0,'o','MarkerSize',12) [DX,DY] = gradient(U); quiver(X,Y,-DX,-DY); surf(X,Y,U); 同号电荷的静电场图像为：
MATLA仿真实验报告 学院：计算机与信息学院 课程：—随机信号分析 姓名： 学号： 班级： 指导老师： 实验一
题目：编写一个产生均值为1,方差为4的高斯随机分布函数程序, 求最大值，最小值，均值和方差，并于理论值比较。 解：具体的文件如下，相应的绘图结果如下图所示 G仁random( 'Normal' ,0,4,1,1024); y=max(G1) x=mi n(G1) m=mea n(G1) d=var(G1) plot(G1);
实验二 题目：编写一个产生协方差函数为CC)=4e":的平稳高斯过程的程序，产生样本函数。估计所产生样本的时间自相关函数和功率谱密度，并求统计自相关函数和功率谱密度，最后将结果与理论值比较。 解：具体的文件如下，相应的绘图结果如下图所示。 N=10000; Ts=0.001; sigma=2; beta=2; a=exp(-beta*Ts); b=sigma*sqrt(1-a*a); w=normrnd(0,1,[1,N]); x=zeros(1,N); x(1)=sigma*w(1); for i=2:N x(i)=a*x(i-1)+b*w(i); end %polt(x); Rxx=xcorr(x0)/N; m=[-N+1:N-1]; Rxx0=(sigma A2)*exp(-beta*abs(m*Ts)); y=filter(b,a,x) plot(m*Ts,RxxO, 'b.' ,m*Ts,Rxx, 'r');
periodogram(y,[],N,1/Ts); 文件旧硯化）插入（1〕 ZMCD 克闻〔D ］窗口曲） Frequency (Hz) 50 100 150 200 250 300 350 400 450 500 NH---.HP)&UO 二 balj/ 」- □歹
《现代机械工程基础实验》之机械工程控制基础综合实验报告 姓名 学号 班级 山东建筑大学机电工程学院 2012.06.04~06
第一部分 Matlab 编程方法及仿真实验 实验1. 三维曲面的绘制（略） 实验2. 系统零极点绘制例：求部分分式展开式和)(t g 一个线性定常系统的传递函数是 1 5422 3)(2 3 ++++= s s s s s G （1） 使用MATLAB 建立传递函数，并确定它的极点和零点，写出)(s G 的部分分式展开式并绘制 系统的脉冲响应。 实验结果：零点-0.6667 极点-0.8796 + 1.1414i -0.8796 - 1.1414i -0.2408 实验3. 系统的阶跃响应 例. )(s G 的阶跃响应 对例2中由（1）式给出的传递函数)(s G ，增加一个0=s 处的极点，使用impulse 命令绘制其拉普拉斯反变换式曲线，得到阶跃响应图。将该响应与对)(s G 使用step 命令所得到的响应比较，确定系统的DC 增益。利用初值定理和终值定理来校验结果。 实验结果：DC 增益= 2
实验4. 双输入反馈系统单位阶跃响应 考虑一个如图1所示的反馈系统，它既有参考输入也有干扰输入，其中对象和传感器的传递函数是 )12)(15.0(4)(++=s s s G p ，105.01 )(+=s s H 控制器是一个增益为80，有一个在3-=s 处的零点，极点/零点比15=α超前控制器。推导 两个独立的MATLAB 模型，其中一个模型的输入为)(s R ，另一个输入为)(s D 。使用这些模型确定闭环零点和极点，并在同一坐标系内绘制它们的阶跃响应。 D (s ) 图1 具有参考和干扰输入的反馈系统方框图 实验结果： 参考输入的CL 极点：-49.3658 -7.3336 + 7.9786i -7.3336 - 7.9786i -3.4670 参考输入的DC 增益：320 干扰输入的CL 零点：-45 干扰输入的CL 极点：-49.3658 -7.3336 + 7.9786i -7.3336 - 7.9786i -3.4670 干扰输入的DC 增益：4 -20
实验三：等量异号点电荷的电势分布 一、实验目的与要求 1.掌握命令窗口中直接输入语句，进行编程绘制等量异号点电荷的电势分布图； 2.掌握二维网格和三维曲面绘图的语句。 二、实验类型 设计 三、实验原理及说明 这里在命令窗口中直接输入简单的语句进行编程设计。MATLAB有几千个通用和专用 五、实验内容和步骤 （一）建立等量异号点电荷的电势方程
物理情景是oxy平面上在x=2,y=0处有一正电荷，x= -2,y=0处有一负电荷，根据 计算两点电荷电场中电势的分布，由于 （二）利用MA TLAB的函数, 绘制等量异号点电荷的电势分布图 首先选定一系列的x和y后，组成了平面上的网络点，再计算对应每一点上的z值。例如-5:0.2:5,-4:0.2:4分别是选取横坐标与纵坐标的一系列数值，meshgrid是生成数据网格的命令，[x,y]是xy平面上的坐标网格点。z是场点(x ,y)的电势，要求写出z的表达式。这里用到MA TLAB的函数mesh（）描绘3D网格图，meshgrid（）描绘在3D图形上加坐标网格，sqrt（）求变量的平方根。mesh（）是三维网格作图命令，mesh(x,y,z)画出了每一个格点(x,y)上对应的z值（电势）。在命令窗口中直接输入简单的语句，如下。 解1
解2
当场点即在电荷处时，会出现分母为零的情况，因此在r里加了一个小量0.01，这样既可以完成计算，又不会对结果的正确性造成太大影响。 另外需要注意的是表达式中的“./ ”、“.^ ”是对数组运算的算符，含义与数值运算中的“./ ”、“.^ ”相同，不同之处是后者只对单个数值变量进行运算，而前者对整个数组变量中的所有元素同时进行运算。 解2为了减少计算量，增加精确度，与先前的示例相比，计算范围由原先的-5 题目: 数字信号处理MATLAB仿真实验 姓名 学院 专业 班级 学号 班内序号 实验一：数字信号的 FFT 分析 1、实验内容及要求 (1) 离散信号的频谱分析： 设信号 此信号的0.3pi 和 0.302pi 两根谱线相距很近，谱线 0.45pi 的幅度很小，请选择合适的序列长度 N 和窗函数，用 DFT 分析其频谱，要求得到清楚的三根谱线。 (2) DTMF 信号频谱分析 用计算机声卡采用一段通信系统中电话双音多频（DTMF ）拨号数字 0～9的数据，采用快速傅立叶变换（FFT ）分析这10个号码DTMF 拨号时的频谱。 2、实验目的 通过本次实验，应该掌握： (a) 用傅立叶变换进行信号分析时基本参数的选择。 (b) 经过离散时间傅立叶变换（DTFT ）和有限长度离散傅立叶变换（DFT ） 后信号频谱上的区别，前者 DTFT 时间域是离散信号，频率域还是连续的，而 DFT 在两个域中都是离散的。 (c) 离散傅立叶变换的基本原理、特性，以及经典的快速算法（基2时间抽选法），体会快速算法的效率。 (d) 获得一个高密度频谱和高分辨率频谱的概念和方法，建立频率分辨率和时间分辨率的概念，为将来进一步进行时频分析（例如小波）的学习和研究打下基础。 (e) 建立 DFT 从整体上可看成是由窄带相邻滤波器组成的滤波器组的概念，此概念的一个典型应用是数字音频压缩中的分析滤波器，例如 DVD AC3 和MPEG Audio 。 3、程序代码 (1) N=5000; n=1:1:N; x=0.001*cos(0.45*pi*n)+sin(0.3*pi*n)-cos(0.302*pi*n-pi/4); y=fft(x,N); magy=abs(y(1:1:N/2+1)); k=0:1:N/2; w=2*pi/N*k; stem(w/pi,magy) axis([0.25,0.5,0,50]) （2） column=[1209,1336,1477,1633]; line=[697,770,852,941]; fs=10000; N=1024; 00010450303024().*cos(.)sin(.)cos(.)x n n n n ππππ=+-- MATLAB的实验仿真 目录 实验一MATLAB在控制系统模型建立与仿真中的应用 (1) 实验二典型系统的时域响应分析 (13) 实验三线性控制系统的根轨迹与频域分析 (17) 实验四线性系统的校正 (22) 附录一 MATLAB6.5 控制系统工具箱函数和结构化的控制语句 (30) 附录二 SIMULINK 基本模块介绍 (34) 实验一 MATLAB 在控制系统模型建立与仿真中的应用 一、 MATLAB 基本操作与使用 1. 实验目的 1) 掌握MATLAB 仿真软件的安装及启动，熟悉 MATLAB 工作环境平台。 2) MATLAB 命令窗口，包括工具条以及菜单选项的使用；MATLAB 语言的基本规定，包括数值的表示、变量命名规定、基本运算符、预定义变量以及表达式等。 3) MATLAB 图形绘制功能、M 文件程序设计和线性控制系统传递函数模型的建立等。 2. 实验仪器 PC 计算机一台，MATLAB 软件1套 3. 实验内容 1) MATLAB 的启动 这里介绍MATLAB 装入硬盘后，如何创建MATLAB 的工作环境。 方法一 MATLAB 的工作环境由matlab.exe 创建，该程序驻留在文件夹matlab\bin\ 中。它的图标是 matlab 。只要从<我的电脑>或<资源管理器>中去找这个程序，然后双击此图标，就会自动创建如图1所示的MATLAB6.5 版的工作平台 。 图1 在英文Windows 平台上的MATLAB6.5 MATLAB 工作平台 方法二 假如经常使用MATLAB ，则可以在Windows 桌面上创建一个MATLAB 快捷方式图标。具体办法为： 把<我的电脑>中的 matlab 图标用鼠标点亮，然后直接把此图标拖到Windows 桌面上即可。此后，直接双击Windows 桌面上的matlab 图标，就可建立图1所示的 MATLAB 工作平台。 2) MATLAB 工作环境平台 桌面平台是各桌面组件的展示平台，默认设置情况下的桌面平台包括 6 个窗口，具体如下： ① MATLAB 窗口 Command Window北邮dsp软件matlab仿真实验报告
Matlab仿真实验教程
相关主题
文本预览