MATLAB仿真及电子信息应用（第2版）.ppt

其中用到两个自定义的子函数：rectcos和rectsin，它们的定义如下：functiony=rectcos(t,n,w)y=(abs(t)=1).*cos(n*w*t);functiony=rectsin(t,n,w)y=(abs(t)=1).*sin(n*w*t);例7-12绘制重复周期为6、脉冲宽度为2、脉冲幅度为1的矩形波的幅度谱和相位谱。MATLAB程序T=6;tao=2;w=2*pi/T;a0=quadl(@rect,-3,3)/T;N=20;an=zeros(1,N);bn=zeros(1,N);fork=1:Nan(k)=quadl(@rectcos,-3,3,[],[],k,w)*2/T;%计算anbn(k)=quadl(@rectsin,-3,3,[],[],k,w)*2/T;%计算bnendc0=a0;%求c0cn=sqrt(an.^2+bn.^2);%求cnphi=-atan(bn./an);%计算相位%根据an的正负修正当bn=0时的相位值fork=1:Nifphi(k)==0ifan(k)0phi(k)=pi;endendend%绘制幅度谱和相位谱n=1:N;subplot(1,2,1);stem(n*w,cn);xlabel(nw);ylabel(cn);axis([0,21,0,0.6]);title(幅度谱);subplot(1,2,2);stem(n*w,phi);xlabel(nw);ylabel(phi);axis([0,21,0,3.5]);title(相位谱);例7-12中周期矩形波的脉冲宽度不变，将重复周期由原来的6变为10，重绘其幅度谱，与上图7-9相比较，可以看出，当脉冲宽度不变、周期增大时，谱线的幅度随之减小，谱线之间的间隔变小，谱线更密。例7-12中周期矩形波的周期不变，脉冲宽度由原来的2变为1，其幅度谱如下，与图7-9相比较，可以看出，当周期不变、脉冲宽度减小时，谱线的幅度随之减小，信号频带宽度增加。7.4.2傅立叶变换MATLAB程序symst;%定义符号变量tf=sym(Heaviside(t+1))-sym(Heaviside(t-1));%输入符号函数fF=fourier(f);%求f的傅立叶变换Fezplot(abs(F));%绘制幅度谱例7-14进行逆傅立叶变换symsxwuv;F1=exp(-w^2/4);%符号函数的变量是wf1=ifourier(F1)%返回函数的参量是默认变量xf1=1/pi^(1/2)*exp(-x^2)F2=exp(-abs(x));%符号函数的变量是xf2=ifourier(F2)%返回函数的参量是变量tf2=1/(1+t^2)/piF3=2*exp(-abs(w))-1;%符号函数的变量是wf3=ifourier(F3,u)%返回函数的参量是指定变量uf3=-dirac(u)+2/pi/(1+u^2)T=4;%时间宽度N=200;%时域采样点数delta=T/N;%时域采样间隔t=linspace(-T/2,T/2-T/N,N);%时间向量f=(t=-1)-(t=1);%f(t)f1=f.*exp(-j*10*t);%f1(t)f2=f.*exp(j*10*t);%f2(t)W=20*pi;%频域范围K=1000;