时频转换分析_matlab源码_bidmc

%用fft对求信号的幅频和相频曲线 function fftfenxi clear;clc;clf; %N=2048; %fft默认计算的信号是从0开始的 fs=12000; load('bidmc01m.mat', 'val') y=val; N=length(y); deta=1/fs; m=0:N-1; f=m*fs/N; %可以查看课本就是这样定义横坐标频率范围的 %下面计算的Y就是x(t)的傅里叶变换数值 Y=fft(y,N); % z=2*abs(Y); z=sqrt(Y.*conj(Y))/fs;%conj(Y)求Y的共轭 figure(1) % plot(f(1:100),z(1:100)); plot(f,z); title('幅频曲线') xiangwei=angle(Y); figure(2) plot(f,xiangwei) title('相频曲线') figure(3) plot(1:N,y,'r') title('原始信号') figure(4) polar(angle(Y),abs(Y)) title('复数')