Matlab编程实例(2) 同期平均

 

%多点同期平均

close all;

clear all;

pi = 3.14159;

Samp2=input('您需要几组信号做同期平均？')

Samp1=1000 %设置采样精度

t = linspace(0,4*pi,Samp1); %创建向量

g=zeros(Samp2,Samp1);

A=10;

p=0;

w=1;

na=1;%噪声放大系数

s=A*sin(w*t+p)%原信号

temp=0;

for m=1:1:Samp2;

g(m,:)=A*sin(w*t+p)+na*randn(1,Samp1);%构造二位矩阵，每行为一加了噪声的信号

end;

for n=1:1:Samp1; 

    temp=0;

    for k=1:1:Samp2;%将每一列求和

    temp=temp+g(k,n);

    end;

y(n)=temp/Samp2;%同期平均，将一列求平均数

end;

subplot(3,1,1);plot(t,s,'r');

title('原信号');

subplot(3,1,2);plot(t,g(1,:),'b');%画出第一行

title('叠加白噪声后信号');

subplot(3,1,3);plot(t,y,'k');

title('同期平均后信号');

 

 

10组同期平均后

 

1000组同期平均后