重庆大学通信工程学院数字信号处理实验报告5.docx

重庆大学通信工程学院数字信号处理实验报告5.docx

《重庆大学通信工程学院数字信号处理实验报告5.docx》由会员分享，可在线阅读，更多相关《重庆大学通信工程学院数字信号处理实验报告5.docx（11页珍藏版）》请在冰豆网上搜索。

重庆大学通信工程学院数字信号处理实验报告5

实验5：

变换域中的线性移不变离散时间系统

1.第一题

b=[0.20-0.3];

a=[1-0.40.8];

freqz（b,a）

2.第二题

b=[0.20-0.3];

a=[0.8-0.41];

freqz（b,a）

3.第三题

b=[0.30-0.2];

a=[1-0.40.8];

freqz（b,a）

4.第四题

b=[0.50.5];

a=1;

freqz（b,a）

5.第五题

b=[10.2-0.40.21];

zplane（b,1）

hn=impz（b,1）

un=ones（1,100）;

yn=conv（un,hn）

6.第六题

clc;

clear;

closeall;

%%稳定非最小相移系统

B=[0.3-0.5];

A=[1-0.70.8];

figure

（1）

zplane（B,A）;

%%全通

B1=[15/3];

A1=[13/5];

figure

（2）

zplane（B1,A1）;

%%最小相移系统

B2=[0.5-0.3];

A2=[1-0.70.8];

figure（3）

zplane（B2,A2）;

%%频率响应

N=1024;

figure（4）

freqz（B,A,N）;

figure（5）

freqz（B2,A2,N）;

9.第九题

clc;

clear;

closeall;

%%产生信号

h=[1234];

fori=0:

18%定义x（n）

x（i+1）=i+2;

end

Y2=conv（x,h）;

figure

（1）

stem（Y2）;

%%重叠相加法

N=5;

Lx=length（x）;%序列长度

M=length（h）;%h（n）长度

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

t=zeros（1,M-1）;

Y1=zeros（1,Lx+M-1）;

a=floor（Lx/N）;%计算出迭代次数

fork=0:

a

A=x（k*N+1:

k*N+N）;

y1=fft（A,N+M-1）;%利用fft进行运算

y2=fft（h,N+M-1）;

y3=y1.*y2;

q=ifft（y3,N+M-1）;

Y1（k*N+1:

k*N+M-1）=q（1:

M-1）+t（1:

M-1）;%重叠相加

Y1（k*N+M:

k*N+N）=q（M:

N）;

t（1:

M-1）=q（N+1:

N+M-1）;

end

Y1（（a+1）*N+1:

（a+1）*N+M-1）=q（N+1:

N+M-1）;

figure

（2）

stem（Y1（1:

Lx+M-1））;