华中科技大学《复变函数与积分变换》积分变换实验课讲义.pptVIP

积分变换实验课 Matlab实现 矩形脉冲信号的傅氏变换 矩形脉冲信号 （例8.2） 频率突变信号 正弦叠加信号 观察正弦叠加信号：频率在 和 的位置对应的幅值特别高。因此标记出这两个频谱峰值对应的频率分量，正好可以验证信号的频率成份。 观察频率突变信号和正弦叠加信号的频谱图可以看出它们虽然有不同的时间过程，但是频谱图却相似，都是在 和 处有突起 脉冲信号 脉冲信号 傅里叶变换的缺陷 时变频率信号在音乐，地震信号，雷达回波，非平稳物理过程中是很常见的，人们希望知道信号在突变时刻对应的频率成份 但是傅里叶变换的积分作用平滑了突变成份，即任一频率点的谱值是由信号在整个时间域上的贡献决定的；反之，信号在某一时刻的状态也是由频谱在整个频率域上的贡献决定的 一个著名的例子就是 函数，时间上的点脉冲具有在频率域上正负无限延伸的均匀频谱 频率突变信号的小波变换 矩形脉冲信号的傅氏变换 编写信号 M=8; tend=1; T=10; N=2^M; dt=T/N; n=0:N-1; t=n*dt; w=zeros(size(t,2),1); Tow=find((tend-t)0); w(Tow,1)=ones(length(Tow),1); plot(t,w,b,LineWidth,2.5);title(矩形脉冲信号);xlabel(t ---); 在命令窗口输入cftbyfft(w,t);得到该信号的傅氏变换频谱图 或者在命令窗口输入exp8_2_;得到该信号的傅氏变换频谱图与理论值的对比图。注意，当M的取值过小时，就会出现混跌引起图形的差别。 频率突变信号 编写信号并绘制图形 clear;M=8;N=2^M; t=linspace(-10,10,N); s1=find(t.0); x(s1)=cos(2*pi*6*t(s1)); s2=find(t=.0); x(s2)=cos(2*pi*3*t(s2)); plot(t,x); title(频率突变信号);xlabel(t---); 保存在sinbreak.mat中 在命令窗口输入cftbyfft(x,t);得到傅氏变换频谱图 正弦叠加信号 编写信号并绘制图形 clear; M=8; N=2^M; t=linspace(-10,10,N); x1=sin(2*pi*6*t); x2=sin(2*pi*3*t); xx=x1+x2; plot(t,xx); title(正弦叠加信号);xlabel(t ---); 在命令窗口输入cftbyfft(xx,t); 得到该信号的傅氏变换频谱图 脉冲信号1 编写信号 clear;M=10; T=10; N=2^M; dt=T/N; n=0:N-1; t=n*dt; w=zeros(size(t)); w(1:6)=100; plot(t,w,b,LineWidth,2.5);title(The Delta Function);xlabel(t---); 在命令窗口输入cftbyfft(w,t);得到该信号的傅氏变换频谱图 脉冲信号2 编写信号 clear;M=10; T=10; N=2^M; dt=T/N; n=0:N-1; t=n*dt; w=zeros(size(t)); w(100:105)=100; plot(t,w,b,LineWidth,2.5);title(The Delta Function);xlabel(t---); 在命令窗口输入cftbyfft(w,t);得到该信号的傅氏变换频谱图 几种信号 山形函数 山形函数 clear; M=10; N=2^M; t=linspace(-10,10,N); mtn=zeros(size(t)); s1=find(t=-1t0); mtn(s1)=1+t(s1); s2=find(t0t=1); mtn(s2)=1-t(s2); plot(t,mtn); title(Mountian Function);xlabel(t ---); 余弦信号 clear; M=10; N=2^M; t=linspace(-10,10,N); xcos=cos(3*t); 高斯函数 clear; M=10; N=2^M; t=linspace(-10,10,N); a=1/4; g=exp(-a*t.^2); 注：高斯函数的傅氏变换仍为高斯函数，即 阶跃信号 clear; M=10; N=2^M; t=li