电脑桌面
添加小米粒文库到电脑桌面
安装后可以在桌面快捷访问

一个简单的单边带电台Matlab仿真VIP免费

一个简单的单边带电台Matlab仿真_第1页
1/4
一个简单的单边带电台Matlab仿真_第2页
2/4
一个简单的单边带电台Matlab仿真_第3页
3/4
一个简单的单边带电台仿真这个仿真以真实的音频信号作为输入,设计一个单边带发信机。将基带信号调制为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,声音仍然是清晰可懂的。

1、当您付费下载文档后,您只拥有了使用权限,并不意味着购买了版权,文档只能用于自身使用,不得用于其他商业用途(如 [转卖]进行直接盈利或[编辑后售卖]进行间接盈利)。
2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。
3、如文档内容存在违规,或者侵犯商业秘密、侵犯著作权等,请点击“违规举报”。

碎片内容

一个简单的单边带电台Matlab仿真

确认删除?
VIP
微信客服
  • 扫码咨询
会员Q群
  • 会员专属群点击这里加入QQ群
客服邮箱
回到顶部