3
回答
亲们 求围观 求火速围观 大神们都来指导指导 Matlab 吧·····
利用AWS快速构建适用于生产的无服务器应用程序,免费试用12个月>>>   
%***************  生成输入信号  *****%
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        
举报
琪乐无穷
发帖于6年前 3回/409阅
顶部