关于matlab的两个小程序的解释说明

微信公开课小程序热门讨论「关于matlab的两个小程序的解释说明」最新回复-好久没做过这类地解答了，希望我说的还算明白_________________________functionpr1()%11差分方程a=[1,1,0.9];%通过a,b构造了一个滤波器，传递函数为H(z)=1/(1z^(1)+0.9*z^(2)),a构造分母，b构造分子b=1;x=impseq(0,20,120)%impseq是用来构造单位冲击函数，时间t的范围为[20,120],当且仅当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)[z,p,g]=tf2zp(b,a)%[Z,P,K]=TF2ZP(NUM,DEN)findsthezeros,poles,andgains:%z为零点，p为极点，g为增益，如下：%(zz1)(zz2)...(zzn)%H(z)=K%(zp1)(zp2)...(zpn)zplane(z,p)%画出极零图function[x,n]=impseq(n0,n1,n2)n=[n1:n2];x=[(nn0)==0];functionpr2()%12experiment__________________________________________b=[0.0181,0.0543,0.0543,0.0181];%其实虚线框中的这部分的作用和上面函数中的FILTER函数的作用一样，就是构造一个滤波器，H(z)还是b为分母，a为分子a=[1.000,1.76,1.1829,0.2781];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([0,1,0,1]);%调节坐标轴范围xlabel(w(pi));ylabel(|H|);title(幅度，相位响应);subplot(2,1,2);plot(w/pi,angH);grid;xlabel(w(pi));ylabel(angle(H));
才悬赏100  哎 慢慢等吧~~ 更多有关「关于matlab的两个小程序的解释说明」的疑问请扫码关注微信公开课+小程序！了解更多