資源簡介
比較簡單的時域平均程序,對仿真信號進行時域平均處理,效果非常明顯。
代碼片段和文件信息
clcclear?allclose?all
%%?時域同步平均仿真?
%?產生一個無噪的正弦信號
w=10;
t=0:(2*pi/w/500):(300*2*pi/w);
x=sin(w*t);
%?產生白噪聲
n=length(t);
ss=888;
randn(‘seed‘ss);
z=randn(1n);
z1=z-mean(z);
z2=z1/max(abs(z1));
z2=z2*5;
%?含噪的正弦信號?
x=x+z2;
%?繪制含噪的正弦信號
figure
plot(x);title(‘仿真含噪信號‘)
N?=?2048;
n1?=?0:N-1;
fs?=?2500/pi;
f=n1*fs/N;
x1=fft(xN);
mag1=abs(x1)*2/N;
figure
plot(f(1:N/2)mag1(1:N/2));
title(‘仿真信號頻譜‘)
%%?Averaging?in?time?domain
M=10000;??%?數據長度
xx=zeros(1M);?%?用于存儲時域平均后的數據
m1=1;m2=M;
nnnn=fix(n/M);?%?確定分段數
%?將數據分成nnnn段,
評論
共有 條評論