本帖最后由 zwj108 于 2012-11-27 18:22 编辑

n=0:3;
xb=impesq(0,0,3);  %用于产生xb(n)=δ(n)信号
hb=impesq(0,0,3)+2.5*impesq(1,0,3)+2.5*impesq(2,0,3)+impesq(3,0,3);  %用于产生hb(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3)
k=-200:200;
w=(pi/100)*k;
Xb=xb*(exp(-1i*pi/100)).^(n'*k);     %信号xb的傅氏变换
Hb=hb*(exp(-1i*pi/100)).^(n'*k);    %信号hb的傅氏变换
n=0:3;
subplot(3,2,1);
%做序列xb的时域和频域的图形
stem(n,xb);
  axis([-2 2 0 2]);
  xlabel('n');
  ylabel('xb(n)');
  title('xb(n)');
subplot(3,2,2);
plot(w/pi,abs(Xb));                     
  axis([-2 2 0 2]);
  xlabel('w/pi');
  ylabel('xb(|(jw)|');
  title('[xb(ejw)]');  
  subplot(3,2,3);
%做序列hb的时域和频域的图形
  stem(n,hb);
  axis([0 4 0 3]);
  xlabel('n');
  ylabel('Hb');
  title('Hb(n)');
subplot(3,2,4);
plot(w/pi,abs(Hb));
  axis([-2 2 0 8]);
  xlabel('w/pi');
  ylabel('Hb(|(jw)|');
  title('[Hb(ejw)]');
  n=0:6;
  y=conv(xb,hb);                        %序列xb,hb的卷积
  Y=y*(exp(-1i*pi/100)).^(n'*k);   %卷积的傅氏变换
  subplot(3,2,5);
%响应y的时域和频域的图形
  stem(n,y);                                %卷积作图
  axis([0 7 0 3]);
  xlabel('n');
  ylabel('y(n)');
  title('xb*hb');
subplot(3,2,6);
plot(w/pi,abs(Y));                       %傅氏变换作图
  axis([-2 2 0 8]);
  xlabel('w/pi');
  ylabel('|Y(jw)|');
  title('[Y(ejw)]');
??? Error using ==> conv2
First and second arguments must be
single or double.
Error in ==> conv at 39
c = conv2(a(:),b(:),shape);