%***************  生成输入信号  *****%
    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

沙发下
一着急把 概率密度那个代码忘了写 求斑竹不要删帖啊  真心求 我这应该不算求作业答案吧 最起码我还是思索了 呵呵呵

楼主要干什么？

上面有注释 谢谢你来看啊

你这个是什么语言？发出来是要解决什么问题吗？

这个是matlab语言 我也是没办法了 所以才拿来让大家看看 虽说这是C的天地  但是我相信还是有人会这个的
我的问题是请哪位帮帮忙 根据上面的整体 然后帮我把中间注释空着的那些代码写一写  就这些 真心感谢大家了 时间真的不够了

看着象matlab,还真是，只能友情帮顶了，上大学时在图书馆看到很多学科都有用matlab的书籍，觉得matlab很牛，却一直连碰都没碰过，遗憾啊，前两天和同事聊起来这个，还被鄙视了，说这个是研究生课程，本科接触不到。

现在本科就要接触了 其实这个上手很快 不难 难的是怎么熟练  我也是大一的时候学习过 好久没有用 所以现在不熟悉了

matlab早忘了，不过最近要用到，正准备学一学

matlab 选修课 已经忘光了