数字信号处理实验1

数字信号处理实验一
第一题：
（1）
（2）
（3）
（4）
（5）
（6）
第二题：
（1）加法程序(seqadd.m)
function [y,n] = seqadd(x1,n1,x2,n2)
% [y,n] = seqadd(x1,n1,x2,n2)
% -----------------------------
% 实现y(n) = x1(n)+x2(n)
% y = 在包含n1和n2的n点上求序列和,
% x1 = 在位置向量n1上的第一序列
% x2 = 在位置向量n2上的第二序列(n2可与n1不同)
n = min(min(n1),min(n2)):max(max(n1),max(n2)); % y(n)的长度
y1 = zeros(1,length(n)); y2 = y1; % 初始化
y1(find((n>=min(n1))&(n<=max(n1))))=x1; % 具有y的长度的x1
y2(find((n>=min(n2))&(n<=max(n2))))=x2; % 具有y的长度的x2
y = y1+y2; % 序列相加
(2)乘法程序（seqmult.m）
function [y,n] = seqmult(x1,n1,x2,n2)
% [y,ny] = seqmult(x1,n1,x2,n2)
% -----------------------------
% 实现y(n) = x1(n)*x2(n)
% y = 在n区间上的乘积序列,ny 包含n1 和n2
% x1 = 在n1上的第一序列
% x2 = 在n2上的第二序列(n2可与n1不同)
n = min(min(n1),min(n2)):max(max(n1),max(n2)); % y(n)的位置向量
y1 = zeros(1,length(n)); y2 = y1; % y1,y2初始置零y1(find((n>=min(n1))&(n<=max(n1))))=x1; % 具有y的长度的x1
y2(find((n>=min(n2))&(n<=max(n2))))=x2; % 具有y的长度的x2
y = y1 .* y2; % 序列相乘
(3)反转程序（seqfold.m）
function [y,yn]=seqfold(x,nx)
y = fliplr(x); % 将序列数值左右翻转
yn = -fliplr(nx); % 将序列位置对零位置左右翻转，故同时改变正负号
(4)移位程序（seqshift.m）
function [y,ny] = seqshift(x,nx,n0)
% [y,ny] = seqshift(x,nx,n0)
% -------------------------
% 实现y(n) = x(n-n0)
% n0为平移样本数
ny = nx + n0; % 位置向量移位
y = x; % 序列的值不变
(5)第二题的程序（ex2.m）
x1=[1:5];
Xn=[x1,zeros(1,5)];
n1=0:9;
x2=[1,2,1,2];
Hn=[x2,zeros(1,6)];
n2=0:9;
%移位
[y,yn] = seqshift(Xn,n1,2);
subplot(3,2,1);stem(yn,y);title('右移2位的Xn'); [y,yn] = seqshift(Hn,n2,2);
subplot(3,2,2);stem(yn,y);title('右移2位的Hn'); %反转
[y,yn]=seqfold(Xn,n1);
subplot(3,2,3);stem(yn,y);title('反转的Xn'); [y,yn]=seqfold(Hn,n2);
subplot(3,2,4);stem(yn,y);title('反转的Hn');
%加法
[y,n] = seqadd(Xn,n1,Hn,n2);
subplot(3,2,5);stem(n,y);title('加法');
%乘法
[y,n] = seqmult(Xn,n1,Hn,n2);
subplot(3,2,6);stem(n,y);title('乘法');
第二题的结果：
第三题：
（1）卷积程序（juanji.m）
function y=juanji(x,nx,h,nh)
[h,nh]=seqfold(h,nh);%将序列hn翻折
nxmin=min(nx);nxmax=max(nx);
nhmin=min(nh);nhmax=max(nh);
A=[nxmin nxmax nhmin nhmax];
Min=min(A);Max=max(A);
N=abs(Max-Min);%确定要使序列hn移位的区间for i=-N:N
[h1,nh1] = seqshift(h,nh,i);%序列移位
[y,n] = seqmult(x,nx,h1,nh1);%序列移位后相乘Y(N+i+1)=sum(y);%两序列相乘后求和
end
X=-N:N;
stem(X,Y)
第三题结果：。