数字信号处理实验课件最终Word文件下载.docx
《数字信号处理实验课件最终Word文件下载.docx》由会员分享，可在线阅读，更多相关《数字信号处理实验课件最终Word文件下载.docx（41页珍藏版）》请在冰豆网上搜索。
该函数库既包含了诸如求和、正弦、余弦、复数运算等简单函数；
也包含了矩阵转置、特征值、贝塞尔函数、快速傅立叶变换等复杂函数。
（5）Matlab应用程序接口（API）：
是一个同Fortran和C等高级语言进行交互的库。
2Simulink
Simulink是用于动态仿真的交互式系统。
允许用户在屏幕上绘制框图来模拟一个系统，并能够动态地控制该系统。
Simulink的模块库为用户提供了多种多样的功能模块用户可以通过简单的操作建立模型，然后选择合适的的输入源模块（如正弦波）作为输入信号，用适当的接收模块（如示波器）观察系统响应、分析系统特性，如仿真结果不符合要求，则可以修改系统模型的参数，继续进行仿真分析。
3Matlab工具箱
工具箱是开放式的，可以应用，也可以根据自己的需要进行扩展。
目前已有20多种工具箱投入使用，包括信号处理（signalprocessing）、图象处理（imageprocessing）、通信（communications）、控制系统（controlsystem）等等。
三命令窗口（commandwindow）的菜单系统
1File菜单
（1）New子菜单
有三个选项：
“M-file”新建一个M文件。
“figure”新建一个图形窗口。
“Model”新建一个Simulink模型窗口
（2）Open子菜单：
打开一个文件
（3）ShowGUILayoutTool:
显示图形用户界面设计向导。
（4）SetPath:
打开路径浏览器，用于设置路径。
2Edit菜单：
复制、粘贴、剪切、撤消等
3View菜单：
用于控制是否显示工具栏。
4Window菜单：
使用户在打开的窗口之间切换
5Help菜单
另外，它还包含演示程序：
demo和tour
四基本数值运算
1简单矩阵的生成
例：
A=[1,2,3,4;
5,6,7,8;
8,7,6,5]
或A=[1234
5678
8765]都生成同样的矩阵。
2常用矩阵的生成
（1）zeros生成全零矩阵
如B=zeros（n）生成n*n的全零矩阵
（2）ones生成全一阵
如C=ones（m,n）生成m*n的全１阵
（3）eye生成单位阵
如D=eye（size（A））生成一个与A大小相同的矩阵
3向量的生成
如x1=1:
5
x2=1:
0.5:
3
x3=5:
-1:
1
则可得x1=
12345
x2=
1.00001.50002.00002.50003.0000
x3=
54321
4加减运算
C=A+B或C=plus（A,B）
C=A-BC=minus（A,B）
5乘法运算
矩阵乘C=A*B或C=mtimes（A，B）
数组乘C=A.*BC=times（A,B）
矩阵乘是两个矩阵的线性代数乘，对于非标量A,B,矩阵A的列数必须等于矩阵B的行数。
标量可以和任意大小的矩阵乘。
数组乘是数组A和B的对应元素相乘。
除非一个是标量，否则A和B必须有相同的大小。
A=[121;
212;
121];
B=[111;
111;
111];
C=A*BD=A.*B
结果分别为
C=D=
444121
555212
6除法运算
（1）矩阵除法
左除C=A\B或C=mldivide（A,B）
右除C=A/BC=mrdivide（A,B）
通常矩阵左除不等于右除。
如A是方阵，左除等效于Ａ的逆阵左乘B矩阵，即inv（A）*B。
B/A大体等效于B*inv（A）（B右乘A的逆阵），精确地B/A=（A'
\B'
）'
A=[123;
010;
321];
B=[121]'
;
求C=A\BD=A/B
因为A为3*3矩阵，inv（A）也为3*3矩阵，所以A\B=inv（A）*BA的列数必须等于B的行数.
结果为
C=而D=A/B
-0.7500D=?
?
Errorusing==>
\
2.0000Matrixdimensionsmustagree.
–0.7500
B=[456;
C=D=
0.6250-0.3750-0.37500.57140-0.4286
01.0000-0.000001.00000
1.12501.12502.1250001.0000
（2）数组乘法
数组左除：
C=A.\B或C=mldivide（A,B）
数组右除：
C=A./BC=mrdivide（A,B）
数组左除是数组B和A中对应元素相除，即B（i,j）/A（i,j），
除非其中一个是标量，否则二者大小必须相同。
数组右除是数组A和B中对应元素相除，即A（i,j）/B（i,j），
除非其中一个是标量，否则二者大小必须相同。
例：
求C=A.\BD=A./B
结果为
c=d=
4.00002.50002.00000.25000.40000.5000
NaN1.0000NaNNaN1.0000NaN
1.00001.00001.00001.00001.00001.0000
注：
NaN不是一个数，当0/0或Inf/Inf时用它作为标记
Inf为表示无穷大的专用变量，当1/0时或浮点溢出用
它作为标记。
7乘方运算
（1）矩阵的乘方运算C=A^p或C=mpower（A,p）
（2）数组的乘方运算 C=A.^BC=power（A,B）
矩阵的乘方运算要求A是一个方阵，且p是一个标量，
A^p的意思是矩阵A的p次方，若p为正整数，则A^p是矩阵A自乘p次，若p为负整数，则首先对A求逆，然后将它自乘p次。
数组的乘方运算是矩阵A和B对应元素的乘方，即A（i,j）的B（i,j）次方。
A=[123;
010;
求A^2
ans=
10106
010
61010
B=eye（3）;
求C=A.^B D=A.^2E=2.^A
c=d=e=
111149248
111010121
111941842
8矩阵旋转
B=fliplr（A）矩阵A关于垂直轴沿左右方向进行列维翻转
B=flipud（A）矩阵A关于水平轴沿上下方向进行行维翻转
B=rot90（A）矩阵A逆时针方向旋转90度
B=rot90（A,k）矩阵A逆时针方向旋转k*90度
456];
B=fliplr（A）;
C=flipud（A）;
D=rot90（A,3）
B=C=D=
32145641
65412352
63
五基本平面做图
1plot函数
plot（x,y）x为横坐标向量，y为纵坐标向量
t=0:
0.1:
4*pi;
y=sin（t）;
plot（t,y）
结果见图1图1
2subplot函数用于在一个图形窗口显示几幅图
格式：
subplot（m,n,i）
功能：
把图形窗口分割为m行n列子窗口，然后选定第i个窗口为当前窗口。
figure
（1）;
clf
n=[-5:
5];
x=2*impseq（-2,-5,5）-impseq（4,-5,5）;
subplot（3,1,1）;
stem（n,x）;
title（'
impluse'
）;
xlabel（'
n'
ylabel（'
x（n）'
axis（[-5,5,-2,3]）;
n=[0:
20];
x1=n.*（stepseq（0,0,20）-stepseq（10,0,20））;
x2=10*exp（-0.3*（n-10））.*（stepseq（10,0,20）-stepseq（20,0,20））;
x=x1+x2;
subplot（3,1,2）;
'
axis（[0,20,-1,11]）;
50];
x=cos（0.04*pi*n）;
subplot（3,1,3）;
axis（[0,50,-1.2,1.2]）;
结果见图2
图2
实验一序列的产生及绘图
一、实验目的
1．熟悉信号处理软件MATLAB的使用。
2．离散信号的基本运算实现。
3．了解基本序列及复杂序列的产生方法。
4．运用卷积方法观察系统的时域特性。
升级会员 