信息与通信Matlab与数字通信系统仿真课件

Matlab基础与数字通信仿真赵婷婷Matlab基础与数字通信仿真赵婷婷内容Matlab简介1Matlab与数字通信仿真4Matlab语言基础2Matlab绘图基础3内容Matlab简介1Matlab与数字通信仿真4Matla内容Matlab简介1Matlab窗口环境Matlab是什么？内容Matlab简介1Matlab窗口环境Matlab是什Matlab是什么？Matlab=MatrixLaboratory=矩阵实验室由MathWorks公司于1984年推出的一套高性能的数值计算和可视化数学软件基本元素是矩阵数学表达式与运算规则符合通常的习惯语法不严格具有强大的数值计算和图示能力具有丰富的工具箱(ToolBox)强大的help功能Matlab是什么？Matlab=MatrixLaboMatlab窗口环境脚本M文件在命令窗口中输入并执行它所用的变量都要在工作空间中获取Matlab窗口环境脚本M文件Matlab窗口环境程序文件以.m格式进行存取程序运行时只需在工作空间中键入其名称即可Matlab窗口环境程序文件MATLAB语句形式变量＝表达式；通过等号将表达式的值赋予变量。当键入回车键时，该语句被执行。语句执行之后，窗口自动显示出语句执行的结果。如果不希望结果回显，只需在语句后加分号即可。Matlab窗口环境MATLAB语句形式Matlab窗口环境MATLAB语句形式变量＝表达式；通过等号将表达式的值赋予变量。当键入回车键时，该语句被执行。语句执行之后，窗口自动显示出语句执行的结果。如果不希望结果回显，只需在语句后加分号即可。Matlab窗口环境MATLAB语句形式Matlab窗口环境命令行编辑器方向键和控制键可以编辑修改已输入的命令↑：回调上一行命令↓：回调下一行命令多行命令（…):如果命令语句超过一行或者太长希望分行输入，可以使用多行命令继续输入

Matlab窗口环境例如：S=1-12+13+4+…9-4-18;命令行编辑器Matlab窗口环境例如：内容Matlab语言基础2基本运算符变量矩阵常用函数基本语句内容Matlab语言基础2基本运算符变量矩阵常用函数Matlab语言基础:变量变量的命名变量无需定义即可使用变量的名字必须以字母开头，之后可以是任意字母、数字或下划线(不能超过31个字符)变量名称区分字母的大小写变量可包含下连字符、数字，但不能包含空格、标点变量操作在命令窗口中，同时存储着输入的命令和创建的所有变量值，它们可以在任何需要的时候被调用采用clear将清除内存内的所有变量采用“clear变量名”清除指定变量Matlab语言基础:变量变量的命名变量操作一些特殊的变量ans：最近生成的无名变量i,j：虚数单位pi：圆周率

nargin：函数的输入变量个数eps：计算机的最小数nargout：函数的输出变量个数inf：无穷大(2/0)

realmin：最小正实数realmax：最大正实数nan：不定量(0/0)flops：浮点运算数Matlab语言基础:变量一些特殊的变量inf：无穷大(2/0)常用的数学运算符运算式中，通常不需要考虑空格多条命令可以放在一行中同行中各条命令可用分号隔开，表示禁止结果显示同行中各条命令也可用逗号隔开，表示需要显示结果＋----加－----减*----乘

/----左除

\----右除

^----幂Matlab语言基础:基本运算符常用的数学运算符＋----加常用的关系运算符

>----大于

>=----大于等于

<----小于

<=

----小于等于

==----等于

~=----不等于运算法则若关系式成立，结果为1若关系式不成立，结果为0Matlab语言基础:基本运算符常用的关系运算符>----大常用的逻辑运算符

&----与

|----或

!----非运算法则若逻辑真，结果为1若逻辑假，结果为0Matlab语言基础:基本运算符常用的逻辑运算符&----与其他常见符号Matlab语言基础:基本运算符

=

----变量赋值%----注释符

‘----共轭转置符

.’----转置符:----冒号运算符举例：

n:s:m产生n~m、步长为s的序列其他常见符号Matlab语言基础:基本运算符=循环语句for-end循环：循环次数事先确定时使用while-end循环：循环次数不能事先确定时使用Matlab语言基础:基本语句fori=n:s:m

语句体ends为步长可以为正数，负数或小数while表达式语句体end表达式为真，执行语句体表达式为假，终止该循环循环语句Matlab语言基础:基本语句fori=n:s:分支语句if-elseif-else-end：具有优先级switch-case-otherwise-end：无优先级Matlab语言基础:基本语句if表达式1

语句体1elseif表达式2

语句体2………..else

语句体elseendswitch变量

case变量值1语句体1case变量值2语句体2……otherwise

语句体otherwiseend分支语句Matlab语言基础:基本语句if表达式1sw常用的数学函数abs----求绝对值sin----正弦cos----余弦tan----正切asin----反正弦acos----反余弦atan----反正切sqrt----求平方根exp----eximag----求虚部real----求实部sign----求符号log----求自然对数log10----以10为底的对数conj----共轭复数Matlab语言基础:常用函数常用的数学函数abs----求绝对值eMatlab语言基础:矩阵矩阵的生成例如：创建一个3×3矩阵输入：a=[123;456;789]

或者a=[123456789]

屏幕输出：a=123456789矩阵的元素直接排列在方括号内行与行之间用分号隔开每行内的元素使用空格或逗号隔开大的矩阵可以用分行输入，回车键代表分号Matlab语言基础:矩阵矩阵的生成例如：创建一个3×3矩矩阵的生成举例：输入a=[1:2:10]屏幕输出：a=13579输入a=[10:-2:1]屏幕输出：a=108642用线性等间距生成向量矩阵start:step:end

start：起始值

step:步长(默认为1，可为负值)

end:终止值Matlab语言基础:矩阵矩阵的生成举例：用线性等间距生成向量矩阵sta一些常用的特殊矩阵eye(2,3)ans=100010eye(2)ans=1001V=[572];A=diag(V)A=500070002Matlab语言基础:矩阵单位矩阵：eye(m,n);eye(m)零矩阵：zeros(m,n);zeros(m)一矩阵：ones(m,n);ones(m)对角矩阵：对角元素向量V=[a1,a2,…,an]A=diag(V)随机矩阵：rand(m,n)产生一个m×n的均匀分布的随机矩阵zeros(2,3)ans=000000zeros(2)ans=0000ones(2,3)ans=111111ones(2)ans=1111rand(2,2)ans=0.41540.87440.01500.9708一些常用的特殊矩阵eye(2,3)eye(2)V=[57矩阵元素的访问通过确认矩阵下标，可以对矩阵的子块进行插入、提取和重排注意：MATLAB中矩阵下标从1开始，而不是从0开始Matlab语言基础:矩阵

A(m,n)：提取第m行，第n列元素A(:,n)：提取第n列元素A(m,:)：提取第m行元素A(m1:m2,n1:n2)：提取m1~m2行和n1~n2列的所有元素A(:)：得到一个长列矢量，其元素按矩阵的列进行排列A=[123456];A(2,1)ans=4A(2,:)ans=456A(:,2)ans=25A(1:2,2:3)ans=2356A(:)ans=142536矩阵元素的访问Matlab语言基础:矩阵A(矩阵的扩展与消除矩阵扩展：如果在原矩阵中一个不存在的地址位置上设定一个数（赋值），则该矩阵会自动扩展行列数，并在该位置上添加这个数，而且在其他没有指定的位置补零。消除子块：如果将矩阵的子块赋值为空矩阵[]，则相当于消除了相应的矩阵子块。Matlab语言基础:矩阵A=[123456];A(3,2)=7A=[123456070]A(2,:)=[]A=[123]矩阵的扩展与消除Matlab语言基础:矩阵A=[12矩阵的大小Matlab语言基础:矩阵[m,n]=size(A,x)：返回矩阵的行列数m与n当x=1，则只返回行数m当x=2，则只返回列数nlength(A)=max(size(A))：返回行数或列数的最大值A=[123456];size(A)ans=23size(A,1)ans=2size(A,2)ans=3矩阵的大小Matlab语言基础:矩阵[m,n]=size(矩阵的运算举例A=[1+i2-2i;-3+3i-4-4i]A=1+i2-2i-3+3i-4-4iB=A’B=1-i-3-3i2+2i-4+4iC=A.’C=1+i-3+3i2-2i-4-4iMatlab语言基础:矩阵转置：转置符号有两种形式

A’----共轭转置A.’----普通转置求逆：inv(A)

求行列式：det(A)矩阵的运算举例B=A’C=A.’Matlab语言基础:矩阵矩阵的四则运算与幂运算

.*----点乘

.\----点右除

./----点左除

.^----点幂点运算是两个维数相同矩阵对应元素间的运算Matlab语言基础:矩阵＋----矩阵加－----矩阵减*----矩阵乘

/----矩阵左除

\----矩阵右除

^----矩阵幂矩阵的四则运算与幂运算.*----内容Matlab绘图基础3内容Matlab绘图基础3Matlab绘图基础MATLAB提供了丰富的绘图功能helpgraph2d可得到所有画二维图形的命令helpgraph3d可得到所有画三维图形的命令Matlab绘图基础MATLAB提供了丰富的绘图功能helpplot(x,y,’属性’)或者stem()：绘图figure：选择图像

figure(1)；figure(2)；…；figure(n)打开不同的图形窗口，以便绘制不同的图形gridon：在所画出的图形坐标中加入栅格

gridoff：除去图形坐标中的栅格

holdon：当前图形保持不变，同时允许在这个坐标内绘制另外一个图形holdoff：使新图覆盖旧的图形Matlab绘图基础plot(x,y,’属性’)或者stem()：绘图Matl举例：

t=[0:pi/20:9*pi];

f=sin(t);

figure(1)

plot(t,f,'r:*')

gridonMatlab绘图基础举例：Matlab绘图基础设定对数坐标Matlab绘图基础semilogx：x轴为对数坐标，y轴为线性坐标semilogy：y轴为对数坐标，x轴为线性坐标设定轴的范围axis（[xminxmaxyminymax]）axis(‘equal’)：将x坐标轴和y坐标轴的单位刻度调整为一样设定对数坐标Matlab绘图基础semilogx：x轴为对数文字标示Matlab绘图基础text(x,y,’字符串’)在图形的指定坐标位置(x,y)处，标示单引号括起来的字符串gtext(‘字符串’)利用鼠标在图形的某一位置标示字符串title(‘字符串’)在所画图形的最上端显示说明该图形标题的字符串xlabel(‘字符串’)，ylabel(‘字符串’)设置x，y坐标轴的名称。legend(‘字符串1’,‘字符串’,…,‘字符串n’)在屏幕上开启一个小视窗，然后依据绘图命令的先后用对应的字符串区分图形上的线文字标示Matlab绘图基础text(x,y,’字符串’)t=[0:pi/20:5*pi];f=sin(t)figure(1)subplot(121)plot(t,f)gridonsubplot(122)plot(t,-f)gridonMatlab绘图基础subplot（mnk）：分割图形显示窗口m:上下分割个数，n:左右分割个数，k:子图编号t=[0:pi/20:5*pi];Matlab绘图基础sub内容Matlab与数字通信仿真4内容Matlab与数字通信仿真4Matlab与数字通信仿真Matlab与数字通信仿真Matlab与数字通信仿真Matlab与数字通信仿真Matlab与数字通信仿真产生随机数randsrc：产生给定符号集合的随机数randn：产生正态分布的随机数randint：产生均匀分布的随机整数

举例：产生0、1分布的随机序列

randsrc(1,N,[0,1;0.5,0.5])randint(1,N)Matlab与数字通信仿真产生随机数randsrc：产生给定Matlab与数字通信仿真产生随机数Matlab与数字通信仿真产生随机数Matlab与数字通信仿真Matlab与数字通信仿真Matlab与数字通信仿真数字调制qammodqamdemodgenqammodgenqamdemodpammod,pamdemod建议自己编写！

举例：进行QPSK调制

输入2比特符号I路Q路00011011-11113-1-11021-1const=[1+j,-1+j,1-j,-1-j];y=genqammod(x,const);Matlab与数字通信仿真数字调制qammod举例：进行Q数字调制Matlab与数字通信仿真数字调制Matlab与数字通信仿真Matlab与数字通信仿真Matlab与数字通信仿真Matlab与数字通信仿真插值upsample建议自己编写！

举例：进行8倍上采样

N=8;y=upsample(x,N);

Matlab与数字通信仿真插值upsample举例：进行8Matlab与数字通信仿真Matlab与数字通信仿真Matlab与数字通信仿真波形成型y=rcosflt(x,Fd,Fs,type_flag,r,delay)filter=rcosine(Fd,Fs,type_flag,r,delay)y=conv(x,filter)建议自己编写！Matlab与数字通信仿真波形成型y=rcosflt(xMatlab与数字通信仿真波形成型Matlab与数字通信仿真波形成型给大家的建议自己去编写课程相关的函数：养成好的习惯：变量、函数的命名有意义，加入注释(%)，增强可读性.m文件（主程序）开头先把环境打扫干净(clc；clear)！注意在子程序中不要用clear尽量使用矩阵运算程序尽量模块化，即采用主程序调用子程序参数值最好集中放在程序的开始部分，以便维护给大家的建议自己去编写课程相关的函数：推荐几本书《现代通信系统－－使用MATLAB》约翰·G·谱罗克斯，马苏德·萨勒赫，刘树棠译西安交通大学出版社《MATLAB与通信仿真》王立宁，乐光新，赞菲人民邮电出版社《精通MATLAB5》张宜华清华大学出版社推荐几本书《现代通信系统－－使用MATLAB》ThankYouforyourattention!ItistimeforQ&A!ThankYouforyourattention!IMatlab基础与数字通信仿真赵婷婷Matlab基础与数字通信仿真赵婷婷内容Matlab简介1Matlab与数字通信仿真4Matlab语言基础2Matlab绘图基础3内容Matlab简介1Matlab与数字通信仿真4Matla内容Matlab简介1Matlab窗口环境Matlab是什么？内容Matlab简介1Matlab窗口环境Matlab是什Matlab是什么？Matlab=MatrixLaboratory=矩阵实验室由MathWorks公司于1984年推出的一套高性能的数值计算和可视化数学软件基本元素是矩阵数学表达式与运算规则符合通常的习惯语法不严格具有强大的数值计算和图示能力具有丰富的工具箱(ToolBox)强大的help功能Matlab是什么？Matlab=MatrixLaboMatlab窗口环境脚本M文件在命令窗口中输入并执行它所用的变量都要在工作空间中获取Matlab窗口环境脚本M文件Matlab窗口环境程序文件以.m格式进行存取程序运行时只需在工作空间中键入其名称即可Matlab窗口环境程序文件MATLAB语句形式变量＝表达式；通过等号将表达式的值赋予变量。当键入回车键时，该语句被执行。语句执行之后，窗口自动显示出语句执行的结果。如果不希望结果回显，只需在语句后加分号即可。Matlab窗口环境MATLAB语句形式Matlab窗口环境MATLAB语句形式变量＝表达式；通过等号将表达式的值赋予变量。当键入回车键时，该语句被执行。语句执行之后，窗口自动显示出语句执行的结果。如果不希望结果回显，只需在语句后加分号即可。Matlab窗口环境MATLAB语句形式Matlab窗口环境命令行编辑器方向键和控制键可以编辑修改已输入的命令↑：回调上一行命令↓：回调下一行命令多行命令（…):如果命令语句超过一行或者太长希望分行输入，可以使用多行命令继续输入

Matlab窗口环境例如：S=1-12+13+4+…9-4-18;命令行编辑器Matlab窗口环境例如：内容Matlab语言基础2基本运算符变量矩阵常用函数基本语句内容Matlab语言基础2基本运算符变量矩阵常用函数Matlab语言基础:变量变量的命名变量无需定义即可使用变量的名字必须以字母开头，之后可以是任意字母、数字或下划线(不能超过31个字符)变量名称区分字母的大小写变量可包含下连字符、数字，但不能包含空格、标点变量操作在命令窗口中，同时存储着输入的命令和创建的所有变量值，它们可以在任何需要的时候被调用采用clear将清除内存内的所有变量采用“clear变量名”清除指定变量Matlab语言基础:变量变量的命名变量操作一些特殊的变量ans：最近生成的无名变量i,j：虚数单位pi：圆周率

nargin：函数的输入变量个数eps：计算机的最小数nargout：函数的输出变量个数inf：无穷大(2/0)

realmin：最小正实数realmax：最大正实数nan：不定量(0/0)flops：浮点运算数Matlab语言基础:变量一些特殊的变量inf：无穷大(2/0)常用的数学运算符运算式中，通常不需要考虑空格多条命令可以放在一行中同行中各条命令可用分号隔开，表示禁止结果显示同行中各条命令也可用逗号隔开，表示需要显示结果＋----加－----减*----乘

/----左除

\----右除

^----幂Matlab语言基础:基本运算符常用的数学运算符＋----加常用的关系运算符

>----大于

>=----大于等于

<----小于

<=

----小于等于

==----等于

~=----不等于运算法则若关系式成立，结果为1若关系式不成立，结果为0Matlab语言基础:基本运算符常用的关系运算符>----大常用的逻辑运算符

&----与

|----或

!----非运算法则若逻辑真，结果为1若逻辑假，结果为0Matlab语言基础:基本运算符常用的逻辑运算符&----与其他常见符号Matlab语言基础:基本运算符

=

----变量赋值%----注释符

‘----共轭转置符

.’----转置符:----冒号运算符举例：

n:s:m产生n~m、步长为s的序列其他常见符号Matlab语言基础:基本运算符=循环语句for-end循环：循环次数事先确定时使用while-end循环：循环次数不能事先确定时使用Matlab语言基础:基本语句fori=n:s:m

语句体ends为步长可以为正数，负数或小数while表达式语句体end表达式为真，执行语句体表达式为假，终止该循环循环语句Matlab语言基础:基本语句fori=n:s:分支语句if-elseif-else-end：具有优先级switch-case-otherwise-end：无优先级Matlab语言基础:基本语句if表达式1

语句体1elseif表达式2

语句体2………..else

语句体elseendswitch变量

case变量值1语句体1case变量值2语句体2……otherwise

语句体otherwiseend分支语句Matlab语言基础:基本语句if表达式1sw常用的数学函数abs----求绝对值sin----正弦cos----余弦tan----正切asin----反正弦acos----反余弦atan----反正切sqrt----求平方根exp----eximag----求虚部real----求实部sign----求符号log----求自然对数log10----以10为底的对数conj----共轭复数Matlab语言基础:常用函数常用的数学函数abs----求绝对值eMatlab语言基础:矩阵矩阵的生成例如：创建一个3×3矩阵输入：a=[123;456;789]

或者a=[123456789]

屏幕输出：a=123456789矩阵的元素直接排列在方括号内行与行之间用分号隔开每行内的元素使用空格或逗号隔开大的矩阵可以用分行输入，回车键代表分号Matlab语言基础:矩阵矩阵的生成例如：创建一个3×3矩矩阵的生成举例：输入a=[1:2:10]屏幕输出：a=13579输入a=[10:-2:1]屏幕输出：a=108642用线性等间距生成向量矩阵start:step:end

start：起始值

step:步长(默认为1，可为负值)

end:终止值Matlab语言基础:矩阵矩阵的生成举例：用线性等间距生成向量矩阵sta一些常用的特殊矩阵eye(2,3)ans=100010eye(2)ans=1001V=[572];A=diag(V)A=500070002Matlab语言基础:矩阵单位矩阵：eye(m,n);eye(m)零矩阵：zeros(m,n);zeros(m)一矩阵：ones(m,n);ones(m)对角矩阵：对角元素向量V=[a1,a2,…,an]A=diag(V)随机矩阵：rand(m,n)产生一个m×n的均匀分布的随机矩阵zeros(2,3)ans=000000zeros(2)ans=0000ones(2,3)ans=111111ones(2)ans=1111rand(2,2)ans=0.41540.87440.01500.9708一些常用的特殊矩阵eye(2,3)eye(2)V=[57矩阵元素的访问通过确认矩阵下标，可以对矩阵的子块进行插入、提取和重排注意：MATLAB中矩阵下标从1开始，而不是从0开始Matlab语言基础:矩阵

A(m,n)：提取第m行，第n列元素A(:,n)：提取第n列元素A(m,:)：提取第m行元素A(m1:m2,n1:n2)：提取m1~m2行和n1~n2列的所有元素A(:)：得到一个长列矢量，其元素按矩阵的列进行排列A=[123456];A(2,1)ans=4A(2,:)ans=456A(:,2)ans=25A(1:2,2:3)ans=2356A(:)ans=142536矩阵元素的访问Matlab语言基础:矩阵A(矩阵的扩展与消除矩阵扩展：如果在原矩阵中一个不存在的地址位置上设定一个数（赋值），则该矩阵会自动扩展行列数，并在该位置上添加这个数，而且在其他没有指定的位置补零。消除子块：如果将矩阵的子块赋值为空矩阵[]，则相当于消除了相应的矩阵子块。Matlab语言基础:矩阵A=[123456];A(3,2)=7A=[123456070]A(2,:)=[]A=[123]矩阵的扩展与消除Matlab语言基础:矩阵A=[12矩阵的大小Matlab语言基础:矩阵[m,n]=size(A,x)：返回矩阵的行列数m与n当x=1，则只返回行数m当x=2，则只返回列数nlength(A)=max(size(A))：返回行数或列数的最大值A=[123456];size(A)ans=23size(A,1)ans=2size(A,2)ans=3矩阵的大小Matlab语言基础:矩阵[m,n]=size(矩阵的运算举例A=[1+i2-2i;-3+3i-4-4i]A=1+i2-2i-3+3i-4-4iB=A’B=1-i-3-3i2+2i-4+4iC=A.’C=1+i-3+3i2-2i-4-4iMatlab语言基础:矩阵转置：转置符号有两种形式

A’----共轭转置A.’----普通转置求逆：inv(A)

求行列式：det(A)矩阵的运算举例B=A’C=A.’Matlab语言基础:矩阵矩阵的四则运算与幂运算

.*----点乘

.\----点右除

./----点左除

.^----点幂点运算是两个维数相同矩阵对应元素间的运算Matlab语言基础:矩阵＋----矩阵加－----矩阵减*----矩阵乘

/----矩阵左除

\----矩阵右除

^----矩阵幂矩阵的四则运算与幂运算.*----内容Matlab绘图基础3内容Matlab绘图基础3Matlab绘图基础MATLAB提供了丰富的绘图功能helpgraph2d可得到所有画二维图形的命令helpgraph3d可得到所有画三维图形的命令Matlab绘图基础MATLAB提供了丰富的绘图功能helpplot(x,y,’属性’)或者stem()：绘图figure：选择图像

figure(1)；figure(2)；…；figure(n)打开不同的图形窗口，以便绘制不同的图形gridon：在所画出的图形坐标中加入栅格

gridoff：除去图形坐标中的栅格

holdon：当前图形保持不变，同时允许在这个坐标内绘制另外一个图形holdoff：使新图覆盖旧的图形Matlab绘图基础plot(x,y,’属性’)或者stem()：绘图Matl举例：

t=[0:pi/20:9*pi];

f=sin(t);

figure(1)

plot(t,f,'r:*')

gridonMatlab绘图基础举例：Matlab绘图基础设定对数坐标Matlab绘图基础semilogx：x轴为对数坐标，y轴为线性坐标semilogy：y轴为对数坐标，x轴为线性坐标设定轴的范围axis（[xminxmaxyminymax]）axis(‘equal’)：将x坐标轴和y坐标轴的单位刻度调整为一样设定对数坐标Matlab绘图基础semilogx：x轴为对数文字标示Matlab绘图基础text(x,y,’字符串’)在图形的指定坐标位置(x,y)处，标示单引号括起来的字符串gtext(‘字符串’)利用鼠标在图形的某一位置标示字符串title(‘字符串’)在所画图形的最上端显示说明该图形标题的字符串xlabel(‘字符串’)，ylabel(‘字符串’)设置x，