直流调速系统的MATLAB仿真(参考程序)

直流调速系统的MATLAB 仿真

一、开环直流速系统的仿真

开环直流调速系统的电气原理如图1所示。直流电动机的电枢由三相晶闸管整流电路经平波电抗器L 供电，通过改变触发器移相控制信号c U 调节晶闸管的控制角α，从而改变整流器的输出电压，实现直流电动机的调速。该系统的仿真模型如图2所示。

图1 开环直流调速系统电气原理图

图2 直流开环调速系统的仿真模型

为了减小整流器谐波对同步信号的影响，宜设三相交流电源电感s 0L =，直流电动机励磁由直流电源直接供电。触发器（6-Pulse ）的控制角（alpha_deg ）由移相控制信号c U 决定，移相特性的数学表达式为

min

c cmax

9090U U αα?-=?-

在本模型中取min 30α=?，cmax 10V U =，所以c 906U α=-。在直流电动机的负载转矩输入端L T 用Step 模块设定加载时刻和加载转矩。

仿真算例1 已知一台四极直流电动机额定参数为N 220V U =，N 136A I =，

N 1460r /min n =，a 0.2R =Ω，2222.5N m GD =?。励磁电压f 220V U =，励磁电

流f 1.5A I =。采用三相桥式整流电路，设整流器内阻rec 0.3R =Ω。平波电抗器

d 20mH L =。仿真该晶闸管-直流电动机开环调速系统，观察电动机在全压起动

和起动后加额定负载时的电机转速n 、电磁转矩e T 、电枢电流d i 及电枢电压d u 的变化情况。N 220V U =

仿真步骤：

1）绘制系统的仿真模型（图2）。 2）设置模块参数（表1） ① 供电电源电压

N rec N 2min 2200.3136

130(V)2.34cos 2.34cos30U R I U α++?=

=≈??

② 电动机参数 励磁电阻：

f f f 220

146.7()1.5

U R I =

==Ω 励磁电感在恒定磁场控制时可取“0”。 电枢电阻：

a 0.2R =Ω

电枢电感由下式估算：

N a N N 0.422019.1

19.10.0021(H)2221460136

CU L pn I ?==?≈???

电枢绕组和励磁绕组间的互感af L ：

N a N e N 2200.2136

0.132(V min/r)1460

U R I K n --?=

=≈?

T e 60600.132 1.262π2πK K =

=?≈ T af f 1.26

0.84(H)1.5

K L I =

== 电机转动惯量

2222.5

0.57(kg m )449.81

GD J g ==≈??

③ 额定负载转矩

L T N 1.26136171.4(N m)T K I ==?≈?

表1 开环直流调速系统主要模型参数

3）设置仿真参数：仿真算法odel5s ，仿真时间5.0s ，直流电动机空载起动，起动2.5s 后加额定负载L 171.4N m T =?。

4）进行仿真并观察、分析结果（图3）。

（可以用语句plot(tout,yout)进行示波器的曲线处理。）

图3 开环直流调速系统的仿真结果

二、转速闭环直流调速系统的仿真

带转速负反馈的有静差直流调速系统的电气原理如图4所示，系统由转速给

定环节*

n U 、转速调节器ASR （放大器p K ）、移相触发器GT 、晶闸管整流器UCR

和直流电动机M 和测速发电机TG 等组成。

图4 转速闭环直流调速系统电气原理图

图5 转速闭环直流调速系统的仿真模型

转速负反馈有静差直流调速系统的仿真模型如图5所示，模型在图2开环调

速系统的基础上，增加了转速给定*

n U ，转速反馈n-feed 、放大器Gain 和反映放

大器输出限幅的饱和特性模块Saturation ，饱和限幅模块的输出是移相触发器的控制电压c U ，转速反馈直接取自电动机的转速输出，没有另加测速发电机，取

转速反馈系数*

nm

n N

U K n =。

仿真算例2 在算例1的基础上观察转速负反馈系统在不同放大器放大倍数时对转速变化的影响。

仿真步骤：

1）绘制系统的仿真模型（图5）。

2）设置模块参数（表2）。

3）设置仿真参数：仿真算法odel5s，仿真时间1.5s，直流电动机空载起动，

起动0.5s后加额定负载

L 171.4N m

T=?。

4）进行仿真并观察、分析结果（图6）：

（用语句plot(tout1,yout1,tout2,yout2,tout3,yout3)进行示波器的曲线处理。）表2 转速闭环直流调速系统主要模型参数

图6 转速闭环直流调速系统的仿真结果

三、转速电流双闭环直流调速系统的仿真

转速电流双闭环直流调速系统的电气原理如图7所示，由于晶闸管整流器不能通过反向电流，因此不能产生反向制动转矩而使电动机快速制动。

图7 转速电流双闭环直流调速系统的电气原理图

双闭环直流调速系统的仿真可以依据系统的动态结构图（图8a）进行，也可以用SIMULINK的Power System模块来组建。两种仿真的不同在于主电路，前者晶闸管和电动机用传递函数来表示，后者晶闸管和电动机使用Power System 模块，而控制部分则是相同的。下面对这两种方法分别进行介绍。

1. 基于动态结构图的双闭环直流调速系统仿真

双闭环直流调速系统的实际动态结构图如图8b所示，它与图8a的不同之处

在于增加了滤波环节，包括电流滤波、转速滤波和两个给定信号的滤波环节。这是因为电流检测信号中常含有交流分量，为了不使它影响到调节器的输入，需加低通滤波。这样的滤波环节的传递函数可用一阶惯性环节来表示，其滤波时间常数oi T 可按需要选定，以滤平电流检测信号为准。然而，在抑制交流分量的同时，滤波环节也延迟了反馈信号的作用，为了平衡这个延迟作用，在给定信号通道上加入一个同等时间常数的惯性环节，称作给定滤波环节。其意义是，让给定信号和反馈信号经过相同的延时，使二者在时间上得到恰当的配合，从而带来设计上的方便。同样，由测速发电机得到的转速反馈电压信号含有换向纹波，因此也需要滤波，滤波时间常数用on T 表示。根据和电流环一样的道理，在转速给定通道上也加入时间常数为on T 的给定滤波环节。

a)

b)

图8 转速电流双闭环直流调速系统的动态结构图

依据系统动态结构图的仿真模型如图9所示，仿真模型与系统动态结构图的各个环节基本上是对应的。需要指出的是，双闭环系统的转速和电流两个调节器都是有饱和特性和带输出限幅的PI 调节器，为了充分反映在饱和和限幅非线性影响下调速系统的工作情况，需要构建考虑饱和和输出限幅的PI 调节器，过程如下：

线性PI 调节器的传递函数为

i pi p p

1()k s

W s k k s s

+=+

=ττ 式中，p k 为比例系数，i k 为积分系数，时间常数p i /k k =τ。

上述PI 调节器的传递函数可以直接调用SIMULINK 中的传递函数或零极点模块，而考虑饱和和输出限幅的PI 调节器模型如图10所示。模型中比例和积分环节分为两个通道，其中积分模块Integrate 的限幅表示调节器的饱和限幅值，而调节器的输出限幅值由饱和模块Saturation 设定。

图9 转速电流双闭环直流调速系统仿真模型

图10 带饱和和输出限幅的PI 调节器

仿真算例3 以算例1的晶闸管-直流电动机系统为基础，设计一个转速电流双闭环控制的调速系统，设计指标为：转速超调量n %10%σ≤，电流超调量

i %10%σ≤，过载倍数 1.5λ=，取电流反馈滤波时间常数oi 0.002s T =，转速反馈滤波时间常数on 0.01s T =，取转速调节器和电流调节器的饱和值为12V ，输出限

幅值为10V 。额定转速时转速给定电压*

nm

10V U =。通过仿真观察系统的转速、电流响应，以及参数变化（主要是调节器参数）对系统响应的影响。

仿真步骤：

1）构建系统的仿真模型（图9）。

2）设置模块参数（调节器参数计算和设定） ① 机电时间常数：m 0.161s T =

电磁时间常数：0.076s l T =

三相晶闸管整流电路平均失控时间：s 0.0017s T =

② 电流调节器ACR 参数的计算*

电流反馈系数：*im i N 10

0.05(V/A)1.5136

U K I λ==≈?

电流环时间常数之和i s oi 0.00170.0020.0037(s)T T T ∑=+=+= ACR 的传递函数为i ACR pi ii

pi i 11

(s)s W K K K s s

ττ+=+=，其中 时间常数i 0.076s l T ==τ 比例系数i pi Σi S i

0.0760.5

2.57220.0037400.05

R

k T K K ?=

=

≈???τ

积分系数pi

ii i

2.57

33.80.076

k k =

=

≈τ ③ 转速调节器ASR 参数的计算*

转速反馈系数：*nm n N 10

0.00667(V min/r)1500

U K n ==≈?

电流环等效时间常数i 220.00370.0074(s)T ∑=?=

转速环时间常数之和n i on 220.00370.010.0174(s)T T T ∑∑=+=?+=。 ASR 的传递函数为n ASR pn pi

pn n 11

(s)s W k k k s s

+τ=+=τ，其中 时间常数n n 50.01740.087(s)hT ∑τ==?= 比例系数i e m pn n n (1)60.050.1320.16110.992250.006670.50.0174

h K K T K hK RT ∑+???=

=≈????

积分系数pn in n

10.99

126.30.087

k k =

=

≈τ （选择中频段宽度5h =）

模型各环节参数如图9所示，其中调节器参数见表3，调节器积分环节限幅

值为12V

±，调节器输出限幅值为10V

±。

表3 转速电流双闭环直流调速系统主要模型参数

3）设置仿真参数：仿真算法odel5s，仿真时间2.0s，电动机空载起动，起

动0.8s后突加额定负载（

dL N 136A

I I

==）。4）进行仿真并观察、分析结果（图11）：

a）

b ）

图11 基于动态结构图的双闭环直流调速系统仿真结果

2. 基于Power System 模块的双闭环直流调速系统仿真

采用SIMULINK 的Power System 模块组成的转速电流双闭环直流调速系统的仿真模型如图12所示，模型由晶闸管-直流电动机组成的主电路和转速、电流调节器组成的控制电路两部分构成。其中的主电路部分，交流电源、晶闸管整流器、触发器、移相控制和电动机等环节使用Power System 模型库中的模块。控制电路的主体是转速和电流两个调节器，以及反馈滤波环节，这部分与前述基于动态结构图的双闭环系统仿真相同。将这两部分拼接起来即组成晶闸管-电动机转速电流双闭环控制的直流调速系统的仿真模型。

模型中转速反馈和电流反馈均直接取自电机测量单元的转速和电流输出端，这样减少了测速和电流检测环节，但并不影响仿真的真实性。电流调节器ACR 的输出端接移相特性模块（Shifter ）的输入端，而ACR 的输出限幅值就决定了控制角的min α（o 30）和max α（o 150）。

图12 基于Power System的双闭环直流调速系统仿真模型应该注意，图12与图9仿真模型的不同在于以晶闸管整流器和电动机模型取代了动态结构图中的晶闸管整流器和电动机传递函数，由于动态结构图中的晶闸管整流器和电动机传动函数是线性的，其电流可以反向（如图11b），因此转速调节过程要快一些，而实际的晶闸管整流器不能通过反向电流，所以仿真的结果略有不同，采用晶闸管整流器和电动机模型的仿真可以更好地反映系统的工作情况（如图13所示）。

改变调节器ASR、ACR参数的仿真结果如图14所示，可见转速、电流的仿真结果基本满足设计指标的要求。

（可用语句plot(tout,yout(:,1))和plot(tout,yout(:,2))进行示波器的曲线处理。）

a）

b）

图13 基于Power System的双闭环直流调速系统仿真结果

a）

b）

图14 改变调节器参数的仿真结果

（

pn 20

k=，

in 60

k=；

pi 5

k=，

ii 15

k=）

直流调速系统的MATLAB仿真(参考程序)汇总.

直流调速系统的MATLAB 仿真 一、开环直流速系统的仿真 开环直流调速系统的电气原理如图1所示。直流电动机的电枢由三相晶闸管整流电路经平波电抗器L 供电，通过改变触发器移相控制信号c U 调节晶闸管的控制角α，从而改变整流器的输出电压，实现直流电动机的调速。该系统的仿真模型如图2所示。 图1 开环直流调速系统电气原理图 图2 直流开环调速系统的仿真模型 为了减小整流器谐波对同步信号的影响，宜设三相交流电源电感s 0L =，直流电动机励磁由直流电源直接供电。触发器（6-Pulse ）的控制角（alpha_deg ）由移相控制信号c U 决定，移相特性的数学表达式为 min c cmax 9090U U αα?-=?-

在本模型中取min 30α=?，cmax 10V U =，所以c 906U α=-。在直流电动机的负载转矩输入端L T 用Step 模块设定加载时刻和加载转矩。 仿真算例1 已知一台四极直流电动机额定参数为N 220V U =，N 136A I =， N 1460r /min n =，a 0.2R =Ω，2222.5N m GD =?。励磁电压f 220V U =，励磁电流f 1.5A I =。采用三相桥式整流电路，设整流器内阻rec 0.3R =Ω。平波电抗器 d 20mH L =。仿真该晶闸管-直流电动机开环调速系统，观察电动机在全压起动和起动后加额定负载时的电机转速n 、电磁转矩 e T 、电枢电流d i 及电枢电压d u 的变化情况。N 220V U = 仿真步骤： 1）绘制系统的仿真模型（图2）。 2）设置模块参数（表1） ① 供电电源电压 N rec N 2min 2200.3136 130(V)2.34cos 2.34cos30U R I U α++?= =≈?? ② 电动机参数 励磁电阻： f f f 220146.7()1.5 U R I = ==Ω 励磁电感在恒定磁场控制时可取“0”。 电枢电阻： a 0.2R =Ω 电枢电感由下式估算： N a N N 0.422019.1 19.10.0021(H)2221460136 CU L pn I ?==?≈??? 电枢绕组和励磁绕组间的互感af L ： N a N e N 2200.2136 0.132(V min/r)1460 U R I K n --?= =≈?

基于MATLAB的电力系统潮流计算

基于MATLAB的电力系统潮流计算 %简单潮流计算的小程序，相关的原始数据数据数据输入格式如下： %B1是支路参数矩阵，第一列和第二列是节点编号。节点编号由小到大编写%对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点%编号，将变压器的串联阻抗置于低压侧处理。 %第三列为支路的串列阻抗参数。 %第四列为支路的对地导纳参数。 %第五烈为含变压器支路的变压器的变比 %第六列为变压器是否是否含有变压器的参数，其中“1”为含有变压器，%“0”为不含有变压器。 %B2为节点参数矩阵，其中第一列为节点注入发电功率参数；第二列为节点%负荷功率参数；第三列为节点电压参数；第六列为节点类型参数，其中 %“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数。 %X为节点号和对地参数矩阵。其中第一列为节点编号，第二列为节点对地%参数。 n=input('请输入节点数:n='); n1=input('请输入支路数:n1='); isb=input('请输入平衡节点号:isb='); pr=input('请输入误差精度:pr='); B1=input('请输入支路参数:B1='); B2=input('请输入节点参数:B2='); X=input('节点号和对地参数:X='); Y=zeros(n); Times=1; %置迭代次数为初始值 %创建节点导纳矩阵 for i=1:n1 if B1(i,6)==0 %不含变压器的支路 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/B1(i,3); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4); Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4); else %含有变压器的支路 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/(B1(i,3)*B1(i,5)); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/B1(i,3);

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

信号与系统的MATLAB仿真

信号与系统的MATLAB 仿真 一、信号生成与运算的实现 1.1 实现)3(sin )()(π±== =t t t t S t f a )(sin )sin()sin(sin )()(t c t t t t t t t S t f a '=' '== ==πππ π ππ m11.m t=-3*pi:0.01*pi:3*pi; % 定义时间范围向量t f=sinc(t/pi); % 计算Sa(t)函数 plot(t,f); % 绘制Sa(t)的波形 运行结果： 1.2 实现)10() sin()(sin )(±== =t t t t c t f ππ m12.m t=-10:0.01:10; % 定义时间范围向量t f=sinc(t); % 计算sinc(t)函数 plot(t,f); % 绘制sinc(t)的波形 运行结果： 1.3 信号相加：t t t f ππ20cos 18cos )(+= m13.m syms t; % 定义符号变量t f=cos(18*pi*t)+cos(20*pi*t); % 计算符号函数f(t)=cos(18*pi*t)+cos(20*pi*t) ezplot(f,[0 pi]); % 绘制f(t)的波形 运行结果：

1.4 信号的调制：t t t f ππ50cos )4sin 22()(+= m14.m syms t; % 定义符号变量t f=(2+2*sin(4*pi*t))*cos(50*pi*t) % 计算符号函数f(t)=(2+2*sin(4*pi*t))*cos(50*pi*t) ezplot(f,[0 pi]); % 绘制f(t)的波形 运行结果： 1.5 信号相乘：)20cos()(sin )(t t c t f π?= m15.m t=-5:0.01:5; % 定义时间范围向量 f=sinc(t).*cos(20*pi*t); % 计算函数f(t)=sinc(t)*cos(20*pi*t) plot(t,f); % 绘制f(t)的波形 title('sinc(t)*cos(20*pi*t)'); % 加注波形标题 运行结果：

基于Matlab的电力系统故障研究仿真

基于Matlab的电力系统故障分析与仿真 摘要：本文介绍了MATLAB软件在电力系统中的应用，以及利用动态仿真工具Simulink和电力系统工具箱PSD进行仿真的基本方法。在仿真平台上，以单机—无穷大系统为建模对象，通过选择模块，参数设置，以及连线，对电力系统的多种故障进行仿真分析。同时，设计一个GUI图形界面，将仿真波形清晰地显示在界面上以便比较和分析。结果表明，仿真波形基本符合理论分析，说明了MATLAB是电力系统仿真研究的有力工具。 关键词：电力系统；仿真；故障；MATLAB；GUI Abstract：This paper introduces the applications of MATLAB in power system analysis, and the basic simulation method of taking use of Simulink and PSD. On MATLAB simulation platform, take a single machine-infinite-bus system as modeling objects, by selecting the module, parameter settings, and connectingmodules to simulate and analysevariousfault of power system. At the same time, in order to facilitate comparison and analysis simulation waveform, design a GUI for showing waveform clearly.The results show that the simulation waveform in line with theoretical analysis, indicates that MATLAB is a powerful tool for researching simulation of power system. Keywords：PowerSystem。 Simulation。 Fault。 Matlab。 GUI 0 前言[1,2] 随着电力工业的发展，电力系统规划、运行和控制的复杂性亦日益增加，电力系统的生产和研究中仿真软件的应用也越来越广泛。现在，我们主要使用的电力系统仿真软件有：EMTP程序，用于电力系统电磁暂态计算，电力系统暂态过电压分析，暂态保护装置的综合选择等。PSCAD／EMTDC程序，典型应用是计算电力系统遭受扰动或参数变化时，参数随时间变化的规律。PSASP，其功能主要有稳态分析、故障分析和机电暂态分析。还有MathWorks公司开发的MATLAB软件。在MATLAB中，电力系统模型可以在Simulink环境下直接搭建，也可以进行封装和自定义模块库，充分显现了其仿真平台的优越性。更重要的是，MATLAB提供了丰富的工具箱资源，以及大量的实用模块，使我们可以更加深入地研究电力系统的行为特性。本篇论文将在熟练掌握MATLAB软件的基础上，对电力系统的故障进行建模、仿真、分析，并且设计一个GUI图形用户界面来反映故障波形。

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.

matlab控制系统仿真.

课程设计报告 题目PID控制器应用 课程名称控制系统仿真院部名称龙蟠学院 专业自动化 班级M10自动化 学生姓名 学号 课程设计地点 C208 课程设计学时一周 指导教师应明峰 金陵科技学院教务处制成绩

一、课程设计应达到的目的 应用所学的自动控制基本知识与工程设计方法，结合生产实际，确定系统的性能指标与实现方案，进行控制系统的初步设计。 应用计算机仿真技术，通过在MATLAB软件上建立控制系统的数学模型，对控制系统进行性能仿真研究，掌握系统参数对系统性能的影响。 二、课程设计题目及要求 1．单回路控制系统的设计及仿真。 2．串级控制系统的设计及仿真。 3．反馈前馈控制系统的设计及仿真。 4．采用Smith 补偿器克服纯滞后的控制系统的设计及仿真。 三、课程设计的内容与步骤 (1)．单回路控制系统的设计及仿真。 （a）已知被控对象传函W(s) = 1 / (s2 +20s + 1)。 （b）画出单回路控制系统的方框图。 （c）用MatLab的Simulink画出该系统。

（d）选PID调节器的参数使系统的控制性能较好，并画出相应的单位阶约响应曲线。注明所用PID调节器公式。PID调节器公式Wc（s）=50(5s+1)/(3s+1) 给定值为单位阶跃响应幅值为3。 有积分作用单回路控制系统PID控制器取参数分别为：50 2 5 有积分作用单回路控制系统PID控制器取参数分别为：50 0 5

大比例作用单回路控制系统PID控制器取参数分别为：50 0 0 （e）修改调节器的参数，观察系统的稳定性或单位阶约响应曲线，理解控制器参数对系统的稳定性及控制性能的影响？ 答：由上图分别可以看出无积分作用和大比例积分作用下的系数响应曲线，这两个PID调节的响应曲线均不如前面的理想。增大比例系数将加快系统的响应，但是过大的比例系数会使系统有比较大的超调，并产生振荡，使稳定性变坏；

直流调速系统的MATLAB仿真参考程序汇总

直流调速系统的MA TLAB仿真 一、开环直流速系统的仿真 开环直流调速系统的电气原理如图1所示。直流电动机的电枢由三相晶闸管整流电路经平波电抗器供电，通过改变触发器移相控制信号调节晶闸管的U L c ，从 而改变整流器的输出电压，实现直流电动机的调速。该系统的仿真控制角模型如图2所示。 L + + GT UCR E d- - 开环直流调速系统电气原理图图1 图2 直流开环调速系统的仿真模型

L?0，直为了减小整流器谐波对同步信号的影响，宜设三相交流电源电感s流电动机励磁由直流电源直接供电。触发器（6-Pulse）的控制角（alpha_deg）由U 决定，移相特性的数学表达式为移相控制信号c???90?min U?90??c U cmax 1 ??。在直流电动机的负载，所以，在本模型中取U?10V6??30?90U?ccmaxmin转矩输 入端用Step模块设定加载时刻和加载转矩。T L仿真算例1已知一台四极直流电动机额定参数为，，136AIU?220V?NN22。励磁电压，励磁电，， 220VminUR?0.2??1460rn?/m?22.5NGD?fNa流。采用三相桥式整流电路，设整流器内阻。平波电抗器??1.5A0.3RI?recf。仿真该晶闸管-直流电动机开环调速系统，观察电动机在全压起动20mHL?d n、电磁转矩、电枢电流和起动后加额定负载时的电机转速及电枢电压的uTi ded变化情况。220V?U N仿真步骤： 1）绘制系统的仿真模型（图2）。 2）设置模块参数（表1） ①供电电源电压 U?RI220?0.3?136NNrec130(V)U??? 2?2.34?cos302.34cos?min②电动机参数 励磁电阻： U220f)146.7(?R???f I1.5f励磁电感在恒定磁场控制时可取“0”。 电枢电阻： ?0.2R?a电枢电感由下式估算： CU0.4?220N?19.1?L?19.1?0.0021(H) a2pnI2?2?1460?136NN L：电枢绕组和励磁绕组间的互感 af U?RI220?0.2?136NNa?K?0.132(V?min/r)?e n1460N 2 6060K??0.132?K?1.26eTπ2π2K1.26T0.84(H)??L?af1.5I f电机转动惯量 222.5GD2 )??0.57(kg?mJ?9.814?4g 额定负载转矩③ 模块参数名参数

基于matlab的电力系统短路电流计算

湖北民族学院 信息工程学院 题目: 基于matlab的电力系统短路电流计算 专业：电气工程及其自动化 班级： 0308407 学号： 030840705 学生姓名： 指导教师： 2011年6 月1 日

信息工程学院课程设计任务书 年月日

信息工程学院课程设计成绩评定表

摘要 随着电力工业的发展，电力系统的规模越来越大，在这种情况下，许多大型的电力科研实验很难进行，尤其是电力系统中对设备和人员等危害最大的事故故障，尤其是短路故障，而在分析解决事故故障时要不断的实验，在现实设备中很难实现，一是实际的条件难以满足；二是从系统的安全角度来讲也是不允许进行实验的。考虑这两种情况，寻求一种最接近于电力系统实际运行状况的数字仿真工具十分重要，而MATLAB软件中的SIMULINK是用来对动态系统进行建模、仿真和分析的集成开发环境，是结合了框图界面和交互仿真能力的非线性动态系统仿真工具，为解决具体的工程问题提供了更为快速、准确和简洁的途径。 关键词：短路电流计算，MATLAB，仿真 Abstract Along with the development of the electric power industry, the scale of the power system is more and more big, in this case, many large power research is difficult to, especially in the power system, equipment and personnel to the harm such as the biggest accident, especially fault fault location, and on the analysis of the accident to solve the fault of the experiment, in the reality constantly in equipment, it is difficult to accomplish a is practical conditions to meet; The security of the system from the perspective is not allowed in the experiment. Consider the two kinds of circumstances, for one of the most close to power system actual the operation condition of digital simulation tool is very important, and MATLAB software SIMULINK is used for the dynamic system modeling, simulation and analysis of the integrated development environment, is combined with the block diagram interface and interactive simulation of nonlinear dynamic system ability of simulation tools, to solve the specific engineering problem, provides a more rapid, accurate and simple way. Keywords: short-circuit current calculation, MATLAB, the simulatio

电力电子技术与电力系统分析matlab仿真

电气2013级卓班电力电子技术与电力系统分析 课程实训报告 专业：电气工程及其自动化 班级： 姓名： 学号： 指导教师：

兰州交通大学自动化与电气工程学院 2016 年 1 月日

电力电子技术与电力系统分析课程实训报告 1 电力电子技术实训报告 1.1 实训题目 1.1.1电力电子技术实训题目一 一．单相半波整流 参考电力电子技术指导书中实验三负载，建立MATLAB/Simulink环境下三相半波整流电路和三相半波有源逆变电路的仿真模型。仿真参数设置如下： (1)交流电压源的参数设置和以前实验相关的参数一样。 (2)晶闸管的参数设置如下： R=0.001Ω，L =0H，V f=0.8V，R s=500Ω，C s=250e-9F on (3)负载的参数设置 RLC串联环节中的R对应R d，L对应L d，其负载根据类型不同做不同的调整。 (4)完成以下任务： ①仿真绘出电阻性负载（RLC串联负载环节中的R d= Ω，电感L d=0，C=inf，反电动势为0）下α=30°，60°，90°，120°，150°时整流电压U d，负载电流L 和晶闸管两端电压U vt1的波形。 d ②仿真绘出阻感性负载下（负载R d=Ω，电感L d为，反电动势E=0）α=30°，60°，90°，120°，150°时整流电压U d，负载电流L d和晶闸管两端电压U vt1的波形。 ③仿真绘出阻感性反电动势负载下α=90°，120°，150°时整流电压U d，负载电流L d和晶闸管两端电压U vt1的波形，注意反电动势E的极性。 (5)结合仿真结果回答以下问题： ①该三项半波可控整流电路在β=60°，90°时输出的电压有何差异?

单闭环直流调速系统的MATLAB计算与仿真

1 绪论 直流电动机具有良好的起、制动性能，宜于在广泛范围内平滑调速，在轧钢机、矿井卷扬机、挖掘机、海洋钻机、金属切割机床、造纸机、高层电梯等需要高性能可控电力拖动的领域中得到了广泛的应用。近年来直流调速系统发展很快，然而直流拖动控制系统毕竟在理论上和实践上都比较成熟，而且从反馈闭环控制的角度来看，它又是交流拖动控制系统的基础，所以首先应该很好的掌握直流系统。我们可以首先从单闭环转速负反馈直流调速系统来研究。由于系统需要观察较多的性能，计算参数较多，而MATLAB中的Simulink实用工具可直接构建其动态模型，省去大量的计算，通过修改动态模型可完善系统性能。 直流调速系统概述 从生产机械要求控制的物理量来看，电力传动自动控制系统有调速系统、位置伺服系统、张力控制系统等其他多种类型，各种系统往往是通过控制转速来实现的，因此调速系统是最基本的驱动控制系统。调速系统目前分为交流和直流调速控制系统，由于直流调速系统的调速范围广，静差率小、稳定性好并且具有良好的动态性能。因此在相当长的时期内，高性能的调速系统几乎都采用了直流调速系统。相比于交流调速系统，直流调速系统在理论上和实践上更加成熟。 直流调速是现代电力拖动自动控制系统中发展较早的自动控制系统。在20世纪60年代发展起来的电力电子技术，使电能可以转换和控制，产生了现代各种高效、节能的新型电源和交直流调速装置，为工业生产，交通运输，建筑、办公、家庭自动化控制设备提供了现代化的高新技术，提高了生产效率和人们的生活质量，因此，人类社会的生产、生活发生了巨大变化。随着新型电力电子器件的研究和开发，先进控制技术的发展，电力电子和电力传动控制装置的性能也不断优化和提高，这一变化的影响将越来越大。 单闭环直流电机调速系统在现代日常生活中的应用越来越广泛，其良好的调速性能、低廉的价格越来越被大众接受。 单闭环直流调速系统由整流变压器、平波电抗器、晶闸管整流调速装置、电动机-发电机、闭环控制系统组成。我们可以通过调节晶闸管的控制角来调节转速，非常方便，高效。

基于MATLAB的电力系统仿真

《电力系统设计》报告题目: 基于MATLAB的电力系统仿 学院：电子信息与电气工程学院 班级: 13级电气 1 班 姓名:田震 学号: 日期：2015年12月6日 基于MATLAB的电力系统仿真 摘要：目前，随着科学技术的发展和电能需求量的日益增长，电力系统规模越来 越庞大，超高压远距离输电、大容量发电机组、各种新型控制装置得到了广泛的应用，这对于合理利用能源，充分挖掘现有的输电潜力和保护环境都有重要意义。另一方面，随着国民经济的高速发展，以城市为中心的区域性用电增长越来越快，大电网负荷中心的用电容量越来越大，长距离重负荷输电的情况日益普遍，电力系统在人们的生活和工作中担任重要角色，电力系统的稳定运行直接影响着人们的日常生活。从技术和安全上考虑直接进行电力试验可能性很小，因此迫切要求运用电力仿真来解决这些问题。 电力系统仿真是将电力系统的模型化、数学化来模拟实际的电力系统的运行，可以帮助人们通过计算机手段分析实际电力系统的各种运行情况，从而有效的了解电力系统概况。本文根据电力系统的特点，利用MATLAB的动态仿真软件Simulink搭建了无穷大电源的系统仿真模型，得到了在该系统主供电线路电源端发生三相短路接地故障并由故障器自动跳闸隔离故障的仿真结果，并分析了这一暂态过程。通过仿真结果说明MATLAB 电力系统工具箱是分析电力系统的有效工具。 关键词：电力系统；三相短路；故障分析；MATLAB仿真 目录 一．前言.............................................. 二．无穷大功率电源供电系统仿真模型构建............... 1.总电路图的设计......................................

基于MATLAB的电力系统仿真

《电力系统设计》报告 题目: 基于MATLAB的电力系统仿学院：电子信息与电气工程学院 班级: 13级电气 1 班 姓名:田震 学号: 20131090124 日期：2015年12月6日

基于MATLAB的电力系统仿真 摘要：目前，随着科学技术的发展和电能需求量的日益增长，电力系统规模越来越庞大，超高压远距离输电、大容量发电机组、各种新型控制装置得到了广泛的应用，这对于合理利用能源，充分挖掘现有的输电潜力和保护环境都有重要意义。另一方面，随着国民经济的高速发展，以城市为中心的区域性用电增长越来越快，大电网负荷中心的用电容量越来越大，长距离重负荷输电的情况日益普遍，电力系统在人们的生活和工作中担任重要角色，电力系统的稳定运行直接影响着人们的日常生活。从技术和安全上考虑直接进行电力试验可能性很小，因此迫切要求运用电力仿真来解决这些问题。 电力系统仿真是将电力系统的模型化、数学化来模拟实际的电力系统的运行，可以帮助人们通过计算机手段分析实际电力系统的各种运行情况，从而有效的了解电力系统概况。本文根据电力系统的特点，利用MATLAB的动态仿真软件Simulink搭建了无穷大电源的系统仿真模型，得到了在该系统主供电线路电源端发生三相短路接地故障并由故障器自动跳闸隔离故障的仿真结果，并分析了这一暂态过程。通过仿真结果说明MATLAB电力系统工具箱是分析电力系统的有效工具。 关键词：电力系统；三相短路；故障分析；MATLAB仿真

目录 一．前言 (4) 二．无穷大功率电源供电系统仿真模型构建 (5) 1.总电路图的设计 (5) 2.各个元件的参数设定 (6) 2.1供电模块的参数设定 (6) 2.2变压器模块的参数设置 (6) 2.3输电线路模块的参数设置 (7) 2.4三相电压电流测量模块 (8) 2.5三相线路故障模块参数设置 (8) 2.6三相并联RLC负荷模块参数设置 (9) 3.仿真结果 (9)

用Matlab计算潮流计算电力系统分析

《电力系统潮流上机》课程设计报告 院系：电气工程学院 班级：电088班 学号： 0812002221 学生姓名：刘东昇 指导教师：张新松 设计周数：两周 日期：2010年 12 月 25 日

一、课程设计的目的与要求 目的：培养学生的电力系统潮流计算机编程能力，掌握计算机潮流计算的相关知识 要求：基本要求： 1.编写潮流计算程序； 2.在计算机上调试通过； 3.运行程序并计算出正确结果； 4.写出课程设计报告 二、设计步骤： 1.根据给定的参数或工程具体要求（如图），收集和查阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 三、设计原理 1．牛顿-拉夫逊原理 牛顿迭代法是取x0 之后，在这个基础上，找到比x0 更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0 的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为

额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。 牛顿—拉夫逊迭代法的一般步骤： （1）形成各节点导纳矩阵Y。 （2）设个节点电压的初始值U和相角初始值e 还有迭代次数初值为0。 （3）计算各个节点的功率不平衡量。 （4）根据收敛条件判断是否满足，若不满足则向下进行。 （5）计算雅可比矩阵中的各元素。 （6）修正方程式个节点电压 （7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。 （8）计算平衡节点输出功率和各线路功率 2．网络节点的优化 １）静态地按最少出线支路数编号 这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。当由n 个节点的出线支路相同时，则可以按任意次序对这n 个节点进行编号。这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也２）动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出来的，在编号过程中认为固定不变的，事实上，在节点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化(增加，减少或保持不变)。因此，如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果，动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。 3．MATLAB编程应用 Matlab 是“Matrix Laboratory”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。由于使用Matlab 编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。 四、设计内容

增量调制MATLAB仿真实验

增量调制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个采样点使客商的编码输出序列以 解码 编码 二电平量化 单位延迟 单位 延迟

09级系统仿真与MATLAB语言实验

系统仿真与MATLAB语言 实验指导书

对参加实验学生的总要求 1、认真复习有关理论知识，明确每次实验目的，了解实验相关软件操作，熟悉实验内容和方法。 2、实验过程中注意仔细观察，认真记录有关数据和图像，并经由指导教师查验后方可结束实验。 3、应严格遵守实验室规章制度，服从实验室教师的安排和管理。 4、对实验仪器的操作使用严格按照实验室要求进行。

实验总要求 1、封面：注明实验名称、实验人员班级、学号（全号）和姓名等。 2、内容方面：注明实验所用设备、仪器及实验步骤方法；记录清楚实验所得的原始数据和图像，并按实验要求绘制相关图表、曲线或计算相关数据；认真分析所得实验结果，得出明确实验结论。并注明该结论所依据的原理和理论；对实验进行反馈回顾，总结出实验方法要领和注意事项，对实验失败的原因进行分析剖解，总结出实验的经验和教训。 3、文字方面，撰写规范，杜绝错别字。 4、杜绝抄袭，杜绝提供不真实的实验内容。

实验一 MATLAB 语言工作环境和基本操作 1 实验目的 1）．熟悉MATLAB 的开发环境； 2）．掌握MATLAB 的一些常用命令； 3）．掌握矩阵、变量、表达式的输入方法及各种基本运算。 2 实验器材 计算机WinXP 、Matlab7.0软件 3 实验内容 (1). 输入 A=[7 1 5;2 5 6;3 1 5]，B=[1 1 1; 2 2 2;3 3 3]， 在命令窗口中执行下列表达式，掌握其含义： A(2, 3) A(:,2) A(3,:) A(:,1:2:3) A(:,3).*B(:,2) A(:,3)*B(2,:) A*B A.*B A^2 A.^2 B/A B./A (2).输入 C=1:2:20，则 C （i ）表示什么？其中 i=1,2,3,…,10； （3）掌握MA TLAB 常用命令 >> who %列出工作空间中变量 >> whos %列出工作空间中变量，同时包括变量详细信息 >>save test %将工作空间中变量存储到test.mat 文件中 >>load test %从test.mat 文件中读取变量到工作空间中 >>clear %清除工作空间中变量 >>help 函数名 %对所选函数的功能、调用格式及相关函数给出说明 >>lookfor %查找具有某种功能的函数但却不知道该函数的准确名称 如： lookfor Lyapunov 可列出与Lyapunov 有关的所有函数。 （4） 在MATLAB 的命令窗口计算: 1) )2sin(π 2) 5.4)4.05589(÷?+ (5). 试用 help 命令理解下面程序各指令的含义： clear t =0:0.001:2*pi; subplot(2,2,1); polar(t, 1+cos(t)) subplot(2,2,2); plot(cos(t).^3,sin(t).^3) subplot(2,2,3); polar(t,abs(sin(t).*cos(t))) subplot(2,2,4); polar(t,(cos(2*t)).^0.5) （6）（选做）设计M 文件计算: x=0:0.1:10 当sum>1000时停止运算，并显示求和结果及计算次数。 i i i x x sum 2100 2 -= ∑ =

基于Matlab计算程序的电力系统运行分析课程设计

课程设计 课程名称：电力系统分析 设计题目：基于Matlab计算程序地电力系统运行分析学院：电力工程学院 专业：电气工程自动化 年级： 学生姓名： 指导教师： 日期： 教务处制

目录 前言 (1) 第一章参数计算 (2) 一、目标电网接线图 (2) 二、电网模型地建立 (3) 第二章潮流计算 (6) 一．系统参数地设置 (6) 二．程序地调试 (7) 三、对运行结果地分析 (13) 第三章短路故障地分析计算 (15) 一、三相短路 (15) 二、不对称短路 (16) 三、由上面表对运行结果地分析及在短路中地一些问题 (21) 心得体会 (26) 参考文献 (27)

前言 电力系统潮流计算是电力系统分析中地一种最基本地计算，是对复杂电力系统正常和故障条件下稳态运行状态地计算.潮流计算地目标是求取电力系统在给定运行状态地计算.即节点电压和功率分布，用以检查系统各元件是否过负荷.各点电压是否满足要求，功率地分布和分配是否合理以及功率损耗等.对现有电力系统地运行和扩建，对新地电力系统进行规划设计以及对电力系统进行静态和暂态稳定分析都是以潮流计算为基础.潮流计算结果可用如电力系统稳态研究，安全估计或最优潮流等对潮流计算地模型和方法有直接影响. 在电力系统中可能发生地各种故障中，危害最大且发生概率较高地首推短路故障.产生短路故障地主要原因是电力设备绝缘损坏.短路故障分为三相短路、两相短路、单相接地短路及两相接地短路.其中三相短路时三相电流仍然对称，其余三类短路统成为不对称短路.短路故障大多数发生在架空输电线路.电力系统设计与运行时，要采取适当地措施降低短路故障地发生概率.短路计算可以为设备地选择提供原始数据.

高频电子线路Matlab仿真实验

高频电子线路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);