[½»Á÷] ÇëÎÊMATLABÀï,ÔõôÇóÒ»¸ö¸ß˹°×ÔëÉùµÄµÄ¾ùÖµ? ÒÑÓÐ1È˲ÎÓë

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
    %*****************  ÔëÉùµÄÊý×ÖÌØÐÔ  ****************%
    %*****************    1.¾ùÖµ       *****************%
    m1 = t.*xn;
    for i=1:N
    E(i) =sum(m1(1:i)); %Çó°×ÔëÉùµÄ¾ùÖµ
    end
    figure(4);
    plot(t,E,'y-');          %Êä³ö°×ÔëÉù¾ùÖµµÄʱÓò²¨ÐÎ
    title('ÔëÉù¾ùÖµ²¨ÐÎ');
    xlabel('t'),ylabel('E');
        hold on;
        grid on;