sinhscy4.m

个人硕士论文的关于数字接收机调制解调方面的matlab程序

%对sin信号进行抽样得到的结论4
clc
clear
fc=1;%sin信号的基准频率
a=-4;%信号的起始时间
b=4;%信号的终止时间
ts=1/(32*fc);%为使信号能较平滑的画出，设置的采样时间间隔
fs=1/ts;%由上面的采样时间间隔得到的采样频率
t=a:ts:b;
f=-fs/2+fs/length(t):fs/length(t):fs/2;
%------
f_s=fc*3;%设置的采样频率，相对于'基准频率'时的sin信号是过采样
t_s=1/f_s;%由上面的采样频率得到的采样时间间隔
n=t_s/ts;
t_c=a:t_s:b;

w1=sin(2*pi*fc*t);%频率为'基准频率'时的sin信号
w1_spect=fftshift(abs(fft(w1)));

x1=sin(2*pi*fc*t_c);%用上面规定的采样频率，对频率为'基准频率'的sin信号进行抽样后得到的抽样值
y1=zeros(1,n*(length(x1)-1)+1);
y1(1:n:n*length(x1))=x1;
y1_spect=fftshift(abs(fft(y1)));

z1=convn(y1,sinc(t/t_s),'same');%重建信号
z1_spect=fftshift(abs(fft(z1)));%重建信号的频谱

%----
w2=sin(2*pi*(fc+f_s)*t);%频率为'基准频率+采样频率'时的sin信号
w2_spect=fftshift(abs(fft(w2)));

x2=sin(2*pi*(fc+f_s)*t_c);%用上面规定的采样频率，对频率为'基准频率+采样频率'的sin信号进行抽样后得到的抽样值
y2=zeros(1,n*(length(x2)-1)+1);
y2(1:n:n*length(x2))=x2;
y2_spect=fftshift(abs(fft(y2)));

z2=convn(y2,sinc(t/t_s),'same');%重建信号
z2_spect=fftshift(abs(fft(z2)));%重建信号的频谱

%------
w3=sin(2*pi*(fc+2*f_s)*t);%频率为'基准频率+2*采样频率'时的sin信号
w3_spect=fftshift(abs(fft(w3)));

x3=sin(2*pi*(fc+2*f_s)*t_c);%用上面规定的采样频率，对频率为'基准频率+2*采样频率'的sin信号进行抽样后得到的抽样值
y3=zeros(1,n*(length(x3)-1)+1);
y3(1:n:n*length(x3))=x3;
y3_spect=fftshift(abs(fft(y3)));

z3=convn(y1,sinc(t/t_s),'same');%重建信号
z3_spect=fftshift(abs(fft(z3)));%重建信号的频谱


figure(1)
subplot(3,2,1)
plot(t,w1);title('频率为基准频率的sin信号');grid
subplot(3,2,2)
plot(f,w1_spect);title('频率为基准频率的sin信号的频谱');grid
subplot(3,2,3)
stem(t_c,x1);title('抽样点');grid
subplot(3,2,4)
plot(f,y1_spect);title('抽样点的频谱');grid
subplot(3,2,5)
plot(t,z1);title('用抽样点重建的信号');grid
subplot(3,2,6)
plot(f,z1_spect);title('重建信号的频谱');grid

figure(2)
subplot(3,2,1)
plot(t,w2);title('频率为基准频率+采样频率的sin信号');grid
subplot(3,2,2)
plot(f,w2_spect);title('频率为基准频率+采样频率的sin信号的频谱');grid
subplot(3,2,3)
stem(t_c,x2);title('抽样点');grid
subplot(3,2,4)
plot(f,y2_spect);title('抽样点的频谱');grid
subplot(3,2,5)
plot(t,z2);title('用抽样点重建的信号');grid
subplot(3,2,6)
plot(f,z2_spect);title('重建信号的频谱');grid

figure(3)
subplot(3,2,1)
plot(t,w3);title('频率为基准频率+2*采样频率的sin信号');grid
subplot(3,2,2)
plot(f,w3_spect);title('频率为基准频率+2*采样频率的sin信号的频谱');grid
subplot(3,2,3)
stem(t_c,x3);title('抽样点');grid
subplot(3,2,4)
plot(f,y3_spect);title('抽样点的频谱');grid
subplot(3,2,5)
plot(t,z3);title('用抽样点重建的信号');grid
subplot(3,2,6)
plot(f,z3_spect);title('重建信号的频谱');grid