如果您无法下载资料，请参考说明：
1、部分资料下载需要金币，请确保您的账户上有足够的金币
2、已购买过的文档，再次下载不重复扣费
3、资料包下载后请先用软件解压，在使用对应软件打开
functionY=dirac(X)%DIRACDeltafunction.%DIRAC(X)iszeroforallX,exceptX==0whereitisinfinite.%DIRAC(X)isnotafunctioninthestrictsense,butrathera%distributionwithint(dirac(x-a)*f(x),-inf,inf)=f(a)and%diff(heaviside(x),x)=dirac(x).%SeealsoHEAVISIDE.%Copyright1993-2003TheMathWorks,Inc.%$Revision:1.1.6.2$$Date:2004/04/1622:21:56$Y=zeros(size(X));Y(X==0)=Inf;functionY=heaviside(X)%HEAVISIDEStepfunction.%HEAVISIDE(X)is0forX<0,1forX>0,andNaNforX==0.%HEAVISIDE(X)isnotafunctioninthestrictsense.%SeealsoDIRAC.%Copyright1993-2003TheMathWorks,Inc.%$Revision:1.1.6.2$$Date:2004/04/1622:23:24$Y=zeros(size(X));Y(X>0)=1;Y(X==0)=NaN;clc,clear;symstwf;f=sym('Heaviside(t)')-sym('Heaviside(t-1)');F=fourier(f)subplot(2,1,1);ezplot(f,[0:2,0:1.2]);subplot(2,1,2);ezplot(abs(F),[-10:10]);title('幅度谱')clc,clear;symstwf;f=exp(-2*abs(t))F=fourier(f)subplot(2,1,1);ezplot(f,[0:2,0:1.2]);subplot(2,1,2);ezplot(abs(F),[-10:10]);title('幅度谱')clc,clear;symstw;F=2*(sin(w)/w);f=ifourier(F,w,t)ezplot(f,[-2:2])clc,clear;b=[1];a=[132];w=linspace(0,5,200);H=freqs(b,a,w);figure(1);subplot(2,1,1);plot(w,abs(H));xlabel('w');ylabel('幅频特性');gridon;subplot(2,1,2);plot(w,angle(H));xlabel('w');ylabel('相频特性');gridon;
立即下载
