南昌丸学卖验掖告.docx

南昌丸学卖验掖告.docx南昌丸学卖验掖告学生姓名:_学号:6100209228班级:电了093班 实验类型:口验证□综合■设计口创新实验口期:2011・04・8实验成绩: 傅里叶变换(一)实验目的1、 掌握对不同的函数进行傅里叶变换的程序编写;2、 熟悉生成联系周期信号的方法;3、 练****matlab编程。(二)=fsana(t,f„N),计算周期信号f的前N个指数形式的傅立叶级数系数,t表示f对应的抽样时间(均为一个周期);再编写函数f=fssyn(F,t),由傅立叶级数系数F合成抽样吋间t对应的函数。设计信号验证这两个是否正确。定义F=fsana(t,f,N)。functionF=fsana(tzf,N)omgl=2*pi/(max(t)-min(t));k=[0:N]*;F=l/length(t)*exp(-j*kron(k*omglzt・1))*f定义f=fssyn(F’t)functionf=fssyn(Fzt)omgl=2*pi/(max(t)-min(t));N=floor(length(F)/2);k=[0:N];f=exp(j*kron(匸,k*omgl))*F;%—个周期时域范围%时域抽样点数运行所定义的函数Tl=2*pi;N1=300;t=linspace(O,T1・T1/N1,N1)';f=cos(t);subplot(2,2,1)plot(t,f);title「原函数TN=IO;FI=fsana(t,f,N);数subplot(2,2,2)stem(abs(Fl),'s*);%生成抽样时间点%生成抽样函数值%调用fsana函数求解前N项傅立叶级数系%绘制离散的幅度曲线title「前N项傅立叶级数系数幅度曲线T;f2=fssyn(Fl,t); %调用fssyn函数求原时域函数subplot(2,2,3)plot(t,f2,fkf);title(,傅立叶逆变换后吋域函数,);运行结果原函数前N项傅立叶级数系数幅度曲线傅立叶逆变换后时域函数F=・・・--・-・--^xl(t)的钱10个傅立叶级数洗漱Fi,0《ivlO,(t),0《tv2,并绘出图形定义F=fsana(t,f,N)functionF=fSana(tAfzN)omgl=2*pi/(max(t)-min(t));k=[0:N]1;F=1/l