西工大软件工程信号与系统实验四.docxVIP

西北工业大学《信号与系统》实验报告学 院：软件与微电子学院学　　号：姓　　名：专 业：实验时间：实验地点： 软件学院实验室310 指导教师：陈勇西北工业大学2017年9月实验目的运用MATLAB软件对理论课程所学的离散信号与连续信号的相关知识进行表示与验证实验要求学会使用MATLAB进行连续时间傅里叶变换的数值近似；了解连续时间傅里叶变换性质；了解系统的时域和频域特性；学会使用MATLAB函数freqs，residue；学会用部分分式展开求微分方程的单位冲激响应。实验设备（环境） 1. Windows 10 64位系统 2. MATLAB R 2014b 软件环境实验内容与步骤4.2 连续时间傅里叶变换的数值近似基本题（a）（b）（c）（d）（e）（f）（g）4.3 连续时间傅里叶变换性质基本题（a）（b）4.4 系统的时域和频域特性基本题（a）（b）（c）（d）4.5 用部分分式展开求微分方程的单位冲激响应基本题（a）（b）（c）实验结果4.2连续时间傅里叶变换的数值近似(A)代码如下：结果如下：(B)代码如下：clear;clc;x1=sym(exp(-2*(t-5))*Heaviside(t-5))x2=sym(exp(2*(t-5))*Heaviside(-t+5))y1=fourier(x1)y2=fourier(x2)y=simple(y1+y2)结果如下： x1 =exp(-2*(t-5))*Heaviside(t-5) x2 =exp(2*(t-5))*Heaviside(-t+5)y1 =1/(2+i*w)*exp(-5*i*w) y2 =1/(2-i*w)*exp(-5*i*w) y =4*exp(-5*i*w)/(4+w^2)(C)代码如下：clear;clc;tau=0.01;T=10;t=[0:tau:T-tau];N=length(t)y=exp(-2*abs(t-5));y1=fft(y)y2=fftshift(tau*fft(y)错误分析：由于N的长度为1000，故计算出的样本Y(jw)值有1000个，由于计算结果太多，因此没有将运行结果保存过来。(D)构造一个频率样本向量w，它按照 w=-(pi/tau)+(0:N-1)*(2*pi/(N*tau));与存在向量Y中的值相对应。(E)代码如下：clear;clc;tau=0.01;T=10;t=[0:tau:T-tau];N=length(t)w=-(pi/tau)+(0:N-1)*(2*pi/(N*tau));y=exp(-2*abs(t-5));y=fftshift(tau*fft(y));for t=1:N x(t)=y(t)*exp(5*j*w(t));end结果分析：设F[ f ( t ) ]傅立叶变换为F(W),由傅立叶变换的时移特性可知，F[ f ( t – t0 ) ]=F ( w )*exp(-j*w*t0),可得?F ( w )= F[ f ( t – t0 ) ]* exp(j*w*t0),因此本题中X(jw)=Y(jw)* exp(5jw), 而Y(jw)可由快速傅立叶变换得到。(F)代码如下：clear;clc;tau=0.01;T=10;t=[0:tau:T-tau];N=length(t)w=-(pi/tau)+(0:N-1)*(2*pi/(N*tau));%通过近似求的X(jw)y=exp(-2*abs(t-5));y=fftshift(tau*fft(y));for i=1:N x(i)=y(i)*exp(5*j*w(i)); F(i)=abs(x(i)); an(i)=angle(x(i));end%直接求出傅立叶变换tt=linspace(-5,5,1000);ww=linspace(-5*pi,5*pi,1000);xx=sym(exp(-2*abs(tt)))yy=fourier(xx);FF=abs(yy);ann=0;%画图subplot(211);plot(w,F,r);hold on;ezplot(FF);title(幅频特性,红线表示近似值 );grid on;subplot(212);plot(w,an,r);hold on;ezplot(ann)title(相频特性，红线表示近似值 );grid on;结果如下：由于红线基本在水平线周围浮动，可得CTFT近似值与解析所得的大概相同，但存在误差。但由图可知，在较高频率上的近似不如较低频率上的好。(G)代码如下：clear;clc;tau=0.01;T=10;t=[0:tau:T-tau];N=length(t)w=-(pi/tau)+(0:N-1)*(2