【導(dǎo)讀】科教學(xué)中通信原理課程實(shí)驗(yàn)的演示系統(tǒng)。到的優(yōu)點(diǎn)，但是在計(jì)算與圖形顯示方面，顯得能力不足。方法，在教學(xué)中也可起到很好的輔助作用，它不僅能完成設(shè)計(jì)分析，還能完成通訊系統(tǒng)的許多實(shí)際條件無(wú)法完成的仿真實(shí)驗(yàn)。論文網(wǎng)用SIMULINK仿真平臺(tái)。庫(kù)，用模塊庫(kù)可以搭建自己的系統(tǒng)，完成系統(tǒng)設(shè)計(jì)分析。充分利用兩個(gè)軟件的優(yōu)點(diǎn)進(jìn)行混合編程，可以編制出界面。限于篇幅，本文僅以雙邊帶抑制載波調(diào)制為例子具體介紹了通。過(guò)MATLAB仿真的過(guò)程，對(duì)其他的章節(jié)實(shí)驗(yàn)都可以做類似分析和仿真，程序都已在附錄中給出。但利用ActiveX自動(dòng)化和DDE法簡(jiǎn)單方便，效率較高，功能強(qiáng)。[9]徐明遠(yuǎn),邵玉斌.MATLAB仿真在通信和電子工程中的應(yīng)用.這篇并不太長(zhǎng)的論文，不僅使。護(hù)的懵懂少年成長(zhǎng)為一個(gè)懂得人間寒暑的高校畢業(yè)生，其間的痛苦，感謝在百忙之中抽出時(shí)間審閱本論文的專家教授。分析和檢驗(yàn)?zāi)康?。這些方法的最大缺點(diǎn)是比較繁雜，而且需要花費(fèi)很多時(shí)間。理進(jìn)行了圖像仿真，在此基礎(chǔ)上集成了一個(gè)仿真的可視化開(kāi)發(fā)軟件。

　　

【正文】 r %直接在命令窗口中調(diào)用 message 函數(shù)，也可以在其他函數(shù)中調(diào)用 %調(diào)用例子 x=[,]。 f1=rand(4,4)。 hf1=hmessage(x,f1,4,4)。 %求平均信息量 hf1=1 hx=message(x,4)。 %此散信息源的平均信息量為 hx=2 c1=hxhf1。 disp(39。信道 1 的容量為 39。) c1 %信道 1 的信道容量為： c1=1 %對(duì)第二個(gè)信道 f2=rand(4,8)。 hf2=hmessage(x,f2,4,8)。 %平均互信息量為 hf2= hx=message(x,4)。 %平均信息量為 hx=2 c2=hxhf2。 %信道的容量為： c2= disp(39。信道 2 的容量為 39。) abs(c2) 其中 函數(shù)如下： %.哈夫曼編碼 function [h,l]=huffman(p)。 if length(find(p0))~=0, %判斷輸入是否概率分布 error(39。Input is not a ,there is negative ponent39。) end if abs(sum(p)1)10e10 error(39。Input is not a ,the sum of the ponents is not equal to 1.39。) end n=length(p)。 %將輸入的概率元素個(gè)數(shù) q=p。 for i=1:n1 [q,l]=sort(q)。 %將輸入的概率元素排序 m(i,: )=[l(1:ni+1),zeros(1,i1)]。 q=[q(1)+q(2),q(3:n),1]。 end for i=1:n1 c(i,: )=blanks(n*n)。 end c(n1,n)=39。039。 c(n1,2*n)=39。139。 for i=2:n1 c(ni,1:n1)=c(ni+1,n*(find(m(ni+1,:)==1))(n2):n*(find(m(ni+1,:)==1)))。 c(ni,n)=39。039。 c(ni,n+1:2*n1)=c(ni,1:n1)。 c(ni,2*n)=39。139。 for j=1:i1 c(ni,(j+1)*n+1:(j+2)*n)=c(ni+1,n*(find(m(ni+1,: )==j+1)1)+1:n*find(m(ni+1, : )==j+1))。 論文網(wǎng) )==i)*n)。 ll(i)=length(find(abs(h(i, : ))~=32))。 end l=sum(p.*ll)。 %常規(guī)雙邊帶調(diào)幅 echo on df= 本文來(lái)自六維論文網(wǎng)原文請(qǐng)找 3249114。 %頻度分辯力 t0=。 %定義 t0 信號(hào)的持續(xù)時(shí)間 的值 tz=。 %定義抽樣時(shí)間 fz=1/tz。 fa=50。 %定義載波頻率 snr=10。 %定義信噪比，用 DB 表示 snr_lin=10^(snr/10)。 %信噪比數(shù)值 a=。 %定義調(diào)拷制系數(shù) t=[0:tz:t0]。 %定義抽樣點(diǎn)數(shù)據(jù) %定義信號(hào) m m=zeros(1,501)。 for i=1:1:125, m(i)= i。 end。 for i=126:1:375, m(i)=m(125)i+125。 end。 for i=376:1:501, m(i)=m(375)i375。 end。 m=m/1000。 c=cos(2*pi*fa.*t)。 %載波信號(hào) m_n=m/max(abs(m))。 [M,m,df1]=fftseq(m,tz,df)。 %付里葉變換 ,fftseq 函數(shù)的程序段見(jiàn)后 M=M/fz。 %頻率縮放 ,便于做圖 f=[0:df1:df1*(length(m)1)]fz/2。 %定義頻率向量 u=(1+a*m_n).*c。 %將調(diào)質(zhì)信號(hào)調(diào)頎在載波信號(hào)上 [U,u,df1]=fftseq(u,tz,df)。 %對(duì)以調(diào)信號(hào)做付里葉變換 U=U/fz。 %h 頻率縮放 signal_power=ampower(u(1:length(t)))。 %計(jì)算信號(hào)的功率 ,ampower 程序見(jiàn)后 pmn=ampower(m(1:length(t)))/(max(abs(m)))^2。 %計(jì)算調(diào)制信號(hào)的功率 eta=(a^2*pmn)/(1+a^2*pmn)。 %計(jì)算調(diào)制效率 noise_power=eta*signal_power/snr_lin。 %計(jì)算噪聲功率 noise_std=sqrt(noise_power)。 %噪聲標(biāo)準(zhǔn)差 noise=noise_std*randn(1,length(u))。 %產(chǎn)生高斯分布噪聲 r=u+noise。 %總 接收信號(hào) [R,r,df1]=fftseq(r,tz,df)。 %對(duì)總信號(hào)進(jìn)行付里葉變換 R=R/fz。 %頻率縮放 %以下為結(jié)果顯示 signal_power % 顯示信號(hào)功率 eta %顯示調(diào)制效率 clf subplot(7,2,1) plot(t,m(1:length(t))) %做出調(diào)制信號(hào)的曲線 axis=([0 ])。 xlabel(39。Time39。) 本文來(lái)自六維論文網(wǎng)原文請(qǐng)找 3249114signal39。)。 subplot(7,2,2) plot(t,u(1:length(t))) %做出以調(diào)信號(hào)的曲線 axis=([0 ])。 xlabel(39。Time39。) title(39。the modulated signal39。)。 subplot(7,2,3) plot(t,c(1:length(t))) %做出載波的曲線 axis=([0 ])。 xlabel(39。Time39。) title(39。The carrier39。)。 subplot(7,1,3) plot(f,abs(fftshift(M))) %作出頻域的調(diào)制信號(hào) xlabel(39。Frequency39。)。 title(39。Spectrum of the message signal39。) subplot(7,1,4) plot(f,abs(fftshift(U))) %作出頻域的以調(diào)信號(hào) title(39。Spectrum of the modulated signal39。)。 xlabel(39。Frequency39。) subplot(7,1,5) plot(t,noise(1:length(t))) %作出噪聲曲線 title(39。noise sample39。)。 xlabel(39。Time39。) subplot(7,1,6) plot(t,r(1:length(t))) %做出總信號(hào)的時(shí)域曲線 title(39。signal and noise39。)。 xlabel(39。title39。) 論文網(wǎng) %.抑制雙邊代調(diào)制 clear echo on t0=2。 %信號(hào)持續(xù)時(shí)間 ts=。 %抽樣時(shí)間間隔 fc=100。 %載波頻率 fs=1/ts。 df=。 %頻率分辨力 t=[t0/2:ts:t0/2]。 %定義時(shí)間序列 %以下三句為定義信號(hào)序列 x=sin(200*t)。 m=x./(200*t)。 m(1001)=1。 %避免產(chǎn)生無(wú)窮大的值 c=cos(2*pi*fc.*t)。 %載波 u=m.*c。 %抑制載波調(diào)制 [M,m,df1]=fftseq(m,ts,df)。 %付里葉變換 M=M/fs。 [U,u,df1]=fftseq(m,ts,df)。 U=U/fs。 %頻率壓縮 f=[0:df1:df1*(length(m)1)]fs/2。 %pause clf subplot(4,2,1) plot(t,m(1:length(t))) %作出為調(diào)信號(hào)的波形 axis([,]) xlabel(39。時(shí)間 39。)。 title(39。未調(diào)信號(hào) 39。) subplot(4,2,3)