SSB(滤波法法)调制解调Matlab.m文件源码

%%

%将已调信号送入信道

%先根据所给信噪比产生高斯白噪声

signal_power = power_x(samu(1:Lt)); %已调信号的平均功率

noise_power = (signal_power * fs)/(snr_lin*4*Bw_eq); %求出噪声方差（噪声均值为0） noise_std = sqrt(noise_power); %噪声标准差

noise = noise_std * randn(1,Lt); %产生噪声

%画出信道高斯白噪声波形及频谱，此时，噪声已实现，为确知信号，可求其频谱

subplot(323);

plot(t,noise);

axis([0,t0,-max(noise),max(noise)]);

xlabel('t');title('噪声信号');

subplot(324);

[noisef,noise,df1,f] = T2F_new(noise,dt,df,fs); %噪声频谱

plot(f,fftshift(abs(noisef))); %画出噪声频谱

xlabel('f');title('噪声频谱');

%%

%信道中的信号 %叠加了噪声的已调信号频谱

sam = samu(1:Lt) + noise(1:Lt);

subplot(325);

plot(t,sam);

axis([0,t0,-max(sam),max(sam)]);

xlabel('t');title('信道中的信号');

subplot(326);

[samf,sam,df1,f] = T2F_new(sam,dt,df,fs); %求出叠加了噪声的已调信号频谱 plot(f,fftshift(abs(samf))); %画出叠加了噪声的已调信号频谱 xlabel('f');title('信道中信号的频谱');

axis([-fc-5*f0,fc+5*f0,0,max(samf)+0.1]);

set(gca,'XTick', -10:10:10);

set(gca,'XGrid','on');

%%

%----------------经过带通滤波器

%经过理想滤波器后的信号及其频谱

DEM = H_high.*samuf; %滤波器输出信号的频谱

[dem] = F2T_new(DEM,fs);%滤波器输出信号的波形

figure(4);