实验4-非周期信号的频域分析.docVIP

实验4 非周期信号的频域分析 一、实验目的 理解非周期信号的频域分析方法，掌握典型信号的幅度谱和相位谱，理解信号的调制特性，掌握傅里叶变换的性质：尺度变换、时移、频移、卷积定理、对称性、微分特性。 二、实验内容 1、利用符号函数fourier和ifourier求傅里叶变换和傅里叶逆变换。 a. 利用符号函数fourier求教材(上册)P114双边指数信号当a=3时的傅里叶变换表达式。 代码: syms t x fourier(exp(-x^2),x,t); ezplot(t,[-3 3]); syms t a f=sym(exp(-a*abs(t))); f=subs(f,a,3); F=fourier(f); F =6/(9+w^2) b. 利用符号函数ifourier求教材(上册)P114公式(3-32)第一式当a=1时的傅里叶逆变换表达式。 代码: ezplot(F); f1=subs(f,a,1); F=fourier(f1); F F = 2/(1+w^2) ezplot(F); f2=ifourier(F,t); f2 f2 = exp(t)*heaviside(-t)+exp(-t)*heaviside(t) ezplot(f2); c. 利用符号函数fourier和ezplot画出及其幅频谱。\ 代码: syms t w f=0.5*exp(-2*t)*heaviside(t); F=fourier(f); ezplot(f); ezplot(abs(F)); 2、矩形脉冲信号和周期矩形脉冲信号的频谱比较 画出教材(上册)P110图3-18右边的离散谱和连续谱。 代码: pinpu(1,4,1); pinpu(1,6,1); syms t1 w F=fourier((heaviside(t1+0.5*1)-heaviside(t1-0.5*1))); ezplot(F); function pinpu(E,T,t) w=-16*pi/T:2*pi/T:16*pi/T; Fn=E*t/T*sinc(0.5*t*w/pi); stem(w,Fn); hold on; w1=-16*pi/T:0.1:16*pi/T; Fn1=E*t/T*sinc(0.5*t*w1/pi); plot(w1,Fn1); hold off; 3、幅度调制信号及其频谱 已知线性调制信号表示式如下： a. ； b. 式中，试分别画出它们的波形图和频谱图 代码: syms t w a=1; f=cos(a*t)*cos(9*a*t); F=fourier(f); F F = 1/2*pi*(dirac(w+10)+dirac(w+8)+dirac(w-10)+dirac(w-8)) ezplot(F); ezplot(f); f1=(1.5+sin(a*t))*cos(9*a*t); F1=fourier(f1); ezplot(f1); ezplot(F1); F1 F1 = 1/2*pi*(3*dirac(w+9)+i*dirac(w+10)-i*dirac(w+8)-i*dirac(w-10)+3*dirac(w-9)+i*dirac(w-8)) 4、傅里叶变换的性质（尺度变换、时移、频移、卷积定理、对称性） a. 设，求的频谱，并与的频谱进行比较。(提示：利用单位阶跃信号的符号函数Heaviside) syms t w f=heaviside(t+1)-heaviside(t-1); y=fourier(f); ezplot(y); f1=heaviside(2*t+1)-heaviside(2*t-1); y1=fourier(f1); ezplot(y1); b. 画出、 和的幅度谱和相位谱，观察信号时移对信号频谱的影响。 代码: ezplot(y); w=[-10 10]; f=1/2*(w*i+2); f1=angle(f); plot(w,f1); syms t w f=0.5*exp(-2*t)*heaviside(t); F=fourier(f); F F = 1/2/(2+i*w) ezplot(abs(F)); w=[-10 10]; F=1/2*(i*w+2); Fa=angle(F); plot(w,Fa); Fb=Fa-0.4*w; plot(w,Fb); Fc=Fa+0.4*w; plot(w,Fc); c. 画出、和的频谱，进行相互比较。 代码: syms t w f=heaviside(t+1)-heaviside(t-1); f1=f*exp(-i*20*