xiao完整版完整版.docx

实验内容： 如下图所示锯齿波信号，分别取一个周期的抽样数据x1(t),0 = t 1和五个周期的数据x(t), 0 = t 5,计算其傅里叶变换X1（w）和X（w），比较有何不同并解释原因。 编写程序如下： T1=1; %单个周期时域范围 N1=10000; %时域抽样点数 t1=linspace(0,T1-T1/N1,N1); %生成抽样时间点 f1=1-2*t1; %生成抽样函数值 OMG=32*pi; %频域范围 K1=100; %频域抽样点数 omg=linspace(-OMG/2,OMG/2-OMG/K1,K1); %生成抽样频率点 X1=T1/N1*exp(-j*kron(omg,t1.))*f1; %傅里叶正变换求解傅里叶系数 fs1=OMG/2/pi/K1*exp(j*kron(t1,omg.))*X1; %傅里叶逆变换还原时域函数 T2=5; %五个周期时域范围 N2=10000; %时域抽样点数 t2=linspace(0,T2-T2/N2,N2); %生成抽样时间点 fs2=0*t2; f2=sawtooth(t2*2*pi,0); %生成五个周期的锯齿波 X2=T2/N2*exp(-j*kron(omg,t2.))*f2; %傅里叶正变换求解傅里叶系数 fs2=fs2+OMG/2/pi/K1*exp(j*kron(t2,omg.))*X2; %傅里叶逆变换还原时域函数 figure; %生成一个2*2矩阵子图 subplot(2,2,1); plot(omg,abs(X1),r); %一个周期时的频谱图 xlabel(Frequency),ylabel(Amplitude) title(单个锯齿周期幅频特性曲线); subplot(2,2,2); plot(t1,fs1,r); %还原的时域函数 xlabel(Time),ylabel(Amplitude) title(Function after recovered); subplot(2,2,3); plot(omg,abs(X2),r); %五个周期时的频谱图 xlabel(Frequency),ylabel(Amplitude) title(五个锯齿周期幅频特性曲线); subplot(2,2,4); plot(t2,fs2,r); %还原的时域函数 xlabel(Time),ylabel(Function after recovered) title(Function after recovered); 相关曲线： 请编写函数F=fsana(t,f,,N),计算周期信号f的前N个指数形式的傅里叶级数系数，t表示f对应的抽样时间（均为一个周期）;再编写函数f=fssyn（F，t）,由傅里叶级数系数F合成抽样时间t对应的函数。设计信号这两个函数是否正确。 3：利用fsana和fssyn计算练习题1中x1(t)的前10个傅里叶级数系数Fi，0=i10，请用这些系数合成周期为0.5的锯齿波y(t),0=t2并绘出图形。 编写程序如下： 编写fsana函数： %函数功能：定义一个功能函数来实现周期函数的傅里叶正变换 t表示f表示的抽样时间 f为一个周期内抽样函数值 N为前几个指数形式的傅立叶级数系数 function F=fsana(t,f,N) omg1=2*pi/(max(t)-min(t)); k=[-N:N]; F=1/length(t)*exp(-j*kron(k*omg1,t.))*f; 编写fssyn函数： %函数功能：定义功能函数fssyn由傅立叶级数系数F合成时间t对应的函数 function f=fssyn(F,t)