...9;n');ylabel('h(n)');figure(2)[z,p,g]=tf2zp(b,a);zplane(z,p)title('零极点');function[x,n]=chongji(n1,n2)n=[n1:n2];x=[n==0];functionshiyan1()a=[1,-1,0.9];b=1;x=chongji(-20,120);n=-20:120;h=filter(b,a,x);figure(1)stem(n,h);title('冲击响应');实验1-2运行结...
时间:2024-04-01 00:04栏目:教育交底
...don;>>xlabel('Time(sec)'),ylabel('f(t)');>>title('频率为100Hz的正弦脉冲信号');>>N=500;>>k=-N:N;>>W=2*pi*k/((2*N+1)*dt);>>Fw=dt*ft*exp(-j*t1'*W);>>subplot(222);>>plot(W,abs(Fw)),gridon;>>xlabel('\o...
时间:2024-04-29 08:35栏目:教案课件
...z(b,a,N);x=[ones(1,5),zeros(1,N-5)];y=filter(b,a,x);subplot(1,2,1);stem(h);title('直线型h(n)');subplot(1,2,2);stem(y);title('直线型y(n)');程序运行成果：2.级联形式①impseq.mfunction[x,n]=impseq(n0,ns,nf)n=[ns:nf];x=[(n-n0)==0];②casfilter.mfunctiony=casfilter(b0...
时间:2024-04-11 01:31栏目:读后感
...pace(-1,1,100);figure(1);figure(1);subplot(311);subplot(311);plot(y)plot(y)title('title('原始信号原始信号')')grid;grid;thr=0.4;thr=0.4;ythard=wthresh(y,'h',thr);ythard=wthresh(y,'h',thr);subplot(312);subplot(312);plot(ythard);plot(ythard);title('t...
时间:2024-03-30 21:35栏目:教案课件
...w=(-ws/2+(0:L-1)*ws/L)/(2*pi);plot(w,abs(X));ylabel('幅度谱');title('矩形窗截断');4-300-200-100010020030024681012141618幅度谱使用矩形窗截断b)使用hamming窗截断N=30;%数据的长度L=512;f1=100;f2=120;fs=600;T=1/fs;ws=2*pi*fs;t=(0:N-1)*T;x=cos(2*pi*f1*t)+cos(...
时间:2024-04-26 15:35栏目:教案课件
...;gridxlabel('\omega/\pi');ylabel('Groupdelay,samples');title('OriginalFilter*);F=0:0.001:0.4;g=grpdelay(n©F,2);%EqualizethepassbandGd=max(g)-g;%DesignthealIpassdelayequalizer[numzden,tau]=iirgrpdelay(6,F,[00.4],Gd);%设计六阶的全通滤波器[bza]=iirgrpdelay(6,F,[00.4]...
时间:2024-04-03 14:11栏目:教案课件
...;gridxlabel('\omega/\pi');ylabel('Groupdelay,samples');title('OriginalFilter*);F=0:0.001:0.4;g=grpdelay(n©F,2);%EqualizethepassbandGd=max(g)-g;%DesignthealIpassdelayequalizer[numzden,tau]=iirgrpdelay(6,F,[00.4],Gd);%设计六阶的全通滤波器[bza]=iirgrpdelay(6,F,[00.4]...
时间:2024-04-03 14:11栏目:教案课件
...变化的振幅xlabel('频率/\itHz');ylabel('振幅');title('图1：fs=100，N=1024');gridon;%两种信号叠加，x=A*sin(2*pi*B*t+C)+2*A*sin(2*pi*1.5*B*t+2.5*C);%信号y=fft(x,N);%对信号进行傅里叶变换yy=abs(y);%求得傅里叶变换后的振幅yy=yy*2/N;%...
时间:2024-04-26 15:35栏目:教案课件
...xlabel('频率/Hz');ylabel('振幅');xlabel('频率/Hz');ylabel('振幅');title('Ndata=136Nfft=128');gridon;n=0:Ndata-1;t=n/fs;%时间序列5绘出Nyquist频率之前的振幅xlabel('频率/Hz');ylabel('振幅');title('Ndata=136Nfft=512');gridon;FFT振幅谱也基本不受影响。相同,这样的...
时间:2024-04-07 00:57栏目:教案课件
...制cos(0t)、cos(30t)、cos(50t)和x(t)的波形图，给图形加title，网格线和x坐标标签，并且程序能够接受从键盘输入的和式中的项数。抄写程序Q3_1如下：clear,%Clearallvariablescloseall,%Closeallfigurewindowsdt=0.00001;%Specifythestepoftimevariablet=-2:d...
时间:2024-04-26 15:35栏目:教案课件
