MATLAB_数字信号处理_03_离散系统中的卷积与z变换

SYSU_SECE_MATLAB_数字信号处理_03

离散系统中的卷积与z变换

%%%题目一
figure(1)
a=[1 -0.7 0];
b=[2 0 -1];
N=16;
n=0:N-1;
x=[n>=3];
subplot(3,1,1);stem(n,x,'filled');
xlabel('n');ylabel('x(n)');title('输入序列x(n)');
hn=impz(b,a,n);
subplot(3,1,2);stem(n,hn,'filled');
xlabel('n');ylabel('h(n)');title('系统冲激响应h(n)');
y=conv(x,hn);
subplot(3,1,3);stem(y,'filled');
xlabel('n');ylabel('y(n)');title('系统输出响应y(n)');

%part1---solution1:
figure(2)
n=0:31;
xn=[n>=3];
subplot(2,1,1);stem(n,xn,'filled');
xlabel('n');ylabel('x(n)');title('T2-(1)输入序列x(n)');
a=[1 -0.5 0 0 0 0 1 -0.5];
b=[1 -1 1 0 0 0 0 0];
y=dlsim(b,a,xn);
subplot(2,1,2);stem(y,'filled');
xlabel('n');ylabel('y(n)');title('T2-(1)响应y(n)-dlsim方法');
%part1---solution2:
figure(3)
n=0:31;
xn=[n>=3];
y01=0;x01=0;
a=[1 -0.5 0 0 0 0 1 -0.5];
b=[1 -1 1 0 0 0 0 0];
xi=filtic(b,a,x01,y01);
subplot(2,1,1);stem(n,xn,'filled');
xlabel('n');ylabel('x(n)');title('T2-(1)输入序列x(n)');
yn=filter(b,a,xn,xi);
subplot(2,1,2);stem(y,'filled');
xlabel('n');ylabel('y(n)');title('T2-(1)响应y(n)-filter方法');
%part2---solution1:
figure(4)
n=0:31;
xn=[n==0]-[n==5];
subplot(2,1,1);stem(n,xn,'filled');
xlabel('n');ylabel('x(n)');title('T2-(2)输入序列x(n)');
a=[1 -0.5 0 0 0 0 1 -0.5];
b=[1 -1 1 0 0 0 0 0];
y=dlsim(b,a,xn);
subplot(2,1,2);stem(y,'filled');
xlabel('n');ylabel('y(n)');title('T2-(3)响应y(n)-dlsim方法');
%part2---solution2:
figure(5)
n=0:31;
xn=[n==0]-[n==5];
y01=0;x01=0;
a=[1 -0.5 0 0 0 0 1 -0.5];
b=[1 -1 1 0 0 0 0 0];
xi=filtic(b,a,x01,y01);
subplot(2,1,1);stem(n,xn,'filled');
xlabel('n');ylabel('x(n)');title('T2-(2)输入序列x(n)');
yn=filter(b,a,xn,xi);
subplot(2,1,2);stem(yn,'filled');
xlabel('n');ylabel('y(n)');title('T2-(2)响应y(n)-filter方法');
%part3---solution1:
figure(6)
n=0:31;
xn=exp(0.1*n).*[n>=3];
subplot(2,1,1);stem(n,xn,'filled');
xlabel('n');ylabel('x(n)');title('T2-(3)输入序列x(n)');
a=[1 -0.5 0 0 0 0 1 -0.5];
b=[1 -1 1 0 0 0 0 0];
y=dlsim(b,a,xn);
subplot(2,1,2);stem(y,'filled');
xlabel('n');ylabel('y(n)');title('T2-(3)响应y(n)-dlsim方法');
%part3---solution2:
figure(7)
n=0:31;
xn=exp(0.1*n).*[n>=3];
y01=0;x01=0;
a=[1 -0.5 0 0 0 0 1 -0.5];
b=[1 -1 1 0 0 0 0 0];
xi=filtic(b,a,x01,y01);
subplot(2,1,1);stem(n,xn,'filled');
xlabel('n');ylabel('x(n)');title('T2-(3)输入序列x(n)');
yn=filter(b,a,xn,xi);
subplot(2,1,2);stem(yn,'filled');
xlabel('n');ylabel('y(n)');title('T2-(3)响应y(n)-filter方法');

%part1:题目三---（1）
figure(8)
a=[1 -0.602012 0.495684 -0.035924];
b=[0.187632 0.241242 0.241242 0.187632];
y0=[5 5];
xi=filtic(b,a,y0);
N=32;n=0:N-1;
x=[n==3];
y_zi=filter(b,a,[n>100],xi);
subplot(3,1,1);stem(n,y_zi,'filled');xlabel('n');ylabel('y_zi(n)');
title('T3-(1)-零输入响应');
y_zs=filter(b,a,x,filtic(b,a,0));
subplot(3,1,2);stem(n,y_zs,'filled');xlabel('n');ylabel('y_zs(n)');
title('T3-(1)-零状态响应');
y=filter(b,a,x,xi);
subplot(3,1,3);stem(n,y,'filled');xlabel('n');ylabel('y(n)');
title('T3-(1)-完全响应');
%part2:题目三---（2）
figure(9)
a=[1 -0.602012 0.495684 -0.035924];
b=[0.187632 0.241242 0.241242 0.187632];
y0=[5 5];
xi=filtic(b,a,y0);
N=32;n=0:N-1;
x=[ones(1,5),zeros(1,N-5)];
y_zi=filter(b,a,[n>100],xi);
subplot(3,1,1);stem(n,y_zi,'filled');xlabel('n');ylabel('y_zi(n)');
title('T3-(2)-零输入响应');
y_zs=filter(b,a,x,filtic(b,a,0));
subplot(3,1,2);stem(n,y_zs,'filled');xlabel('n');ylabel('y_zs(n)');
title('T3-(2)-零状态响应');
y=filter(b,a,x,xi);
subplot(3,1,3);stem(n,y,'filled');xlabel('n');ylabel('y(n)');
title('T3-(3)-完全响应');
%part3:题目三---（3）
figure(10)
a=[1 -0.602012 0.495684 -0.035924];
b=[0.187632 0.241242 0.241242 0.187632];
y0=[5 5];
xi=filtic(b,a,y0);
N=32;n=0:N-1;
x=cos(2*pi*n/3)+sin(3*pi*n/10);
y_zi=filter(b,a,[n>100],xi);
subplot(3,1,1);stem(n,y_zi,'filled');xlabel('n');ylabel('y_zi(n)');
title('T3-(3)-零输入响应');
y_zs=filter(b,a,x,filtic(b,a,0));
subplot(3,1,2);stem(n,y_zs,'filled');xlabel('n');ylabel('y_zs(n)');
title('T3-(2)-零状态响应');
y=filter(b,a,x,xi);
subplot(3,1,3);stem(n,y,'filled');xlabel('n');ylabel('y(n)');
title('T3-(3)-完全响应');

%%%题目四=============================
syms n a w0

%part1:
x1=n*a^n;
X1=ztrans(x1)
%part2:
x2=exp(-a*n)*sin(n*w0);
X2=ztrans(x2)

%%%题目五
%part1:
figure(11)
syms z
b=[10 20 0 0];
a=[1 8 19 12];
[r,p,c]=residuez(b,a);
N=10;n=0:N-1;
xn=r(1)*p(1).^n+r(2)*p(2).^n+r(3)*p(3).^n;
subplot(2,1,1);stem(n,xn,'filled');title('')
h_=impz(b,a,N);
subplot(2,1,2);stem(n,h_,'filled');title('')
%part2:
figure(12)
syms z
b=[1 0 0 0];
a=[1 -0.9 -0.81 0.729];
[r,p,c]=residuez(b,a);
N=10;n=0:N-1;
xn=5*r(1)*p(1).^n+r(2)*p(2).^n.*(n-1)+r(3)*p(3).^n;
subplot(2,1,1);stem(n,xn,'filled');title('')
h_=impz(b,a,N);
subplot(2,1,2);stem(n,h_,'filled');title('')

% 以上用到的MATLAB函数
% impz 求解数字系统的冲激响应
% filter(b,a,x,zi) 在zi中指定x的初始状态
% filtic(b,a,y,x) 求给定输入x和y时的初始状态
% dlsim(b,a,x) 求输入信号为x时系统的响应
% conv(x,n) 求取两个有限长序列x和h的卷积

上述题目基本全部改编/出自自西安电子科技大学出版社出版的《数字信号处理实验（MATLAB版）》，可以在其中找到原理部分，我们理论课使用的教材也是西电的：《数字信号处理》（第四版）高西全 丁玉美

声明：不保证代码的正确性，只是提供一种思路……