实验十五信号的可视化和信号的时域运算.ppt

Matlab仿真单元 * 第二部分 上机实验 实验目的： (1)了解Matlab的数值计算、符号运算、图形控制等基本功能及使用方法； (2)掌握Matlab的程序设计方法，熟练掌握Matlab用于信号与系统分析的专用函数及使用方法。 Matlab的文件有两类：脚本文件和函数文件。 1．脚本文件 在Matlab环境下直接输入的语句，如果放在一个以 .m 为后缀的文件中，就称为脚本文件。有了脚本文件，可直接在命令窗口中输入脚本文件名以执行脚本文件中的每一条语句，这与在Matlab命令窗口中直接逐条输入语句的结果完全一致。 2．函数文件 Matlab的函数文件至少由两部分构成： ·函数定义行； 　·函数体。 另外在函数的任何位置都可以插入注释行，注释行以%开头。 在command windows下，输入： t=1:10回车后，显示t=1,2,……t=10 如换成：t=1:10; 后，则不显示 (多分号) 若换成：t=1:0.01:10; (0.01为步长) 如画余弦波形： t=1:0.01:10; plot(t,cos(t)); %波形感觉平滑 若间隔较大， t=1:0.1:10; plot(t,cos(t)); %波形失真 以上命令不能保存 要保存，需用m文件。 （1）脚本文件 如：test1.m N=1:10; Stem(a,sin(n)); 此m文件称为脚本文件，无输入输出参数，像C语言中的main。 输入方式： 选file—new—mfile后，即可。 （2）函数文件 如例1： function f=u(t); f=(t0); 定义了u(t)函数（阶跃函数），存为u.m至work目录中。 其中：function为关键字；f为函数；u为函数名；t 为参数； f=(t0); 为功能，当t0为真，返回“1”；当t0为假，返回“0”； 以u.m文件存至work目录中，退出。以后就可用。 例2：求卷积和的函数文件 function [f, n]=dconv(f1, f2, n1, n2) %f1：序列f1(n)的非零样值向量； %f2：序列f2(n)的非零样值向量； %n1：序列f1(n)的对应时间序号向量； %n2：序列f2(n)的对应时间序号向量； %f：卷积和序列f (n)= f1(n)* f2(n)的非零样值向量； %n：卷积和序列f (n)的对应时间序号向量； f=conv(f1, f2); N0=n1(1)+n2(1); L=length(n1)+length(n2)-1; n=N0: N0+L-1; 以上是，卷积和序列f (n) = f1(n)* f2(n)的M函数文件。 我们看到，当函数具有多个输出变量时，则以方括号括起，多个输入变量和多个输出变量之间均以逗号分割。 实验十五 信号的可视化 和信号的时域运算 一、实验目的 1.了解Matlab的基本功能及使用方法； 2.掌握用Matlab画信号波形的基本方法； 3.借助的Matlab符号运算功能实现信号的反褶、移位、尺度变换等运算。 二、实验原理 例1 设f1(t)=sin(Ωt)、 f2(t)=sin(8Ωt) ，画出f3(t)= f1(t)+f2(t)和f4(t)= f1(t).f2(t)的时间波形。 t=-5:0.01:5; %在-5~+5间取值，间隔0.01 f1=sin(pi/2*t); f2=sin(8*pi/2*t); %令Ω=π/2 f3=f1+f2; f4=f1.*f2; %点乘表示逐点相乘 % 直接相乘表示矩阵相乘 subplot(4,1,1), plot(t, f1) ；%subplot将图形进行空间分割; plot为连续函数画图命令，plot（自变量，函数） axis([-5 5 -1.2 1.2]), grid on； % axis([横轴范围 纵轴取值]), grid on为打开网格； set(gca, ‘xtick’, [-4 -2 0 2 4]), set(gca, ‘ytick’, [-1 0 1]) %在水平或垂直方向上，画虚线； subplot(4,1,2), plot(t, f2)； axis([-5 5 -1.2 1.2]), grid on