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;
E = int(m1,t,0,Inf); %求白噪声的均值
figure(4);
plot(t,E,'y-'); %输出白噪声均值的时域波形
title('噪声均值波形');
xlabel('t'),ylabel('E');
hold on;
grid on;
目前程序是这样,但是均值部分有错,??? Error using ==> mtimes
Inner matrix dimensions must agree.如果想求出xn的均值应该怎么写?