数字信号的处理实验问题解释Word文件下载.docx

数字信号的处理实验问题解释Word文件下载.docx

《数字信号的处理实验问题解释Word文件下载.docx》由会员分享，可在线阅读，更多相关《数字信号的处理实验问题解释Word文件下载.docx（60页珍藏版）》请在冰豆网上搜索。

subplot（4,2,1）;

stem（n,a）;

xlabel（'

n'

）;

xlim（[05]）;

ylabel（'

A'

subplot（4,2,2）;

stem（n,b）;

B'

subplot（4,2,3）;

stem（n,c）;

C'

subplot（4,2,4）;

stem（n,d）;

D'

subplot（4,2,5）;

stem（n,e）;

E'

subplot（4,2,6）;

stem（n,f）;

F'

subplot（4,2,7）;

stem（n,g）;

G'

（2）用MATLAB实现下列序列：

a）x（n）=0.8n0≤n≤15

b）x（n）=e（0.2+3j）n0≤n≤15

c）x（n）=3cos（0.125πn+0.2π）+2sin（0.25πn+0.1π）0≤n≤15

d）将c）中的x（n）扩展为以16为周期的函数x16（n）=x（n+16）,绘出四个周期。

e）将c）中的x（n）扩展为以10为周期的函数x10（n）=x（n+10）,绘出四个周期。

N=0:

15;

%a）x（n）=0.8n0≤n≤15

xa=0.8.^N;

figure;

subplot（2,1,1）;

stem（N,xa）;

xlabel（'

xlim（[016]）;

xa'

%b）x（n）=e（0.2+3j）n0≤n≤15

xb=exp（（0.2+3*j）*N）;

subplot（2,1,2）;

stem（N,xb）;

xb'

%c）x（n）=3cos（0.125πn+0.2π）+2sin（0.25πn+0.1π）0≤n≤15

xc=3*cos（0.125*pi*N+0.2*pi）+2*sin（0.25*pi*N+0.1*pi）;

subplot（3,1,1）;

stem（N,xc）;

xc'

%d）将c）中的x（n）扩展为以16为周期的函数x16（n）=x（n+16）,绘出四个周期。

k=0:

3;

m=0;

fori=1:

4

forj=1:

16

m=m+1;

n（m）=N（j）+16*k（i）;

x16（m）=3*cos（0.125*pi*n（m）+0.2*pi）+2*sin（0.25*pi*n（m）+0.1*pi）;

end

end

subplot（3,1,2）;

stem（n,x16）;

x16'

%e）将c）中的x（n）扩展为以10为周期的函数x10（n）=x（n+10）,绘出四个周期。

10

x10（j）=x16（j）;

fori=1:

3

form=1:

x10（i*10+m）=x10（m）;

end

40;

subplot（3,1,3）;

stem（n,x10）;

x10'

（3）x（n）=[1,-1,3,5],产生并绘出下列序列的样本:

a）x1（n）=2x（n+2）-x（n-1）-2x（n）

b）

clearall

T=4;

x=[1-135];

x（5:

8）=x（1:

4）;

stem（1:

8,x）;

grid;

ifi-1<

x1（i）=2*x（i+2）-x（i-1）-2*x（i）;

else

x1（i）=2*x（i+2）-x（i-1+T）-2*x（i）;

end

x1（5:

8）=x1（1:

8,x1）;

（4）绘出下列时间函数的图形，对x轴、y轴以及图形上方均须加上适当的标注：

a）x（t）=sin（2πt）0≤t≤10s

b）x（t）=cos（100πt）sin（πt）0≤t≤4s

ta=0:

0.05:

10;

xa=sin（2*pi*ta）;

plot（ta,xa）;

t'

幅度'

tb=0:

0.01:

xb=cos（100*pi*tb）.*sin（pi*tb）;

plot（tb,xb）;

（5）编写函数stepshift（n0,n1,n2）实现u（n-n0）,n1<

n0<

n2,绘出该函数的图形，起点为n1，终点为n2。

n0=5;

ns=1;

nf=10;

%ns为起点;

nf为终点；

在=n=n0处生成单位阶跃序列

n=[ns:

nf];

x=[（n-n0）>

=0];

stem（n,x）;

（6）给一定因果系统

求出并绘制H（z）的幅频响应与相频响应。

clearall;

b=[1,sqrt

（2）,1];

a=[1,-0.67,0.9];

[h,w]=freqz（b,a）;

am=20*log10（abs（h））;

plot（w,am）;

ph=angle（h）;

plot（w,ph）;

（7）计算序列{8-2-123}和序列{23-1-3}的离散卷积，并作图表示卷积结果。

a=[8-2-123];

b=[23-1-3];

c=conv（a,b）;

%计算卷积

M=length（c）-1;

n=0:

1:

M;

（8）求以下差分方程所描述系统的单位脉冲响应h（n）,0≤n≤50

y（n）+0.1y（n-1）-0.06y（n-2）=x（n）-2x（n-1）

N=50;

a=[1-2];

b=[10.1-0.06];

x=[1zeros（1,N-1）];

N-1;

y=filter（a,b,x）;

stem（k,y）;

幅度'

实验二信号的采样与重建

一，实验目的

（1）通过观察采样信号的混叠现象，进一步理解奈奎斯特采样频率的意义。

（2）通过实验，了解数字信号采样转换过程中的频率特征。

（3）对实际的音频文件作插和抽取操作，体会低通滤波器在插和抽取中的作用。

二，实验容

（1）采样混叠，对一个模拟信号Va（t）进行等间采样，采样频率为200HZ，得到离散时间信号V（n）.Va（t）由频率为30Hz,150Hz,170Hz,250Hz,330Hz的5个正弦信号的加权和构成。

Va（t）=6cos（60pi*t）+3sin（300pi*t）+2cos（340pi*t）+4cos（500pi*t）+10sin（660pi*t）观察采样后信号的混叠效应。

程序：

clear,

closeall,

t=0:

0.1:

20;

Ts=1/2;

Ts:

V=8*cos（0.3*pi*t）+5*cos（0.5*pi*t+0.6435）-10*sin（0.7*pi*t）;

Vn=8*cos（0.3*pi*n）+5*cos（0.5*pi*n+0.6435）-10*sin（0.7*pi*n）;

subplot（221）

plot（t,V）,

gridon,

subplot（222）

stem（n,Vn,'

.'

）,

（2）输入信号X（n）为归一化频率f1=0.043，f2=0.31的两个正弦信号相加而成，N=100，按因子M=2作抽取：

（1）不适用低通滤波器；

（2）使用低通滤波器。

分别显示输入输出序列在时域和频域中的特性。

clear;

N=100;

M=2;

f1=0.043;

f2=0.31;

x=sin（2*pi*f1*n）+sin（2*pi*f2*n）;

y1=x（1:

2:

100）;

y2=decimate（x,M,'

fir'

figure

（1）;

stem（n,x（1:

N））;

title（'

inputsequence'

fudu'

figure

（2）;

N/2-1;

stem（n,y1）;

outputsequencewithoutLP'

figure（3）;

m=0:

N/M-1;

stem（m,y2（1:

N/M））;

outputsequencewithLP'

figure（4）;

[h,w]=freqz（x）;

plot（w（1:

512）,abs（h（1:

512）））;

frequencyspectrumoftheinputsequence'

w'

figure（5）;

[h,w]=freqz（y1）;

frequencyspectrumoftheoutputsequencewithoutLP'

figure（6）;

[h,w]=freqz（y2）;

（3）输入信号X（n）为归一化频率f1=0.043，f2=0.31的两个正弦信号相加而成，长度N=50，插因