matlab冲激响应函数(关于matlab的两个小程序的解释说明)

本文目录
关于matlab的两个小程序的解释说明
好久没做过这类地解答了,希望我说的还算明白
_________________________
function pr1()
%11差分方程
a=; %通过a,b构造了一个滤波器,传递函数为H(z)=1/(1-z^(-1)+0.9*z^(-2)),a构造分母,b构造分子
b=1;
x=impseq(0,-20,120) %impseq是用来构造单位冲击函数,时间t的范围为,当且仅当t=0处为1,即x(21)=1,其余均为零
n=-20:120;
h=filter(b,a,x); %以a,b构造的的滤波器,当输入为x时,输出为h,因为x是单位冲击函数,所以h就是该滤波器的单位冲击响应了
figure(1)
stem(n,h); %以n为横轴,h为纵轴画出火柴杆图
title(’冲激响应’);
xlabel(’n’);
ylabel(’h(n)’);
%
figure(2)
= TF2ZP(NUM,DEN) finds the zeros, poles, and gains:
% z为零点,p为极点,g为增益,如下:
% (z-z1)(z-z2)...(z-zn)
% H(z) = K ------------------------
% (z-p1)(z-p2)...(z-pn)
zplane(z,p) %画出极零图
function =impseq(n0,n1,n2)
n=;
x=;
function pr2()
%12 experiment
__________________________________________
b=; %其实虚线框中的这部分的作用和上面函数中的FILTER函数的作用一样,就是构造一个滤波器,H(z)还是b为分母,a为分子
a=;
m=0:length(b)-1;
l=0:length(a)-1;
K=500;
k=1:K;
w=pi*k/K;
H=(b*exp(-j*m’*w))./(a*exp(-j*l’*w)); %求出单位冲击响应函数
__________________________________________________________________________________
magH=abs(H); %求出单位冲击响应函数的幅值
angH=angle(H); %求出单位冲击响应函数的相位
figure(1)
subplot(2,1,1); %做出子图
plot(w/pi,magH);
grid; %添加栅格
axis(); %调节坐标轴范围
xlabel(’w(pi)’);
ylabel(’|H|’);
title(’幅度,相位响应’);
subplot(2,1,2);
plot(w/pi,angH);
grid;
xlabel(’w(pi)’);
ylabel(’angle(H)’);
matlab如何求一个函数的冲击响应
1、g=finverse(f):返回符号函数f的反函数g。其中,f是一个符号函数表达式,其变量为x。求得的反函数g是一个满足g(f(x))=x的符号函数。
》》
syms
x;
》》
f=sym(2/sin(x));
》》
finverse(f)
ans
=
asin(2/x)
2、g=finverse(f,v):返回自变量v的符号函数f的反函数。求得的反函数g是一个满足g(f(v))=v的符号函数。当f包含不止一个符号变量时,往往调用这个格式。
当finverse求得的解不惟一时,matlab会给出警告。
》》
syms
x;
》》
f=sym(x^2+1);
》》
finverse(f)
warning:
finverse(x^2+1)
is
not
unique.
》
in
c:\matlab6p5\toolbox\symbolic\@sym\finverse.m
at
line
43
ans
=
(-1+x)^(1/2)
怎样用MATLAB画出窗函数低通滤波器的单位脉冲响应
hd=sin(wd*n)./(pi*n); % 理想冲激响应
w=hamming(N)’; % 海明窗
h=hd.*w; % 实际冲激响应
H=20*log10(abs(fft(h,1024))); % 实际滤波器的幅频特性
HH=;
subplot(221),stem(nn,hd,’k’);
xlabel(’n’);title(’理想冲激响应’);axis();
subplot(222),stem(nn,w,’k’);axis();
title(’海明窗’);xlabel(’n’);
subplot(223),stem(nn,h,’k’);
axis();xlabel(’n’);title(’实际冲激响应’);
w=(-512:511)/511;
subplot(224),plot(w,HH,’k’);
axis();xlabel(’\omega/\pi’);title(’滤波器幅频特性’);
set(gcf,’color’,’w’);参考这段程序
matlab求冲激响应
clc
clear;
n=-20:100;
x=zeros(1,length(n));
x(21)=1;
b=;
a=;
filter(b,a,x)
matlab中如何画出冲激函数
主程序如下:
n=;
x1=impseq(0,-5,5);
subplot(2,2,1);stem(n,x1);title(’单位取样序列’);
xlabel(’n’);ylabel(’x(n)’);
在这之前,要定义一个函数,就是一个.m文件,如下:
function=impseq(n0,n1,n2)
if((n0《n1)|(n0》n2)|(n1》n2))
error(’参数必须满足n1《=n0《=n2)
end
n=;
x=;
注意.m文件名要与函数名一致,即使用impseq!

更多文章:
全球新冠肺炎疫情背景下航运发展(盐田港复苏日志:半年历劫从“低谷”到“爆仓” 疫情之后巨轮如何越洋航行)
2026年9月7日 17:10
matlab求解带字母参数方程组(我想matlab求一个关于x,y的方程组 ab c d f e h m n 都是参数)
2026年9月7日 16:30
oracle中的循环语句(下面哪个不是oracle程序设计中的循环语句 a for)
2026年9月7日 15:30
电脑里2个系统怎么删除一个(电脑开机显示有两个系统,如何删除一个)
2026年9月7日 12:20
scrollthrough意思(“scroll”是什么意思)
2026年9月7日 08:00
怎么激活keygen(注册机如何激活cad2008一个简单激活cad2008的方法)
2026年9月7日 06:30



