二、卷积的实现

1．（1）编写卷积积分程序，实现两个矩形脉冲信号f1(t)、f2(t)的卷积运算，其中f1(t)?g10(t?5)，f2(t)?g24(t?12)并调试运行出结果。（必做）

y1=[ones(1,10),zeros(1,20)]; y2=[ones(1,24),zeros(1,20)]; y=conv(y1,y2); n1=1:length(y1); n2=1:length(y2); L=length(y);

subplot(3,1,1);plot(n1,y1);axis([1,L,0,2]);

subplot(3,1,2);plot(n2,y2);axis([1,L,0,2]); n=1:L;

subplot(3,1,3);plot(n,y);axis([1,L,0,15]);

（2）改变矩形脉冲信号的参数，如f1(t)?g10(t?5)，f2(t)?g20(t?5)重新计算卷积结果。（必做）

自定义函数：

function [f,k]=sconv(f1,f2,k1,k2,p) f=conv(f1,f2); f=f*p;

k0=k1(1)+k2(1);

k3=length(f1)+length(f2)-2; k=k0:p:k0+k3*p; subplot(2,2,1) plot(k1,f1); title('f1(t)'); xlabel('t'); ylabel('f1(t)'); subplot(2,2,2); plot(k2,f2); title('f2(t)') subplot(2,2,3) plot(k,f);

h=get(gca,'position'); h(3)=2.5*h(3); set(gca,'position',h); title('f(t)=f1(t)*f2(t)')

xlabel('t'); ylabel('f(t)')

程序： clear p=0.01;

k1=0:p:10;

f1=ones(1,length(k1)); k2=-5:p:15;

f2=ones(1,length(k2));

[f,k]=sconv(f1,f2,k1,k2,p)

2．编写卷积积分程序，实现矩形脉冲信号f1(t)?g10(t?5)，与f2(t)?e?(t)的卷积运算，并调试运行出结果。（必做）（要求坐标必须是正确的，f2范围取0-5。）

clear p=0.01; k1=0:p:10;

f1=ones(1,length(k1)); k2=0:p:10; f2=exp(-k2);

[f,k]=sconv(f1,f2,k1,k2,p);

?t

3．尝试不调用库函数conv，实现上题中（1）信号的卷积。（选做）

clear p=0.01; k1=0:p:10;

f1=ones(1,length(k1)); k2=0:p:24;

f2=ones(1,length(k2));

[f,k]=sconv(f1,f2,k1,k2,p)

三、连续信号的频域特性分析

1．求门函数g4(t)的傅立叶变换，并画出其频率特性曲线图。（必做） 自定义函数：

function y=sf1(t,w);

y=(t>=-2 & t<=2).*1.*exp(-j*w*t);

程序：

w=linspace(-6*pi,6*pi,512); N=length(w); F=zeros(1,N); for k=1:N

F(k)=quadl('sf3',-2,2,[],[],w(k)); end figure(1); plot(w,real(F)); xlabel('\\omega'); ylabel('F(j\\omega)');

title('门函数g4(t)频率特性曲线');