Matlab源程序代碼.doc

正弦波的源程序：(一)，用到的函數(shù)1,f2t函數(shù)function x=f2t(X)global dt df t f T N%x=f2t(X)%x為時(shí)域的取樣值矢量%X為x的傅氏變換%X與x長(zhǎng)度相同并為2的整冪%本函數(shù)需要一個(gè)全局變量dt(時(shí)域取樣間隔)X=X(N/2+1:N),X(1:N/2);x=ifft(X)/dt;end 2,t2f函數(shù)。function X=t2f(x)global dt df N t f T%X=t2f(x)%x為時(shí)域的取樣值矢量%X為x的傅氏變換%X與x長(zhǎng)度相同,并為2的整冪。%本函數(shù)需要一個(gè)全局變量dt(時(shí)域取樣間隔)H=fft(x);X=H(N/2+1:N),H(1:N/2)*dt;end(二),主程序。1，%（1）繪出正弦信號(hào)波形及頻譜global dt df t f Nclose allk=input(取樣點(diǎn)數(shù)2k, k取10左右);if isempty(k), k=10; endf0=input(f0=取1(kz)左右);if isempty(f0), f0=1; endN=2k;dt=0.01; %msdf=1/(N*dt); %KHzT=N*dt; %截短時(shí)間Bs=N*df/2; %系統(tǒng)帶寬f=-Bs+df/2:df:Bs; %頻域橫坐標(biāo)t=-T/2+dt/2:dt:T/2; %時(shí)域橫坐標(biāo)s=sin(2*pi*f0*t); %輸入的正弦信號(hào)S=t2f(s); %S是s的傅氏變換a=f2t(S); %a是S的傅氏反變換a=real(a);as=abs(S);subplot(2,1,1) %輸出的頻譜 plot(f,as,b);gridaxis(-2*f0,+2*f0,min(as),max(as) xlabel(f (KHz)ylabel(|S(f)| (V/KHz) %figure(2)subplot(2,1,2)plot(t,a,black) %輸出信號(hào)波形畫圖gridaxis(-2/f0,+2/f0,-1.5,1.5)xlabel(t(ms)ylabel(a(t)(V)gtext(頻譜圖)最佳基帶系統(tǒng)的源程序：(一)，用到的函數(shù)f2t函數(shù)和t2f函數(shù)。代碼(二),主程序globaldt t f df N Tclose allclear Eb_N0 Pek=input(取樣點(diǎn)數(shù)2k, k取13左右);if isempty(k), k=13; endz=input(每個(gè)信號(hào)取樣點(diǎn)數(shù)2z, z(二),主程序global dt t f df N Tclose allclear Eb_N0 Pek=input(取樣點(diǎn)數(shù)2k, k取13左右);if isempty(k), k=13; endz=input(每個(gè)信號(hào)取樣點(diǎn)數(shù)2z, z(二),主程序global dt t f df N Tclose allk=input(取樣點(diǎn)數(shù)2k, k取10左右);if isempty(k), k=13; endz=input(每個(gè)信號(hào)取樣點(diǎn)數(shù)2z,zk);if isempty(z), z=5; endf0=input(f0=取1(kz)左右);if isempty(f0), f0=1; endN=2kL=2z;M=N/L;Rb=2;Ts=0.5 %碼元寬度是0.5usdt=Ts/L; %時(shí)域采樣間隔df=1/(N*dt) %MHzT=N*dt %截短時(shí)間Bs=N*df/2 %系統(tǒng)帶寬f=-Bs+df/2:df:Bs; %頻域橫坐標(biāo)t=-T/2+dt/2:dt:T/2; %時(shí)域橫坐標(biāo)EP=zeros(size(f)+eps;for ii1=1:30b=round(rand(1,M); %產(chǎn)生隨機(jī)碼s=zeros(1,N);for ii=1:L/2;s(ii+0:M-1*L)=b;endS=t2f(s); %S是s的傅氏變換a=f2t(S); %a是S的傅氏反變換a=abs(a);P=S.*conj(S)/T; %功率譜EP=(EP*(ii1-1)+P+eps)/ii1;figure(1) %輸出的功率畫圖plot(f,10*log10(EP),b);gridaxis(-15,15,-100,20);xlabel(f (KHz);ylabel(P (V/KHz) ;title(功率譜圖)figure(2) %輸出信號(hào)波形畫圖plot(t,a,b)gridaxis(-3,3,0,1.3);xlabel(t (ms);ylabel(s(t) (V);title(單極性RZ的時(shí)域圖)end8psk的源程序：functionpb,ps=cm_sm32(snr_in_dB) N=10000;E=1;snr=10(snr_in_dB/10);sgma=sqrt(E/(3*2*snr);s000=1 0;s001=cos(pi/4) sin(pi/4);s011=0 1;s010=cos(3*pi/4) sin(3*pi/4);s110=-1 0;s111=cos(5*pi/4) sin(5*pi/4);s101=0 -1;s100=cos(7*pi/4) sin(7*pi/4);for i=1:N,temp=rand;if(temp0.125),dsource1(i)=0;dsource2(i)=0;dsource3(i)=0;elseif(temp0.25),dsource1(i)=0;dsource2(i)=0;dsource3(i)=1;elseif(temp0.375),dsource1(i)=0;dsource2(i)=1;dsource3(i)=0;elseif(temp0.5),dsource1(i)=0;dsource2(i)=1;dsource3(i)=1;elseif(temp0.625), dsource1(i)=1;dsource2(i)=0;dsource3(i)=0;elseif(temp0.75),dsource1(i)=1;dsource2(i)=0;dsource3(i)=1;elseif(temp0.875),dsource1(i)=1;dsource2(i)=1;dsource3(i)=0;elsedsource1(i)=1;dsource2(i)=1;dsource3(i)=1;end;end;numofsymbolerror=0;numofbiterror=0;for i=1:N,n(1)=gngauss(sgma);n(2)=gngauss(sgma);if(dsource1(i)=0)&(dsource2(i)=0)&(dsource3(i)=0),r=s000+n;elseif(dsource1(i)=0)&(dsource2(i)=0)&(dsource3(i)=1),r=s001+n;elseif(dsource1(i)=0)&(dsource2(i)=1)&(dsource3(i)=0),r=s010+n;elseif(dsource1(i)=0)&(dsource2(i)=1)&(dsource3(i)=1),r=s011+n;elseif(dsource1(i)=1)&(dsource2(i)=0)&(dsource3(i)=0),r=s100+n;elseif(dsource1(i)=1)&(dsource2(i)=0)&(dsource3(i)=1),r=s101+n;elseif(dsource1(i)=1)&(dsource2(i)=1)&(dsource3(i)=0),r=s110+n;elser=s111+n;end;c000=dot(r,s000);c001=dot(r,s001);c010=dot(r,s010);c011=dot(r,s011);c100=dot(r,s100);c101=dot(r,s101);c110=dot(r,s110);c111=dot(r,s111);c_max=max(c000 c001 c010 c011 c100 c101 c110 c111);if(c000=c_max),decis1=0;decis2=0;decis3=0;elseif(c001=c_max),decis1=0;decis2=0;decis3=1;elseif(c010=c_max),decis1=0;decis2=1;decis3=0;elseif(c011=c_max),decis1=0;decis2=1;decis3=1;elseif(c100=c_max),decis1=1;decis2=0;decis3=0;elseif(c101=c_max),decis1=1;decis2=0;decis3=1;elseif(c110=c_max),decis1=1;decis2=1;decis3=0;elsedecis1=1;decis2=1;decis3=1;end;symbolerror=0;if(decis1=dsource1(i),numofbiterror=numofbiterror+1;symbolerror=1;end;if(decis2=dsource2(i),numofbiterror=numofbiterror+1;symbolerror=1;end;if(decis3=dsource3(i),numofbiterror=numofbiterror+1;symbolerror=1;end;if(symbolerror=1),numofsymbolerror=numofsymbolerror+1;end;end;ps=numofsymbolerror/N;pb=numofbiterror/(3*N);4PSK的源程序：% 對(duì)M=4的PSK通信系統(tǒng)進(jìn)行蒙特卡羅仿真 %echo onSNRindB1=0:2:10;SNRindB2=0:0.1:10;for i=1:length (SNRindB1),pb , ps=cm_sm32 (SNRindB1 (i); %仿真比特和符號(hào)誤碼率smld_bit_err_prb (i)=pb;smld_symbol_err_prb (i)=ps;end;for i=1:length (SNRindB2),SNR=exp(SNRindB2(i)*log(10)/10); %信噪比theo_err_prb(i)=Qfunct(sqrt(2*SNR); %理論比特誤碼率end;%隨