数字信号出来实验指导.doc

数字信号出来实验指导 实验一 用于信号处理的MATLAB基本操作 = 1 \* Arabic 1．实验目的： 学习和掌握MATLAB最基本的矩阵运算与绘图工具。 学习用MATLAB实现信号处理的基本操作。 = 2 \* Arabic 2．实验步骤： 熟悉MATLAB的使用环境和方法。 练习使用基本的向量生成、矩阵运算、绘图等语句。 利用冒号（：）生成向量 X1=[1 2 3 4 5] X2=[1.000 1.500 2.000 2.500] X3=[5 4 3 2 1] 分别生成3*3，3*4的全0矩阵，全1矩阵和随机矩阵 分别输入矩阵 1 2 3 1.0 1.1 1.2 A= 4 5 6 B= 2.0 2.1 2.2 7 8 9 3.0 3.1 3.2 分别计算A+B，A-B，A+3，A-4，A*B，A.*B，C=inv（A）,A/B,A./B 分别计算sin(x1),cos(x1),exp(x1),lg(x2),sqrt(x2) = 3 \* Arabic 3.实验说明 向量的生成 利用冒号（：）生成向量 X=j:k (jk) 生成X=[j,j+1,j+2,…,k-1,k] X=j:I:k 如果I0 且jk或I0且jk,则生成向量 X=[j,j+I,j+2I,…,k-I,k] 如果I0 且jk,则生成向量 X=[j,j+I,j+2I,…,k] 例 X1=1：5 X1=[1 2 3 4 5] X2=1：0.5：3 X2=[1.000 1.500 2.000 2.500 3.000] X3=5：-1：1 X3=[5 4 3 2 1] 矩阵的生成 Zeros 生成全0阵 B= Zeros(n) 生成n*n的全0矩阵 B= Zeros(m,n) 生成m*n的全0矩阵 B= Zeros(Size(A)) 生成与与矩阵A大小相同的全0矩阵 注：m,n必须是一个非负数 ones生成全1阵 rand生成均匀分布的随机阵 二 矩阵的算术运算 加法和减法 对于同维矩阵指令的 A+B B 对于矩阵和标量（一个数）的加减运算，指令为： A+3 9 乘法和除法运算 A*B 是数学中的矩阵乘法，遵循矩阵乘法规则 A．*B 是同维矩阵对应位置元素做乘法 B=inv（A）是求矩阵的逆 A/B 是数学中的矩阵除法，遵循矩阵除法规则 A．/B 是同维矩阵对应位置元素相除 另 表示矩阵的转置运算 三．数组函数 下面列举一些基本函数，他们的用法和格式都相同。 Sin(A),cos(A),exp(A),log(A)(相当于ln) Sgrt(A)开平方 abs(A)求模 real(A)求实部 imag(A)求虚部 式中A可以是标量也可以是矩阵 例： 利用等差向量产生一个正弦值向量 t=0:0.1:10 A=sin(t) PLot(A) 这时候即可看到一个绘有正弦曲线的窗口弹出 另：每条语句后面加“；”表示不要显示当前语句的执行结果 不加“；”表示要显示当前语句的执行结果。 练习：t=0:0.01:1; y=sin(2*pi*t); plot(t,y); xlabel(0-1); ylabel(-1-1); text(0.6,0,y=sin(2*pi*t)); H=legend(sin); 四．绘图（二维绘图） 利用前例 如果要使向量的横纵坐标一一对应，则应写为plot(t,A) 如不要横坐标对应时，可只写plot(A) 如果A是一个矩阵不是向量，则在同一窗口中可绘出与矩阵行数相同的曲线且颜色不同。 同时也可指定曲线的颜色和格式。 如plot(t,A,’b.’) 则原来的连续曲线就变成了蓝色的点线。 下面给出常用选项。 r红 g绿 b蓝 y 黄 m 洋红 c青 w白 k黑 ━实线 ┅┅虚线 *用星号标数据点 。用点标 o（字母‘o‘）用圆圈标 x（字母’x‘）用叉标 -.点划线 下面给出窗口分割语句 12 1 2 Y Y+1====+1+1 N Y+Y 将一个窗口分割成x行，y列 共x*y个窗口