Matlab编程实例(4) 相位角与相关系数曲线

 

%相位角与相关系数曲线

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)));

 

包含噪声

不包含噪声