matlab仿真自由落体

通信仿真技术与实践上机作业一

（实例 1.1）试对空气中在重力作用下不同质量物体的下落过程进行建模和仿真。已知重力加速度g=9.8m/s^2，在初始时刻t0=0s时物体由静止开始坠落。考虑空气阻力的影响。

（1）建立数学模型

质量为m的物体在自由坠落过程中受到竖直向下的恒定重力和向上的空气阻力f 的作用，由牛顿第二定律，我们知道，重力mg，加速度a以及物体质量m之间的关系是：

mg-f=ma

f=k*(v^2) k=空气阻力系数，为一恒定值

a=g-k(v^2/m)

（2）数学模型的解析分析

v(t)=at

s(t)=21at^2

（3）根据数学模型建立计算机仿真模型（编程）

将方程转换为一种在自变量（时间）上的“递推”表达式

v(t+dt)=v(t)+dv=v(t)+adt

s(t+dt)=s(t)+ds=s(t)+v(t)dt

（4）执行仿真和结果分析

% 自由落体.m

% 模拟受到空气阻力的小球

g=9.8; % 重力加速度

a=g;

m=10;

k=0.5; % 空气阻力系数

v=0; % 设定初始速度条件

s=0; % 设定初始位移条件

t=0; % 设定起始时间

dt=0.1; % 设置计算步长

N=20; % 设置仿真递推次数. 仿真时间等于N与dt的乘积

for f=1:N

v=v+a*dt; % 计算新时刻的速度

a=a-k.*(v^2)./m; % 空气阻力f=k*(v^2)/m

s(f+1)=s(f)+v*dt; % 新位移

t(f+1)=t(f)+dt; % 时间更新

end

% 作图: 受空气阻力落体结果与自由落体结果对比

t=0:dt:N*dt;

subplot(1,2,2)

plot(t,s,'o');

xlabel('时间 t'); ylabel('位移 s');

legend('受空气阻力的落体');

运行得到的结果:

（5）仿真程序的功能扩展---以动态方式来观察物体坠落的过程

受到空气阻力落体动画.m

〔实例1.2〕对乒乓球的弹跳过程进行仿真。忽略空气对球的影响，乒乓球垂直下落，落点为光滑的水平面，乒乓球接触落点立即反弹。如果不考虑弹跳中的能量损耗，则反弹前后的瞬时速率不变，但方向相反。如果考虑撞击损耗，则反弹速率有所降低。我们希望通过仿真得出乒乓球位移随时间变化的关系曲线，并进行弹跳过程的“实时”动画显示。

（1）数学模型

首先对乒乓球弹跳过程进行一些理想化假设。设球是刚性的，质量为m，垂直下落。碰击面为水平光滑平面。在理想情况下碰击无能量损耗。如果考虑碰击面损耗，则碰击前后速度方向相反，大小按比例系0<=K<=1下降。在t时刻的速度设为v=v(t)，位移设为y=y(t)，并以碰击点为坐标原点，水平方向为坐标横轴建立直角坐标系。球体的速度以竖直向上方向为正方向。重力加速度为g=9.8m/s^2。

初始条件假设：设初始时刻t0=0球体的初始速度为v0=v(t0)，初始位移为y0=y(t0)。

受力分析：在空中时小球受重力F=mg作用，其中，g=-(dt/dv) 。则在t+dt时刻小球的速度为

v(t+dt)=v(t)-gdt 在t+dt时刻小球的位移为

y(t+dt)=y(t)+v(t)dt

在小球撞击水平面的瞬间，即y(t)=0的时刻，它的速度方向改变，大小按比例K衰减。

当K=1时，就是无损耗弹跳情况。因此，小球反弹瞬间（t+dt时刻）的速度为

v(t+dt)=-Kv(t)-gdt; 0<=K<=1

反弹瞬间的位移为

y(t+dt)=y(t)-Kv(t)dt=-Kv(t)dt

（2）仿真模型设计（程序）

从数学模型中可见，小球在空中自由运动时刻与撞击时刻的动力方程不同。通过小球所处位置（位移）是否为零可判定小球处于何种状态。程序文件代码如下。

乒乓球弹跳曲线.m

% 乒乓球弹跳曲线.m

g=9.8; % 重力加速度

v0=0; % 初始速度

y0=1; % 初始位置

m=1; % 小球质量

t0=0; % 起始时间

K=1; % 弹跳的损耗系数

N=5000; % 仿真的总步进数

dt=0.001; % 仿真步长

v=v0; % 初状态

y=y0;

for k=1:N

if (y>0)|(v>0) % 小球在空中的(含刚刚弹起瞬间)动力方程计算

v =v -g*dt;

y =y +v*dt;

else % 碰击瞬间的计算

y =y-K.*v*dt;

v =-K.*v-g*dt;

end

s(k)=y; % 将当前位移记录到s数组中以便作图

end

t=t0:dt:dt*(N-1);% 仿真时间长度

subplot(1,2,1);plot(t,s);

title('无损耗时的弹跳');

xlabel('时间 t');

ylabel('位移 y(t)');

axis([0 5 0 1.1]);

g=9.8; % 重力加速度

v0=0; % 初始速度

y0=1; % 初始位置

m=1; % 小球质量

t0=0; % 起始时间

K=0.85; % 弹跳的损耗系数

N=5000; % 仿真的总步进数

dt=0.001; % 仿真步长

v=v0; % 初状态

y=y0;

for k=1:N

if (y>0)|(v>0) % 小球在空中的(含刚刚弹起瞬间)动力方程计算 v =v -g*dt;

y =y +v*dt;

else % 碰击瞬间的计算

y =y-K.*v*dt;

v =-K.*v-g*dt;

end

s(k)=y; % 将当前位移记录到s数组中以便作图

end

t=t0:dt:dt*(N-1);% 仿真时间长度

subplot(1,2,2);

plot(t,s);

title('有损耗的弹跳曲线');

xlabel('时间 t');

ylabel('位移 y(t)');

axis([0 5 0 1.1]);

% 乒乓球仿真动画.m

实例1.3〕试用蒙特卡罗方法求出半径为1的圆的面积。并与理论值对比。（1）数学模型

设有两个相互独立的随机变量x,y，服从[0,2]上的均匀分布。那么由它们所确

定的坐标点(x;y)是均匀分布于边长为2的一个正方形区域中。如下图所示。

该正方形的内接圆的半径为1。显然，坐标点(x,y)落入圆中的概率p等于该圆面积Sc与正方形面积S之比，即 Sc=pS

因此，只要通过随机试验统计出落入圆中点的频度，即可计算出圆的近似面积来。当随机试验的次数充分大的时候，计算结果就趋近于理论真值。

（2）仿真试验

求圆周率1.m

% ch1example3prg1.m

sita=0:0.01:2*pi;

x=sin(sita);

y=cos(sita);% 计算半径为1的圆周上的点，以便作出圆周观察

m=0; % 在圆内在落点计数器

x1=2*rand(1000,1)-1;% 产生均匀分布于[-1, +1]直接的两个独立随机数x1,y1 y1=2*rand(1000,1)-1;

N=1000; % 设置试验次数

for n=1:N % 循环进行重复试验并统计

p1=x1(1:n);

q1=y1(1:n);

if (x1(n)*x1(n)+y1(n)*y1(n))<1 % 计算落点到坐标原点的距离,判别落点是否在圆内

m=m+1; % 如果落入圆中,计数器加1

end

plot(p1,q1,'.',x,y,'-k',[-1 -1 1 1 -1],[-1 1 1 -1 -1],'-k');

axis equal; % 坐标纵横比例相同

axis([-2 2 -2 2]); % 固定坐标范围

text(-1,-1.2,['试验总次数 n=',num2str(n)]);% 显示试验结果

text(-1,-1.4,['落入圆中数 m=',num2str(m)]);

text(-1,-1.6,['近似圆面积 S_c=',num2str(m/n*4)]);

set(gcf,'DoubleBuffer','on'); % 双缓冲避免作图闪烁

drawnow; % 显示结果

end

程序执行中，将动态显示随机落点情况和当前的统计计算结果。随着试验次数增加，计算结果将趋近于半径为1的圆面积的真值pi。

蒙特卡罗法求圆周率2.m

% ch1example3prg2.m

tic % 计时器启动

n=10000; % 每次随机落点10000个

for k=1:1000 % 重复试验1000次

x1=2*rand(n,1)-1;

y1=2*rand(n,1)-1; % 随机落点产生

m(k)=sum((x1.*x1+y1.*y1)<1);% 求落入圆中的点数和

end

S_c=mean(m).*4./n % 计算并显示结果

time=toc

S_c =

3.1414

time =

0.6886

这种算法计算速度更快，仿真的结果更准确

〔实例1.4〕实际物理试验中，当我们让一个乒乓球垂直下落到一个完全水平的玻璃板上后，乒乓球不断弹跳，直到能量耗尽。我们要建立更加接近真实的物理环境的弹跳模型，就必须考虑被忽视的微小的扰动因素, 根据大数定理，在数学上我们就可以将水平作用力建模为一个高斯随机变量。

（1）分析

在三维空间内仿真小球下落的过程：

ax(t)=Fx(t)/m

dvx(t)=ax(t)dt

dsx(t)=vx(t)dt

Y、Z平面也如此

（2）仿真代码

三维弹跳动画.m

%三维弹跳动画

g=9.8; % 重力加速度

v0=0; % 初始速度

y0=1; % 初始位置

m=0.1; % 小球质量

t0=0; % 起始时间

K=0.8; % 弹跳的损耗系数

N=5000; % 仿真的总步进数

dt=0.005;% 仿真步长

v=v0; % 初状态

y=y0;

vx=0;

vz=0;

sx=0;

sz=0;

for k=1:N

if y > 0 % 小球在空中的动力方程计算

v=v -g*dt;

y=y +v*dt;

else % 碰击瞬间的计算

y=-K.*v*dt;

v=-K.*v-g*dt;

end

Fx=randn; % x水平方向的随机力,方差为1

ax=Fx./m; % Fx导致的x水平方向的加速度

vx=vx+ax*dt; % 小球在x水平方向的瞬时速度

sx=sx+vx*dt; % 小球在x水平方向的位移

Fz=randn; % z水平方向的随机力,方差为1

az=Fz./m; % Fz导致的z水平方向的加速度

vz=vz+az*dt; % 小球在z水平方向的瞬时速度

sz=sz+vz*dt; % 小球在z水平方向的位移

plot3(sx,sz,y,'.'); grid on;hold on;

axis([-2 2 -2 2 0 1]); % 坐标范围固定

set(gcf,'DoubleBuffer','on'); % 双缓冲避免作图闪烁

xlabel('水平方向x');ylabel('水平方向z');zlabel('垂直方向y'); drawnow; % 作图

end

multisim buck电路仿真

第一章概述 1、1 直流―直流变换的分类 直流—直流变换器(DC-DC)就是一种将直流基础电源转变为其她电压种类的直流变换装置。目前通信设备的直流基础电源电压规定为?48V,由于在通信系统中仍存在?24V(通信设备)及+12V、+5V(集成电路)的工作电源,因此,有必要将?48V基础电源通过直流—直流变换器变换到相应电压种类的直流电源,以供实际使用。D C/DC变换就是将固定的直流电压变换成可变的直流电压,也称为直流斩波。主要有 (1)Buck电路——降压斩波,其输出平均电压小于输入电压,极性相同。 (2)Boost电路——升压斩波,其输出平均电压大于输入电压,极性相同。 (3)Buck－Boost电路——降压―升压斩波,其输出平均电压大于或小于输入电压,极性相反,电感传输。 (4)Cuk电路——降压或升压斩波,其输出平均电压大于或小于输入电压,极性相反,电容传输。 此外还有Sepic、Zeta电路。 1、2 直流—直流变换器的发展 当今软开关技术的发展使得DC/DC发生了质的飞跃,美国VICOR公司(美国怀格公司,国际知名的电源模块生产厂家)设计制造的多种ECI软开关DC/DC变换器,其最大输出功率有300W、600W、800W等,相应的功率密度为(6、2、10、17)W/cm3,效率为(80～90)%。日本NEMIC—LAMBDA(联美兰达,日本的开关电源厂商、2012年兰达被TDK收购,名称也改为TDK-LAMBDA)公司最新推出的一种采用软开关技术的高频开关电源模块RM系列,其开关频率为(200～300)kHz,功率密度已达到27W/cm3,采用同步整流器(MOSFET代替肖特基二极管),使整个电路效率提高到90%。

Matlab Simulink 仿真步骤

MATLAB基础与应用简明教程 张明等编著 北京航空航天大学出版社(2001.01) MATLAB软件环境是美国New Mexico大学的Cleve Moler博士首创的，全名为MATrix LABoratory（矩阵实验室）。它建立在20世纪七八十年代流行的LINPACK（线性代数计算）和ESPACK（特征值计算）软件包的基础上。LINPACK和ESPACK软件包是从Fortran语言开始编写的，后来改写为C语言，改造过程中较为复杂，使用不便。MA TLAB是随着Windows环境的发展而迅速发展起来的。它充分利用了Windows环境下的交互性、多任务功能语言，使得矩阵计算、数值运算变得极为简单。MA TLAB语言是一种更为抽象的高级计算机语言，既有与C语言等同的一面，又更为接近人的抽象思维，便于学习和编程。同时，它具有很好的开放性，用户可以根据自己的需求，利用MA TLAB提供的基本工具，灵活地编制和开发自己的程序，开创新的应用。 本书重点介绍了MA TLAB的矩阵运算、符号运算、图形功能、控制系统分析与设计、SimuLink仿真等方面的内容。 Chap1 MATLAB入门与基本运算 本章介绍MATLAB的基本概念，包括工作空间；目录、路径和文件的管理方式；帮助和例题演示功能等。重点介绍矩阵、数组和函数的运算规则、命令形式，并列举了可能得到的结果。由于MA TLAB的符号工具箱是一个重要分支，其强大的运算功能在科技领域有特殊的帮助作用。 1.1 MATLAB环境与文件管理 1.2 工作空间与变量管理 1.2.1 建立数据 x1=[0.2 1.11 3]; y1=[1 2 3;4 5 6]建立一维数组x1和二维矩阵y1。分号“；”表示不显示定义的数据。 MATLAB还提供了一些简洁方式，能有规律地产生数组： xx=1:10 %xx从1到10，间隔为1 xx=-2:0.5:1 %xx从-2到1，间隔为0.5 linespace命令等距离产生数组，logspace在对数空间中等距离产生数组。对于这一类命令，只要给出数组的两端数据和维数就可以了。 xx=linespace(d1,d2,n) %表示xx从d1到d2等距离取n个点 xx=logspace(d1,d2,n) %表明xx从10d1到10d2等距离取n个点 1.2.2 who和whos命令 who: 查看工作空间中有哪些变量名 whos: 了解这些变量的具体细节 1.2.3 exist命令 查询当前的工作空间内是否存在一个变量，可以调用exist()函数来完成。 调用格式：i=exist(…A?); 式中，A为要查询的变量名。返回的值i表示A存在的形式： i=1 表示当前工作空间内存在一个变量名为A的矩阵； i=2 表示存在一个名为A.m的文件； i=3 表示MATLAB的工作路径下存在一个名为A.mex的文件；

MATLAB仿真实验报告

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

Buck电路的设计与仿真

uck 电路的设计与仿真 1、Buck 电路设计： 设计一降压变换器，输入电压为 20V ，输出电压5V ，要求纹波电压为输出 电压的0.5%，负载电阻10欧姆，求工作频率分别为10kHz 和50kHz 时所需的 电感、电容。比较说明不同开关频率下，无源器件的选择。 解：(1)工作频率为10kHz 时， A.主开关管可使用MOSFET ,开关频率为10kHz ; B 输入20V ，输出5V ，可确定占空比 Dc=25% ; C.根据如下公式选择电感 这个值是电感电流连续与否的临界值，L>L c 则电感电流连续，实际电感值 可选为1.1~1.2倍的临界电感，可选择为4 10?H ； D.根据纹波的要求和如下公式计算电容值 C=^^T s2 J =4.17 10 牛 8L^U 。 8 沃 4.5 沃 10 X0.0055 1 0000 (2)工作频率为50kHz 时， A.主开关管可使用MOSFET ,开关频率为50kHz ; B 输入20V ，输出5V ，可确定占空比 Dc=25% ; C.根据如下公式选择电感 . (1—DJR T (1 —0.25)汇10,. 1 L c (1 _DJR T 2 s (1-0.25)1° 亠 2 10000 = 3.75 10* H 5 (1-0.25) 0.75 10, H 50000 这个值是电感电流连续与否的临界值, L>Lc 则电感电流连续，实际电感值

L c T s 2

可选为1.2倍的临界电感，可选择为0.9 10" H ； D.根据纹波的要求和如下公式计算电容值 分析：在其他条件不变的情况下，若开关频率提高 n 倍，则电感值减小为 1/n ，电容值也减小到1/n 。从上面推导中也得出这个结论 2、Buck 电路仿真: 利用sim power systems 中的模块建立所设计降压变换器的仿真电路。输 入电压为20V 的直流电压源，开关管选 MOSFET 模块（参数默认），用Pulse Gen erator 模块产生脉冲驱动开关管 建模: 分别做两种开关频率下的仿真 工作频率为10kHz 时 U o (1-D c ) 8L U o T s 2 5 (1-0.25) 1 8 0.9 10J 0.005 5 500002 = 0.833 10*F matlab20120510 ?

电力电子技术MATLAB仿真报告模板

《电气专业核心课综合课程设计》 题目：基于MATLAB的电力电子技术 仿真分析 学校： 院（系）： 专业班级： 学生姓名： 学号： 指导教师： 目录

1.整流电路仿真………………………………………………………………………………页码 1.1单相半波可控整流系统………………………………………………………………页码 1.1.1晶闸管的仿真…………………………………………………………………页码 1.1.2单相半波可控整流电路的仿真………………………………………………页码 1.2晶闸管三相桥式整流系统的仿真…………………………………………………页码 1.3相位控制的晶闸管单相交流调压器带系统的仿真………………………………页码 2.斩波电路仿真………………………………………………………………………………页码 2.1降压斩波电路（Buck变换器）………………………………………………………页码 2.1.1可关断晶闸管（GTO）的仿真…………………………………………………页码 2.1.2 Buck变换器的仿真………………………………………………………页码 2.2升压斩波电路（Boost变换器）………………………………………………………页 码 2.2.1绝缘栅双极型晶体管（IGBT）的仿真…………………………………………页码 2.2.2 Boost变换器的仿真……………………………………………………………页码4.逆变电路仿真………………………………………………………………………………页码 4.1晶闸管三相半波有源逆变器的仿真………………………………………………页码 5.课程设计总结………………………………………………………………………………页码参考文献……………………………………………………………………………………页码 电气专业核心课综合课程设计任务书

题目Buck电路的设计与仿真

题目：Buck 电路的设计与仿真 1、Buck 电路设计： 设计一降压变换器，输入电压为20V ，输出电压5V ，要求纹波电压为输出电压的0.5％，负载电阻10欧姆，求工作频率分别为10kHz 和50kHz 时所需的电感、电容。比较说明不同开关频率下，无源器件的选择。 解：（1）工作频率为10kHz 时， A.主开关管可使用MOSFET ，开关频率为10kHz ； B.输入20V ，输出5V ，可确定占空比Dc=25%； C.根据如下公式选择电感 H T R D L s c c 41075.310000 1210)25.01(2)1(-?=??-=-= 这个值是电感电流连续与否的临界值，L>c L 则电感电流连续，实际电感值可选为1.2倍的临界电感，可选择为H 4105.4-?； D.根据纹波的要求和如下公式计算电容值 =?-=2008)1(s c T U L D U C 2410000 15005.0105.48)25.01(5?????-?-=F 41017.4-? （2）工作频率为50kHz 时， A.主开关管可使用MOSFET ，开关频率为50kHz ； B.输入20V ，输出5V ，可确定占空比Dc=25%； C.根据如下公式选择电感 H T R D L s c c 41075.050000 1210)25.01(2)1(-?=??-=-= 这个值是电感电流连续与否的临界值，L>Lc 则电感电流连续，实际电感值可选为1.2倍的临界电感，可选择为H 4109.0-?； D.根据纹波的要求和如下公式计算电容值 =?-=2008)1(s c T U L D U C 2450000 15005.0109.08)25.01(5?????-?-=F 410833.0-? 分析： 在其他条件不变的情况下，若开关频率提高n 倍，则电感值减小为1/n ，电容值也减小到1/n 。从上面推导中也得出这个结论。 2、Buck 电路仿真： 利用simpowersystems 中的模块建立所设计降压变换器的仿真电路。输入电压为20V 的直流电压源，开关管选MOSFET 模块(参数默认)，用Pulse Generator 模块产生脉冲驱动开关管。分别做两种开关频率下的仿真。 （一）开关频率为10Hz 时； (1)使用理论计算的占空比，记录直流电压波形，计算稳态直流电压值，计算稳态直流纹波电压，并与理论公式比较，验证设计指标。 由第一步理论计算得占空比Dc=25%； 实验仿真模型如下所示（稳态直流电压值为4.299V ）：

MATLAB电路仿真报告

课程设计 题目电路仿真 学院自动化学院 专业自动化专业 班级自动化0806班 姓名孙功武 指导教师徐腊梅 年月日

摘要 《电路原理》是电类专业必修的一门重要的技术基础课，它具有基础科学和技术科学的二重性，不仅是电类学生学习后续课程的基础，也直接为解决电工电子工程中的一些实际问题服务。大一下学期开始，通过对本课程的学习，我初步掌握了近代电路理论的一些基本知识和概念，能分析计算一些常见的，比较简单的基本电路，初步具有了解决实际问题的能力，并为后续课程的学习准备了必要的电路理论知识。其分析电路的常见方法有：节点电压法，网孔电流法，叠加原理分析法，戴维宁定理和诺顿定理等等。本文主要讨论用网口电流法来分析直流电路中关于电阻电路的计算方法。在这个分析解决问题的过程中需要运用到MATLAB软件。MATLAB是矩阵实验室（Matrix Laboratory）的简称，是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。本文就是通过对MATLAB编程计算出的结果和Simulink仿真出的结果进行对比，来的出所要的结论和效果。 关键词：电路原理，网孔电流法，MATLAB，SIMULINK，

Abstract The circuit principle of electricity class specialized is compulsory course is an important technology, it has the basic science and technology is not only the scientific duality, electricity class student learning courses, and subsequent direct solution for electrical and electronic engineering of some actual problems. A semester began, this course of study, I have mastered some modern circuit theory, the elementary knowledge and the concept of some common to analysis and calculation, the basic circuit is simple, is the ability to solve practical problems, and for subsequent course of study prepared necessary circuit theory knowledge. The analysis of the common method: circuit node voltage, current, mesh superposition principle analysis method, DaiWeiNing theorem and NORTON's theorem, etc. This paper discusses how to use the WangKou current method to analyze the dc resistance circuit in the calculation. In the analysis and problem solving process needs to apply MATLA B software. MATLAB is Matrix lab (Matrix of Laboratory), is the MathWorks company business mathematics software is used to develop, data visualization algorithm, data analysis and numerical calculation of senior technical calculation language and interactive environment, including MATLAB and Simulink two most. This paper is based on MATLAB calculation results and Simulink results, comparing to the conclusion of the effect.

multisimbuck电路仿真设计

第一章概述 1.1 直流―直流变换的分类 直流—直流变换器（DC-DC）是一种将直流基础电源转变为其他电压种类的直流变换装置。目前通信设备的直流基础电源电压规定为?48V，由于在通信系统中仍存在?24V（通信设备）及+12V、+5V（集成电路）的工作电源，因此，有必要将?48V基础电源通过直流—直流变换器变换到相应电压种类的直流电源，以供实际使用。D C/DC变换是将固定的直流电压变换成可变的直流电压，也称为直流斩波。主要有 （1）Buck电路——降压斩波，其输出平均电压小于输入电压，极性相同。 （2）Boost电路——升压斩波，其输出平均电压大于输入电压，极性相同。 （3）Buck－Boost电路——降压―升压斩波，其输出平均电压大于或小于输入电压，极性相反,电感传输。 （4）Cuk电路——降压或升压斩波，其输出平均电压大于或小于输入电压，极性相反，电容传输。 此外还有Sepic、Zeta电路。 1.2 直流—直流变换器的发展 当今软开关技术的发展使得DC/DC发生了质的飞跃，美国VICOR公司(美国怀格公司,国际知名的电源模块生产厂家)设计制造的多种ECI软开关DC/DC变换器，其最大输出功率有300W、600W、800W等，相应的功率密度为(6.2、10、17)W/cm3，效率为（80～90）%。日本NEMIC—LAMBDA(联美兰达,日本的开关电源厂商.2012年兰达被TDK收购,名称也改为TDK-LAMBDA)公司最新推出的一种采用软开关技术的高频开关电源模块RM系列，其开关频率为（200～300)kHz，功率密度已达到27W/cm3，采用同步整流器（MOSFET代替肖特基二极管），使整个电路效率提高到90%。

MATLAB通信系统仿真实验报告1

MATLAB通信系统仿真实验报告

实验一、MATLAB的基本使用与数学运算 目的：学习MATLAB的基本操作，实现简单的数学运算程序。 内容： 1-1要求在闭区间[0，2π]上产生具有10个等间距采样点的一维数组。试用两种不同的指令实现。 运行代码：x=[0:2*pi/9:2*pi] 运行结果： 1-2用M文件建立大矩阵x x=[0.10.20.30.40.50.60.70.80.9 1.11.21.31.41.51.61.71.81.9 2.12.22.32.42.52.62.72.82.9 3.13.23.33.43.53.63.73.83.9] 代码：x=[0.10.20.30.40.50.60.70.80.9 1.11.21.31.41.51.61.71.81.9 2.12.22.32.42.52.62.72.82.9 3.13.23.33.43.53.63.73.83.9] m_mat 运行结果： 1-3已知A=[5，6;7，8]，B=[9，10;11,12]，试用MATLAB分别计算 A+B,A*B,A.*B,A^3，A.^3，A/B,A\B. 代码：A=[56;78]B=[910;1112]x1=A+B X2=A-B X3=A*B X4=A.*B X5=A^3 X6=A.^3X7=A/B X8=A\B

运行结果： 1-4任意建立矩阵A，然后找出在[10,20]区间的元素位置。 程序代码及运行结果： 代码：A=[1252221417;111024030;552315865]c=A>=10&A<=20运行结果： 1-5总结：实验过程中，因为对软件太过生疏遇到了些许困难，不过最后通过查书与同学交流都解决了。例如第二题中，将文件保存在了D盘，而导致频频出错，最后发现必须保存在MATLAB文件之下才可以。第四题中，逻辑语言运用到了ij，也出现问题，虽然自己纠正了问题，却也不明白错在哪了，在老师的讲解下知道位置定位上不能用ij而应该用具体的整数。总之第一节实验收获颇多。

Buck电路设计与MATLAB仿真

Buck电路设计与仿真 姓名：朱龙胜 班级：电气1102 学号：11291065 日期：2014年5月10日 指导老师：郭希铮 北京交通大学

计算机仿真技术作业四 题目：Buck 电路的设计与仿真 1、Buck 电路设计： 设计一降压变换器，输入电压为20V ，输出电压5V ，要求纹波电压为输出电压的0.5％，负载电阻10欧姆，求工作频率分别为10kHz 和50kHz 时所需的电感、电容。比较说明不同开关频率下，无源器件的选择。 2、Buck 电路理论计算: 由以下公式计算: 20.252.0.5A (1) 3.5% 8() 4.2o d o o o s o s d o LB OB V D V V I R V T D V LC DT V V I I L = == =?-==-== 1.占空比: 负载电流: 纹波电压: 电流连续条件： 得到下列计算结果 3、Buck 电路仿真: 利用simpowersystems 中的模块建立所设计降压变换器的仿真电路。输入电压为20V 的直流电压源，开关管选MOSFET 模块(参数默认)，用Pulse Generator 模块产生脉冲驱动开关管。分别做两种开关频率下的仿真。 (1)使用理论计算的占空比(D=0.25)，记录直流电压波形，计算稳态直流电压值，计算稳态直流纹波电压，并与理论公式比较，验证设计指标。 4、仿真过程:: A ．建立模型: 建立仿真模型如下如所示 :

B. 记录数据: 仿真算法选择ode23tb,最大步长为0.1s ，占空比D=0.25进行仿真，记录数据如下表所 C ．仿真过程: 当f s =10KHz,L=0.375mH C=500μF, 占空比D=0.25,电流连续的临界状态时，记录稳态直流电压值V o =4.736V ，稳态直流电压理论值5V 计算稳态直流纹波电压的理论值 2(1D)0.025V 8s o o T V V CL -?==,通过图中得到直流纹波电压为0.0267V 当fs=10KHz,L=0.375mH, C=500μF,占空比D=0.25,电流连续的临界状态时, 由(1)o S L V D T I L -?= ，得电感电流波动理论值是1A ，由图像得到电感电流波动值是 1A ，与理论计算相符合

simulink-matlab仿真教程

simulink matlab 仿真环境教程 Simulink 是面向框图的仿真软件。 演示一个Simulink 的简单程序 【例1.1】创建一个正弦信号的仿真模型。 步骤如下： (1) 在MATLAB 的命令窗口运行simulink 命令，或单击工具栏中的图标，就可以打开Simulink 模块库浏览器 (Simulink Library Browser) 窗口，如图1.1所示。 (2) 单击工具栏上的图标或选择菜单“File ”——“New ”——“Model ”，新建一个名为“untitled ”的空白 模型窗口。 (3) 在上图的右侧子模块窗口中，单击“Source ”子模块库前的“+”(或双击Source)，或者直接在左侧模块和工具箱栏单击Simulink 下的Source 子模块库，便可看到各种输入源模块。 (4) 用鼠标单击所需要的输入信号源模块“Sine Wave ”(正弦信号)，将其拖放到的空白模型窗口“untitled ”，则“Sine Wave ”模块就被添加到untitled 窗口；也可以用鼠标选中“Sine Wave ”模块，单击鼠标右键，在快捷菜单中选择“add to 'untitled'”命令，就可以将“Sine Wave ”模块添加到untitled 窗口，如图1.2 所示。 图7.1 Simulink 界面

(5) 用同样的方法打开接收模块库“Sinks”，选择其中的“Scope ”模块(示波器)拖放到“untitled”窗口中。 (6) 在“untitled”窗口中，用鼠标指向“Sine Wave”右侧的输出端，当光标变为十字符时，按住鼠标拖向“Scope”模块的输入端，松开鼠标按键，就完成了两个模块间的信号线连接，一个简单模型已经建成。如图1.3所示。 (7) 开始仿真，单击“untitled”模型窗口中“开始仿真”图标，或者选择菜单“Simulink”——“Start”，则仿真开始。双击“Scope”模块出现示波器显示屏，可以看到黄色的正弦波形。如图1.4所示。 (8) 保存模型，单击工具栏的图标，将该模型保存为“Ex0701.mdl”文件。 1.2 Simulink的文件操作和模型窗口 1.2.1 Simulink的文件操作 1. 新建文件 新建仿真模型文件有几种操作： ?在MATLAB的命令窗口选择菜单“File”“New”“Model”。 图7.2 Simulink界面 图7.3 Simulink模型窗口 图7.4 示波器窗口

课程设计之matlab仿真报告

西安邮电大学 专业课程设计报告书 院系名称：电子工程学院学生姓名：李群学号05113096 专业名称：光信息科学与技术班级：光信1103 实习时间：2014年4月8日至2014年4月 18日

一、课程设计题目： 用matlab 仿真光束的传输特性。 二、任务和要求 1、用matlab 仿真光束通过光学元件的变换。 ① 设透镜材料为k9玻璃，对1064nm 波长的折射率为1.5062，镜片中心厚度为3mm ，凸面曲 率半径，设为100mm ，初始光线距离透镜平面20mm 。用matlab 仿真近轴光线（至少10条）经过平凸透镜的焦距，与理论焦距值进行对比，得出误差大小。 ② 已知透镜的结构参数为101=r ，0.11=n ，51=d ，5163.121==' n n （K9玻璃）， 502-=r ，0.12=' n ，物点A 距第一面顶点的距离为100，由A 点计算三条沿光轴夹角分别为10、20、 30的光线的成像。试用Matlab 对以上三条光线光路和近轴光线光路进行仿真，并得出实际光线的球差大小。 ③ 设半径为1mm 的平面波经凸面曲率半径为25mm ，中心厚度3mm 的平凸透镜。用matlab 仿 真平面波在透镜几何焦平面上的聚焦光斑强度分布，计算光斑半径。并与理论光斑半径值进行对比，得出误差大小。（方法：采用波动理论，利用基尔霍夫—菲涅尔衍射积分公式。） 2、用MATLAB 仿真平行光束的衍射强度分布图样。（夫朗和费矩形孔衍射、夫朗和费圆孔衍射、夫朗和费单缝和多缝衍射。） 3、用MATLAB 仿真厄米—高斯光束在真空中的传输过程。（包括三维强度分布和平面的灰度图。） 4、（补充题）查找文献，掌握各类空心光束的表达式，采用费更斯-菲涅尔原理推导各类空心光束在真空中传输的光强表达式。用matlab 对不同传输距离处的光强进行仿真。 三、理论推导部分 第一大题 （1）十条近轴光线透过透镜时，理想情况下光线汇聚透镜的焦点上，焦点到像方主平面的距离为途径的焦距F ，但由于透镜的折射率和厚度会影响光在传输过程中所走的路径（即光程差Δ）。在用MATLAB 仿真以前先计算平行光线的传输路径。，R 为透镜凸面的曲率半径，h 为入射光线的高度，θ1为入射光线与出射面法线的夹角，θ2为出射光线与法线的夹角，n 为透镜材料的折射率。设透镜的中心厚度为d ，则入射光线经过透镜的实际厚度为：L=(R-d) 光线的入射角为：sinq1=h/R 折射角度满足：sinq2=nsinq1 而实际的光束偏折角度为：θ2-θ1。 由此可以看出，当平行光线照射透镜时，在凸面之前光线平行于光轴，在凸面之后发生了偏折，于光轴交汇一点，这一点成为焦点f ，折线的斜率为（-tan(θ2-θ1)）。 （2）根据题意可得，本题所讨论的是与光轴夹角不同的三条光线，经过透镜的两次反射后的成像问题。利用转面公式计算。

BUCK电路闭环控制系统的MATLAB仿真

BUCK 电路闭环PID 控制系统 的MATLAB 仿真 一、课题简介 BUCK 电路是一种降压斩波器，降压变换器输出电压平均值Uo 总是小于输入电压U i 。通常电感中的电流是否连续，取决于开关频率、滤波电感L 和电容C 的数值。 简单的BUCK 电路输出的电压不稳定，会受到负载和外部的干扰，当加入PID 控制器，实现闭环控制。可通过采样环节得到PWM 调制波，再与基准电压进行比较，通过PID 控制器得到反馈信号，与三角波进行比较，得到调制后的开关波形，将其作为开关信号，从而实现BUCK 电路闭环PID 控制系统。 二、BUCK 变换器主电路参数设计 2.1设计及内容及要求 1、 输入直流电压(VIN)：15V 2、 输出电压(VO)：5V 3、 输出电流(IN)：10A 4、 输出电压纹波峰-峰值 Vpp ≤50mV 5、 锯齿波幅值Um=1.5V 6、开关频率(fs)：100kHz 7、采样网络传函H(s)=0.3 8、BUCK 主电路二极管的通态压降VD=0.5V ，电感中的电阻压降 VL=0.1V ，开关管导通压降 VON=0.5V,滤波电容C 与电解电容 RC 的乘积为 F *Ωμ75

2.2主电路设计 根据以上的对课题的分析设计主电路如下： 图2-1 主电路图 1、滤波电容的设计 因为输出纹波电压只与电容的容量以及ESR 有关， rr rr C L N 0.2V V R i I == ? （1） 电解电容生产厂商很少给出ESR ，但C 与R C 的乘积趋于常数，约为50~80μ*ΩF [3]。在本课题中取为75μΩ*F ，由式(1)可得R C =25mΩ，C =3000μF 。 2、滤波电感设计 开关管闭合与导通状态的基尔霍夫电压方程分别如式(2)、(3)所示： IN O L ON L ON /V V V V L i T ---=?（2） O L D L OFF /V V V L i T ++=? （3） off 1/on s T T f += (4) 由上得： L in o L D on V V V V L T i ---=? (5) 假设二极管的通态压降V D =0.5V ，电感中的电阻压降V L =0.1V ，开关管导通压降V ON =0.5V 。利用ON OFF S 1T T f +=，可得T ON =3.73μS ，将此值回代式(5)，可得L =17.5μH

MATLAB仿真实验报告

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/ 」- □歹

buck电路设计

Buck变换器设计——作业 一．Buck主电路设计 1.占空比D计算 2.电感L计算 3.电容C计算 4.开关元件Q的选取 二. Buck变换器开环分析 三. Buck闭环控制设计 1.闭环控制原理 2.补偿环节Gc(s)的设计——K因子法 3.PSIM仿真 4. 补偿环节Gc(s)的修正——应用sisotool 5．修正后的PSIM仿真 四．标称值电路PSIM仿真 五．设计体会 Buck变换器性能指标: 输入电压：标准直流电压48V，变化范围：43V~53V

输出电压：直流电压24V ，5A 输出电压纹波：100mv 电流纹波：0.25A 开关频率：fs=250kHz 相位裕度：60 幅值裕度：10dB 一. Buck 主电路设计： 1.占空比D 计算 根据Buck 变换器输入输出电压之间的关系求出占空比D 的变化范围。 .50V 48V 24U U D .4530V 53V 24U U D 0.558 V 43V 24U U D innom o nom max in o min min in o max ========= 2.电感L 计算 uH 105f i 2)D U -(U i 2)T U -(U L s L min o inmax L on(min) o inmax =?=?= 3.电容C 计算 uF 25.1250000 *1.0*825 .0vf 8i C s L ==??= 电容耐压值：由于最大输出电压为24.1V ，则电容耐压值应大于24.1V 。 考虑到能量储存以及伏在变化的影响，要留有一定的裕度，故电容选取120uf/50V 电容。 4.开关元件Q 的选取

北邮dsp软件matlab仿真实验报告

题目: 数字信号处理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仿真技术报告

专业课程报告 题目：仿真技术与应用课程报告 学院电气工程学院 专业班级 学生姓名 指导教师 提交日期 2013年11月 5日

评语 课程总评成绩： 指导老师： 2013年 12 月20 日

目录 一、负荷预测技术发展情况 (2) 1.1.神经网络理论 (2) 1.2.灰色数学理论 (2) 1.3.组合预测法 (3) 1.4. 线性回归负荷预测方法 (3) 二、算法实现 (4) 三、编程代码 (6) 四、算例测试 (6) 五、心得体会 (9) 六、参考文献 (9)

一、负荷预测技术发展情况 负荷预测技术包括神经网络技术、灰色数学理论、组合预测法、线性回归负荷预测方法 1.1.神经网络理论 神经网络理论是利用神经网络的学习功能，让计算机学习包含在历史负荷数据中的映射关系，再利用这种映射关系预测未来负荷。由于该方法具有很强的鲁棒性、记忆能力、非线性映射能力以及强大的自学习能力，因此有很大的应用市场，但其缺点是学习收敛速度慢，可能收敛到局部最小点；并且知识表达困难，难以充分利用调度人员经验中存在的模糊知识。 神经网络技术进行电力负荷预测，其优点是可以模仿人脑的智能化处理，对大量非结构性、非精确性规律具有自适应功能，具有信息记忆、自主学习、知识推理和优化计算的特点，特别的，其自学习和自适应功能是常规算法和专家系统技术所不具备的。因此，预测被当作人工神经网络（简记为ANN）最有潜力的应用领域之一，许多人都试图应用反传学习算法训练ANN。以用作时间序列预测。误差反向传播算法又称为BP法，提出一个简单的三层人工神经网络模型，就能实现从输入到输出间非线性映射任何复杂函数关系。因此，我们可以将对电力负荷影响最大的几种因素作为输入，即当天的天气温度、天气晴朗度（又称为能见度）、风向风力、峰谷负荷及相关负荷等，争取获得较好的预测结果。 1.2.灰色数学理论 灰色系统理论是中国学者邓聚龙教授1982年3月在国际上首先提出来的，在国际期刊《SYSTEMS AND CONTROL LETTER》刊物上发表，题为“Control Problems of Grey Systems”,引起了国际上的充分重视。 灰色数学理论是把负荷序列看作一真实的系统输出，它是众多影响因子的综合作用结果。这些众多因子的未知性和不确定性，成为系统的灰色特性。灰色系统理论把负荷序列通过生成变换，使其变化为有规律的生成数列再建模，用于负荷预测。 灰色系统理论的形成是有过程的。早年邓教授从事控制理论和模糊系统的研究，取得了许多成果。后来，他接受了全国粮食预测的课题，为了搞好预测工作，他研究了概率统计追求大样本量，必须先知道分布规律、发展趋势，而时间序列法只致力于数据的拟合，不注重规律的发展。邓教授希望在可利用数据不多的情况下，找到了较长时期起

BUCK电路的Saber仿真

功率变换器计算机仿真与设计题目BUCK变换器电路设计 学生姓名 学号 学院 专业电气工程及自动化 班级 指导教师 2013年 10月 20日

一、设计要求 1.1 设计指标： 设计一个BUCK直流变换器，主电路拓扑如图1.1（参数需重新设置），使得其满足以下性能要求： 高压侧蓄电池输入电压V in：30-60V(额定电压48V) 低压侧直流母线输出电压V out：24V 输出电压纹波V out（p-p）：25mV 输出电流I out：2A 开关频率f s：200kHz 电感电流临界连续时I G：0.1A 12 图1.1

二、开环参数计算及仿真 2.1 主电路参数计算： (1)高压侧输入电压V in 变化范围为30-60V ，低压侧输出电压V out 为24V ，则占空比： 8.030 24 min max === in out V V D 4.06024 max min === in out V V D 5.048 24 === innom out nom V V D (2)由于输出电流I out 为2A ，故负载电阻：12out out V R I = =Ω (3)根据电感电流临界连续时I G ：0.1A ，可由下式计算得滤波电感感值： H T I U L U T I L OFF o o CCM o μμ3605)4.01(2 .024 2max min )min(=?-?--==?=?? (4)根据输出电压纹波V out （p-p ）为25mV ，可由下式计算得滤波电容容值： uF f V I C I T C idt C V p p out ripple o p p out T 510 200102582 .082211133)(0) (2=????==?==---? 取F C f μ10=，其中开关频率f 为200KHZ 。 在实际器件中，电容存在寄生电阻，因此实际器件仿真时，电容的选取如下： Ω ====???=??=?-m 125ESR ,600C ,u 520C 25,10652.0min max pp 6 uF F mV V C ESR I V 取而 2.2 开关管及二极管应力计算： （1）开关管的选取 功率管承受的最大电压为60V ，流过开关管电流最大值为2A ，开关管电压电流降额系数均为0.5，则开关管电压要大于或等于120V ，电流最大值要大于4A 。粗略以最大占空比计算电流的有效值为3.2A ，则最大功率为384W ，取400W 。根据仿真,可选irf460作为开关管。 （2）二极管的选取