主頁 >  其他 > 實驗 9 音頻信號采集及處理

實驗 9 音頻信號采集及處理

2021-01-03 11:33:07 其他

實驗 9 音頻信號采集及處理

  • 代碼剖析
    • 聲音回放
    • 顯示音頻時域波形
    • 顯示音頻頻域波形
    • 聲音主要頻率1
    • 聲音主要頻率2
    • 合成聲音播放
    • 播放加噪信號
    • 加噪之后的時域和頻譜圖
    • 顯示低通濾波器特性
    • 播放濾波后的聲音
    • 濾波后的時域和頻譜圖
    • 關閉所有figure
  • 鳴謝
  • GUI代碼
  • 代碼及ding.wav下載

代碼剖析

聲音回放

[voice,fs]=audioread('ding.wav'); 
sound(voice,fs); %聲音回放

顯示音頻時域波形

figure(1);
[voice,fs]=audioread('ding.wav'); 
plot(voice); %繪出時域波
xlabel('t'); ylabel('A'); %坐標名稱
title(' 初始音頻信號時域波形 ','FontName','宋體');
grid on;

顯示音頻頻域波形

figure(2);
[voice,fs]=audioread('ding.wav'); 
n=length(voice); %計算長度
voice1=fft(voice,n);%快速傅里葉變換
plot(abs(voice1)); %繪出原始音頻信號頻譜
title(' 初始音頻信號頻域波形 ','FontName','宋體');
xlabel('f'); ylabel('A');
grid on;

聲音主要頻率1

[signal,fs] = audioread('ding.wav' ); 
len_signal = length( signal );                                      %信號的長度
T = 1/fs;                                                           %采樣周期
t = T:T:len_signal/fs;                                              %時間t序列
fft_N = 2^nextpow2( len_signal );                              %計算fft變換的點數
SIGNAL = fft( signal,fft_N );                                  %快速傅里葉變換
SIGNAL_AMP = abs( SIGNAL );         %只做幅頻分析,對fft變換后序列取模
SIGNAL_AMP_HALF = SIGNAL_AMP(1:(fft_N/2));%由于對稱性,只取前一半
[SIGNAL_pks,locs] = findpeaks( SIGNAL_AMP_HALF );         %尋找頻譜峰值點
sort_pks = sort( SIGNAL_pks,'descend' );                  %對頻譜峰值點進行排序
sort_loc = find( SIGNAL_pks>=sort_pks(2) ); %找到峰值最大的2個點排序后在pks中的位置
n1 = locs( sort_loc(1) );                                 %最高峰值點在序列中的位置
n2 = locs( sort_loc(2) );         %第二高峰值點在序列中的位置

Y1 = SIGNAL_AMP_HALF(n1);        %最高峰值點在頻譜中的幅度
Y2 = SIGNAL_AMP_HALF(n2);         %第二高峰值點在頻譜中的幅度

f1 = (n1-1)*fs/fft_N                                     %最高峰值點頻率
y1 = (Y1/(Y1+Y2))*cos(2*pi*f1*t);       %第一個頻率信號
figure(9);
plot( t,y1 );                    %繪制原信號時域波形
axis( [T len_signal*T -0.08 0.08] );                              %限定顯示范圍
title( '分解信號時域波形圖1' ,'FontName','宋體');                                      %繪制標題
xlabel('時間(s)','FontName','宋體'); ylabel('幅值','FontName','宋體');  

聲音主要頻率2

[signal,fs] = audioread('ding.wav' ); 

len_signal = length( signal );                                      %信號的長度
T = 1/fs;                                                           %采樣周期
t = T:T:len_signal/fs;                                              %時間t序列
fft_N = 2^nextpow2( len_signal );                              %計算fft變換的點數
SIGNAL = fft( signal,fft_N );                                  %快速傅里葉變換
SIGNAL_AMP = abs( SIGNAL );         %只做幅頻分析,對fft變換后序列取模
SIGNAL_AMP_HALF = SIGNAL_AMP(1:(fft_N/2));%由于對稱性,只取前一半
[SIGNAL_pks,locs] = findpeaks( SIGNAL_AMP_HALF );         %尋找頻譜峰值點
sort_pks = sort( SIGNAL_pks,'descend' );                  %對頻譜峰值點進行排序
sort_loc = find( SIGNAL_pks>=sort_pks(2) ); %找到峰值最大的2個點排序后在pks中的位置
n1 = locs( sort_loc(1) );                                 %最高峰值點在序列中的位置
n2 = locs( sort_loc(2) );         %第二高峰值點在序列中的位置
Y1 = SIGNAL_AMP_HALF(n1);        %最高峰值點在頻譜中的幅度
Y2 = SIGNAL_AMP_HALF(n2);         %第二高峰值點在頻譜中的幅度
f1 = (n1-1)*fs/fft_N;                                     %最高峰值點頻率
f2 = (n2-1)*fs/fft_N                                     %第二高峰值點頻率
y2 = (Y2/(Y1+Y2))*cos(2*pi*f2*t);      %第二個頻率信號
figure(8)
plot( t,y2 );                    %繪制原信號時域波形
axis( [T len_signal*T -0.08 0.08] );                              %限定顯示范圍
title( '分解信號時域波形圖2' ,'FontName','宋體');                                      %繪制標題
xlabel('時間(s)','FontName','宋體'); ylabel('幅值','FontName','宋體');  

合成聲音播放

[signal,fs] = audioread('ding.wav' ); 

len_signal = length( signal );                                      %信號的長度
T = 1/fs;                                                           %采樣周期
t = T:T:len_signal/fs;                                              %時間t序列
fft_N = 2^nextpow2( len_signal );                              %計算fft變換的點數
SIGNAL = fft( signal,fft_N );                                  %快速傅里葉變換
SIGNAL_AMP = abs( SIGNAL );         %只做幅頻分析,對fft變換后序列取模
SIGNAL_AMP_HALF = SIGNAL_AMP(1:(fft_N/2));%由于對稱性,只取前一半
[SIGNAL_pks,locs] = findpeaks( SIGNAL_AMP_HALF );         %尋找頻譜峰值點
sort_pks = sort( SIGNAL_pks,'descend' );                  %對頻譜峰值點進行排序
sort_loc = find( SIGNAL_pks>=sort_pks(2) ); %找到峰值最大的2個點排序后在pks中的位置
n1 = locs( sort_loc(1) );                                 %最高峰值點在序列中的位置
n2 = locs( sort_loc(2) );         %第二高峰值點在序列中的位置
Y1 = SIGNAL_AMP_HALF(n1);        %最高峰值點在頻譜中的幅度
Y2 = SIGNAL_AMP_HALF(n2);         %第二高峰值點在頻譜中的幅度
f1 = (n1-1)*fs/fft_N;                                     %最高峰值點頻率
f2 = (n2-1)*fs/fft_N;                                     %第二高峰值點頻率

y1 = (Y1/(Y1+Y2))*cos(2*pi*f1*t);       %第一個頻率信號
y2 = (Y2/(Y1+Y2))*cos(2*pi*f2*t);      %第二個頻率信號
A = max( signal );      %原信號的最大幅度
hecheng = (y1+y2)/max( y1+y2 );        %合成信號歸一化
sound( hecheng,fs );                                      %合成音頻信號播放

播放加噪信號

[voice,fs]=audioread('ding.wav'); 
n=length(voice); %計算長度
t=0:1/fs:(n-1)/fs;
noise=0.03*sin(2*pi*100000*t'); %100kHz 正弦波噪聲
s=voice+noise; %加噪后的音頻信號
sound(s,fs); %播放加噪的語音

加噪之后的時域和頻譜圖

[voice,fs]=audioread('ding.wav'); 
n=length(voice); %計算長度
t=0:1/fs:(n-1)/fs;
noise=0.03*sin(2*pi*100000*t'); %100kHz 正弦波噪聲

s=voice+noise; %加噪后的音頻信號
figure(3);
n=length(s); %計算長度
voice1=fft(s,n);%快速傅里葉變換
plot(abs(voice1)); %繪出原始音頻信號頻譜
title(' 加噪音頻信號頻域波形 ','FontName','宋體');
xlabel('f'); ylabel('A');
grid on;

figure(6);
plot(s); %繪出原始音頻信號頻譜
title(' 加噪音頻信號時域波形 ','FontName','宋體');
xlabel('f'); ylabel('A');
grid on;

顯示低通濾波器特性

[voice,fs]=audioread('ding.wav'); 
n=length(voice); %計算長度
t=0:1/fs:(n-1)/fs;
noise=0.03*sin(2*pi*100000*t'); %100kHz 正弦波噪聲

s=voice+noise; %加噪后的音頻信號

rp=0.1;     %通帶最大衰減
rs=60;      %阻帶最小衰減
Ft=1000;
Fp=60;      %通帶截止頻率fp
Fs=100;      %阻帶截止頻率fs
wp=2*pi*Fp/Ft;%WP是通帶截止頻率
ws=2*pi*Fs/Ft; %WS是阻帶截止頻率
[n,wn]=ellipord(wp,ws,rp,rs,'s');% 呼叫 ellipod 計算橢圓 DF階數 N和通帶截止頻率 wp,
[bz,az]=ellip(n,rp,rs,wn);% 呼叫 ellip 計算橢圓帶通 DF系統函式系數向量B和 A
[h,w]=freqz(bz,az);

figure(4);
plot(w*fs/(2*pi),abs(h)); %繪制 IIR 低通濾波器特性曲線
title('IIR 低通濾波器特性曲線 ','FontName','宋體');
grid on;
xlabel('f/Hz');ylabel(' 幅度')

播放濾波后的聲音

[voice,fs]=audioread('ding.wav'); 
n=length(voice); %計算長度
t=0:1/fs:(n-1)/fs;
noise=0.03*sin(2*pi*100000*t'); %100kHz 正弦波噪聲

s=voice+noise; %加噪后的音頻信號

rp=0.1;     %通帶最大衰減
rs=60;      %阻帶最小衰減
Ft=1000;
Fp=60;      %通帶截止頻率fp
Fs=100;      %阻帶截止頻率fs
wp=2*pi*Fp/Ft;%WP是通帶截止頻率
ws=2*pi*Fs/Ft; %WS是阻帶截止頻率
[n,wn]=ellipord(wp,ws,rp,rs,'s');% 呼叫 ellipod 計算橢圓 DF階數 N和通帶截止頻率 wp,
[bz,az]=ellip(n,rp,rs,wn);% 呼叫 ellip 計算橢圓帶通 DF系統函式系數向量B和 A
[h,w]=freqz(bz,az);
z=filter(bz,az,s); %濾波
sound(z,fs); %播放加噪的語音

濾波后的時域和頻譜圖

[voice,fs]=audioread('ding.wav'); 
n=length(voice); %計算長度
t=0:1/fs:(n-1)/fs;
noise=0.03*sin(2*pi*100000*t'); %100kHz 正弦波噪聲

s=voice+noise; %加噪后的音頻信號

rp=0.1;     %通帶最大衰減
rs=60;      %阻帶最小衰減
Ft=1000;
Fp=60;      %通帶截止頻率fp
Fs=100;      %阻帶截止頻率fs
wp=2*pi*Fp/Ft;%WP是通帶截止頻率
ws=2*pi*Fs/Ft; %WS是阻帶截止頻率
[n,wn]=ellipord(wp,ws,rp,rs,'s');% 呼叫 ellipod 計算橢圓 DF階數 N和通帶截止頻率 wp,
[bz,az]=ellip(n,rp,rs,wn);% 呼叫 ellip 計算橢圓帶通 DF系統函式系數向量B和 A
[h,w]=freqz(bz,az);

figure(5);
z=filter(bz,az,s); %濾波

figure(5);
n=length(z); %計算長度
voice1=fft(z,n);%快速傅里葉變換
plot(abs(voice1)); %繪出原始音頻信號頻譜
title('濾波后音頻信號頻域波形','FontName','宋體');
xlabel('f'); ylabel('A');
grid on;

figure(7);
plot(z); %繪出原始音頻信號頻譜
title(' 濾波后音頻信號時域波形 ','FontName','宋體');
xlabel('f'); ylabel('A');
grid on;

關閉所有figure

clc;
close all;

鳴謝

感謝我的朋友趙思文軒、陳勇輝為我提供的想法和代碼支持

GUI代碼

classdef mydsp < matlab.apps.AppBase

    % Properties that correspond to app components
    properties (Access = public)
        UIFigure       matlab.ui.Figure
        dingwavButton  matlab.ui.control.Button
        Button_4       matlab.ui.control.Button
        Button_5       matlab.ui.control.Button
        Button_7       matlab.ui.control.Button
        Button_8       matlab.ui.control.Button
        Button_9       matlab.ui.control.Button
        Button_10      matlab.ui.control.Button
        Button_11      matlab.ui.control.Button
        figureButton   matlab.ui.control.Button
        Button_12      matlab.ui.control.Button
        Button_13      matlab.ui.control.Button
        Button_14      matlab.ui.control.Button
    end

    
    properties (Access = public)
         % Description
    end
    

    % Callbacks that handle component events
    methods (Access = private)

        % Callback function
        function ButtonPushed(app, event)
            clc;
            close all;
            global R;
            R= audiorecorder( 44100, 16 ,1);
            record(R)
            R
        end

        % Callback function
        function Button_2Pushed(app, event)
            global R;
            pause(R);
            R
        end

        % Button pushed function: dingwavButton
        function dingwavButtonPushed(app, event)
            [voice,fs]=audioread('ding.wav'); 
            sound(voice,fs); %聲音回放
        end

        % Button pushed function: Button_4
        function Button_4Pushed(app, event)
            figure(1);
            [voice,fs]=audioread('ding.wav'); 
            plot(voice); %繪出時域波
            xlabel('t'); ylabel('A'); %坐標名稱
            title(' 初始音頻信號時域波形 ','FontName','宋體');
            grid on;
        end

        % Button pushed function: Button_5
        function Button_5Pushed(app, event)
            figure(2);
            [voice,fs]=audioread('ding.wav'); 
            n=length(voice); %計算長度
            voice1=fft(voice,n);%快速傅里葉變換
            plot(abs(voice1)); %繪出原始音頻信號頻譜
            title(' 初始音頻信號頻域波形 ','FontName','宋體');
            xlabel('f'); ylabel('A');
            grid on;
        end

        % Callback function
        function Button_6Pushed(app, event)
            wp=2*pi*20000;ws=2*pi*24000;Rp=0.1;As=60;
            [N,wc]=buttord(wp,ws,Rp,As,'s');
            [B,A]=butter(N,wc,'s');
            [voice,fs]=audioread('ding.wav'); 
            
%             n=length(voice); %計算長度
%             voice1=filter(B,A,voice);
%             voice2=fft(voice1,n);%快速傅里葉變換
%             plot(abs(voice2)); %繪出原始音頻信號頻譜
            
        end

        % Button pushed function: Button_7
        function Button_7Pushed(app, event)
            [voice,fs]=audioread('ding.wav'); 
            n=length(voice); %計算長度
            t=0:1/fs:(n-1)/fs;
            noise=0.03*sin(2*pi*100000*t'); %100kHz 正弦波噪聲
            s=voice+noise; %加噪后的音頻信號
            sound(s,fs); %播放加噪的語音
        end

        % Button pushed function: Button_8
        function Button_8Pushed(app, event)
            [voice,fs]=audioread('ding.wav'); 
            n=length(voice); %計算長度
            t=0:1/fs:(n-1)/fs;
            noise=0.03*sin(2*pi*100000*t'); %100kHz 正弦波噪聲
            
            s=voice+noise; %加噪后的音頻信號
            figure(3);
            n=length(s); %計算長度
            voice1=fft(s,n);%快速傅里葉變換
            plot(abs(voice1)); %繪出原始音頻信號頻譜
            title(' 加噪音頻信號頻域波形 ','FontName','宋體');
            xlabel('f'); ylabel('A');
            grid on;
            
            figure(6);
            plot(s); %繪出原始音頻信號頻譜
            title(' 加噪音頻信號時域波形 ','FontName','宋體');
            xlabel('f'); ylabel('A');
            grid on;
        end

        % Button pushed function: Button_9
        function Button_9Pushed(app, event)
            [voice,fs]=audioread('ding.wav'); 
            n=length(voice); %計算長度
            t=0:1/fs:(n-1)/fs;
            noise=0.03*sin(2*pi*100000*t'); %100kHz 正弦波噪聲
            
            s=voice+noise; %加噪后的音頻信號
            
            rp=0.1;     %通帶最大衰減
            rs=60;      %阻帶最小衰減
            Ft=1000;
            Fp=60;      %通帶截止頻率fp
            Fs=100;      %阻帶截止頻率fs
            wp=2*pi*Fp/Ft;%WP是通帶截止頻率
            ws=2*pi*Fs/Ft; %WS是阻帶截止頻率
            [n,wn]=ellipord(wp,ws,rp,rs,'s');% 呼叫 ellipod 計算橢圓 DF階數 N和通帶截止頻率 wp,
            [bz,az]=ellip(n,rp,rs,wn);% 呼叫 ellip 計算橢圓帶通 DF系統函式系數向量B和 A
            [h,w]=freqz(bz,az);
            
            figure(4);
            plot(w*fs/(2*pi),abs(h)); %繪制 IIR 低通濾波器特性曲線
            title('IIR 低通濾波器特性曲線 ','FontName','宋體');
            grid on;
            xlabel('f/Hz');ylabel(' 幅度')
        end

        % Button pushed function: Button_10
        function Button_10Pushed(app, event)
            [voice,fs]=audioread('ding.wav'); 
            n=length(voice); %計算長度
            t=0:1/fs:(n-1)/fs;
            noise=0.03*sin(2*pi*100000*t'); %100kHz 正弦波噪聲
            
            s=voice+noise; %加噪后的音頻信號
            
            rp=0.1;     %通帶最大衰減
            rs=60;      %阻帶最小衰減
            Ft=1000;
            Fp=60;      %通帶截止頻率fp
            Fs=100;      %阻帶截止頻率fs
            wp=2*pi*Fp/Ft;%WP是通帶截止頻率
            ws=2*pi*Fs/Ft; %WS是阻帶截止頻率
            [n,wn]=ellipord(wp,ws,rp,rs,'s');% 呼叫 ellipod 計算橢圓 DF階數 N和通帶截止頻率 wp,
            [bz,az]=ellip(n,rp,rs,wn);% 呼叫 ellip 計算橢圓帶通 DF系統函式系數向量B和 A
            [h,w]=freqz(bz,az);
            z=filter(bz,az,s); %濾波
            sound(z,fs); %播放加噪的語音
        end

        % Button pushed function: Button_11
        function Button_11Pushed(app, event)
            [voice,fs]=audioread('ding.wav'); 
            n=length(voice); %計算長度
            t=0:1/fs:(n-1)/fs;
            noise=0.03*sin(2*pi*100000*t'); %100kHz 正弦波噪聲
            
            s=voice+noise; %加噪后的音頻信號
            
            rp=0.1;     %通帶最大衰減
            rs=60;      %阻帶最小衰減
            Ft=1000;
            Fp=60;      %通帶截止頻率fp
            Fs=100;      %阻帶截止頻率fs
            wp=2*pi*Fp/Ft;%WP是通帶截止頻率
            ws=2*pi*Fs/Ft; %WS是阻帶截止頻率
            [n,wn]=ellipord(wp,ws,rp,rs,'s');% 呼叫 ellipod 計算橢圓 DF階數 N和通帶截止頻率 wp,
            [bz,az]=ellip(n,rp,rs,wn);% 呼叫 ellip 計算橢圓帶通 DF系統函式系數向量B和 A
            [h,w]=freqz(bz,az);
            
            figure(5);
            z=filter(bz,az,s); %濾波

            figure(5);
            n=length(z); %計算長度
            voice1=fft(z,n);%快速傅里葉變換
            plot(abs(voice1)); %繪出原始音頻信號頻譜
            title('濾波后音頻信號頻域波形','FontName','宋體');
            xlabel('f'); ylabel('A');
            grid on;
            
            figure(7);
            plot(z); %繪出原始音頻信號頻譜
            title(' 濾波后音頻信號時域波形 ','FontName','宋體');
            xlabel('f'); ylabel('A');
            grid on;
        end

        % Button pushed function: figureButton
        function figureButtonPushed(app, event)
            clc;
            close all;
        end

        % Button pushed function: Button_12
        function Button_12Pushed(app, event)
            [signal,fs] = audioread('ding.wav' ); 
            len_signal = length( signal );                                      %信號的長度
            T = 1/fs;                                                           %采樣周期
            t = T:T:len_signal/fs;                                              %時間t序列
            fft_N = 2^nextpow2( len_signal );                              %計算fft變換的點數
            SIGNAL = fft( signal,fft_N );                                  %快速傅里葉變換
            SIGNAL_AMP = abs( SIGNAL );         %只做幅頻分析,對fft變換后序列取模
            SIGNAL_AMP_HALF = SIGNAL_AMP(1:(fft_N/2));%由于對稱性,只取前一半
            [SIGNAL_pks,locs] = findpeaks( SIGNAL_AMP_HALF );         %尋找頻譜峰值點
            sort_pks = sort( SIGNAL_pks,'descend' );                  %對頻譜峰值點進行排序
            sort_loc = find( SIGNAL_pks>=sort_pks(2) ); %找到峰值最大的2個點排序后在pks中的位置
            n1 = locs( sort_loc(1) );                                 %最高峰值點在序列中的位置
            n2 = locs( sort_loc(2) );         %第二高峰值點在序列中的位置
            
            Y1 = SIGNAL_AMP_HALF(n1);        %最高峰值點在頻譜中的幅度
            Y2 = SIGNAL_AMP_HALF(n2);         %第二高峰值點在頻譜中的幅度
            
            f1 = (n1-1)*fs/fft_N                                     %最高峰值點頻率
            y1 = (Y1/(Y1+Y2))*cos(2*pi*f1*t);       %第一個頻率信號
            figure(9);
            plot( t,y1 );                    %繪制原信號時域波形
            axis( [T len_signal*T -0.08 0.08] );                              %限定顯示范圍
            title( '分解信號時域波形圖1' ,'FontName','宋體');                                      %繪制標題
            xlabel('時間(s)','FontName','宋體'); ylabel('幅值','FontName','宋體');                                %繪制橫縱坐標的標簽
        end

        % Button pushed function: Button_14
        function Button_14Pushed(app, event)
            [signal,fs] = audioread('ding.wav' ); 
            
            len_signal = length( signal );                                      %信號的長度
            T = 1/fs;                                                           %采樣周期
            t = T:T:len_signal/fs;                                              %時間t序列
            fft_N = 2^nextpow2( len_signal );                              %計算fft變換的點數
            SIGNAL = fft( signal,fft_N );                                  %快速傅里葉變換
            SIGNAL_AMP = abs( SIGNAL );         %只做幅頻分析,對fft變換后序列取模
            SIGNAL_AMP_HALF = SIGNAL_AMP(1:(fft_N/2));%由于對稱性,只取前一半
            [SIGNAL_pks,locs] = findpeaks( SIGNAL_AMP_HALF );         %尋找頻譜峰值點
            sort_pks = sort( SIGNAL_pks,'descend' );                  %對頻譜峰值點進行排序
            sort_loc = find( SIGNAL_pks>=sort_pks(2) ); %找到峰值最大的2個點排序后在pks中的位置
            n1 = locs( sort_loc(1) );                                 %最高峰值點在序列中的位置
            n2 = locs( sort_loc(2) );         %第二高峰值點在序列中的位置
            Y1 = SIGNAL_AMP_HALF(n1);        %最高峰值點在頻譜中的幅度
            Y2 = SIGNAL_AMP_HALF(n2);         %第二高峰值點在頻譜中的幅度
            f1 = (n1-1)*fs/fft_N;                                     %最高峰值點頻率
            f2 = (n2-1)*fs/fft_N;                                     %第二高峰值點頻率
            
            y1 = (Y1/(Y1+Y2))*cos(2*pi*f1*t);       %第一個頻率信號
            y2 = (Y2/(Y1+Y2))*cos(2*pi*f2*t);      %第二個頻率信號
            A = max( signal );      %原信號的最大幅度
            hecheng = (y1+y2)/max( y1+y2 );        %合成信號歸一化
            sound( hecheng,fs );                                      %合成音頻信號播放
        end

        % Button pushed function: Button_13
        function Button_13Pushed(app, event)
             [signal,fs] = audioread('ding.wav' ); 
            
            len_signal = length( signal );                                      %信號的長度
            T = 1/fs;                                                           %采樣周期
            t = T:T:len_signal/fs;                                              %時間t序列
            fft_N = 2^nextpow2( len_signal );                              %計算fft變換的點數
            SIGNAL = fft( signal,fft_N );                                  %快速傅里葉變換
            SIGNAL_AMP = abs( SIGNAL );         %只做幅頻分析,對fft變換后序列取模
            SIGNAL_AMP_HALF = SIGNAL_AMP(1:(fft_N/2));%由于對稱性,只取前一半
            [SIGNAL_pks,locs] = findpeaks( SIGNAL_AMP_HALF );         %尋找頻譜峰值點
            sort_pks = sort( SIGNAL_pks,'descend' );                  %對頻譜峰值點進行排序
            sort_loc = find( SIGNAL_pks>=sort_pks(2) ); %找到峰值最大的2個點排序后在pks中的位置
            n1 = locs( sort_loc(1) );                                 %最高峰值點在序列中的位置
            n2 = locs( sort_loc(2) );         %第二高峰值點在序列中的位置
            Y1 = SIGNAL_AMP_HALF(n1);        %最高峰值點在頻譜中的幅度
            Y2 = SIGNAL_AMP_HALF(n2);         %第二高峰值點在頻譜中的幅度
            f1 = (n1-1)*fs/fft_N;                                     %最高峰值點頻率
            f2 = (n2-1)*fs/fft_N                                     %第二高峰值點頻率
            y2 = (Y2/(Y1+Y2))*cos(2*pi*f2*t);      %第二個頻率信號
            figure(8)
            plot( t,y2 );                    %繪制原信號時域波形
            axis( [T len_signal*T -0.08 0.08] );                              %限定顯示范圍
            title( '分解信號時域波形圖2' ,'FontName','宋體');                                      %繪制標題
            xlabel('時間(s)','FontName','宋體'); ylabel('幅值','FontName','宋體');    
        end
    end

    % Component initialization
    methods (Access = private)

        % Create UIFigure and components
        function createComponents(app)

            % Create UIFigure and hide until all components are created
            app.UIFigure = uifigure('Visible', 'off');
            app.UIFigure.Position = [100 100 640 480];
            app.UIFigure.Name = 'UI Figure';

            % Create dingwavButton
            app.dingwavButton = uibutton(app.UIFigure, 'push');
            app.dingwavButton.ButtonPushedFcn = createCallbackFcn(app, @dingwavButtonPushed, true);
            app.dingwavButton.Position = [76 430 100 24];
            app.dingwavButton.Text = '播放ding.wav';

            % Create Button_4
            app.Button_4 = uibutton(app.UIFigure, 'push');
            app.Button_4.ButtonPushedFcn = createCallbackFcn(app, @Button_4Pushed, true);
            app.Button_4.Position = [76 374 135 24];
            app.Button_4.Text = '顯示音頻時域波形';

            % Create Button_5
            app.Button_5 = uibutton(app.UIFigure, 'push');
            app.Button_5.ButtonPushedFcn = createCallbackFcn(app, @Button_5Pushed, true);
            app.Button_5.Position = [243 373 123 25];
            app.Button_5.Text = '顯示音頻頻域波形';

            % Create Button_7
            app.Button_7 = uibutton(app.UIFigure, 'push');
            app.Button_7.ButtonPushedFcn = createCallbackFcn(app, @Button_7Pushed, true);
            app.Button_7.Position = [76 229 100 24];
            app.Button_7.Text = '播放加噪信號';

            % Create Button_8
            app.Button_8 = uibutton(app.UIFigure, 'push');
            app.Button_8.ButtonPushedFcn = createCallbackFcn(app, @Button_8Pushed, true);
            app.Button_8.Position = [231 229 147 24];
            app.Button_8.Text = '加噪之后的時域和頻譜圖';

            % Create Button_9
            app.Button_9 = uibutton(app.UIFigure, 'push');
            app.Button_9.ButtonPushedFcn = createCallbackFcn(app, @Button_9Pushed, true);
            app.Button_9.Position = [76 179 123 24];
            app.Button_9.Text = '顯示低通濾波器特性';

            % Create Button_10
            app.Button_10 = uibutton(app.UIFigure, 'push');
            app.Button_10.ButtonPushedFcn = createCallbackFcn(app, @Button_10Pushed, true);
            app.Button_10.Position = [249 179 111 24];
            app.Button_10.Text = '播放濾波后的聲音';

            % Create Button_11
            app.Button_11 = uibutton(app.UIFigure, 'push');
            app.Button_11.ButtonPushedFcn = createCallbackFcn(app, @Button_11Pushed, true);
            app.Button_11.Position = [400 179 135 24];
            app.Button_11.Text = '濾波后的時域和頻譜圖';

            % Create figureButton
            app.figureButton = uibutton(app.UIFigure, 'push');
            app.figureButton.ButtonPushedFcn = createCallbackFcn(app, @figureButtonPushed, true);
            app.figureButton.Position = [260 96 100 24];
            app.figureButton.Text = '關閉所有figure';

            % Create Button_12
            app.Button_12 = uibutton(app.UIFigure, 'push');
            app.Button_12.ButtonPushedFcn = createCallbackFcn(app, @Button_12Pushed, true);
            app.Button_12.Position = [76 291 119 26];
            app.Button_12.Text = '聲音主要頻率1';

            % Create Button_13
            app.Button_13 = uibutton(app.UIFigure, 'push');
            app.Button_13.ButtonPushedFcn = createCallbackFcn(app, @Button_13Pushed, true);
            app.Button_13.Position = [241 291 119 26];
            app.Button_13.Text = '聲音主要頻率2';

            % Create Button_14
            app.Button_14 = uibutton(app.UIFigure, 'push');
            app.Button_14.ButtonPushedFcn = createCallbackFcn(app, @Button_14Pushed, true);
            app.Button_14.Position = [400 291 100 26];
            app.Button_14.Text = '合成聲音播放';

            % Show the figure after all components are created
            app.UIFigure.Visible = 'on';
        end
    end

    % App creation and deletion
    methods (Access = public)

        % Construct app
        function app = mydsp

            % Create UIFigure and components
            createComponents(app)

            % Register the app with App Designer
            registerApp(app, app.UIFigure)

            if nargout == 0
                clear app
            end
        end

        % Code that executes before app deletion
        function delete(app)

            % Delete UIFigure when app is deleted
            delete(app.UIFigure)
        end
    end
end

代碼及ding.wav下載

鏈接:https://pan.baidu.com/s/1UMIm8A_ThtiTailmHmYnSg
提取碼:6666
復制這段內容后打開百度網盤手機App,操作更方便哦

轉載請註明出處,本文鏈接:https://www.uj5u.com/qita/243942.html

標籤:其他

上一篇:2020我與csdn的相遇

下一篇:docker實戰部署SRS4.0流媒體服務器對接GB28181監控視頻控制臺

標籤雲
其他(157675) Python(38076) JavaScript(25376) Java(17977) C(15215) 區塊鏈(8255) C#(7972) AI(7469) 爪哇(7425) MySQL(7132) html(6777) 基礎類(6313) sql(6102) 熊猫(6058) PHP(5869) 数组(5741) R(5409) Linux(5327) 反应(5209) 腳本語言(PerlPython)(5129) 非技術區(4971) Android(4554) 数据框(4311) css(4259) 节点.js(4032) C語言(3288) json(3245) 列表(3129) 扑(3119) C++語言(3117) 安卓(2998) 打字稿(2995) VBA(2789) Java相關(2746) 疑難問題(2699) 细绳(2522) 單片機工控(2479) iOS(2429) ASP.NET(2402) MongoDB(2323) 麻木的(2285) 正则表达式(2254) 字典(2211) 循环(2198) 迅速(2185) 擅长(2169) 镖(2155) 功能(1967) .NET技术(1958) Web開發(1951) python-3.x(1918) HtmlCss(1915) 弹簧靴(1913) C++(1909) xml(1889) PostgreSQL(1872) .NETCore(1853) 谷歌表格(1846) Unity3D(1843) for循环(1842)

熱門瀏覽
  • 網閘典型架構簡述

    網閘架構一般分為兩種:三主機的三系統架構網閘和雙主機的2+1架構網閘。 三主機架構分別為內端機、外端機和仲裁機。三機無論從軟體和硬體上均各自獨立。首先從硬體上來看,三機都用各自獨立的主板、記憶體及存盤設備。從軟體上來看,三機有各自獨立的作業系統。這樣能達到完全的三機獨立。對于“2+1”系統,“2”分為 ......

    uj5u.com 2020-09-10 02:00:44 more
  • 如何從xshell上傳檔案到centos linux虛擬機里

    如何從xshell上傳檔案到centos linux虛擬機里及:虛擬機CentOs下執行 yum -y install lrzsz命令,出現錯誤:鏡像無法找到軟體包 前言 一、安裝lrzsz步驟 二、上傳檔案 三、遇到的問題及解決方案 總結 前言 提示:其實很簡單,往虛擬機上安裝一個上傳檔案的工具 ......

    uj5u.com 2020-09-10 02:00:47 more
  • 一、SQLMAP入門

    一、SQLMAP入門 1、判斷是否存在注入 sqlmap.py -u 網址/id=1 id=1不可缺少。當注入點后面的引數大于兩個時。需要加雙引號, sqlmap.py -u "網址/id=1&uid=1" 2、判斷文本中的請求是否存在注入 從文本中加載http請求,SQLMAP可以從一個文本檔案中 ......

    uj5u.com 2020-09-10 02:00:50 more
  • Metasploit 簡單使用教程

    metasploit 簡單使用教程 浩先生, 2020-08-28 16:18:25 分類專欄: kail 網路安全 linux 文章標簽: linux資訊安全 編輯 著作權 metasploit 使用教程 前言 一、Metasploit是什么? 二、準備作業 三、具體步驟 前言 Msfconsole ......

    uj5u.com 2020-09-10 02:00:53 more
  • 游戲逆向之驅動層與用戶層通訊

    驅動層代碼: #pragma once #include <ntifs.h> #define add_code CTL_CODE(FILE_DEVICE_UNKNOWN,0x800,METHOD_BUFFERED,FILE_ANY_ACCESS) /* 更多游戲逆向視頻www.yxfzedu.com ......

    uj5u.com 2020-09-10 02:00:56 more
  • 北斗電力時鐘(北斗授時服務器)讓網路資料更精準

    北斗電力時鐘(北斗授時服務器)讓網路資料更精準 北斗電力時鐘(北斗授時服務器)讓網路資料更精準 京準電子科技官微——ahjzsz 近幾年,資訊技術的得了快速發展,互聯網在逐漸普及,其在人們生活和生產中都得到了廣泛應用,并且取得了不錯的應用效果。計算機網路資訊在電力系統中的應用,一方面使電力系統的運行 ......

    uj5u.com 2020-09-10 02:01:03 more
  • 【CTF】CTFHub 技能樹 彩蛋 writeup

    ?碎碎念 CTFHub:https://www.ctfhub.com/ 筆者入門CTF時時剛開始刷的是bugku的舊平臺,后來才有了CTFHub。 感覺不論是網頁UI設計,還是題目質量,賽事跟蹤,工具軟體都做得很不錯。 而且因為獨到的金幣制度的確讓人有一種想去刷題賺金幣的感覺。 個人還是非常喜歡這個 ......

    uj5u.com 2020-09-10 02:04:05 more
  • 02windows基礎操作

    我學到了一下幾點 Windows系統目錄結構與滲透的作用 常見Windows的服務詳解 Windows埠詳解 常用的Windows注冊表詳解 hacker DOS命令詳解(net user / type /md /rd/ dir /cd /net use copy、批處理 等) 利用dos命令制作 ......

    uj5u.com 2020-09-10 02:04:18 more
  • 03.Linux基礎操作

    我學到了以下幾點 01Linux系統介紹02系統安裝,密碼啊破解03Linux常用命令04LAMP 01LINUX windows: win03 8 12 16 19 配置不繁瑣 Linux:redhat,centos(紅帽社區版),Ubuntu server,suse unix:金融機構,證券,銀 ......

    uj5u.com 2020-09-10 02:04:30 more
  • 05HTML

    01HTML介紹 02頭部標簽講解03基礎標簽講解04表單標簽講解 HTML前段語言 js1.了解代碼2.根據代碼 懂得挖掘漏洞 (POST注入/XSS漏洞上傳)3.黑帽seo 白帽seo 客戶網站被黑帽植入劫持代碼如何處理4.熟悉html表單 <html><head><title>TDK標題,描述 ......

    uj5u.com 2020-09-10 02:04:36 more
最新发布
  • 2023年最新微信小程式抓包教程

    01 開門見山 隔一個月發一篇文章,不過分。 首先回顧一下《微信系結手機號資料庫被脫庫事件》,我也是第一時間得知了這個訊息,然后跟蹤了整件事情的經過。下面是這起事件的相關截圖以及近日流出的一萬條資料樣本: 個人認為這件事也沒什么,還不如關注一下之前45億快遞資料查詢渠道疑似在近日復活的訊息。 訊息是 ......

    uj5u.com 2023-04-20 08:48:24 more
  • web3 產品介紹:metamask 錢包 使用最多的瀏覽器插件錢包

    Metamask錢包是一種基于區塊鏈技術的數字貨幣錢包,它允許用戶在安全、便捷的環境下管理自己的加密資產。Metamask錢包是以太坊生態系統中最流行的錢包之一,它具有易于使用、安全性高和功能強大等優點。 本文將詳細介紹Metamask錢包的功能和使用方法。 一、 Metamask錢包的功能 數字資 ......

    uj5u.com 2023-04-20 08:47:46 more
  • vulnhub_Earth

    前言 靶機地址->>>vulnhub_Earth 攻擊機ip:192.168.20.121 靶機ip:192.168.20.122 參考文章 https://www.cnblogs.com/Jing-X/archive/2022/04/03/16097695.html https://www.cnb ......

    uj5u.com 2023-04-20 07:46:20 more
  • 從4k到42k,軟體測驗工程師的漲薪史,給我看哭了

    清明節一過,盲猜大家已經無心上班,在數著日子準備過五一,但一想到銀行卡里的余額……瞬間心情就不美麗了。最近,2023年高校畢業生就業調查顯示,本科畢業月平均起薪為5825元。調查一出,便有很多同學表示自己又被平均了。看著這一資料,不免讓人想到前不久中國青年報的一項調查:近六成大學生認為畢業10年內會 ......

    uj5u.com 2023-04-20 07:44:00 more
  • 最新版本 Stable Diffusion 開源 AI 繪畫工具之中文自動提詞篇

    🎈 標簽生成器 由于輸入正向提示詞 prompt 和反向提示詞 negative prompt 都是使用英文,所以對學習母語的我們非常不友好 使用網址:https://tinygeeker.github.io/p/ai-prompt-generator 這個網址是為了讓大家在使用 AI 繪畫的時候 ......

    uj5u.com 2023-04-20 07:43:36 more
  • 漫談前端自動化測驗演進之路及測驗工具分析

    隨著前端技術的不斷發展和應用程式的日益復雜,前端自動化測驗也在不斷演進。隨著 Web 應用程式變得越來越復雜,自動化測驗的需求也越來越高。如今,自動化測驗已經成為 Web 應用程式開發程序中不可或缺的一部分,它們可以幫助開發人員更快地發現和修復錯誤,提高應用程式的性能和可靠性。 ......

    uj5u.com 2023-04-20 07:43:16 more
  • CANN開發實踐:4個DVPP記憶體問題的典型案例解讀

    摘要:由于DVPP媒體資料處理功能對存放輸入、輸出資料的記憶體有更高的要求(例如,記憶體首地址128位元組對齊),因此需呼叫專用的記憶體申請介面,那么本期就分享幾個關于DVPP記憶體問題的典型案例,并給出原因分析及解決方法。 本文分享自華為云社區《FAQ_DVPP記憶體問題案例》,作者:昇騰CANN。 DVPP ......

    uj5u.com 2023-04-20 07:43:03 more
  • msf學習

    msf學習 以kali自帶的msf為例 一、msf核心模塊與功能 msf模塊都放在/usr/share/metasploit-framework/modules目錄下 1、auxiliary 輔助模塊,輔助滲透(埠掃描、登錄密碼爆破、漏洞驗證等) 2、encoders 編碼器模塊,主要包含各種編碼 ......

    uj5u.com 2023-04-20 07:42:59 more
  • Halcon軟體安裝與界面簡介

    1. 下載Halcon17版本到到本地 2. 雙擊安裝包后 3. 步驟如下 1.2 Halcon軟體安裝 界面分為四大塊 1. Halcon的五個助手 1) 影像采集助手:與相機連接,設定相機引數,采集影像 2) 標定助手:九點標定或是其它的標定,生成標定檔案及內參外參,可以將像素單位轉換為長度單位 ......

    uj5u.com 2023-04-20 07:42:17 more
  • 在MacOS下使用Unity3D開發游戲

    第一次發博客,先發一下我的游戲開發環境吧。 去年2月份買了一臺MacBookPro2021 M1pro(以下簡稱mbp),這一年來一直在用mbp開發游戲。我大致分享一下我的開發工具以及使用體驗。 1、Unity 官網鏈接: https://unity.cn/releases 我一般使用的Apple ......

    uj5u.com 2023-04-20 07:40:19 more