产生信号和噪声
- 格式:doc
- 大小:13.00 KB
- 文档页数:1
%产生信号和噪声
clear all
n=0:1:499; %0‐499,间隔为1
x=sin(0.1*n*pi);
u=2*(rand(1,500)-0.5); % rand(1,500) 1到500之间随机分布
s=x+u;
figure;%图像的设置
subplot(311),plot(n,x,'r'); %subplot(mnp) m表示是图排成行m,n表示图排成n列,p表示图所在的位置,p=1表示从左到右从上到下的第一个位置。
n--x轴,x--y轴,r--red。
title('原始信号波形');
axis([0 100 -1.2 1.2]);%表示横纵轴坐标的上下限
xlabel('n');
ylabel('x[n]');
subplot(312),plot(n,u,'r');
title('产生的随机噪声');
axis([0 100 -2.2 2.2]);
xlabel('n');
ylabel('噪声幅度');
subplot(313),plot(n,s,'r');
title('叠加噪声后的信号波形');
axis([0 100 -4.2 4.2]);
xlabel('n');
ylabel('s[n]');
%计算x的自相关
figure;
[r,n]=xcorr(x,x);%x的自相关函数
subplot(211),plot(n,r,'g');
title('原始信号的自相关图形');
[r,n]=xcorr(s,s);%s的自相关函数
subplot(212),plot(n,r,'g');
title('叠加噪声后信号的自相关图形');
%产生冲击串
figure;
i=(n/20==floor(n/20)); %floor(n/2)表示对(n/2)向下取整。
re=xcorr(i,s);%i和s的互相关函数
subplot(211),stem(i);%显示i的信号波形
title('冲击串图象');
N=length(re);
n=1:N;
subplot(212),plot(n,re);
title('经过躁声的周期性函数与冲击串的互相关图象');。