离散时间系统分析及MATLAB实现

离散时间系统分析及MATLAB实现
离散时间系统分析及MATLAB实现

基于matlab实现OFDM的编码.

clc; clear all; close all; fprintf('OFDM系统仿真\n'); carrier_count=input('输入系统仿真的子载波数: \n');%子载波数128,64,32,16 symbols_per_carrier=30;%每子载波含符号数 bits_per_symbol=4;%每符号含比特数,16QAM调制 IFFT_bin_length=1024;%FFT点数 PrefixRatio=1/4;%保护间隔与OFDM数据的比例1/6~1/4 GI=PrefixRatio*IFFT_bin_length ;%每一个OFDM符号添加的循环前缀长度为1/4*IFFT_bin_length ,即256 beta=1/32;%窗函数滚降系数 GIP=beta*(IFFT_bin_length+GI);%循环后缀的长度40 SNR=10; %信噪比dB %================信号产生=================================== baseband_out_length=carrier_count*symbols_per_carrier*bits_per_symbol;%所输入的比特数目 carriers=(1:carrier_count)+(floor(IFFT_bin_length/4)-floor(carrier_count/2));%共轭对称子载波映射复数数据对应的IFFT点坐标 conjugate_carriers = IFFT_bin_length - carriers + 2;%共轭对称子载波映射共轭复数对应的IFFT点坐标 rand( 'twister',0); %每次产生不相同得伪随机序列 baseband_out=round(rand(1,baseband_out_length));%产生待调制的二进制比特流figure(1); stem(baseband_out(1:50)); title('二进制比特流') axis([0, 50, 0, 1]); %==============16QAM调制==================================== complex_carrier_matrix=qam16(baseband_out);%列向量 complex_carrier_matrix=reshape(complex_carrier_matrix',carrier_count,symbols_per

离散时间系统特性分析

实验五实验报告 实验名称:离散时间系统特性分析

一、实验目的: 1 。深入理解单位样值响应,离散系统的频率响应的概念; 2。 掌握通过计算机进行求得离散系统的单位样值响应,以及离散系统的频率 响应的方法。 二、实验原理: 对于离散系统的单位样值而言,在实际处理过程中,不可能选取无穷多项的取值。往往是选取有限项的取值,当然这里会产生一个截尾误差,但只要这个误差在相对小一个范围里,可以忽略不计。 另外,在一些实际的离散系统中,往往不是事先就能得到描述系统的差分方程的,而是通过得到系统的某些相应值,则此时系统的分析就需借助计算机的数值处理来进行,得到描述系统的某些特征,甚至进而得到描述系统的数学模型。 本实验首先给出描述系统的差分方程,通过迭代的方法求得系统的单位样值响应,进而求得该离散系统的频率响应。限于试验条件,虽然给出了系统方程,但处理的方法依然具有同样的实际意义。 具体的方法是: 1 在给定系统方程的条件下,选取激励信号为δ(n),系统的起始状态为零 状态,通过迭代法,求得系统的单位样值响应h(n)(n=0,…,N )。 2 利用公式 其中Ω的取值范围为0~2π 。计算系统的频率响应。 三、实验内容 1 已知系统的差分方程为 利用迭代法求得系统的单位样值响应,取N =10。 2 利用公式 其中

#include #include #define N 10 #define M 20 #define pi 3.1415926 struct pinlv{ double fu; double xiang;}; double h[N+1],x[N+1]; struct pinlv PL(double w) { double a=0, b=0,fu,xiang; int k; struct pinlv FX; for(k=0;k<=N;k++){ a=a+h[k]*cos(-k*w); b=b+h[k]*sin(-k*w);} fu=sqrt(a*a+b*b); xiang=atan(b/a); if((a<0)&(b>0)) xiang=xiang+pi; if((a<0)&(b<0)) xiang=xiang-pi; FX.fu=fu; FX.xiang=xiang; return(FX); } main() { int i,j; double w0; struct pinlv FX[M+1]; FILE *fp1,*fp2; fp1=fopen("H:\\单位样值响应.txt","w"); fp2=fopen("H:\\频率特性.txt","w"); h[-1]=0;h[-2]=0; for(i=-1;i<=N;i++) x[i]=0; x[0]=1; for(i=0;i<=N;i++) h[i]=1.3*h[i-1]-0.4*h[i-2]+x[i-1]; printf("系统的单位样值响应为\n"); fprintf(fp1,"系统的单位样值响应(从x[0]开始)为\n"); fprintf(fp1,"激励x[i] 响应y[i]\n"); for(i=0;i<=N;i++)

环境系统分析小结

环境系统分析小结 环工卓越班章雷1302031005 摘要:系统分析是对研究对象进行有目的、有步骤的探索和研究过程,它运用科学的方法和工具,确定一个系统所应具备的功能和相应的环境条件,以确定实现系统目标的最佳方案。 关键词:环境系统分析;环境生态;环境质量模型. 环境系统分析是以环境质量的变化规律、污染物对人体和生态的影响、环境自净能力以及有关环境工程技术原理为依据,运用系统工程学的理论和方法,研究如何建立起一个合理的环境污染预防控制系统的数学模型,并研究如何利用它来分析各种污染控制过程可调因素(或各种可替换方案)对环境目标或费用、能耗等的影响,以及寻求最优决策方案。 环境系统分析的理论基础和专门技术基础。理论基础:环境科学、环境经济学、环境工程学和系统工程学的基本理论(如运筹学)。专门技术基础:数学建模、计算科学、环境影响评估方法、生命周期评估、系统化的图与网络分析方法。 《环境系统分析》是我国高等学校环境工程、环境科学专业的一门专业基础课,课程任务和教学目标包括:1.使

学生了解污染物在水体和大气中的迁移、扩散和变化规律,建立相应的环境系统模型;2.使学生掌握建立环境数学模型的一般知识;3.使学生了解湖泊、水库水体富营养化的原因和水体富营养化的控制技术;4.使学生掌握区域性环境污染控制系统规划的基本原理和方法;5.使学生建立采用最优化技术求解水污染控制系统规划问题的概念,并有能力解决一般性问题。本门课程一共有十一个章节,主要内容有:环境系统分析概论、数学模型概述、环境质量基本模型、水体水质模型、流域非点源模型、大气质量模型、环境质量评价方法与模型以及环境规划,还有环境决策分析。其中,水体水质模型主要指内陆水体模型,包括湖泊水库水质模型和河流水质模型;环境规划包括水环境规划和大气环境规划。 通过对这门课程的学习,我们对环境系统的分析方法有了一定的了解,它的最大特征是追求环境系统的最优化。 环境系统分析的最优化方法的选用主要有对确定性问题,可采用线性规划、动态规划、非线性规划、整体规划等。对非确定性问题,可用马尔可夫过程,排队论,对策论等方法进行最优化。有的系统优化问题还应用网络理论、图论和模糊数学等进行最优化。 根据本书前言部分介绍这门课程的学科基础包括数学、运筹学、环境科学与环境工程学等,内容较为丰富,通过选

OFDM技术仿真(MATLAB代码)

第一章绪论 1.1简述 OFDM是一种特殊的多载波传输方案,它可以被看作是一种调制技术,也可以被当作一种复用技术。多载波传输把数据流分解成若干子比特流,这样每个子数据流将具有低得多的比特速率,用这样的低比特率形成的低速率多状态符号再去调制相应的子载波,就构成多个低速率符号并行发送的传输系统。正交频分复用是对多载波调制(MCM,Multi-Carrier Modulation)的一种改进。它的特点是各子载波相互正交,所以扩频调制后的频谱可以相互重叠,不但减小了子载波间的干扰,还大大提高了频谱利用率。 符号间干扰是多径衰落信道宽带传输的主要问题,多载波调制技术包括正交频分复用(OFDM)是解决这一难题中最具前景的方法和技术。利用OFDM技术和IFFT方式的数字实现更适宜于多径影响较为显著的环境,如高速WLAN 和数字视频广播DVB等。OFDM作为一种高效传输技术备受关注,并已成为第4代移动通信的核心技术。如果进行OFDM系统的研究,建立一个完整的OFDM 系统是必要的。本文在简要介绍了OFDM 基本原理后,基于MATLAB构建了一个完整的OFDM动态仿真系统。 1.2 OFDM基本原理概述 1.2.1 OFDM的产生和发展 OFDM的思想早在20世纪60年代就已经提出,由于使用模拟滤波器实现起来的系统复杂度较高,所以一直没有发展起来。在20世纪70年代,提出用离散傅里叶变换(DFT)实现多载波调制,为OFDM的实用化奠定了理论基础;从此以后,OFDM在移动通信中的应用得到了迅猛的发展。 OFDM系统收发机的典型框图如图1.1所示,发送端将被传输的数字信号转换成子载波幅度和相位的映射,并进行离散傅里叶变换(IDFT)将数据的频谱表达式变换到时域上。IFFT变换与IDFT变换的作用相同,只是有更高的计算效

实验三___离散时间系统的时域分析

实验三 离散时间系统的时域分析 1.实验目的 (1)理解离散时间信号的系统及其特性。 (2)对简单的离散时间系统进行分析,研究其时域特性。 (3)利用MATLAB对离散时间系统进行仿真,观察结果,理解其时域特性。 2.实验原理 离散时间系统,主要是用于处理离散时间信号的系统,即是将输入信号映射成的输出的某种运算,系统的框图如图所示: (1)线性系统 线性系统就是满足叠加原理的系统。如果对于一个离散系统输入信号为时,输出信号分别为,即:。 而且当该系统的输入信号为时,其中a,b为任意常数,输出为,则该系统就是一个线性离散时间系统。 (2)时不变系统 如果系统的响应与激励加于系统的时刻无关,则该系统是时不变系统。对于一个离散时间系统,若输入,产生输出为,则输入为,产生输出为,即: 若,则。 通常我们研究的是线性时不变离散系统。 3.实验内容及其步骤 (1)复习离散时间系统的主要性质,掌握其原理和意义。 (2)一个简单的非线性离散时间系统的仿真 系统方程为: x = cos(2*pi*0.05*n); x1[n] = x[n+1] x2[n] = x[n] x3[n] = x[n-1] y = x2.*x2-x1.*x3; 或者:y=x*x- x[n+1]* x[n-1] 是非线性。 参考:% Generate a sinusoidal input signal clf; n = 0:200; x = cos(2*pi*0.05*n); % Compute the output signal x1 = [x 0 0]; % x1[n] = x[n+1] x2 = [0 x 0]; % x2[n] = x[n] x3 = [0 0 x]; % x3[n] = x[n-1]

OFDM系统设计及其Matlab实现

课程设计 。 课程设计名称:嵌入式系统课程设计 专业班级: 07级电信1-1 学生姓名:__王红__________ 学号:_____107_____ 指导教师:李国平,陈涛,金广峰,韩琳 课程设计时间:— |

1 需求分析 运用模拟角度调制系统的分析进行频分复用通信系统设计。从OFDM系统的实现模型可以看出,输入已经过调制的复信号经过串/并变换后,进行IDFT或IFFT和并/串变换,然后插入保护间隔,再经过数/模变换后形成OFDM调制后的信号s(t)。该信号经过信道后,接收到的信号r(t)经过模/数变换,去掉保护间隔,以恢复子载波之间的正交性,再经过串/并变换和DFT或FFT后,恢复出OFDM的调制信号,再经过并/串变换后还原出输入符号 2 概要设计 1.简述OFDM通信系统的基本原理 2.简述OFDM的调制和解调方法 3.概述OFDM系统的优点和缺点 4.基于MATLAB的OFDM系统的实现代码和波形 : 3 运行环境 硬件:Windows XP 软件:MATLAB 4 详细设计 OFDM基本原理 一个完整的OFDM系统原理如图1所示。OFDM的基本思想是将串行数据,并行地调制在多个正交的子载波上,这样可以降低每个子载波的码元速率,增大码元的符号周期,提高系统的抗衰落和干扰能力,同时由于每个子载波的正交性,大大提高了频谱的利用率,所以非常适合移动场合中的高速传输。

在发送端,输入的高比特流通过调制映射产生调制信号,经过串并转换变成N条并行的低速子数据流,每N个并行数据构成一个OFDM符号。插入导频信号后经快速傅里叶反变换(IFFT)对每个OFDM符号的N个数据进行调制,变成时域信号为: [ 式 式1中:m为频域上的离散点;n为时域上的离散点;N为载波数目。为了在接收端有效抑制码间干扰(InterSymbol Interference,ISI),通常要在每一时域OFDM符号前加上保护间隔(Guard Interval,GI)。加保护间隔后的信号可表示为式,最后信号经并/串变换及D/A转换,由发送天线发送出去。 式 接收端将接收的信号进行处理,完成定时同步和载波同步。经A/D转换,串并转换后的信号可表示为:

离散控制系统分析方法

实验二 离散控制系统分析方法 一、实验目的 利用MATLAB 对各种离散控制系统进行时域分析。 二、实验指导 1.控制系统的稳定性分析 由前面章节学习的内容可知,对线性系统而言,如果一个连续系统的所有极点都位于s 平面的左半平面,则该系统是一个稳定系统。对离散系统而言,如果一个系统的全部极点都位于z 平面的单位圆内部,则该系统是一个稳定系统。一个连续的稳定系统,如果所有的零点都位于s 平面的左半平面,即所有零点的实部小于零,则该系统是一个最小相位系统。一个离散的稳定系统,如果所有零点都位于z 平面的单位圆内,则称该系统是一个最小相位系统。由于Matlab 提供了函数可以直接求出控制系统的零极点,所以使用Matlab 判断一个系统是否为最小相位系统的工作就变得十分简单。 2.控制系统的时域分析 时域分析是直接在时间域对系统进行分析。它是在一定输入作用下,求得输出量的时域表达式,从而分析系统的稳定性、动态性能和稳态误差。这是一种既直观又准确的方法。 Matlab 提供了大量对控制系统的时域特征进行分析的函数,适用于用传递函数表示的模型。其中常用的函数列入表1,供学生参考。 例1.z z z H 5.05 .1)(2+= 试绘出其单位阶跃响应及单位斜波输入响应。 解:为求其单位阶跃响应及单位斜波输入响应,编制程序如下: num=[1.5]; den=[1 0.5 0];sysd=tf(num,den,0.1) [y,t,x]=step(sysd);

subplot(1,2,1) plot(t,y); xlabel('Time-Sec'); ylabel('y(t)'); gtext('单位阶跃响应') grid; u=0:0.1:1; subplot(1,2,2) [y1,x]=dlsim(num,den,u); plot(u,y1) xlabel('Time-Sec'); ylabel('y(t)'); gtext('单位速度响应') grid 二、实验内容 1、MATLAB在离散系统的分析应用 对于下图所示的计算机控制系统结构图1,已知系统采样周期为T=0.1s,被 控对象的传递函数为 2 () s(0.11)(0.05s1) G s s = ++ ,数字控制器 0.36 () 0.98 z D z z - = + ,试 求该系统的闭环脉冲传递函数和单位阶跃响应。 图1 计算机控制系统结构图 实验步骤: 1).求解开环脉冲传递函数,运用下面的matlab语句实现:>> T=0.1; >> sys=tf([2],[0.005 0.15 1 0]); %将传函分母展开>> sys1=c2d(sys,T,'zoh'); >> sys2=tf([1 -0.36],[1 0.98],0.1); >> sys3=series(sys2,sys1) 执行语句后,屏幕上显示系统的开环脉冲传递函数为: sys3 = 0.03362 z^3 + 0.05605 z^2 - 0.01699 z - 0.002717 --------------------------------------------------

无线通信原理 基于matlab的ofdm系统设计与仿真..

基于matlab的ofdm系统设计与仿真

摘要 OFDM即正交频分复用技术,实际上是多载波调制中的一种。其主要思想是将信道分成若干正交子信道,将高速数据信号转换成并行的低速子数据流,调制到相互正交且重叠的多个子载波上同时传输。该技术的应用大幅度提高无线通信系统的信道容量和传输速率,并能有效地抵抗多径衰落、抑制干扰和窄带噪声,如此良好的性能从而引起了通信界的广泛关注。 本文设计了一个基于IFFT/FFT算法与802.11a标准的OFDM系统,并在计算机上进行了仿真和结果分析。重点在OFDM系统设计与仿真,在这部分详细介绍了系统各个环节所使用的技术对系统性能的影响。在仿真过程中对OFDM信号使用QPSK调制,并在AWGN信道下传输,最后解调后得出误码率。整个过程都是在MATLAB环境下仿真实现,对ODFM系统的仿真结果及性能进行分析,通过仿真得到信噪比与误码率之间的关系,为该系统的具体实现提供了大量有用数据。

第一章 ODMF 系统基本原理 1.1多载波传输系统 多载波传输通过把数据流分解为若干个子比特流,这样每个子数据流将具有较低的比特速率。用这样的低比特率形成的低速率多状态符号去调制相应的子载波,构成了多个低速率符号并行发送的传输系统。在单载波系统中,一次衰落或者干扰就会导致整个链路失效,但是在多载波系统中,某一时刻只会有少部分的子信道会受到衰落或者干扰的影响。图1-1中给出了多载波系统的基本结构示意图。 图1-1多载波系统的基本结构 多载波传输技术有许多种提法,比如正交频分复用(OFDM)、离散多音调制(DMT)和多载波调制(MCM),这3种方法在一般情况下可视为一样,但是在OFDM 中,各子载波必须保持相互正交,而在MCM 则不一定。 1.2正交频分复用 OFDM 就是在FDM 的原理的基础上,子载波集采用两两正交的正弦或余弦函数集。函数集{t n ωcos }, {t m ωsin } (n,m=0,1,2…)的正交性是指在区间(T t t +00,)内有正弦函数同理:)0()()(2/0cos *cos 00===≠?? ???=? +m n m n m n T T tdt m t n T t t ωω 其中ωπ2=T (1-1)

离散时间系统的时域分析

第七章离散时间系统的时域分析 §7-1 概述 一、离散时间信号与离散时间系统 离散时间信号:只在某些离散的时间点上有值的 信号。 离散时间系统:处理离散时间信号的系统。 混合时间系统:既处理离散时间信号,又处理连 续时间信号的系统。 二、连续信号与离散信号 连续信号可以转换成离散信号,从而可以用离散时间系统(或数字信号处理系统)进行处理: 三、离散信号的表示方法:

1、 时间函数:f(k)<——f(kT),其中k 为序号,相当于时间。 例如:)1.0sin()(k k f = 2、 (有序)数列:将离散信号的数值按顺序排列起来。例如: f(k)={1,0.5,0.25,0.125,……,} 时间函数可以表达任意长(可能是无限长)的离散信号,可以表达单边或双边信号,但是在很多情况下难于得到;数列的方法表示比较简单,直观,但是只能表示有始、有限长度的信号。 四、典型的离散时间信号 1、 单位样值函数:? ??==其它001)(k k δ 下图表示了)(n k ?δ的波形。

这个函数与连续时间信号中的冲激函数 )(t δ相似,也有着与其相似的性质。例如: )()0()()(k f k k f δδ=, )()()()(000k k k f k k k f ?=?δδ。 2、 单位阶跃函数:? ??≥=其它001)(k k ε 这个函数与连续时间信号中的阶跃函数)(t ε相似。用它可以产生(或表示)单边信号(这里称为单边序列)。 3、 单边指数序列:)(k a k ε

比较:单边连续指数信号:)()()(t e t e t a at εε=,其 底一定大于零,不会出现负数。 (a) 0.9a = (d) 0.9a =? (b) 1a = (e) 1a =? (c) 1.1a = (f) 1.1a =?

离散控制系统分析方法

实验二离散控制系统分析方法 一、实验目的 利用MATLAB对各种离散控制系统进行时域分析。 二、实验指导 1.控制系统的稳定性分析 由前面章节学习的内容可知,对线性系统而言,如果一个连续系统的所有极点都位于s平面的左半平面,则该系统是一个稳定系统。对离散系统而言,如果一个系统的全部极点都位于z平面的单位圆内部,则该系统是一个稳定系统。一个连续的稳定系统,如果所有的零点都位于s平面的左半平面,即所有零点的实部小于零,则该系统是一个最小相位系统。一个离散的稳定系统,如果所有零点都位于z平面的单位圆内,则称该系统是一个最小相位系统。由于Matlab提供了函数可以直接求出控制系统的零极点,所以使用Matlab判断一个系统是否为最小相位系统的工作就变得十分简单。 2.控制系统的时域分析 时域分析是直接在时间域对系统进行分析。它是在一定输入作用下,求得输出量的时域表达式,从而分析系统的稳定性、动态性能和稳态误差。这是一种既直观又准确的方法。 Matlab提供了大量对控制系统的时域特征进行分析的函数,适用于用传递函数表示的模型。其中常用的函数列入表1,供学生参考。

例1.z z z H 5.05 .1)(2+= 试绘出其单位阶跃响应及单位斜波输入响应。 解:为求其单位阶跃响应及单位斜波输入响应,编制程序如下: num=[1.5]; den=[1 0.5 0];sysd=tf(num,den,0.1) [y,t,x]=step(sysd); subplot(1,2,1) plot(t,y); xlabel('Time-Sec'); ylabel('y(t)'); gtext('单位阶跃响应') grid; u=0:0.1:1; subplot(1,2,2) [y1,x]=dlsim(num,den,u); plot(u,y1) xlabel('Time-Sec'); ylabel('y(t)'); gtext('单位速度响应') grid 二、 实验内容 1、MATLAB 在离散系统的分析应用 对于下图所示的计算机控制系统结构图1,已知系统采样周期为T=0.1s ,被

用MATLAB实现OFDM仿真分析

3.1 计算机仿真 仿真实验是掌握系统性能的一种手段。它通过对仿真模型的实验结果来确定实际系统的性能。从而为新系统的建立或系统的改进提供可靠的参考。通过仿真,可以降低新系统失败的可能性,消除系统中潜在的瓶颈。优化系统的整体性能,衡量方案的可行性。从中选择最后合理的系统配置和参数配置。然后再应用于实际系统中。因此,仿真是科学研究和工程建设中不可缺少的方法。 3.1.1 仿真平台 ●硬件 CPU:Pentium III 600MHz 内存:128M SDRAM ●软件 操作系统:Microsoft Windows2000 版本5.0 仿真软件:The Math Works Inc. Matlab 版本6.5 包括MATLAB 6.5的M文件仿真系统。 Matlab是一种强大的工程计算软件。目前最新的6.x版本 (windows环境)是一种功能强、效率高、便于进行科学和工程计算的交互式软件包。其工具箱中包括:数值分析、矩阵运算、通信、数字信号处理、建模和系统控制等应用工具程序,并集应用程序和图形于一便于使用的集成环境中。在此环境下所解问题的Matlab语言表述形式和其数学表达形式相同,不需要按传统的方法编程。Matlab的特点是编程效率高,用户使用方便,扩充能力强,语句简单,内涵丰富,高效方便的矩阵和数组运算,方便的绘图功能。 3.1.2 基于MATLAB的OFDM系统仿真链路 根据OFDM 基本原理,本文给出利用MATLAB编写OFDM系统的仿真链路流程。串行数据经串并变换后进行QDPSK数字调制,调制后的复信号通过N点IFFT变换,完成多载波调制,使信号能够在N个子载波上并行传输,中间插入10训练序列符号用于信道估计,加入循环前缀后经并串转换、D /A后进入信道,接收端经过N点FFT变换后进行信道估计,将QDPSK解调后的数据并串变换后得到原始信息比特。 本文采用MATLAB语言编写M文件来实现上述系统。M文件包括脚本M文件和函数M文件,M文件的强大功能为MATLAB的可扩展性提供了基础和保障,使MATLAB能不断完善和壮大,成为一个开放的、功能强大的实用工具。M文件通过input命令可以轻松实现用户和程序的交互,通过循环向量化、数组维数预定义等提高M文件执行速度,优化内存管理,此外,还可以通过类似C++语言的面向对象编程方法等等。

离散时间系统的分析

课程设计报告 课程设计题目:离散时间系统分析学号:201420130206 学生姓名:董晓勇 专业:通信工程 班级:1421301 指导教师:涂其远 2015年12月18日

离散时间系统的分析 一、设计目的和意义 1 . 目的: (1)深刻理解卷积和、相加、相乘运算,掌握求离散序列卷积和、相加相乘的计算方法;(2)加深理解和掌握求离散序列Z变换的方法; (3)加深和掌握离散系统的系统函数零点、函数极点和系统时域特性、系统稳定性的关系。 2 . 意义: 在对《信号与系统》一书的学习中,进行信号与系统的分析是具有十分重要的意义,同时也是必不可少的。利用matlab函数,只需要简单的编程,就可以实现系统的时域、频域分析,对系统特性进行分析,为实际的系统设计奠定了基础。本设计在离散系统Z域分析理论的基础上,利用matlab对离散系统的稳定性和频域响应进行了分析。 二、设计原理

第一部分:对离散时间系统的时域进行分析呈 对离散时间信号的代数运算(相加、相乘、卷积和),是在时域进行分析。相加用“+”来完成,相乘用“·*”来完成,卷积和则用conv 函数来实现,具体形式为y=conv(x1,x2,….),其中x1,x2,…..为输入的离散序列 ,y 为输出变量。 在零初始状态下,matlab 控制工具箱提供了一个filter 函数,可以计算差分方程描述的系统的响应,其调用形式为: y=filter(b,a,f) 其中,a=[a0,a1,a2,…]、b=[b0,b1,b2,….]分别是系统方程左、右边的系数向量,f 表示输入向量,y 表示输出向量。 第二部分:对离散时间系统的Z 域进行分析 matlab 工具箱提供了计算Z 正变换的函数ztrans,其调用形式为: F=zrtans(f) %求符号函数f 的Z 变换,返回函数的自变量为z 。 Matlab 的zplane 函数用于系统函数的零极点图的绘制,调用方式为: zplane(b,a)其中,b 、a 分别为系统函数分子、分母多项式的系数向量。 matlab 中,利用freqz() 函数可方便地求得系统的频率响应,调用格式为: freqz(b,a,N) 该调用方式将绘制系统在0~PI 范围内N 个频率等分点的幅频特性和相频特性图。 三、 详细设计步骤 1.自己设计两个离散时间序列x1、x2,对其进行相加,相乘,卷积运算,并显示出图形。 2.根据已知的LTI 系统:y[n]-0.7y[n-1]-0.6y[n-2]+y[n-3]=x[n]+0.5[n-1],得其在Z 域输 入输出的传递函数为: 1 12310.5()10.70.6z H z z z z ----+= --+ 利用matlab 求:(1)系统函数的零点和极点,并在z 平面显示他们的分布;(2)画出幅频响应和相频响应的特性曲线。 四、 设计结果及分析 (1).自行设计产生两个离散序列信号,对其进行相加、乘及卷积运算

离散LSI系统分析

信号与系统分析实验报告 实验项目名称:离散线性时不变系统分析; 连续时间系统分析 所属课程名称:信号与系统实验教程 实验类型:验证型 指导教师: 实验日期:2013.06.04 班级: 学号: 姓名:

离散线性时不变系统分析 一、实验目的 1. 掌握离散线性时不变系统的单位序列响应、单位阶跃响应和任意激励下响应的MATLAB 求解方法。 2. 掌握离散线性时不变系统的频域分析方法; 3. 掌握离散线性时不变系统的复频域分析方法; 4. 掌握离散线性时不变系统的零极点分布与系统特性的关系。 二、实验原理及方法 1. 离散线性时不变系统的时域分析 描述一个N 阶线性时不变离散时间系统的数学模型是线性常系统差分方程,N 阶线性时不变离散系统的差分方程一般形式为 ) ()(0 i n x b k n y a M i i N k k -=-∑∑== (2.1) 也可用系统函数来表示 12001212120 () ()()() ()1M i M i i M N N k N k k b z b b z b z b z Y z b z H z X z a z a z a z a z a z ----=----=++++== == ++++∑∑ (2.2) 系统函数()H z 反映了系统响应和激励间的关系。一旦上式中k a ,i b 的数据确定了,系统的性质也就确定了。特别注意0a 必须进行归一化处理,即01a =。 对于复杂信号激励下的线性系统,可以将激励信号在时域中分解为单位序列或单位阶跃 序列的线性叠加,把这些单元激励信号分别加于系统求其响应,然后把这些响应叠加, 即可得到复杂信号作用于系统的零状态响应。因此,求解系统的单位序列响应和单位阶跃响应尤为重要。由图2-1可以看出一个离散LSI 系统响应与激励的关系。 ()()() z X z H z =()()*() n x n h n 图2-1 离散LSI 系统响应与激励的关系 (1) 单位序列响应(单位响应) 单位响应()h n 是指离散线性时不变系统在单位序列()n δ激励下的零状态响应,因此 ()h n 满足线性常系数差分方程(2.1)及零初始状态,即 ()() N M k i k i a h n k b n i δ==-=-∑∑, (1)(2)0h h -=-== (2.3) 按照定义,它也可表示为 ()()()h n h n n δ=* (2.4) 对于离散线性时不变系统,若其输入信号为()x n ,单位响应为()h n ,则其零状态响应

连续和离散系统分析

实验一连续与离散系统分析 一、实验目得 学习连续系统与离散系统响应得matlab求解方法; 二、实验主要仪器设备与材料 计算机 三、实验方法、步骤及结果测试 实验方法:编程,上机调试,分析实验结果; 步骤: 编程实现上述各实验内容 四、实验结果 1、某系统得传递函数为: 试求系统得冲激响应与阶跃响应。 2、编制程序求解下列两个系统得单位冲激响应与阶跃响应,并绘出其图形。要求

分别用filter、conv、impz三种函数完成。给出理论计算结果与程序计算结果并讨论。 (I) 理论计算结果: 程序计算结果: A:单位冲激响应 (1)用Filter函数(2)用Conv函数 (3)用impz函数 单位冲激响应: n 0 1 2 3 4 5 h(n) 1 -1、75 1、19 -0、67 0、355 -0、18 单位阶跃响应: n 0 1 2 3 4 5 y(n) 1 -0、75 0、44 -0、234 0、12 -0、06

B:单位阶跃响应(1)用Fil ter 函数 (2)用Conv 函数 (3)用Imp z函数 (II ) 理论计算结果: 程序计算结果: A:单位冲激响应(1)用f ilter 函数 单位冲激响应: n 0 1 2 3 4 5 h(n) 0 0、25 0、25 0、25 0、25 单位阶跃响应: N 0 1 2 3 4 5 y(n) 0 0、25 0、5 0、75 1 1

(2)用Conv函数 (3)用Impz函数 B:单位阶跃响应 (1)用filter函数 (2)用Conv函数 (3)用Impz函数

连续系统离散化处理基本方法

在数字计算机上对连续系统进行仿真时,首先遇到的问题是如何解决数字计算机在数值及时间上的离散性与被仿真系统数值及时间上的连续性这一基本问题。 从根本意义上讲,数字计算机所进行的数值计算仅仅是“数字”计算,它表示数值的精度受限于字长,这将引入舍入误差;另一方面,这种计算是按指令一步一步进行的,因而,还必须将时间离散化,这样就只能得到离散时间点上系统性能。用数字仿真的方法对微分方程的数值积分是通过某种数值计算方法来实现的。任何一种计算方法都只能是原积分的一种近似。因此,连续系统仿真,从本质上是对原连续系统从时间、数值两个方面对原系统进行离散化,并选择合适的数值计算方法来近似积分运算,由此得到的离散模型来近似原连续模型。如何保证离散模型的计算结果从原理上确能代表原系统的行为,这是连续系统数字仿真首先必须解决的问题。 设系统模型为:),,(t u y f y =&,其中u (t )为输入变量,y (t )为系统变量;令仿真时间间隔为h ,离散化后的输入变量为)(?k t u ,系统变量为)(?k t y ,其中k t 表示t=kh 。如果)()(?k k t u t u ≈,)()(?k k t y t y ≈,即0)()(?)(≈-=k k k u t u t u t e ,0)()(?)(≈-=k k k y t y t y t e (对所有k=0,1,2,…),则可认为两模型等价,这称为相似 原理(参见图)。 实际上,要完全保证0)(,0)(==k y k u t e t e 是很困难的。进一步分析离散化引的误差,随着计算机技术的发展,由计算机字长引入的舍入误差可以忽略,关键是数值积分算法,也称为仿真建模方法。相似原理用于仿真时,对仿真建模方法有三个基本要求: (1)稳定性:若原连续系统是稳定的,则离散化后得到的仿真模型也应是稳定的。关于稳定性的详细讨论将在节中进行。 (2)准确性:有不同的准确性评价准则,最基本的准则是: 绝对误差准则:δ≤-=)()(?)(k k k y t y t y t e 相对误差准则:δ≤-= )(?)()(?)(k k k k y t y t y t y t e 其中 规定精度的误差量。 原连续模型 仿真模型 )(≈k y t e 图 相

2010年本科毕业设计:基于MATLAB的OFDM系统仿真及分析

2010年本科毕业设计:基于MATLAB的OFDM系统仿真及分 析 MATLABOFDM 正交频分复用(OFDM) 是第四代移动通信的核心技术。该文首先简要介绍了OFDM的发展状况及基本原理, 文章对OFDM 系统调制与解调技术进行了解析,得 到了OFDM 符号的一般表达式,给出了OFDM 系统参数设计公式和加窗技术的原理 及基于IFFT/FFT 实现的OFDM 系统模型,阐述了运用IDFT 和DFT 实现OFDM 系统的根源所在,重点研究了理想同步情况下,保护时隙(CP)、加循环前缀前后和不同的信道内插方法在高斯信道和多径瑞利衰落信道下对OFDM系统性能的影响。在给出OFDM系统模型的基础上,用MATLAB语言实现了传输系统中的计算机仿真并给出 参考设计程序。最后给出在不同的信道条件下,研究保护时隙、循环前缀、信道 采用LS估计方法对OFDM系统误码率影响的比较曲线,得出了较理想的结论。 : 正交频分复用;仿真;循环前缀;信道估计 I Title: MATLAB Simulation and Performance Analysis of OFDM System ABSTRACT OFDM is the key technology of 4G in the field of mobile communication. In this

article OFDM basic principle is briefly introduced. This paper analyzes the modulation and demodulation of OFDM system, obtaining a general expression of OFDM mark, and giving the design formulas of system parameters, principle of windowing technique, OFDM system model based on IFFT/FFT, the origin which achieves the OFDM system by using IDFT and DFT. Then, the influence of CP and different channel estimation on the system performance is emphatically analyzed respectively in Gauss and Rayleigh fading channels in the condition of ideal synchronization. Besides, based on the given system model OFDM system is computer simulated with MATLAB language and the referential design procedure is given. Finally, the BER curves of CP and channel estimation are given and compared. The conclusion is satisfactory. KEYWORDS:OFDM; Simulation; CP; Channel estimation II

实验七--离散系统分析的MATLAB实现讲解学习

实验七 离散系统分析的MATLAB 实现 一、实验目的 1、掌握利用MATLAB 绘制系统零极点图的方法; 2、掌握离散时间系统的零极点分析方法; 3、学习离散系统响应的MATLAB 求解方法; 4、掌握用MATALB 实现离散系统频率特性分析的方法; 5、深刻理解离散系统的系统函数零极点对系统频响的影响,可以根据 零极点知识设计简单的滤波器。 二、基本原理 (一)离散系统零极点 线性时不变离散系统可用线性常系数差分方程描述,即 ()()N M i j i j a y n i b x n j ==-=-∑∑ (1) 其中()y k 为系统的输出序列,()x k 为输入序列。 将式(1)两边进行Z 变换, 00 () () ()() () M j j j N i i i b z Y z B z H z X z A z a z -=-== == ∑∑ (2) 将式(2)因式分解后有: 11 () ()() M j j N i i z q H z C z p ==-=-∏∏ (3) 其中C 为常数,(1,2,,)j q j M =L 为()H z 的M 个零点,(1,2,,)i p i N =L 为()H z 的 N 个极点。 系统函数()H z 的零极点分布完全决定了系统的特性,若某系统函数的零极点已知,则系统函数便可确定下来。因此,系统函数的零极点分布对离散系统特性的分析具有非常重要意义。 (二)离散系统零极点图及零极点分析 1、零极点图的绘制 设离散系统的系统函数为 () ()() B z H z A z =

则系统的零极点可用MATLAB 的多项式求根函数roots()来实现,调用格式为: p=roots(A) 其中A 为待求根多项式的系数构成的行矩阵,返回向量p 则是包含多项式所有根的列向量。如多项式为231 ()48 B z z z =+ +,则求该多项式根的MATLAB 命令为为: A=[1 3/4 1/8]; P=roots(A) 运行结果为: P = -0.5000 -0.2500 需注意的是,在求系统函数零极点时,系统函数可能有两种形式:一种是分子、分母多项式均按z 的降幂次序排列;另一种是分子、分母多项式均按1z -的升幂次序排列。这两种方式在构造多项式系数向量时稍有不同。 (1)()H z 按z 的降幂次序排列:系数向量一定要由多项式最高次幂开始,一直到常数项,缺项要用0补齐。如 34322()3221 z z H z z z z z +=++++ 其分子、分母多项式系数向量分别为A=[1 0 2 0]、B=[1 3 2 2 1]。 (2)()H z 按1z -的升幂次序排列:分子和分母多项式系数向量的维数一定要相同,不足的要用0补齐,否则0z =的零点或极点就可能被漏掉。如 1 1212()11124 z H z z z ---+=++ 其分子、分母多项式系数向量分别为A=[1 2 0]、B=[1 1/2 1/4]。 用roots()求得()H z 的零极点后,就可以用plot()函数绘制出系统的零极点图。下面是求系统零极点,并绘制其零极点图的MATLAB 实用函数ljdt(),同时还绘 制出了单位圆。函数ljdt()的程序如下: function ljdt(A,B) % The function to draw the pole-zero diagram for discrete system p=roots(A) %求系统极点 q=roots(B) %求系统零点 p=p'; %将极点列向量转置为行向量 q=q'; %将零点列向量转置为行向量 x=max(abs([p q 1])); %确定纵坐标范围

连续离散系统频域分析

课程实验报告 学年学期2015-2016学年第二学期 课程名称信号与系统 实验名称连续和离散系统的频域分析实验室北校区5号楼计算机房 专业年级电气141 学生姓名宋天绍 学生学号2014011595 提交时间 成绩 任课教师吴凤娇 水利与建筑工程学院

实验二:连续和离散系统的频域分析 一:实验目的 1:学习傅里叶正变换和逆变换,理解频谱图形的物理含义 2:了解连续和离散时间系统的单位脉冲响应 3:掌握连续时间系统的频率特性 二:实验原理 1. 傅里叶正变换和逆变换公式 正变换:()()j t F f t e dt ωω∞ --∞ =? 逆变换:1()()2j t f t F e d ωωωπ ∞ -∞ = ? 2. 频域分析 t j t j e d d e t e ωωωπ ωωωπ??∞∞-∞∞-E =E =)(21)(21)(将激励信号分解为无穷多个正弦分量的和。 ?∞∞-H E =ωωωπωd e t r t j zs )()(21)(,R(ω)为)(t r zs 傅里叶变换;π ωωd )(E 各频率分量的复数振幅 激励单位冲激响应时的零状态响应→ )(t δ)(t h 单位阶跃响应时的零状态响应激励→)(t u )(t g 3 各函数说明: (1)impulse 冲激响应函数:[Y,X,T]=impulse(num,den); ) 1()2()1() 1()2()1()()()(1 1++++++++==--n a s a s a m b s b s b s A s B s H n n m m num 分子多项式系数; num=[b(1) b(2) … b(n+1)]; den 分母多项式系数; den=[a(1) a(2) … a(n+1)]; Y,X,T 分别表示输出响应,中间状态变量和时间变量; 如:3 52 )(2 +++= s s s s H ,等价于)(2)()(3)(5)(t e t e t r t r t r +=++ 定义den=[1 5 3];num=[1 2]; [Y,X,T]=impulse(num,den); (2)step 阶跃响应函数:[Y,X,T]=step(num,den);num 分子多项式;den 分母多项式 Y,X,T 分别表示输出响应,中间状态变量和时间变量; 如:3 52 )(2+++= s s s s H ,den=[1 5 3];num=[1 2]; [Y,X,T]= step (num,den); (3)impz 数字滤波器的冲激响应 [h,t]=impz(b,a,n) b 分子多项式系数;a 分母多项式系数;n 采样样本 h 离散系统冲激响应;t 冲激时间,其中t=[0:n-1]', n=length(t)时间样本数

相关文档
最新文档