Matlab源程序代码.doc

正弦波的源程序：(一)，用到的函数1,f2t函数function x=f2t(X)global dt df t f T N%x=f2t(X)%x为时域的取样值矢量%X为x的傅氏变换%X与x长度相同并为2的整幂%本函数需要一个全局变量dt(时域取样间隔)X=X(N/2+1:N),X(1:N/2);x=ifft(X)/dt;end 2,t2f函数。function X=t2f(x)global dt df N t f T%X=t2f(x)%x为时域的取样值矢量%X为x的傅氏变换%X与x长度相同,并为2的整幂。%本函数需要一个全局变量dt(时域取样间隔)H=fft(x);X=H(N/2+1:N),H(1:N/2)*dt;end(二),主程序。1，%（1）绘出正弦信号波形及频谱global dt df t f Nclose allk=input(取样点数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; %截短时间Bs=N*df/2; %系统带宽f=-Bs+df/2:df:Bs; %频域横坐标t=-T/2+dt/2:dt:T/2; %时域横坐标s=sin(2*pi*f0*t); %输入的正弦信号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) %输出信号波形画图gridaxis(-2/f0,+2/f0,-1.5,1.5)xlabel(t(ms)ylabel(a(t)(V)gtext(频谱图)最佳基带系统的源程序：(一)，用到的函数f2t函数和t2f函数。代码(二),主程序globaldt t f df N Tclose allclear Eb_N0 Pek=input(取样点数2k, k取13左右);if isempty(k), k=13; endz=input(每个信号取样点数2z, z(二),主程序global dt t f df N Tclose allclear Eb_N0 Pek=input(取样点数2k, k取13左右);if isempty(k), k=13; endz=input(每个信号取样点数2z, z(二),主程序global dt t f df N Tclose allk=input(取样点数2k, k取10左右);if isempty(k), k=13; endz=input(每个信号取样点数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; %时域采样间隔df=1/(N*dt) %MHzT=N*dt %截短时间Bs=N*df/2 %系统带宽f=-Bs+df/2:df:Bs; %频域横坐标t=-T/2+dt/2:dt:T/2; %时域横坐标EP=zeros(size(f)+eps;for ii1=1:30b=round(rand(1,M); %产生随机码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) %输出信号波形画图plot(t,a,b)gridaxis(-3,3,0,1.3);xlabel(t (ms);ylabel(s(t) (V);title(单极性RZ的时域图)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的源程序：% 对M=4的PSK通信系统进行蒙特卡罗仿真 %echo onSNRindB1=0:2:10;SNRindB2=0:0.1:10;for i=1:length (SNRindB1),pb , ps=cm_sm32 (SNRindB1 (i); %仿真比特和符号误码率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;%随