Matlab信号处理——噪声信号与原信号之间夹角的余弦

1、构造带噪声信号，将信号循环平移后与原先信号求余弦1、流程图显示信号y初始化定义带噪声信号ynlen=1截取第一段波y1形截取第二段波y2合并y1、y2得Y求Y与y的夹角余弦nlen=nlen+1nlen=1000?NY显示结果2、程序代码%程序说明：构造带噪声信号，将信号循环平移后与原先信号求余弦%时间2012年3月12日，韩宝安,学号10011311477clc; %清空 close all; %关闭所有窗口clear all; %清空所有变量t=linspace(0,2*pi,300); %定义变量y=10*sin(t)+randn(1,length(t); %构造带噪声信号subplot

2、(2,1,1); %两行一列第一个图plot(t,y); %显示带噪声信号grid on; %开启网格线xlabel(t,fontsize,15); %横坐标显示为t，大小为15号字ylabel(y,fontsize,15); %纵坐标显示为y，大小为15号字title(figure1:带噪声信号的波形,fontsize,15); %显示标题，15号字set(gca,xtick,pi/2,pi,3*pi/2,2*pi,fontsize,15); %设置横坐标标注位置xticklabels=pi/2,pi,3pi/2,2pi; %标注位置显示内容set(gca,xticklabel,xtickl

3、abels); %显示for nlen=1:length(t) %在空间长度上循环 index1=1:length(t)-nlen; %定义区间1 y1=y(index1); %截取区间1上波形 index2=length(t)-nlen+1:length(t); %定义区间2 y2=y(index2); %截取区间2上波形 Y=y2,y1; %合共两波形 dotYy=dot(y,Y); %求点积 mo_Y=norm(Y); %求模 mo_y=norm(y); %求模 cosYy(nlen)=dotYy/mo_Y/mo_y; %求夹角余弦endsubplot(2,1,2); %两行一列第二个图

4、plot(t(1:length(t),cosYy); %显示夹角余弦grid on; %开启网格线xlabel(nlen,fontsize,15); %横坐标显示为t，大小为15号字ylabel(cosYy,fontsize,15); %纵坐标显示为y，大小为15号字title(figure2:带噪声信号平移后与原带噪声信号夹角的余弦,fontsize,15); %显示标题，15号字set(gca,xtick,pi/2,pi,3*pi/2,2*pi,fontsize,15); %设置横坐标标注位置xticklabels=pi/2,pi,3pi/2,2pi; %标注位置显示内容set(gca,x

5、ticklabel,xticklabels); %显示%程序验证，已知创建了一个长度为300的线性空间，那么如果将此信号向右移动pi/2，则信号将与原信号正交，即夹角余弦接近0nlen=input(程序验证：请输入移动距离nlen（0=nlen=2）:); %确定移动距离ind=int32(300/2/pi*nlen); %将移动距离转化成移动点数index1=1:length(t)-ind; %定义区间1y1=y(index1); %截取区间1上波形index2=length(t)-ind+1:length(t); %定义区间2y2=y(index2); %截取区间2上波形Y=y2,y1;

6、%合共两波形figure; %在另一窗口显示图像subplot(2,1,1); %两行一列第一个图plot(t,y); %显示带噪声信号grid on; %开启网格线xlabel(nlen,fontsize,15); %横坐标显示为t，大小为15号字ylabel(y,fontsize,15); %纵坐标显示为y，大小为15号字title(figure1:带噪声信号的波形,fontsize,15); %显示标题，15号字set(gca,xtick,pi/2,pi,3*pi/2,2*pi,fontsize,15); %设置横坐标标注位置xticklabels=pi/2,pi,3pi/2,2pi;

7、%标注位置显示内容set(gca,xticklabel,xticklabels); %显示subplot(2,1,2); %两行一列第二个图plot(t,Y); %显示移动后的波形grid on; %开启网格线xlabel(nlen,fontsize,15); %横坐标显示为t，大小为15号字ylabel(y,fontsize,15); %纵坐标显示为y，大小为15号字title(figure2:带噪声信号平移后的波形,fontsize,15); %显示标题，15号字set(gca,xtick,pi/2,pi,3*pi/2,2*pi,fontsize,15); %设置横坐标标注位置xtickl

8、abels=pi/2,pi,3pi/2,2pi; %标注位置显示内容set(gca,xticklabel,xticklabels); %显示Y=y2,y1; %合共两波形dotYy=dot(y,Y); %求点积mo_Y=norm(Y); %求模mo_y=norm(y); %求模cosYy=dotYy/mo_Y/mo_y; %求夹角余弦fprintf(移动上述距离后的夹角余弦为：%gn,cosYy); %输出夹角余弦3、程序运行结果主程序部分：Fig 1. 主程序部分图形程序验证部分：（1）、从Fig1中可以看出当移动距离为/2时，夹角余弦接近为0，移动后的波形见Fig2，经过程序计算，此时的夹角余弦为0.000271333，已很接近0。Fig 2.移动距离为/2时