lab1 信号与系统

1、Lab 1: Introduction to MATLAB目录n1. Warm-upn2. Arithmetic Operationsn3. Complex numbersn4. Array indexingn5. Allocating memoryn6. Special characters and functionsn7. Control flown8. Plotting目录n9. Programming in MATLAB (M-files)n10. MATLAB soundn11. Loading and saving data1.Warm-up:n启动Matlabl双击桌面快捷方式启

2、动 文件保存路径为:MATLABR2012aworkl选择路径MATLABR2012abin中的matlab.exe启动 文件保存路径为:MATLABR2012abin操作页面简介标题栏菜单栏常用工具栏命令窗口历史命令窗口工作空间当前目录工作空间的保存工作空间的保存操作页面简介编辑窗口图形窗口nintro Matlab简介nhelp帮助功能qhelp plot画图qhelp ops运算符qhelp arith四则运算n计算功能试着计算： a=1.2, b=2.3, c=4.5 and d=4.注意每行后面的分号！注意每行后面的分号！命令窗口 Command Window34abdc2.算数运算

3、 Arithmetic Operationsn常用运算符号：+ - * /n其他运算符号：（用于向量）q.* 向量间的乘法q./向量间的除法q.向量内所有元素的幂n例如：1212x ,.,andy,.,nnx xxy yy算数运算 Arithmetic Operationsn对比说明：3.复数 Complex numbersn z=3 + 4i % note that you do not need the * after 4n conj(z) % computes the conjugate of z共轭n angle(z) % computes the phase of z相位n real

4、(z) % computes the real part of z实部n imag(z) % computes the imaginary part of z虚部n abs(z) % computes the magnitude of z绝对值n也可以用其他变量来定义一个虚数：也可以用其他变量来定义一个虚数： img=sqrt(-1) z=3+4*img exp(pi*img)4.数组 Array indexingn创造一个由整数1到6的数组: （步长为1） x=1,2,3,4,5,6 or x=1:6n创造一个由整数1到6的数组: （步长为2） x=1:2:6Ans = 1 3 5n提取或插

5、入一个向量:例如： x=1:3 4 6 100:110 x(3:7) length(x) % gives the size of the array or vector x(2:2:length(x)5.内存分配 Allocating memoryn可用zeros为一维向量分配内存： Y=zeros(100,1); 为1001维向量分配内存 Y(30)ans = 0n可用ones定义一个全1矩阵： Y=ones(1,5)ans= 1 1 1 1 16.Matlab中的预定义变量常见的命令和功能nclf -清除图形窗nclc -清除命令窗口中的显示内容nclear - 清除Matlab工作空间中

6、保存的变量nexit -退出Matlabnlength - 给出数组x的长度nfind -返回数组x中非零的元素 例： x=1:10; find(x 4)ans = 5 6 7 8 9 107.控制命令 Control flowqif statementsqswitch statementsqfor loops都以都以 end 结束结束qwhile loopsqbreak statements下面看几个例子下面看几个例子Examples关系和逻辑运算nRelational OperatorsSymbol Meaning=Less than equal= Greater than equal G

7、reater than=Equal =Not equalnLogical OperatorsSymbol Meaning& AND| OR NOT8.画图 Plottingnplot(.) 二维图形的基本函数nstem(.) 杆图nbar(.) 条形图Example: x=1:20; plot(x) %see Figure 1 stem(x) % see Figure 28.画图 Plottingnplot(X,Y,S) S用来表示线的颜色、类型等nsubplot(m,n,p) 图形窗口的分割m:行数n:列数 p:在第p个小窗口中建立坐标轴例： subplot(2,1,1), plot

8、(x) subplot(2,1,2), stem(x)图形编辑函数函数作用作用函数函数作用作用hold on/off控制是保持图形还是刷新图形axis(xmin xmax ymin ymax zmin zmax)分别设置X、Y、Z轴的取值范围title(图形名称)标注图形名称axis equal横、纵坐标轴采用等长刻度xlabel(X轴说明)X轴坐标说明axis on/off显示/取消坐标轴xlabel(Y轴说明)Y轴坐标说明grid on/off控制画/不画网格线text (x,y,图形说明) 在指定的位置添加图形说明box on/off控制加/不加边框线9. Programming in

9、MATLAB (M-files)n用M-file编程File New M-fileq例如：写一个计算x向量平均值的程序n1.建立一个新文件，命名“average.m”，输入以下代码function y=average(x)L=length(x);sum=0;for i=1:Lsum=sum+x(i);endy=sum/L; % the average of xn2.保存文件，FileSave，输入文件名“average.m”n3.在命令窗口执行程序，输入以下代码： x=1:100; y=average(x)ans = 50.500010. 声音 MATLAB soundnsoundsc 播放音频Soundsc(Y,FS)Y- 信号FS-频率q听一下400-Hz的声音 t=0:1/8192:1; x=cos(2*pi*400*t); soundsc(x,8192);q听一下噪音 noise=randn(8192,1); % generate 8192 samples of noise soundsc(noise,8192);11. 载入和保存数据Loading and savin