信号与系统实验报告

信号与系统实验报告

第三次课：四、五章上机实验

1．连续时间傅里叶变换性质：4.3节(b)。

clc

clear

load splat

y=y(1:8192);

N=8192

fs=8192;

Y=fftshift(fft(y));

Y1=conj(Y);

y1=real(Y1);

sound(y1,fs);

2．求由微分方程描述的单位冲激响应：4.5节(b)。

clc

clear

A1=[1 1.5 0.5];

B1=[1 -2];

[R1 P1]=residue(B1,A1)

运行结果：

R1 =

6

-5

P1 =

-1.0000

-0.5000

3．计算离散时间傅里叶变换：5.1节(a),(b),(c)。ａ．

clc

clear

Xn=[ones(1,10),zeros(1,6)];

N=0:15;

stem(N,Xn)

运行结果：

ｂ．

clc

clear

Xn=[ones(1,10),zeros(1,6)];

X=fft(Xn,100);

N=100;

K=0:N-1;

W=2*pi*K/N;

subplot(2,1,1)

stem(K,abs(X))

title('|X(ejw)|')

subplot(2,1,2)

stem(K,angle(X))

title('

运行结果：

ｃ．

clc

clear

N=100;

K=0:N-1;

Xn=[ones(1,10),zeros(1,6)]; X=fft(Xn,100);

W=2*pi*K/N-pi;

subplot(2,1,1)

stem(W,abs(fftshift(X))); title('|X(ejw)|');

subplot(2,1,2)

title('

stem(W,angle(fftshift(X)));

运行结果：

４．吉布斯现象：根据英文教材Example 3.5验证Fig3.9的吉布斯现象(a)~(d)

clc

clear

x=1/2;

k=0;

t=-2:0.01:2;

x1=[-2 -1 -1 1 1 2];y1=[0 0 1 1 0 0];

for i=1:2:49

k=i-k;x=x+(-1)^(k+1)*2/(i*pi)*cos(i*pi*t/2);

if(i==1)

figure;

plot(x1,y1,t,x);title('N=1')

elseif(i==3)

figure;

plot(x1,y1,t,x);title('N=3')

elseif(i==7)

figure;

plot(x1,y1,t,x);title('N=7') elseif(i==9)

figure;

plot(x1,y1,t,x);title('N=9') elseif(i==49)

figure; plot(x1,y1,t,x);title('N=49') end

end

运行结果：

第四次课：七、九、十章上机实验

１．由欠采样引起的混叠：7.1节(a),(b),(c),(d)。ａ．

clc

clear

W=2000*pi;

T=1/8192;

N=[0:8191];

Tn=N*T;

X=sin(W*Tn);

stem(Tn,X)

运行结果：

ｂ．

clc

clear

W=2000*pi;

T=1/8192;

N=0:50

Tn=N*T;

X1=sin(W*Tn); X2=sin(W*N*T) subplot(2,1,1) plot(Tn,X1) subplot(2,1,2) stem(N,X2)

运行结果：

ｃ．

clc

clear

W=2000*pi;

T=1/8192;

N=[0:8191];

X=sin(W*N*T);

N=length(X);

X=fftshift(fft(X,N))*(2*pi/N); W=linspace(-1,1-1/N,N)/(2*T); stem(W,abs(X));

运行结果：

ｄ．

clc

clear

W=3000*pi;

W=4000*pi;

T=1/8192;

N=[0:8191];

Tn=N*T;

X=sin(W*Tn); subplot(2,2,1) stem(Tn,X)

N=0:50

Tn=N*T;

X1=sin(W*Tn); X2=sin(W*N*T) subplot(2,2,2) plot(Tn,X1) subplot(2,2,3) stem(N,X2)

T=1/8192;

N=[0:8191];