亲们 求围观 求火速围观 大神们都来指导指导 Matlab 吧·····
%*************** 生成输入信号 *****%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
请大家帮帮忙吧 帮我把后边打注释的完善下 还有就是看看那些有疑问的地方 真心谢谢大家了