%相位角与相关系数曲线
close all;
clear all;
Samp1=200; %设置信号的采样精度
Samp2=200; %设置相位角p分割精度
A=10;%信号幅值
w=1;%信号角频率
na=5;%噪声放大系数
g=zeros(Samp2,Samp1);%产生Samp2 x Samp1的二维零矩阵,每一行对应一个p值时的函数
t = linspace(0,6*pi/w,Samp1); %创建函数向量
p= linspace(0,pi/3,Samp2);%创建相位角向量
f=A*sin(w*t)+na*randn(1,Samp1);%构造一个标准函数,不含相位角
for m=1:1:Samp2;%m代表不同的相位角p
g(m,:)=A*sin(w*t+p(m))+na*randn(1,Samp1);%构造二位矩阵
r(m)=corr(f,g(m,:));%计算相关系数
end;
plot(p,r);
function r=corr(f,g);
%计算函数f与g的相关系数
r=dot(f,g)/(sqrt(dot(f,f))*sqrt(dot(g,g)));
包含噪声
![](https://images2015.cnblogs.com/blog/299683/201603/299683-20160329131600551-1474874938.jpg)
不包含噪声
![](https://images2015.cnblogs.com/blog/299683/201603/299683-20160329131619738-398329529.jpg)