信号相关性matlab仿真代码示例

clear all;

close all;

Fs = 10000;  %信号采样率        

dt = 1/Fs; %信号采样时间间隔，计算得到          

f=100;%信号频率

N=5;%一共采集了N个周期

L =round(N/dt/f);  %总长度，由采样率，信号频率，一共几个周期计算而来      

ti= (0:L-1)*dt; % 构建t轴，间隔是dt，数字从零到最后一个时刻

T = dt*L%整个采样序列的总时间

p=0;%用角度表示的相位；

P=p*pi/180;%角度变化为弧度；

% P=pi/2;%弧度表示的相位；

A=1.5;%信号幅度，峰峰值的一半

% X=A*sin(2*pi*f*ti)-0.5*A*sin(2*pi*2*f*ti);%原始信号是两个正弦波，一个频率是f幅度A，一个是2f幅度2A

S=A*sin(2*pi*f*ti+P)

subplot(2,1,1);

plot(S,'linewidth',2);

a1=1;

Y=zeros(1,200);

for m=10:200%比较一个频率变化的序列和原始信号的关系。

x1=a1*sin(2*pi*m*ti);%构建一系列正弦信号，频率从10-1000hz，

Y(m)=abs(dot(S,x1));%每一个信号和原始信号做点积，结果放到Y序列。

subplot(211);

plot(S,'r','linewidth',3);

hold on;

subplot(211);

plot(x1,'linewidth',3);

hold off;

subplot(212);

bar(Y);

pause(0.001);

end

S=A*sin(2*pi*f*ti+p);%原始信号是正弦波，一个频率是f幅度A相位p

figure;

y=zeros(1,360);

for m=1:360%比较一个频率变化的序列和原始信号的关系。

M=m*pi/180;%角度变化为弧度；

x1=a1*sin(2*pi*f*ti+M);%构建一系列正弦信号，相位从1-360，

y(m)=dot(S,x1);%每一个信号和原始信号做点积，结果放到Y序列。

subplot(211)

plot(S,'linewidth',3);

hold on;

subplot(211)

plot(x1,'linewidth',3);

hold off;

subplot(212)

bar(y);

pause(0.0001);

end

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删