HMM应用实例.rar_hmm建模_与小波包程序_隐马尔可夫_隐马尔可夫编码程序_马尔可夫

load('s3.txt'); s=s3(:,2) fs=10000;%设定采样频率 N=12800; n=0:N-1; t=n/fs; y=fft(s,N);%进行fft变换 mag=abs(y)/N;%求幅值 m=length(y); f=(0:m/2-1)'*fs/m;%进行对应的频率转换 figure(1); subplot(131) plot(f,mag(1:m/2));%做频谱图 xlabel('Frequency(Hz)'); ylabel('Amplitude'); load('t11.txt'); s=t11(:,2) fs=10000;%设定采样频率 N=12800;%有误，由于个数为2273个，应该改为2273 n=0:N-1; t=n/fs; y=fft(s,N);%进行fft变换 mag=abs(y)/N;%求幅值 m=length(y); f=(0:m/2-1)'*fs/m;%进行对应的频率转换 figure(1); subplot(132) plot(f,mag(1:m/2));%做频谱图 xlabel('Frequency(Hz)'); ylabel('Amplitude'); load('c1.txt'); s=c1(:,2) fs=10000;%设定采样频率 N=12800; n=0:N-1; t=n/fs; y=fft(s,N);%进行fft变换 mag=abs(y)/N;%求幅值 m=length(y); f=(0:m/2-1)'*fs/m;%进行对应的频率转换 figure(1); subplot(133) plot(f,mag(1:m/2));%做频谱图 xlabel('Frequency(Hz)'); ylabel('Amplitude');