傅里叶变换性质研究.docVIP

傅里叶变换的性质研究实验目的学习傅里叶变换的数值实现方法。学习程序优化方法实验内容学习MATLAB的三个知识点：程序优化技巧生成连续周期信号的方法自定义函数中的参数有效性检验涉及MATLAB函数熟悉下列函数的应用Fourier——傅里叶变换Ifourier——傅里叶逆变换Kron ——张量积（叉乘）Sqrt——平方根MATLAB提供了符号fourier和ifourier实现傅里叶变换和逆变换。 实验内容1、计算tu（t）和sin（t）的傅里叶变换。在matlab命令窗口中编写实验程序：syms t %定义符号tF1=fourier(t*heaviside(t)) %计算tu（t）的傅里叶变换F1F2=fourier(sin(t)） %计算sin（t）的傅里叶变换F2编写好程序后，按下enter键，执行程序得到如下结果：F1 = i*(pi*dirac(1,w)*w^2+i)/w^2 F2 = i*pi*(-dirac(w-1)+dirac(w+1))2 请绘制矩形脉冲F(t)={10其他 |t|1/2的波形（解 直接计算法 T=2;N=200;t=linspace(-T/2,T/2-T/N,N);f=0*t;f(t-1/2t1/2)=1;OMG=16*pi;K=100;omg=linspace(-OMG/2,OMG/2-OMG/K,K);F=0*omg;for k=1:K F(k)=F(k)+T/N*f(n)*exp(-j*omg(k)*t(n));endfs=0*t;for n=1:N fs(n)=fs(n)+OMG/2/pi/K*F(k)*exp(j*omg*t(n));endfigure;subplot(121);plot(t,f,t,fs);xlabel(t);ylabel(input);legend(f,fs);subplot(122);plot(omg,F);xlabel(omg);ylabel(output);legend(F,fs) 矢量计算法T=2;N=200;t=linspace(-T/2,T/2-T/N,N);f=0*t;f(t-1/2t1/2)=1;OMG=16*pi;K=100;omg=linspace(-OMG/2,OMG/2-OMG/K,K);F=0*omg;for k=1:K F(k)=T/N*exp(-j*omg(k)*t).*f;endfs=0*t;for n=1:N fs(n)=OMG/2/pi/K*exp(j*omg*t(n)).*F;endfigure;subplot(121);plot(t,f,t,fs);xlabel(t);ylabel(input);legend(f,fs);subplot(122);plot(omg,F);xlabel(omg);ylabel(output);legend(F,fs矩阵计算法T=2;N=200;t=linspace(-T/2,T/2-T/N,N);f=0*t;f(t-1/2t1/2)=1;OMG=16*pi;K=100;omg=linspace(-OMG/2,OMG/2-OMG/K,K);F=0*omg;for k=1:K F(k)=T/N*exp(-j*omg(k)*t).*f;endfs=0*t;for n=1:N fs(n)=OMG/2/pi/K*exp(j*omg*t(n)).*F;endU=exp(-j*kron(omg,t.));F=T/N*U*f;V=exp(j*kron(t,omg.));fs=OMG/2/pi/K*V*F;figure;subplot(121);plot(t,f,t,fs);xlabel(t);ylabel(input);legend(f,fs);subplot(122);plot(omg,F);xlabel(omg);ylabel(output);legend(F,fs)实验总结本实验学习了三种实现傅里叶变换和逆变换的方法，