基于matlab的2PSK的系统仿真

基于matlab的2PSK的系统仿真

摘要：Simulink 是Mathworks 公司推出的基于Matlab 平台的著名仿真环境Simulin 作为一种专业和功能强大且操作简单的仿真工具，目前已被越来越多的工程技术人员所青睐，它搭建积木式的建模仿真方式既简单又直观，而且已经在各个领域得到了广泛的应用。本次课程设计是基于M A T L A B 的2P S K 和2D P S K 仿真，通过系统分析，步骤来完成本次设计任务。通过课程设计从理论学习的轨道逐步引向实际应用，把理论上熟悉的定性分析、定量计算逐步和工程估算、实验调整等手段结合起来，掌握工程设计的步骤和方法，了解科学实验的程序和实施方法，为以后毕业设计和从事信息处理技术的实际工作打下基础。 关键词：MATLAB ；2PSK ，2DPSK ；仿真

1.二进制相移键控

设计流程图如图2-1所示。

c ()

s t c

图 2-1 2PSK 调制解调框图

1.1 PSK 调制原理

在二进制数字调制中，当正弦载波的相位随二进制数字基带信号离散变化时，则产生二进制移相键控(2PSK)信号。2PSK 信号调制有两种方法，即模拟调制法和键控法。通常用已调信号载波的 0°和 180°分别表示二进制数字基带信号的 1 和 0，模拟调制法用两个反相的载波信号进行调制。2PSK 以载波的相位变化作为参考基准的，当基带信号为0时相位相对于初始相位为0°，当基带信号为1时相对于初始相位为180°。

键控法，是用载波的相位来携带二进制信息的调制方式。通常用0°和180°来分别代表0和1。其时域表达式为：

t nT t g a e c n s n PSK ωcos )(2??

?

???-=∑

其中，2PSK 的调制中an 必须为双极性码。本次设计中采用模拟调制法。两种方法原理图分别如图2-2和图2-3所示。

图 2-2 模拟调制法原理图

图 2-3 键控法原理图

2PSK信号的时间波形

1.2 PSK解调原理

由于2PSK的幅度是恒定的，必须进行相干解调。经过带通滤波的信号在相乘器中与本地载波相乘，然后用低通滤波器滤除高频分量，在进行抽样判决。判决器是按极性来判决的。即正抽样值判为1，负抽样值判为0。2PSK信号的相干解调原理图如图2-4所示，各点的波形如图2-5所示。

由于2PSK信号的载波回复过程中存在着180°的相位模糊，即恢复的本地载波与所需相干载波可能相同，也可能相反，这种相位关系的不确定性将会造成解调出的数字基带信号与发送的基带信号正好相反，即“1”变成“0”吗“0”变成“1”，判决器输出数字信号全部出错。这种现象称为2PSK方式的“倒π”现象或“反相工作”。但在本次仿真中是直接给其同频同相的载波信号，所以不存在此问题。

图 2-4 2PSK 的相干解调原理图

图 2-5 相干解调中各点波形图

图中，假设相干载波的基准相位与2PSK 信号的基准一致（通常默认为0相位）。但是由于2PSK 信号的载波回复过程中存在着180°的相位模糊，即恢复的本地载波与所需相干载波可能相同，也可能相反，这种相位关系的不确定性将会造成解调出的数字基带信号与发送的基带信号正好相反，即“1”变成“0”吗“0”变成“1”，判决器输出数字信号全部出错。这种现象称为2PSK 方式的“倒π”现象或“反相工作”。

2PSK 信号在一个码元的持续时间Ts 内可以表示为 u 1T (t) 发送“1”时 S T (t)= u oT (t)=- u 1T (t) 发送“0”时 期中 Acos ωc t 0< t < Ts u 1T (t)= 0 其他

2e

a

b

c

d

e

设发送端发出的信号如上式所示，则接收端带通滤波器输出波形y(t)为

[a+n

c (t)]cosω

c

t-n

s

(t)sinω

c

t 发送“1”时

y(t)=

[-a+n

c (t)]cosω

c

t-n

s

(t)sinω

c

t 发送“0”时

y(t)经过想干解调（相乘—低通）后，送入抽样判决器的输入波形为

a+n

c

(t) 发送“1”时

x(t)=

-a+n

c

(t) 发送“0”时

由于n

c

(t)是均值为0，方差为σ2的高斯噪声，所以x(t)的一维概率密度函数为

由最佳判决门限分析可知，在发送“1”和“0”概率相等时，即P(1)=P(0)时，最佳门限b*=0.此时，发“1”而错判为“0”的概率为

P(0/1)=P(x≦0)=∫0

-∞f

1

(x)dx=1/2erfc(r)

式中：r=a2/2σ2

n

同理，发“0而错判为“1”的概率为

P(1/0)=P(x＞0)=∫0

-∞f

(x)dx=1/2erfc(r)

故2PSK信号相干解调系统的总误码率为

P

e

=P(1)P(0/1)+P(0)P(0/1)= 1/2erfc(r) 在大信噪比(r>>1)的条件下,上式可近似为

P

e ≈e-r/2r

π

2.基于matlab的2PSK的系统仿真

通过编写M文件程序：

2PSK调制解调程序及注释

clear all

close all

i=10;

j=5000;

fc=4; %载波频率

fm=i/5; %码元速率

B=2*fm;

t=linspace(0,5,j);

a=round(rand(1,i)); %随机序列,基带信号figure(3);

stem(a);

st1=t;

for n=1:10

if a(n)<1;

for m=j/i*(n-1)+1:j/i*n

st1(m)=0;

end

else

for m=j/i*(n-1)+1:j/i*n

st1(m)=1;

end

end

end

figure(1);

subplot(411);

plot(t,st1);

title('基带信号st1');

axis([0,5,-1,2]);

%由于PSK中的是双极性信号，因此对上面所求单极性信号取反来与之一起构成双极性码

st2=t;

for k=1:j;

if st1(k)>=1;

st2(k)=0;

else

st2(k)=1;

end

end;

subplot(412);

plot(t,st2);

title('基带信号反码st2');

axis([0,5,-1,2]);

st3=st1-st2;

subplot(413);

plot(t,st3);

title('双极性基带信号st3');

axis([0,5,-2,2]);

s1=sin(2*pi*fc*t);

subplot(414);

plot(s1);

title('载波信号s1');

e_psk=st3.*s1;

figure(2);

subplot(511);

plot(t,e_psk);

title('e_2psk');

noise=rand(1,j);

psk=e_psk+noise; %加入噪声

subplot(512);

plot(t,psk);

title('加噪后波形');

psk=psk.*s1; %与载波相乘

subplot(513);

plot(t,psk);

title('与载波s1相乘后波形');

[f,af] = T2F(t,psk); %通过低通滤波器[t,psk] = lpf(f,af,B);

subplot(514);

plot(t,psk);

title('低通滤波后波形');

for m=0:i-1;

if psk(1,m*500+250)<0;

for j=m*500+1:(m+1)*500;

psk(1,j)=0;

end

else

for j=m*500+1:(m+1)*500;

psk(1,j)=1;

end

end

end

subplot(515);

plot(t,psk);

axis([0,5,-1,2]);

title('抽样判决后波形')

2DPSK调制解调程序及注释

clear all

close all

i=10;

j=5000;

fc=4; %载波频率

fm=i/5; %码元速率

t=linspace(0,5,j);

a=round(rand(1,i));

figure(4);

stem(a);

st1=t;

for n=1:10

if a(n)<1;

for m=j/i*(n-1)+1:j/i*n

st1(m)=0;

end

else

for m=j/i*(n-1)+1:j/i*n

st1(m)=1;

end

end

end

figure(1);

subplot(321);

plot(t,st1);

title('绝对码');

axis([0,5,-1,2]);

b=zeros(1,i); %全零矩阵b(1)=a(1);

for n=2:10

if a(n)>=1;

if b(n-1)>=1

b(n)=0;

else

b(n)=1;

end

else

b(n)=b(n-1);

end

end

for n=1:10

if b(n)<1;

for m=j/i*(n-1)+1:j/i*n st1(m)=0;

end

else

for m=j/i*(n-1)+1:j/i*n st1(m)=1;

end

end

end

subplot(323);

plot(t,st1);

title('相对码st1');

axis([0,5,-1,2]);

st2=t;

for k=1:j;

if st1(k)>=1;

st2(k)=0;

else

st2(k)=1;

end

end;

subplot(324);

plot(t,st2);

title('相对码反码st2');

axis([0,5,-1,2]);

s1=sin(2*pi*fc*t);

subplot(325);

plot(s1);

title('载波信号s1');

s2=sin(2*pi*fc*t+pi);

subplot(326);

plot(s2);

title('载波信号s2');

d1=st1.*s1;

d2=st2.*s2;

figure(2);

subplot(411);

plot(t,d1);

title('st1*s1');

subplot(412);

plot(t,d2);

title('st2*s2');

e_dpsk=d1+d2;

subplot(413);

plot(t,e_dpsk);

title('调制后波形');

noise=rand(1,j);

dpsk=e_dpsk+noise; %加入噪声

subplot(414);

plot(t,dpsk);

title('加噪声后信号');

dpsk=dpsk.*s1; %与载波s1相乘figure(3);

subplot(411);

plot(t,dpsk);

title('与载波相乘后波形');

[f,af]=T2F(t,dpsk); %通过低通滤波器[t,dpsk]=lpf(f,af,B);

subplot(412);

plot(t,dpsk);

title('低通滤波后波形');

st=zeros(1,i); %全零矩阵

for m=0:i-1;

if dpsk(1,m*500+250)<0;

st(m+1)=0;

for j=m*500+1:(m+1)*500;

dpsk(1,j)=0;

end

else

for j=m*500+1:(m+1)*500;

st(m+1)=1;

dpsk(1,j)=1;

end

end

end

subplot(413);

plot(t,dpsk);

axis([0,5,-1,2]);

title('抽样判决后波形')

dt=zeros(1,i); %全零矩阵

dt(1)=st(1);

for n=2:10;

if (st(n)-st(n-1))<=0&&(st(n)-st(n-1))>-1; dt(n)=0;

else

dt(n)=1;

end

end

st=t;

for n=1:10

if dt(n)<1;

for m=j/i*(n-1)+1:j/i*n

st(m)=0;

end

else

for m=j/i*(n-1)+1:j/i*n

st(m)=1;

end

end

end

subplot(414);

plot(t,st);

axis([0,5,-1,2]);

title('码反变换后波形');

%利用FFT计算信号的频谱并与信号的真实频谱的抽样比较。

%脚本文件T2F.m定义了函数T2F，计算信号的傅立叶变换。

function [f,sf]= T2F(t,st)

%This is a function using the FFT function to calculate a signal's Fourier %Translation

%Input is the time and the signal vectors,the length of time must greater %than 2

%Output is the frequency and the signal spectrum

dt = t(2)-t(1);

T=t(end);

df = 1/T;

N = length(st);

f=-N/2*df:df:N/2*df-df;

sf = fft(st);

sf = T/N*fftshift(s）；

用到的低通滤波器函数

function [t,st]=lpf(f,sf,B)

%This function filter an input data using a lowpass filter

%Inputs: f: frequency samples

% sf: input data spectrum samples

% B: lowpass's bandwidth with a rectangle lowpass

%Outputs: t: time samples

% st: output data's time samples

df = f(2)-f(1);

T = 1/df;

hf = zeros(1,length(f)); %全零矩阵

bf = [-floor( B/df ): floor( B/df )] + floor( length(f)/2 );

hf(bf)=1;

yf=hf.*sf;

[t,st]=F2T(f,yf);

st = real(st);

用到的反傅立叶函数

%脚本文件F2T.m定义了函数F2T，计算信号的反傅立叶变换。

function [t,st]=F2T(f,sf)

%This function calculate the time signal using ifft function for the input %signal's spectrum

df = f(2)-f(1);

Fmx = ( f(end)-f(1) +df);

dt = 1/Fmx;

N = length(sf);

T = dt*N;

%t=-T/2:dt:T/2-dt;

t = 0:dt:T-dt;

sff = fftshift(sf);

st = Fmx*ifft(sff);

产生随机信号，按流程图2-1所示顺序对每一模块编程后。程序中注有需注意语句及解释。运行程序，实现2PSK的调制与解调过程。本次设计采用模拟调制法和相干解调法。

参考文献

[1] 李白萍，吴冬梅.通信原理与技术[M].北京：人民邮电出版社,2003

[2] 樊昌信，曹丽娜.通信原理[M].北京：国防工业出版社，2001

[3] 曹志刚，钱亚生.现代通信原理[M].北京：清华大学出版社，1992

[4] 李明明，李白萍.电子信息类专业MATLAB实验教程[M].北京：北京大学出版社，2011

[5] 刘学勇.详解MATLAB/Simulink通信系统建模与仿真[M].北京：电子工业出版社，2011

基于 MATLAB 的QPSK系统仿真设计与实现

通信系统仿真设计实训报告1.课题名称：基于MATLAB 的QPSK系统仿真设计与实现 学生学号： 学生姓名： 所在班级： 任课教师： 2016年10月25日

目录 1.1QPSK系统的应用背景简介 (3) 1.2 QPSK实验仿真的意义 (3) 1.3 实验平台和实验内容 (3) 1.3.1实验平台 (3) 1.3.2实验内容 (3) 二、系统实现框图和分析 (4) 2.1、QPSK调制部分， (4) 2.2、QPSK解调部分 (5) 三、实验结果及分析 (6) 3.1、理想信道下的仿真 (6) 3.2、高斯信道下的仿真 (7) 3.3、先通过瑞利衰落信道再通过高斯信道的仿真 (8) 总结： (10) 参考文献： (11) 附录 (12)

1.1QPSK系统的应用背景简介 QPSK是英文Quadrature Phase Shift Keying的缩略语简称，意为正交相移键控，是一种数字调制方式。在19世纪80年代初期,人们选用恒定包络数字调制。这类数字调制技术的优点是已调信号具有相对窄的功率谱和对放大设备没有线性要求,不足之处是其频谱利用率低于线性调制技术。19世纪80年代中期以后,四相绝对移相键控(QPSK)技术以其抗干扰性能强、误码性能好、频谱利用率高等优点,广泛应用于数字微波通信系统、数字卫星通信系统、宽带接入、移动通信及有线电视系统之中。 1.2 QPSK实验仿真的意义 通过完成设计内容，复习QPSK调制解调的基本原理，同时也要复习通信系统的主要组成部分，了解调制解调方式中最基础的方法。了解QPSK的实现方法及数学原理。并对“通信”这个概念有个整体的理解，学习数字调制中误码率测试的标准及计算方法。同时还要复习随机信号中时域用自相关函数，频域用功率谱密度来描述平稳随机过程的特性等基础知识，来理解高斯信道中噪声的表示方法，以便在编程中使用。 理解QPSK调制解调的基本原理，并使用MATLAB编程实现QPSK信号在高斯信道和瑞利衰落信道下传输，以及该方式的误码率测试。复习MATLAB编程的基础知识和编程的常用算法以及使用MATLAB仿真系统的注意事项，并锻炼自己的编程能力，通过编程完成QPSK调制解调系统的仿真，以及误码率测试，并得出响应波形。在完成要求任务的条件下，尝试优化程序。 通过本次实验，除了和队友培养了默契学到了知识之外，还可以将次实验作为一种推广，让更多的学生来深入一层的了解QPSK以至其他调制方式的原理和实现方法。可以方便学生进行测试和对比。足不出户便可以做实验。 1.3 实验平台和实验内容 1.3.1实验平台 本实验是基于Matlab的软件仿真，只需PC机上安装MATLAB 6.0或者以上版本即可。 （本实验附带基于Matlab Simulink （模块化）仿真，如需使用必须安装simulink 模块） 1.3.2实验内容 1.构建一个理想信道基本QPSK仿真系统,要求仿真结果有 a.基带输入波形及其功率谱 b.QPSK信号及其功率谱

信号与系统的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通信原理仿真 学号： 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仿真 期末大作业 姓名：班级：学号：指导教师：

2012春期末大作业 题目：设单位负反馈控制系统前向通道传递函数由)()(21s G s G 和串联，其中： ) 1(1)()(21++==s A s G s K s G A 表示自己学号最后一位数（可以是零），K 为开环增益。要求： （1）设K=1时，建立控制系统模型，并绘制阶跃响应曲线（用红色虚线，并标注坐标和标题）；求取时域性能指标，包括上升时间、超调量、调节时间、峰值时间； （2）在第（1）问中，如果是在命令窗口绘制阶跃响应曲线，用in1或者from workspace 模块将命令窗口的阶跃响应数据导入Simulink 模型窗口，用示波器显示阶跃响应曲线；如果是在Simulink 模型窗口绘制阶跃响应曲线，用out1或者to workspace 模块将Simulink 模型窗口的阶跃响应数据导入命令窗口并绘制阶跃响应曲线。 （3）用编程法或者rltool 法设计串联超前校正网络，要求系统在单位斜坡输入信号作用时，速度误差系数小于等于0.1rad ，开环系统截止频率s rad c /4.4''≥ω，相角裕度大于等于45度，幅值裕度大于等于10dB 。

仿真结果及分析： （1）、（2）、将Simulink模型窗口的阶跃响应数据导入命令窗口并绘制阶跃响应曲线 通过在Matlab中输入命令： >> plot(tout,yout,'r*-') >> title('阶跃响应曲线') 即可得出系统阶跃响应曲线，如下： 求取该控制系统的常用性能指标：超调量、上升时间、调节时间、峰值时间的程序如下： G=zpk([],[0,-1],5)。 S=feedback(G,1)。

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仿真大作业报告

. 基于MATLAB的电机学计算机辅助分析与仿真 实验报告

一、实验内容及目的 1.1 单相变压器的效率和外特性曲线 1.1.1 实验内容 一台单相变压器，N S =2000kVA, kV kV U U N N 11/127/21=，50Hz ，变压器的参数 和损耗为008.0* ) 75(=C k o R ，0725.0*=k X ，kW P 470=，kW P C KN o 160)75(=。 (1)求此变压器带上额定负载、)(8.0cos 2滞后=?时的额定电压调整率和额定效率。 (2)分别求出当0.1,8.0,6.0,4.0,2.0cos 2=?时变压器的效率曲线，并确定最大效率和达到负载效率时的负载电流。 (3)分析不同性质的负载（),(8.0cos 0.1cos ),(8.0cos 222超前，滞后===???）对变压器输出特性的影响。 1.1.2 实验目的 （1）计算此变压器在已知负载下的额定电压调整率和额定效率 （2）了解变压器效率曲线的变化规律 （3）了解负载功率因数对效率曲线的影响 （4）了解变压器电压变化率的变化规律 （5）了解负载性质对电压变化率特性的影响 1.1.3 实验用到的基本知识和理论 （1）标幺值、效率区间、空载损耗、短路损耗等概念 （2）效率和效率特性的知识 （3）电压调整率的相关知识 1.2串励直流电动机的运行特性 1.2.1实验内容 一台16kw 、220V 的串励直流电动机，串励绕组电阻为0.12Ω，电枢总电阻为0.2Ω。电动势常数为.电机的磁化曲线近似的为直线。其中为比例常数。假设电枢电流85A 时，磁路饱和(为比较不同饱和电流对应的效果，饱和电流可以自己改变）。

基于MATLAB的QPSK通信系统仿真设计毕业设计论文

基于MATLAB的QPSK通信系统仿真设计 摘要 随着移动通信技术的发展，以前在数字通信系统中采用FSK、ASK、PSK 等调制方式，逐渐被许多优秀的调制技术所替代。本文主要介绍了QPSK调制与解调的实现原理框图，用MATLAB软件中的SIMULINK仿真功能对QPSK调制与解调这一过程如何建立仿真模型，通过对仿真模型的运行，得到信号在QPSK 调制与解调过程中的信号时域变化图。通过该软件实现方式，可以大大提高设计的灵活性，节约设计时间，提高设计效率，从而缩小硬件电路设计的工作量，缩短开发周期。 关键词 QPSK，数字通信，调制，解调，SIMULINK -I-

Abstract As mobile communications technology, and previously in the adoption of digital cellular system, ASK, FSK PSK modulation, etc. Gradually been many excellent mod ulation technology substitution, where four phase-shift keying QPSK technology is a wireless communications technology in a binary modulation method. This article prim arily describes QPSK modulation and demodulation of the implementation of the prin ciple of block diagrams, focuses on the MATLAB SIMULINK software emulation in on QPSK modulation and demodulation the process how to build a simulation model, through the operation of simulation model, I get signal in QPSK modulation and dem odulation adjustment process domain change figure. The software implementation, ca n dramatically improve the design flexibility, saving design time, increase efficiency, design to reduce the workload of hardware circuit design, and shorten the developmen t cycle. Keywords QPSK, Digital Communication,modulation,demodulation,SIMULINK -II-

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 大作业

上海电力学院 通信原理Matlab仿真 实验报告 实验名称： 8QAM误码率仿真 试验日期： 2014年 6月3日 专业：通信工程 姓名：罗侃鸣 班级： 2011112班 学号： 20112272

一、实验要求 写MATLAB程序，对图示的信号星座图完成M=8的QAM通信系统Monte Carlo仿真，在不同SNRindB=0：15时，对N=10000（3比特）个符号进行仿真。画出该QAM系统的符号误码率。 二、实验原理 1 QAM调制原理 QAM(Quadrature Amplitude Modulation)正交幅度调制技术，是用两路独立的基带信号对两个相互正交的同频载波进行抑制载波双边带调幅，利用这种已调信号的频谱在同一带宽内的正交性，实现两路并行的数字信息的传输。该调制方式通常有8QAM，16QAM，64QAM。 QAM调制实际上就是幅度调制和相位调制的组合，相位+ 幅度状态定义了一个数字或数字的组合。QAM的优点是具有更大的符号率，从而可获得更高的系统效率。通常由符号率确定占用带宽。因此每个符号的比特（基本信息单位）越多，频带效率就越高。 调制时，将输入信息分成两部分：一部分进行幅度调制；另一部分进行相位调制。对于星型8QAM信号，每个码元由3个比特组成，可将它分成第一个比特和后两个个比特两部分。前者用于改变信号矢量的振幅，后者用于差分相位调制，通过格雷编码来改变当前码元信号矢量相位与前一码元信号矢量相位之间的相位差。 QAM是一种高效的线性调制方式，常用的是8QAM，16QAM，64QAM等。当随着M 的增大，相应的误码率增高，抗干扰性能下降。 2 QAM星座图 QAM调制技术对应的空间信号矢量端点分布图称为星座图。QAM的星座图呈现星状分

MATLAB Simulink系统建模与仿真 实验报告

MATLAB/Simulink 电力系统建模与仿真 实验报告 姓名：****** 专业：电气工程及其自动化 班级：******************* 学号：*******************

实验一无穷大功率电源供电系统三相短路仿真 1.1 无穷大功率电源供电系统仿真模型构建 运行MATLAB软件，点击Simulink模型构建，根据电路原理图，添加下列模块： （1）无穷大功率电源模块（Three-phase source） （2）三相并联RLC负荷模块（Three-Phase Parallel RLC Load） （3）三相串联RLC支路模块（Three-Phase Series RLC Branch） （4）三相双绕组变压器模块（Three-Phase Transformer (Two Windings)） （5）三相电压电流测量模块（Three-Phase V-I Measurement） （6）三相故障设置模块（Three-Phase Fault） （7）示波器模块（Scope） （8）电力系统图形用户界面（Powergui） 按电路原理图连接线路得到仿真图如下： 1.2 无穷大功率电源供电系统仿真参数设置 1.2.1 电源模块 设置三相电压110kV，相角0°，频率50Hz，接线方式为中性点接地的Y形接法，电源电阻0.00529Ω，电源电感0.000140H，参数设置如下图：

1.2.2 变压器模块 变压器模块参数采用标幺值设置，功率20MVA，频率50Hz，一次测采用Y型连接，一次测电压110kV，二次侧采用Y型连接，二次侧电压11kV，经过标幺值折算后的绕组电阻为0.0033，绕组漏感为0.052，励磁电阻为909.09，励磁电感为106.3，参数设置如下图： 1.2.3 输电线路模块 根据给定参数计算输电线路参数为：电阻8.5Ω，电感0.064L，参数设置如下图： 1.2.4 三相电压电流测量模块 此模块将在变压器低压侧测量得到的电压、电流信号转变成Simulink信号，相当于电压、电流互感器的作用，勾选“使用标签（Use a label）”以便于示波器观察波形，设置电压标签“Vabc”，电流标签“Iabc”，参数设置如下图：

matlab机电系统仿真大作业

一曲柄滑块机构运动学仿真 1、设计任务描述 通过分析求解曲柄滑块机构动力学方程，编写matlab程序并建立Simulink 模型，由已知的连杆长度和曲柄输入角速度或角加速度求解滑块位移与时间的关系，滑块速度和时间的关系，连杆转角和时间的关系以及滑块位移和滑块速度与加速度之间的关系，从而实现运动学仿真目的。 2、系统结构简图与矢量模型 下图所示是只有一个自由度的曲柄滑块机构，连杆与长度已知。 图2-1 曲柄滑块机构简图 设每一连杆（包括固定杆件）均由一位移矢量表示，下图给出了该机构各个杆件之间的矢量关系 图2-2 曲柄滑块机构的矢量环

3．匀角速度输入时系统仿真 3.1 系统动力学方程 系统为匀角速度输入的时候，其输入为输出为；。 (1) 曲柄滑块机构闭环位移矢量方程为: (2)曲柄滑块机构的位置方程 (3)曲柄滑块机构的运动学方程 通过对位置方程进行求导，可得 由于系统的输出是与，为了便于建立A*x=B形式的矩阵，使x=[], 将运动学方程两边进行整理，得到 将上述方程的v1与w3提取出来，即可建立运动学方程的矩阵形式 3.2 M函数编写与Simulink仿真模型建立 3.2.1 滑块速度与时间的变化情况以及滑块位移与时间的变化情况 仿真的基本思路：已知输入w2与，由运动学方程求出w3和v1，再通过积分，即可求出与r1。 （1）编写Matlab函数求解运动学方程 将该机构的运动学方程的矩阵形式用M函数compv（u）来表示。 设r2=15mm，r3=55mm，r1（0）=70mm，。 其中各个零时刻的初始值可以在Simulink模型的积分器初始值里设置

M函数如下： function[x]=compv(u) %u(1)=w2 %u(2)=sita2 %u(3)=sita3 r2=15; r3=55; a=[r3*sin(u(3)) 1;-r3*cos(u(3)) 0]; b=[-r2*u(1)*sin(u(2));r2*u(1)*cos(u(2))]; x=inv(a)*b; （2）建立Simulink模型 M函数创建完毕后，根据之前的运动学方程建立Simulink模型，如下图： 图3-1 Simulink模型 同时不要忘记设置r1初始值70，如下图： 图3-2 r1初始值设置

QPSK通信系统性能分析与MATLAB仿真报告

淮海工学院课程设计报告书 课程名称：通信系统的计算机仿真设计 题目：QPSK通信系统性能分析 与MATLAB仿真 学院：电子工程学院 学期：2013-2014-2 专业班级： 姓名： 学号： 评语： 成绩： 签名： 日期：

QPSK通信系统性能分析与MATLAB仿真 1 绪论 1.1 研究背景与研究意义 数字信号传输系统分为基带传输系统和频带传输系统,频带传输系统也叫数字调制系统,该系统对基带信号进行调制,使其频谱搬移到适合在信道(一般为带通信道)上传输的频带上。数字调制和模拟调制一样都是正弦波调制,即被调制信号都为高频正弦波。数字调制信号又称为键控信号,数字调制过程中处理的是数字信号,而载波有振幅、频率和相位3个变量,且二进制的信号只有高低电平两个逻辑量即1和0,所以调制的过程可用键控的方法由基带信号对载频信号的振幅、频率及相位进行调制,最基本的方法有3种:正交幅度调制(QAM) 、频移键控( FSK) 、相移键控( PSK) 。根据所处理的基带信号的进制不同分为二进制和多进制调制(M进制) 。 本实验采用QPSK。QPSK是英文Quadrature Phase Shift Keying的缩略语简称，意为正交相移键控，是一种数字调制方式。在19世纪80年代初期人们选用恒定包络数字调制。这类数字调制技术的优点是已调信号具有相对窄的功率谱和对放大设备没有线性要求不足之处是其频谱利用率低于线性调制技术。19世纪80年代中期以后四相绝对移相键控(QPSK)技术以其抗干扰性能强、误码性能好、频谱利用率高等优点广泛应用于数字微波通信系统、数字卫星通信系统、宽带接入、移动通信及有线电视系统之中。 1.2 课程设计的目的和任务 目的在于使学生在课程设计过程中能够理论联系实际，在实践中充分利用所学理论知识分析和研究设计过程中出现的各类技术问题，巩固和扩大所学知识面，为以后走向工作岗位进行设计打下一定的基础。 课程设计的任务是： (1)掌握一般通信系统设计的过程，步骤，要求，工作内容及设计方法，掌握用计算机仿真通信系统的方法。 (2)训练学生网络设计能力。 (3)训练学生综合运用专业知识的能力，提高学生进行通信工程设计的能力。1.3 可行性分析 QPSK是英文Quadrature Phase Shift Keying的缩略语简称，意为正交相移键控，是一种数字调制方式。在19世纪80年代初期,人们选用恒定包络数字调制。这类数字调制技术的优点是已调信号具有相对窄的功率谱和对放大设备没有线性要求,不足之处是其频谱利用率低于线性调制技术。19世纪80年代中期以后,

实验一 典型环节的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仿真-电机特性曲线)

电机大作业 专业班级：电气XXXX 姓名：XXX 学号：XXX 指导老师：张威

一、研究课题（来源：教材习题 4-18 ） 1. 74 、R 2 0.416 、X 2 3.03 、R m 6. 2 X m 75 。电动机的机械损耗p 139W,额定负载时杂散损耗p 320W, 试求额定负载时的转差率、定子电流、定子功率因数、电磁转矩、输出转矩和效 率。 二、编程仿真 根据T 形等效电路: 3D - R Q 运用MATLAB 进行绘图。MATLAB 文本中，P N PN ，U N UN ，尺 R 1， X 1 X1 , R 2 R 2，X 2 X 2，R m Rm, X m Xm ，p pjixiesunh ao ， p pzasansunhao 。定子电流I11，定子功率因数 Cosangle1，电磁转矩Te ， 效率 Xiaolv 。 1.工作特性曲线绘制 MATLA 文本： R1=0.715;X 仁1.74;Rm=6.2;Xm=75;R2=0.416;X2=3.03;pjixiesu nhao=139; pzasa nsu nhao=320;p=2;m 仁 3; ns=1500;PN=17000;UN=380;fN=50; Z1=R1+j*X1; Zm=Rm+j*Xm; for i=1:2500 s=i/2500; nO=n s*(1-s); Z2=R2/s+j*X2; Z=Z1+Zm*Z2/(Zm+Z2); 有一台三相四极的笼形感应电动机， 参数为P N 17kW 、U N 380V (△联 Rm 结）、尺 0. 715 、X j lcr S

U1=UN; I1=U1/Z; l110=abs(l1); An gle 仁an gle(ll); Cosa ngle10=cos(A ngle1); P仁3*U1*l110*Cosa ngle10; l2=l1*Zm/(Zm+Z2); Pjixie=m1*(abs(I2))A2*(1-s)/s*R2; V=(1-s)*pi*fN; Te0=Pjixie/V; P20=Pjixie-pjixies un hao-pzasa nsun hao; Xiaolv0=P20/P1; P2(i)=P20; n (i)=n0; l11(i)=l110; Cosa ngle1(i)=Cosa ngle10; Te(i)=Te0; Xiaolv(i)=Xiaolv0; hold on; end figure(1) plot(P2, n); xlabel('P2[W]');ylabel(' n[rpm]'); figure(2) plot(P2,l11); xlabel('P2[W]');ylabel('l1[A]'); figure(3) plot(P2,Cosa nglel); xlabel('P2[W]');ylabel('go nglvyi nshu'); figure(4) plot(P2,Te); xlabel('P2[W]');ylabel('Te[Nm]'); figure(5) plot(P2,Xiaolv); xlabel('P2[W]');ylabel('xiaolv');

基于MATLAB的QPSK系统仿真设计与实现

通信系统仿真设计实训报告1.课题名称：基于 MATLAB 的QPSK系统仿真设计与实现 学生学号： 学生： 所在班级： 任课教师： 2016年 10月25日

目录 1.1QPSK系统的应用背景简介 (3) 1.2 QPSK实验仿真的意义 (3) 1.3 实验平台和实验容 (3) 1.3.1实验平台 (3) 1.3.2实验容 (3) 二、系统实现框图和分析 (4) 2.1、QPSK调制部分， (4) 2.2、QPSK解调部分 (5) 三、实验结果及分析 (6) 3.1、理想信道下的仿真 (6) 3.2、高斯信道下的仿真 (7) 3.3、先通过瑞利衰落信道再通过高斯信道的仿真 (8) 总结： (10) 参考文献： (11) 附录 (12)

1.1QPSK系统的应用背景简介 QPSK是英文Quadrature Phase Shift Keying的缩略语简称，意为正交相移键控，是一种数字调制方式。在19世纪80年代初期,人们选用恒定包络数字调制。这类数字调制技术的优点是已调信号具有相对窄的功率谱和对放大设备没有线性要求,不足之处是其频谱利用率低于线性调制技术。19世纪80年代中期以后,四相绝对移相键控(QPSK)技术以其抗干扰性能强、误码性能好、频谱利用率高等优点,广泛应用于数字微波通信系统、数字卫星通信系统、宽带接入、移动通信及有线电视系统之中。 1.2 QPSK实验仿真的意义 通过完成设计容，复习QPSK调制解调的基本原理，同时也要复习通信系统的主要组成部分，了解调制解调方式中最基础的方法。了解QPSK的实现方法及数学原理。并对“通信”这个概念有个整体的理解，学习数字调制中误码率测试的标准及计算方法。同时还要复习随机信号中时域用自相关函数，频域用功率谱密度来描述平稳随机过程的特性等基础知识，来理解高斯信道中噪声的表示方法，以便在编程中使用。 理解QPSK调制解调的基本原理，并使用MATLAB编程实现QPSK信号在高斯信道和瑞利衰落信道下传输，以及该方式的误码率测试。复习MATLAB编程的基础知识和编程的常用算法以及使用MATLAB仿真系统的注意事项，并锻炼自己的编程能力，通过编程完成QPSK调制解调系统的仿真，以及误码率测试，并得出响应波形。在完成要求任务的条件下，尝试优化程序。 通过本次实验，除了和队友培养了默契学到了知识之外，还可以将次实验作为一种推广，让更多的学生来深入一层的了解QPSK以至其他调制方式的原理和实现方法。可以方便学生进行测试和对比。足不出户便可以做实验。 1.3 实验平台和实验容 1.3.1实验平台 本实验是基于Matlab的软件仿真，只需PC机上安装MATLAB 6.0或者以上版本即可。 （本实验附带基于Matlab Simulink （模块化）仿真，如需使用必须安装simulink 模块） 1.3.2实验容 1.构建一个理想信道基本QPSK仿真系统,要求仿真结果有 a.基带输入波形及其功率谱 b.QPSK信号及其功率谱

MATLAB与系统仿真实验指导书

实验指导书（标准格式） 《MATLAB与系统仿真》 实验指导书 课程编号：ME2121025课程名称：MATLAB与系统仿真适应专业：电气工程及其自动化课程类别：限选 实验教学种类：上机课程总学时：30＋6 实验学时：6 执笔人：勾燕洁 西安电子科技大学机电工程学院 2006 年9 月

目录 实验一MA TLAB的基本使用及矩阵操作------------------- 3 实验二MA TLAB编程与图形处理---------------------------- 9 实验三MA TLAB中的数值运算与系统仿真---------------- 15

实验一MA TLAB的基本使用及矩阵操作 一、实验目的和要求 熟悉MA TLAB的界面和基本操作，掌握矩阵的建立方法及各种运算。 1．熟悉MA TLAB软件的界面和帮助系统。 2．掌握MA TLAB软件中关于矩阵建立、矩阵初等变换以及矩阵算术、关系、逻辑运算的各种命令。 3．掌握MA TLAB软件中M函数和M文件的编写，以及程序结构与控制，学会编写一般程序。 二、实验内容 1．启动与退出 2．数、数组、矩阵的输入 3．MA TLAB的基本命令 4．矩阵大小的测试 5．矩阵元素的操作 6．特殊矩阵的产生 7．矩阵的算术运算 8．矩阵的关系运算 9．矩阵的逻辑运算 10．常用函数 三、实验仪器、设备（软、硬件）及仪器使用说明 PC机一台，注意正确开、关机及打开软件。 四、实验原理 无 五、实验方法与步骤 1．启动与退出 双击MA TLAB图标，进入MA TLAB命令窗口，即可输入命令，开始运算； 观察各个窗口以及菜单； 单击File菜单中的Exit，或使用MATLAB的Exit命令退出。

运动控制MATLAB仿真

大作业: 直流双闭环调速MATLAB仿真 运动控制技术课程名称： 名：姓电气学院院：学 自动化业：专 号：学 孟濬指导教师： 2012年6月2日

------------------------------------- -------------学浙大江 李超 一、Matlab仿真截图及模块功能描述 Matlab仿真截图如下，使用Matlab自带的直流电机模型： 模块功能描述： ⑴电机模块（Discrete DC_Machine）：模拟直流电机 ⑵负载转矩给定（Load Torque）：为直流电机添加负载转矩 ⑶Demux：将向量信号分离出输出信号 ⑷转速给定（Speed Reference）：给定转速 ⑸转速PI调节（Speed Controller）：转速PI调节器，对输入给定信号与实际信号

的差值进行比例和积分运算，得到的输出值作为电流给定信号。改变比例和积分运算系数可以得到不同的PI控制效果。 ⑹电流采样环节（1/z）：对电流进行采样，并保持一个采样周期 ⑺电流滞环调节（Current Controller）：规定一个滞环宽度，将电流采样值与给定值进行对比，若：采样值＞给定值+0.5*滞环宽度，则输出0； 若：采样值＜给定值—0.5*滞环宽度，则输出1； 若：给定值—0.5*滞环宽度＜采样值＜给定值+0.5*滞环宽度，则输出不变 输出值作为移相电压输入晶闸管斩波器控制晶闸管触发角 ：根据输入电压改变晶闸管触发角，从而改变电机端电压。GTO⑻晶闸管斩波．⑼续流二极管D1：在晶闸管关断时为电机续流。 ⑽电压传感器Vd：测量电机端电压 ⑾示波器scope：观察电压、电流、转速波形 系统功能概括如下：直流电源通过带GTO的斩波器对直流电机进行供电，输出量电枢电流ia和转速wm通过电流环和转速环对GTO的通断进行控制，从而达到对整个电机较为精确的控制。 下面对各个部分的功能加以详细说明： （1）直流电机 双击电动机模块，察看其参数：

一阶、二阶系统时域和频域仿真

西安交通大学 基于MATLAB/Simulink 的一阶、二阶系统的时域和频 域仿真 ——以单位阶跃信号为输入信号 日期：2013年4月 一阶系统时域和频域仿真 1、建立一阶系统典型数学模型 ()1 1 G s Ts =+ 2、建立simulink 仿真方框图

1T.s+1 Transfer Fcn Step Scope ① 时间常数T=1时，一阶系统时域响应为 12345678 910 00.5 1 一阶系统时域相应（T=1） Matlab 程序： %一阶系统仿真编程 num=[1]; den=[1 1]; bode(num,den); grid on ; gtext('低频段频率-20dB/dec'); 运行程序，有时间常数T=1时，一阶系统的频域响应为

10 -210 -1 10 10 1 10 2 -90-45 一阶系统频域响应 P h a s e (d e g ) Bode Di a gram Frequency (rad/s) -40-30-20-100 低频段斜率-20dB/dec System: sys Frequency (rad/s): 1.01Magni t ude (dB): -3.07 M a g n i t u d e (d B ) ② 时间常数T=3时，一阶系统单位阶跃时域响应 12345678910 00.5 1 一阶系统单位阶跃响应（T=3） Matlab 程序： %一阶系统仿真编程 num=[1]; den=[3 1]; bode(num,den);

grid on ; gtext('低频段频率-20dB/dec'); 运行程序，有时间常数T=3时，一阶系统的频域响应为 10 -210 -1 10 10 1 -90-45 P h a s e (d e g ) Bode Di a gram Frequency (rad/s) -30-20-100 低频段频率-20dB/dec System: sys Frequency (rad/s): 0.334Magni t ude (dB): -3.03 M a g n i t u d e (d B ) 3、分析以上一阶系统在不同时间常数下的单位阶跃响应，可以看出时间常数越小，系统响应越快；而且一阶系统的转角频率为1/T ，在转角频率以上时，幅频特性曲线以-20dB/dec 下降，而相频特性以0°和90°为渐近线。

基于MATLAB的系统仿真实验

基于MATLAB的系统仿真实验 实验指导书 新乡学院 二○○八年三月

说明 MATLAB是MATrix LABoratory的缩写，是一种基于矩阵的数学与工程计算系统，可以用作动态系统的建模与仿真。研究系统的结构和参数的变化对系统性能的影响可利用MATLAB强大的计算和作图功能，因此本实验采用MATLAB仿真研究连续控制系统和离散控制系统的性能分析过程。通过该实验提高学生对控制系统的分析与设计能力，加深对《自动控制原理》课程内容的理解。 由于学时有限，该实验由同学们参考有关MATLAB应用的书籍利用课后时间完成。

第一部分 基于MATLAB 连续系统的仿真 一、 实验目的 在研究系统的结构和参数的变化对系统性能的影响时，采用解析和作图的方法比较麻烦，而且误差也大，用MATLAB 仿真实现则简单方便，精度高。本实验采用MATLAB 实现控制系统的数学描述、控制系统的时域分析及根轨迹和频率特性分析。通过该实验，加深学生对系统阶次，型号，参数与系统性能的关系的理解。 二、实验环境 在计算机Windows 环境下安装好MATLAB6.3以上版本后，双击MATLAB 图标或成“开始”菜单打开MATLAB ，即可进入MATLAB 集成环境。 三、MATLAB 应用实例 1．拉氏变换和反变换 例 求22)(2++=t t t f 的拉氏变换 解 键入 syms s t; ft=t^2+2*t+2; st=laplace(ft,t,s) 运行结果为 st= 2/s^3+2/s^2+2/s 例 求) 2)(34(6 )(2 ++++= s s s s s F 的拉氏反变换 解 键入 syms s t; Fs=(s+6)/(s^2+4*s+3)/(s+2); ft=ilaplace(Fs,s,t) 运行结果为 ft= 3/2*exp(-3*t)+5/2*exp(-t)-4*exp(-2*t) 2。求根运算 例 求多项式 43)(2 3 ++=s s s p 的根，再由根建多项式。 解 键入 p=[1 3 0 4]; r=root(p) 运行结果为 r= -3.3553 0.1777+1.0773i 1.7777-1.0773i