小波应用范例
【例 1-10】已知序列试绘制及
其Fourier变换的幅值图,并用变换后的数值求解逆变换,其中采样频率为。
我们知道,假定采样频率为 (即的采样间隔进行采样),则该信号包含有和两种频率的波,其振幅均为1。
程序如下:
clf %清除图形框的内容
N=100;dt=1; %设置最大点数
n=0:N-1; t=n*dt; %给出时间序列
xn=cos(2*pi*0.24*t)+cos(2*pi*0.26*t); %给出原始信号的值序列
Xk=fft(xn,N); %对原始信号进行Fourier变换
magXk=abs(Xk);phaXk=angle(Xk); %求出Fourier变换的振幅和相位
subplot(2,2,1),plot(t,xn); xlabel('时间/s') %绘出原始信号
title('原始信号(N=100)');
xx=ifft(Xk,N); %Fourier逆变换
x=real(xx); %取变换后的实部,如做实验可以验证其虚部为零
subplot(2,2,2),plot(
1