实指数、复指数序列实验打印.docxVIP

%%实指数、复指数序列 clc ; clear all; close all; n0=0; nf=10; ns=3; n3=n0:nf; x3=(0.75).^n3; n4=n0:nf; x4=exp((-0.2+0.7j)*n4); stem(n3,x3); title(实指数序列) figure subplot(2,2,1),stem(n4,real(x4)); title(复指数序列实部) subplot(2,2,2),stem(n4,imag(x4)); title(复指数序列虚部) subplot(2,2,3),stem(n4,abs(x4)); title(复指数序列幅度) subplot(2,2,4),stem(n4,angle(x4)); title(复指数序列相位) 图1 信号时域实指数图示 图2 信号时域复指数图示 clc; clear all; close all; fs=500; N=200; t=0:1/fs:(N-1)/fs; %时间序列 n=0:N-1; %采样点序列 f=n*fs/N; %采样频率序列 y=cos(2*pi*100*t); %信号函数 Y=abs(fft(y)); %傅里叶变换 figure plot(t,y); title(信号波形图); Y=Y-mean(Y); figure plot(f(1:N/2),Y(1:N/2)) title(信号频率图) 图3 信号傅里叶变换波形图 图4 信号傅里叶变换频率图 clc; clear all; close all; N=1000; %采样点数 n=[0:N-1]; %采样点数序列 fs=500; %采样频率 t=n/fs; %时间序列 lag=100; %延迟矢量 x=sin(2*pi*t+0.5*pi); %信号函数 [c,lags]=xcorr(x,lag,unbiased); %自相关计算 figure plot(t,x); xlabel(t); ylabel(x(t)); grid; figure plot(lags/fs,c); xlabel(t); ylabel(Rx(t)); title(x的自相关函数); grid; 图5 原始信号图示 图6 自相关函数图示 clc; clear all; close all; N=1000; n=[0:N-1]; fs=500; t=n/fs; w1=2*pi; w2=1*pi; o=0.5*pi; a=0.25*pi; x=sin(w1*t+o); y=sin(w2*t+o-a); [c1,rxx]=xcorr(x); [c2,rxy]=xcorr(x,y); [c3,ryy]=xcorr(y); figure subplot(3,1,1) plot(x); title(信号x的波形图); grid; subplot(3,1,2) plot(y); title(信号y的波形图); grid; subplot(3,1,3) plot(rxy,c2); title(信号x和y的互相关函数波形图); grid; 图7 互相关函数图示