%***************  生成输入信号  *****%
    f = 2000                  %设定基频率为2000HZ
    Fs=44100;                 %设定采样频率Fs=44.1kHz
    N=256;                    %取的样本点数N
    n=0:N-1;                  %建立矩阵
    t=n/Fs;                   %采样时间t
    w = 2*pi*f;
    T = 2*pi/w;                   %求出正弦函数的周期
    x1 = sin(w*t);
    x2 =awgn(x1,5,'measured');      %加入SNR为5dB的高斯白噪声得到输入信号x2
    xn =x2-x1;                    %间接获得白噪声xn

    %*****************   时域波形   *************%
    figure(1);                   
    plot(t,x1,'k-');          %正弦波信号时域波形
    title('正弦信号时域波形');
    xlabel('t'),ylabel('x1(t)');
    hold on;
    grid on;

    figure(2);                 
    plot(t,xn,'r-');        %高斯白噪声时域波形
    title('高斯白噪声信号时域波形');
    xlabel('t'),ylabel('xn(t)');
    hold on;
    grid on;

    figure(3)               
    plot(t,x2,'b-');       %输入信号时域波形
    title('输入信号时域波形');
    xlabel('t'),ylabel('x2(t)');
    hold on;
    grid on;
    
    %*****************  噪声的数字特性  ****************%
    
    %*****************    1.均值       *****************% 
    
    E = mean(xn);      %求白噪声的均值
     figure(4);
    plot(t,E,'b-');          %白噪声均值的时域波形
    title('噪声均值波形');
    xlabel('t'),ylabel('E');
    hold on;
    grid on;

   

    %****************      2.方差       *****************%
    D = var(xn);
     figure(5);
    plot(t,D,'b-');          %白噪声方差的时域波形
    title('噪声方差波形');
    xlabel('t'),ylabel('D');
    hold on;
    grid on;  
    
    %****************     3.均方值      ****************%
    
    E2 = E+D;
    figure(6);
    plot(t,E2,'b-');          %白噪声均方值的时域波形
    title('噪声均方值波形');
    xlabel('t'),ylabel('E2');
    hold on;
    grid on;

     %***************   4.输入信号自相关函数    ****************%
     R=xcorr(xn);     %*******自相关函数*******%
     t1=(-length(xn)+1:length(xn)-1)/Fs;
     figure(7);
     plot(t1,R,'r-');          %白噪声自相关函数的时域波形
     title('噪声自相关波形');
     xlabel('t1'),ylabel('R');
     hold on;
     grid on;  
    

     
    %****************      5.概率密度      *****************%  这个多次尝试都做不出来
    figure()
    

    
    
    

    %****************       6.频谱       ****************%    这里编译之后显示有错误 可能要改一下
     figure(9)
     H_arr=fft(xn,N);      
     cm=abs(H_arr);
     fl=(0:length(H_arr)-1)'*44100/length(H_arr);
     plot(f,H_arr,'b-')
     title('线性信号的频谱')
     xlabel('f'),ylabel('S_x_i(f)');
     hold on;
     grid on


    %***************        7.功率谱密度      *****************%  还有这里 出图像之后是直线 不知道对不
    figure(10)
    X_arr=fft(R);      
    cm=abs(X_arr);
    fl=(0:length(X_arr)-1)'*44100/length(X_arr);
    plot(f,X_arr,'r-')
    title('线性信号的功率谱')
    xlabel('f'),ylabel('S_x_i(f)');
    hold on;
    grid on