DSP用m程序(参考).docVIP

数字信号处理课程用m程序 DSP_1.m %简单线性卷积的计算 clear;close all;clc; xn = [0 0.5 1 1.5], % 输入序列 hn = [1 1 1], % 脉冲响应 yn = conv(hn,xn), % 输出序列 L = length(xn)+length(hn)-1,% 输出长度 n = 0:L-1; % 输出序列时间 subplot(3,1,1); %绘输入序列 xnL=[xn,zeros(1,length(hn)-1)];%输入补齐 stem(n,xnL,r*); title(输入序列); ylabel(序列值xn); grid on; subplot(3,1,2);%绘脉冲响应 hnL=[hn,zeros(1,length(xn)-1)]; stem(n,hnL,r*); title(脉冲响应); ylabel(序列值hn); grid on; subplot(3,1,3);%绘输出序列 stem(n,yn,r*); title(输出序列); ylabel(序列值yn); grid on; DSP_2.m %任意序列线性卷积的计算 clear;close all;%在命令窗口中按Ctrl+C中止程序 % xn = input(输入序列xn=), % 输入序列 % xns= input(输入序列起点xns=),% 输入序列的起点 % hn = input(脉冲响应hn=), % 脉冲响应 % hns= input(脉冲响应起点hns=),% 脉冲响应的起点 xn = [ 0.5 1 1.5], % 输入序列 xns= 1; % 输入序列的起点 hn = [1 1 1], % 脉冲响应 hns=0; %脉冲响应的起点 yn = conv(hn,xn), % 输出序列 yns=xns+hns, % 输出序列的起点 subplot(3,1,1); %绘输入序列 xnL=[xn,zeros(1,length(hn)-1)];%输入补齐 nx = xns:length(yn)+xns-1;% 输入序列时间 stem(nx,xnL,r*); title(输入序列); ylabel(序列值xn); grid on; subplot(3,1,2);%绘脉冲响应 hnL=[hn,zeros(1,length(xn)-1)];%脉冲补齐 nh = hns:length(yn)+hns-1;% 脉冲序列时间 stem(nh,hnL,r*); title(脉冲响应); ylabel(序列值hn); grid on; subplot(3,1,3);%绘输出序列 ny = yns:length(yn)+yns-1;% 输出序列时间 stem(ny,yn,r*); title(输出序列); ylabel(序列值yn); grid on; DSP_3.m % z变换的应用--频率响应 %求教材P67题1-36的系统频率响应 clear; close all;clc; b=[1,0]; %降幂输入系统函数的分子 a=[1,-1/2]; %降幂输入系统函数的分母 sys=tf(b,a,-1), %按Z显示系统函数 freqz(b,a), %以512点绘制频率响应曲线 DSP_4.m % z变换的应用--求零极点增益 % 根据系统函数求零极点增益形式 % 计算零-极点-增益 绘零极点图 clear; close all; b=[1,1,0]; %降幂输入系统函数的分子不足补0 a=[1,-1,1/2]; %降幂输入系统函数的分母不足补0 sys=tf(b,a,-1), %按Z显示系统函数 [r,p,k]=tf2zp(b,a),% 根据系统函数求零极点增益形式 zs=roots(b), % 计算零点 ps=roots(a), % 计算极点 ks=b(1)/a(1), % 计算增益 zplane(b,a); % 绘零极点图 DSP_5.m % z变换的应用--生成二阶因式形式 % 根据系统函数生成二阶因式形式 % 参考孙洪P39--Q3_46 clear; close all; b=[2,5,9,5,3]; %降幂输入系统函数分子不足补0 a=[5,45,2,1,1];%降幂输入系统函数分母不足补0 sys=tf(b,a,-1), %按Z显示系统函数 sos=tf2sos(b,a), %生成二阶因式形式 DSP_6.m % z变换的应用--求单位脉冲响应 % 根据系统函数求单位脉冲响应 %求教