一个简单的单边带电台仿真这个仿真以真实的音频信号作为输入，设计一个单边带发信机。将基带信号调制为SSB信号后送入带通型高斯噪声信道，加入给定功率的噪声之后，再送入单边带接收机。单边带接收机将型号解调下来，通过计算机声卡将解调信号播放出来试听效果，从而对信道信噪比与解调音质之间的关系进行主观测试。现设计一个单边带发信机、带通信道和相应的接收机，参数定位如下值。(1)输入信号为一个话音信号，采样率为8000Hz。话音输入后首先进行预滤波，滤波器是一个频率范围在[300,3400]Hz的带通滤波器，其目的是将话音频频谱限制在3400Hz以下。单边带调制的载波频率设计为10kHz，调制输出上边带。要求观测单边带调制前后的信号功率谱。(2)信道是一个带限高斯噪声信道，其通带频率范围是[10000,13500]Hz。能够根据信噪比SNR的要求加入高斯噪声。(3)接收机采用相干解调方式。为了模拟载波频率误差对解调话音音质的影响，设本地载波频率为，与发信机载波平率相差200Hz。解调滤波器设计为300～3400Hz的带通滤波器。1.对音频输入与调制的仿真根据设置的参数，系统中信号最高频率约为14kHz。为了较好地显示调制波形，系统仿真采样率设为50kHz，满足取样定理。由于话音信号的采样率为8000Hz，与系统仿真采样率不等，因此，在进行信号处理之前，必须将话音的采样率提高到50kHz，用插值函数来做这一任务。先编写程序将基带音频信号读入，进行[300，3400]Hz的带通滤波，并将信号采样率提高到50kHz，进行单边带调制之后，将调制输出结果保存为wav文件，文件名为。程序如下：%FileName:%功能：音频录入与调制clc;clearall;%录音，采样点数为1000，采样率为8000%jilu=wavrecord(5*8000,8000,'double');%wavwrite(jilu,'');[wav,fs]=wavread('');%计算声音的时间长度t_end=1/fs*length(wav);%仿真系统采样率Fs=50000;%仿真系统采样时间点t=1/Fs:1/Fs:t_end;%设计300～3400Hz的带通滤波器H(z)[fenzi,fenmu]=butter(3,[3003400]/(fs/2));%对音频信号进行滤波wav=filter(fenzi,fenmu,wav);%利用插值函数将音频信号的采样率提升为Fs=50kHzwav=interp1([1/fs:1/fs:t_end],wav,t,'spline');%音频信号的希尔伯特变换wav_hilbert=imag(hilbert(wav));%载波频率fc=10000;%单边带调制SSB_OUT=wav.*cos(2*pi*fc*t)-wav_hilbert.*sin(2*pi*fc*t);figure(1);subplot(2,2,1);plot(wav(53550:53750));axis([0200]);subplot(2,2,2);psd(wav,10000,Fs);axis([02500-2010]);subplot(2,2,3);plot(SSB_OUT(53550:53750));axis([0200]);subplot(2,2,4);psd(SSB_OUT,10000,Fs);%将SSB调制输出存盘备用wavwrite*SSB_OUT,'');程序仿真结果：2.对指定信噪比信道的仿真仿真指定信噪比信道，仿真函数如下。%FileName:%功能：计算实际信噪比clear;[in,Fs]=wavread('');SNRdB=20;out=ChanelSimulink(in,SNRdB);wavwrite(out,Fs,'');计算信噪比为20db时的信道输出，将结果保存为文件。输出仿真的实际测量信噪比为db。%%功能：仿真指定信噪比的信道functionout=ChanelSimulink(in,SNRdb)%SNR_db设定信噪比%in输入信号序列%out信道输出序列%系统采样率Fs=50000;Power_of_in=var(in);Power_of_noise=Power_of_in/(10.^(SNRdb/10));%信道带宽bandwidth=13500-10000;%噪声功率谱密度值W/HzNO=Power_of_noise/bandwidth;Gause_noise=sqrt(NO*Fs/2).*randn(size(in));%噪声通道10～[num,den]=butter(4,[1000013500]/(Fs/2));signal_of_filter_out=filter(num,den,in);noise_of_filter_out=filter(num,den,Gause_noise);SNR_dB=10*log10(var(signal_of_filter_out)/var(noise_of_filter_out));%测量得出信噪比%信道输出out=signal_of_filter_out+noise_of_filter_out;SNR_dBclear;Fs=50000;%读入信道输出信号数据[recvsignal,Fs]=wavread('');t=(1/Fs:1/Fs:length(recvsignal)/Fs)';%本地载波频率fc_local=10000-200;%本地载波local_carrier=cos(2*pi*fc_local.*t);%相干解调xianggan_out=recvsignal.*local_carrier;%设计300～3400Hz的带通滤波器[fenzi,fenmu]=butter(3,[3003400]/(Fs/2));demod_out=filter(fenzi,fenmu,xianggan_out);sound(demod_out/max(demod_out),Fs);wavwrite(demod_out,Fs,'');subplot(1,2,1);plot(5*demod_out(53550:53750));axis([0200]);subplot(1,2,2);psd(5*demod_out,10000,Fs);axis([025000-2010]);解调后图像，可见仿真结果还行！解调输出信号被保存为音频文件，并由sound函数播放。聆听播放解调输出信号的声音可知，在20dB信道信噪比条件下，即使解调本地载波频率误差达到200Hz，声音仍然是清晰可懂的。
