信号与系统matlab实验与答案

频率越高，图像更加密集。

用 MATLAB 产生音阶信号 1 2 3 4 5 6 7 1，并播放，抽样频率可设为 8000 Hz。

利用 MATLAB 产生信号 x1 t

cost 和 x1 t cos20t ，请画出信号 x1 t +x2 t

和信号 x1 t x2 t 的波形图。

all ;clear;clc;close

f1=262;fs=8000;f2=293;f3=329;f4=349;f5=392;f6=440;f7=

493;f8=523; N_zeros=300;

n=0:(fs/2-N_zeros-1);

x1=[sin(n*2*pi*f1/fs) zeros(1,N_zeros)]; x2=[sin(n*2*pi*f2/fs) zeros(1,N_zeros)]; x3=[sin(n*2*pi*f3/fs) zeros(1,N_zeros)]; x4=[sin(n*2*pi*f4/fs) zeros(1,N_zeros)]; x5=[sin(n*2*pi*f5/fs) zeros(1,N_zeros)]; x6=[sin(n*2*pi*f6/fs) zeros(1,N_zeros)]; x7=[sin(n*2*pi*f7/fs) zeros(1,N_zeros)]; x8=[sin(n*2*pi*f8/fs) zeros(1,N_zeros)]; notes=[x1 x2 x3 x4 x5 x6 x7 x8 ]; sound(notes,fs);

习题：计算下列信号 x n 和 h n 的卷积和。

1. x n h n u n u n 4 ;

nx=0:9;x=ones(1,length(nx)); nh=0:4;h=ones(1,length(nh)); y=conv(x,h);

% 下限 =下限 1+下限 2 ny_min=min(nx)+min(nh); % 上限 =上限 1+上限 2 ny_max=max(nx)+max(nh); ny=ny_min:ny_max;

subplot(3,1,1);stem(nx,x);

xlabel('n');ylabel('x(n)');axis([ny_min ny_max 0 max(x)]); subplot(3,1,2);stem(nh,h);

xlabel('n');ylabel('h(n)');axis([ny_min ny_max 0 max(h)]); subplot(3,1,3);stem(ny,y);

xlabel('n');ylabel('x(n)*h(n)');axis([ny_min ny_max 0 max(y)]);

2. x n

u n u n 4 , h n

0.8n u n

u n 10

nx=0:3;x=ones(1,length(nx));

nh=0:9;h=0.8.^nh.*ones(1,length(nh));

y=conv(x,h);ny_min=min(nx)+min(nh);ny_max=max(nx)+max (nh);

ny=ny_min:ny_max; subplot(3,1,1);stem(nx,x); xlabel( 'n' );ylabel( 'x(n)' );axis([ny_min ny_max 0 max(x)]); subplot(3,1,2);stem(nh,h); xlabel( 'n' );ylabel( 'h(n)' );axis([ny_min ny_max 0 max(h)]); subplot(3,1,3);stem(ny,y); xlabel( 'n' );ylabel( 'x(n)*h(n)' );axis([ny_min ny_max 0 max(y)]);