西工大信号上机实验4非周期信号的频谱分析.docxVIP

上机实验4 非周期信号的频谱分析实验目的掌握用MATLAB编程，分析门信号的频谱；掌握用MATLAB编程，分析冲激信号的频谱；掌握用MATLAB编程，分析直流信号的频谱；掌握用MATLAB编程，分析阶跃信号的频谱；掌握用MATLAB编程，分析单边指数信号的频谱。实验内容与方法1 验证性实验（直接利用符号法进行编程）门信号的傅里叶变换MATLAB程序：clear all;Syms t w; %定义符号变量ut=sym(Heaviside(t+0.5)-Heaviside(t-0.5)); %借助阶跃函数定义门函数subplot(2,1,1);ezplot(ut);Hold on;axis([-1 1 0 1.1]); %给定坐标系范围plot([-0.5 -0.5],[0 1]); %绘制左跳变沿plot([0.5 0.5],[0 1]); %绘制右跳变沿Fw=fourier(ut,t,w); %对门函数进行傅立叶变换FFw=maple(convert,Fwpiecewise);%将傅立叶变换的符号表达式转换为样条曲线FFp=abs(FFw); %求幅度频谱subplot(2,1,2);ezplot(FFp,[-10*pi 10*pi]);%绘制幅度频谱axis([-10*pi 10*pi 0 1.1]);%end运行结果图如下：直流信号MATLAB程序：%直流信号的傅立叶变换clear all;clc;display(please input the value of a);a=input(a=);%键盘输入a的值syms t; %定义符号变量f=exp(-a*abs(t));subplot(1,2,1);ezplot(t); %绘制双边指数信号的波形axis([-2*pi 2*pi 0 1]);ylabel(时域波形);F=fourier(f);subplot(1,2,2);ezplot(abs(F));%绘制幅度频谱axis([-3 3 0 2/a]);ylabel(幅度频谱);%end运行结果如下：3）阶跃信号MATLAB程序：%阶跃信号的傅立叶变换clear all;clc;syms w; xw=1/(j*w); %0.5*sgn(t)的频谱ezplot(imag(xw)); %绘制其虚部axis([-3 3 -1.5*pi 1.5*pi]);hold on;y=0:0.01:pi;plot(0,y); %绘制直流部分y=pi:piplot(0,y,^);title(阶跃信号频谱);xlabel(\omega);axis([-pi pi -6 6]);x=-pi:0.001:pi;plot(x,0); %绘制横坐标y=-6:0.01:6plot(0,y); %绘制纵坐标%end运行结果如下：单边指数信号MATLAB程序：%单边指数信号的傅立叶变换clear all;clc;syms t v w x phae im re; %定义符号变量f=exp(-2*t)*sym(Heaviside(t));%f(t)=exp(-2*t)*u(t)FW=fourier(f);subplot(3,1,1);ezplot(f); %绘制时域波形axis([-1 2.5 0 1.1]);subplot(3,1,2);ezplot(abs(FW)); %绘制幅度谱im=imag(FW);re=real(FW);phase=atan(im/re);subplot(3,1,3);ezplot(phase); %绘制相位谱%end运行结果如下：2 程序设计题试确定下列信号的傅里叶变换的数学表达式f(t)=MATLAB程序：clear all;syms t w;f=exp(-2*t)*sym(Heaviside(t))+1;F=fourier(f,t,w)运行结果如下：f(t)=MATLAB程序：clear all;syms t w;f=exp(-1*t)*sym(Heaviside(t-1))+sym(Heaviside(t+1)-Heaviside(t-1));F=fourier(f,t,w)运行结果如下：f(t)=2U(t)+δ(t-4)MATLAB程序：clear all;syms t w;f=2*sym(Heaviside(t))+dirac(t-4)F=fourier(f)运行结果如下：