subplot(2,2,2),plot(t,abs(f),title(模)subplot(2,2,4),plot(t,angle(f),title(相角)例矩形脉冲信号矩形脉冲信号可用rectpuls函数产生,调用格式为y=rectpuls(t,width),幅度是1,宽度是width,以t=0为对称中心。{=-2:0.01:2;width=1;ff=-2*rectpuls(t,width);plot(t,ft)grid on;,例单位阶跃信号单位阶跃信号u(t)用"t>=0"产生,调用格式为ft=(t>=0)t=-1:0.01:5;ft=(t>=0);plot(tft),grid on,axis([-1,5,-0.5,1.5]),例正弦信号符号算法%定义符号变量tsymsty=sin(pi/4*t)%符号函数表达式%符号函数画图命令ezplot(y,[-16,16])或者%定义符号函数表达式f=sym('sin(pi/4*t))ezplot(f,[-16,16])例单位阶跃信号MATTLAB符号数学函数Heaviside表示阶跃信号,但要画图需在工作目录创建Heaviside的M文件Matlab中自带heaviside函数function f=Heaviside(t)f-(t>0);保存,文件名是Heaviside,调用该函数即可画图,例t=-1:0.01:3;f-heaviside(t);plot(t,f)axis([-1,3,-0.2,1.2])或者y=sym(Heaviside(t);ezplot(y,[-1,5]);grid on例正弦序列信号正弦序列信号可直接调用MATLAB函数cos,例cos(のk+β),当2元/の是整数或分数时,才是周期信号。画cos(k元/8+?),cos(2k)波形程序是:
subplot(2,2,2),plot(t,abs(f)),title('模') subplot(2,2,4),plot(t,angle(f)),title('相角') 例 矩形脉冲信号 矩形脉冲信号可用 rectpuls 函数产生, 调用格式为 y=rectpuls(t,width),幅度是 1,宽度是 width,以 t=0 为对称中心。 t=-2:0.01:2; width=1; ft=2*rectpuls(t,width); plot(t,ft) grid on; 例 单位阶跃信号 单位阶跃信号 u(t)用“t>=0”产生,调用格式为 ft=(t>=0) t=-1:0.01:5; ft=(t>=0); plot(t,ft); grid on; axis([-1,5,-0.5,1.5]); 例 正弦信号符号算法 syms t %定义符号变量 t y=sin(pi/4*t) %符号函数表达式 ezplot(y,[-16,16]) %符号函数画图命令 或者 f=sym('sin(pi/4*t)') %定义符号函数表达式 ezplot(f,[-16,16]) 例 单位阶跃信号 MATTLAB 符号数学函数 Heaviside 表示阶跃信号,但要画图需在工作 目录创建 Heaviside 的 M 文件 Matlab 中自带 heaviside 函数 function f=Heaviside(t) f=(t>0); 保存,文件名是 Heaviside ,调用该函数即可画图,例 t=-1:0.01:3; f=heaviside(t); plot(t,f) axis([-1,3,-0.2,1.2]) 或者 y=sym('Heaviside(t)'); ezplot(y,[-1,5]);grid on 例 正弦序列信号 正弦序列信号可直接调用 MATLAB 函数 cos,例 cos(k ) ,当 2 / 是整数或分数时,才是周期信号。画cos(k / 8 ) ,cos(2k)波形程序是:
k=0:40;subplot(2,1,1)stem(k,cos(k*pi/8),filled')title(cos(k*pi/8))subplot(2,1,2)stem(k,cos(2*k),filled)title('cos(2*k))[1 k=0例单位序列信号(k)=10k+0本题先建立一个画单位序列(k+k.)的M函数文件,画图时调用。M文件建立方法file/new/m-file在文件编辑窗输入程序,保存文件名用函数名。%klk2是画图时间范围,ko是脉冲位置functiondwxulie(k1,k2,ko)k=kl:k2;n=length(k),f=zeros(1,n);f(1,-k0-k1+1)=1;stem(k,f,'filled')axis([k1,k2,0,1.5])title(单位序列8(k))保存文件名dwxulie.m画图时在命令窗口调用,例:dwxulie(-5,5,0)(1 k≥0例单位阶跃序列信号k)=10k<0本题也可先建立一个画单位阶跃序列ε(k+k。)的M函数文件,画图时调用。functionjyxulie(k1,k2,ko)k=k1:-k0-1;kk=-ko:k2;n=length(k);nn=length(kk) ;u=zeros(1,n);uu=ones(1,nn);stem(kk,uu,filled')hold onstem(k,u,'filled')hold offtitle(单位阶跃序列)axis([k1 k2 0 1.5])保存文件名jyxulie.m
k=0:40; subplot(2,1,1) stem(k,cos(k*pi/8),'filled') title('cos(k*pi/8)') subplot(2,1,2) stem(k,cos(2*k),'filled') title('cos(2*k)') 例 单位序列信号 0 0 1 0 ( ) k k k 本题先建立一个画单位序列 ( ) 0 k k 的 M 函数文件,画图时调用。M 文件建立方法: file / new / m-file 在文件编辑窗输入程序,保存文件名用函数名。 function dwxulie(k1,k2,k0) % k1 , k2 是画图时间范围,k0 是脉冲位置 k=k1:k2; n=length(k); f=zeros(1,n); f(1,-k0-k1+1)=1; stem(k,f,'filled') axis([k1,k2,0,1.5]) title('单位序列 δ(k)') 保存文件名 dwxulie.m 画图时在命令窗口调用,例:dwxulie(-5,5,0) 例 单位阶跃序列信号 0 0 1 0 ( ) k k k 本题也可先建立一个画单位阶跃序列 ( ) 0 k k 的 M 函数文件,画图时调用。 function jyxulie(k1,k2,k0) k=k1:-k0-1; kk=-k0:k2; n=length(k); nn=length(kk) ; u=zeros(1,n); uu=ones(1,nn); stem(kk,uu,'filled') hold on stem(k,u,'filled') hold off title('单位阶跃序列') axis([k1 k2 0 1.5]) 保存文件名 jyxulie.m