matlab打阶跃函数和冲击函数,阶跃函数与冲激函数的MATLAB的实现与求解.ppt

阶跃函数与冲激函数的MATLAB的实现与求解

PPT模板下载：/moban/ Z域中阶跃函数与冲激函数的求解 及MATLAB的实现 1 2 3 4 何为冲激响应、阶跃响应 冲激响应与阶跃响应的关系 MATLAB的实现 MATLAB仿真结果 目录 何为冲激响应、阶跃响应？ 冲激响应 冲激响应也叫单位样本响应是指输入单位样本序列 时数字滤波器的输出 记为 h[n] 阶跃响应 阶跃响应也叫单位阶跃响应是指输入单位阶跃序列 时离散时间系统的输 出记为 s[n] 数字滤波器 h[n] 冲激响应和阶跃响应的关系 h[n] H[z] MATLAB的实现 1.impz　　功能：求解数字系统的冲激响应。　　调用格式：　　［h，t］＝impz(b，a，n)；求解数字系统的冲激响应h，取样点数为n。 2.dstep　　功能：求解数字系统的阶跃应。 　　调用格式：　　［s，t］＝dstep(b，a，n)；求解数字系统的阶跃响应h，取样点数为n。 MATLAB实现 3. residuez　　功能：求逆Z变换。　　调用格式：　　[r,p,k] = residuez(num,den); 例题 4.72：一个因果LTI离散时间系统由差分方程 y[n]=0.4y[n-1]+0.05y[n-2]+3x[n] 描述，其中x[n]和y[n]分别是系统的输入和输出序列。1.求该系统的传输函数H(z)2.求该系统的冲激响应h[n]3.求该系统的阶跃响应s[n] MATLAB求解 a=[1,-0.4,0.05];b=[3];n=20; hn=impz(b,a,n);%冲激响应? sn=dstep(b,a,n);%阶跃响应? num = b;den = a;[r,p,k] = residuez(num,den);%逆Z变换 disp('Residues');disp(r') disp('Poles');disp(p') disp('Constants');disp(k) subplot(2,1,1),stem(hn,'filled');%显示冲激响应曲线? title('单位样本响应');ylabel('h(n)');xlabel('n'); axis([0,n,1.1*min(hn),1.1*max(hn)]); subplot(2,1,2),stem(sn,'filled');%显示阶跃响应曲线? title('单位阶跃响应');ylabel('s(n)');xlabel('n'); axis([0,n,1.1*min(sn),1.1*max(sn)]); 结果显示 输入文字 在此录入上述图表的综合分析结论 在此录入上述图表的综合分析结论 在此录入上述图表的综合分析结论 在此录入上述图表的综合分析结论 在此输入标题 输入文字 在此录入上述图表的综合分析结论 在此录入上述图表的综合分析结论 在此录入上述图表的综合分析结论 在此录入上述图表的综合分析结论 “ “ 谢谢您的耐心阅读！ a=[1,-0.4,0.05];b=[3];n=20; hn=impz(b,a,n);%3??¤?ìó|? sn=dstep(b,a,n);%?×???ìó|? num = b;den = a;[r,p,k] = residuez(num,den); disp('Residues');disp(r') disp('Poles');disp(p') disp('Constants');disp(k) subplot(2,1,1),stem(hn,'filled');%??ê?3??¤?ìó|?ú??? title('?μí3μ?μ￥???ù±??ìó|');ylabel('h(n)');xlabel('n'); axis([0,n,1.1*min(hn),1.1*max(hn)]); subplot(2,1,2),stem(sn,'filled');%??ê??×???ìó|?ú??? title('?μí3μ?μ￥???×???ìó|');ylabel('s(n)');xlabel('n'); axis([0,n,1.1*min(sn),1.1*max(sn)]); * * a=[1,-0.4,0.05];b=[3];n=20; hn=impz(b,a,n);%3??¤?ìó|? sn=dstep(b,a,n);%?×???ìó|? num = b;den = a;[r,p,k] = residuez(num,den); disp('Residues');disp(r') disp('Poles');disp(p') disp('Constants