statmoments.m

function [v,unv]=statmoments(p,n) % STATMOMENTS Computes statistical central moments of image histogram. % [v,unv]=statmoments(p,n) computes up to the Nth statistical % central moment of a histogram whose components are in vector % P. The length of must equal 256 or 65536. % ****from ggbondg**** % T

clear                              %清空内存及清屏
clc
%% 原始信号
Fs=44101;                          %音频信号采样率
signal = wavread（‘test.wav‘）;      %读入信号
figure                              
plot（signal）;                      %画出原始音频信号的波形
title（‘原始音频信号‘）
%msgbox（‘正在播放原始音频信号！‘）;  %弹出提示框
%wavplay（signalFs）                 %播放原始音频信号
FFT_NUM = 65536;
freq = （-FFT_NUM/2:1:FFT_NUM/2-1）/FFT_NUM*Fs/1e3;
freq2_val = fft（ signal  FFT_NUM）/Fs;          %做FFT求出频谱
freq2_val=freq2_val.*conj（freq2_val）/FFT_NUM*Fs;
SpecVal（1:FFT_NUM/2） = freq2_val（FFT_NUM/2+1:FFT_NUM）;
SpecVal（FFT_NUM/2+1:FFT_NUM） = freq2_val（1:FFT_NUM/2）; 
figure
plot（ freq 10*log10（ abs（SpecVal）） ）;grid
%% 带噪信号
SNR=30;                            %设置信噪比
noise=10^（-SNR/10）*（randn（size（signal）））; 
b=firlvbo（2565e3FsFs‘hpf‘）;