300字范文 > 【数学信号处理】基于matlab数字信号频谱分析【含Matlab源码 1544期】

【数学信号处理】基于matlab数字信号频谱分析【含Matlab源码 1544期】

时间：2019-08-05 16:23:34

一、获取代码方式

获取代码方式1：

完整代码已上传我的资源：【数学信号处理】基于matlab数字信号频谱分析【含Matlab源码 1544期】

获取代码方式2：

通过订阅紫极神光博客付费专栏，凭支付凭证，私信博主，可获得此代码。

备注：订阅紫极神光博客付费专栏，可免费获得1份代码（有效期为订阅日起，三天内有效）；

二、部分源代码

%% 声音的采集clc;clear;close all;filename = '三狮.wav';[y,Fs] = audioread(filename); y = y(:,1); % 取单声道L = length(y); % 采样点T = 1/Fs; % 周期t = (0:L-1)*T; info = audioinfo(filename); % wav详细信息disp(info); % 输出音频信号采样频率，采样位数等信息sound(y,Fs); %以采样率 Fs 向扬声器发送音频信号 yi = 0;%% 绘图i = i+1;s = '原音频信号时域波形图';plot_time(t,y,s,i); % 绘制时域波形图i = i+1;s = '原音频信号频域波形图';plot_freq(Fs,L,y,s,i); % 绘制频域波形图%% 窗函数对频谱分析的影响%% 矩形窗win_r = boxcar(length(y));y_win_r = y.*win_r; % 时域乘积sound(y_win_r,Fs);i = i+1;s = '矩形窗时域波形图';i = i+1;s = '矩形窗频率域波形图';plot_freq(Fs,L,y_win_r,s,i); % 矩形窗 1/2win_r_2 = boxcar(floor(length(y)/2)); % 取语音长度1/2为窗口长度win_r_2 = [win_r_2;zeros(ceil(length(y)/2),1)]; % 对齐补0y_win_r_w = y.*win_r_2;sound(y_win_r_w,Fs);i = i+1;s = '矩形窗(N = L/2)时域波形图';plot_time(t,y_win_r_w,s,i);i = i+1;s = '矩形窗(N = L/2)频率域波形图';plot_freq(Fs,L,y_win_r_w,s,i); %% 汉宁窗win_h = hann(length(y));y_win_h = y.*win_h;sound(y_win_h,Fs);i = i+1;s = '汉宁窗时域波形图';plot_time(t,y_win_h,s,i);i = i+1;s = '汉宁窗频率域波形图';plot_freq(Fs,L,y_win_h,s,i);% 汉宁窗 1/2win_h_2 = hann(floor(length(y)/2));win_h_2 = [win_h_2;zeros(ceil(length(y)/2),1)]; % 对齐补0y_win_h_2 = y.*win_h_2;sound(y_win_h_2,Fs);i = i+1;s = '汉宁窗(N = L/2)时域波形图';plot_time(t,y_win_h_2,s,i);i = i+1;s = '汉宁窗(N = L/2)频率域波形图';plot_freq(Fs,L,y_win_h_2,s,i);%% 哈明窗win_hm = hamming(length(y));y_win_hm = y.*win_hm;sound(y_win_hm,Fs);i = i+1;s = '哈明窗时域波形图';plot_time(t,y_win_hm,s,i);i = i+1;s = '哈明窗频率域波形图';plot_freq(Fs,L,y_win_hm,s,i);% 哈明窗 1/2win_hm_2 = hamming(floor(length(y)/2));win_hm_2 = [win_hm_2;zeros(ceil(length(y)/2),1)]; % 对齐补0y_win_hm_2 = y.*win_hm_2;sound(y_win_hm_2,Fs);i = i+1;s = '哈明窗(N = L/2)频率域波形图';plot_freq(Fs,L,y_win_hm_2,s,i);%% 布莱克曼窗win_blk = blackman(length(y));y_win_blk = y.*win_blk;sound(y_win_blk,Fs);i = i+1;s = '布莱克曼窗时域波形图';plot_time(t,y_win_blk,s,i);i = i+1;s = '布莱克曼窗频率域波形图';plot_freq(Fs,L,y_win_blk,s,i);% 布莱克曼窗 1/2win_blk_2 = blackman(floor(length(y)/2));win_blk_2 = [win_blk_2;zeros(ceil(length(y)/2),1)]; % 对齐补0y_win_blk_2 = y.*win_blk_2;sound(y_win_blk_2,Fs);i = i+1;s = '布莱克曼窗(N = L/2)时域波形图';plot_time(t,y_win_blk_2,s,i);i = i+1;s = '布莱克曼窗(N = L/2)频率域波形图';plot_freq(Fs,L,y_win_blk_2,s,i);%% 截取2stime = 2;samples = [1,floor(L/2)];[y_cut,Fs_cut] = audioread(filename,samples);y_cut = y_cut(:,1); % 取单声道L_cut = length(y_cut); % 采样点T_cut = 1/Fs_cut; % 周期t_cut = (0:L_cut-1)*T_cut; sound(y_cut,Fs_cut);i = i+1;s = '音频信号截取2s时域波形图';plot_time(t_cut,y_cut,s,i); % 绘制时域波形图i = i+1;s = '音频信号截取2s频域波形图';plot_freq(Fs_cut,L_cut,y_cut,s,i); % 绘制频域波形图%% 降采样对频谱分析的影响Fs_res = Fs/2;y_res=resample(y,Fs_res,Fs);%信号降采样处理% y_res = downsample(y,2);sound(y_res,Fs_res);t_res = 1/(Fs_res) * (0:L_res-1);i = i+1;s = '降采样时域波形图';plot_time(t_res,y_res,s,i);i = i+1;s = '降采样频率域波形图';plot_freq(Fs_res,L_res,y_res,s,i);