正弦波的源程序.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函数。代码(二),主程序global dt t df Nclose allk=input(取样点数2k, k取13左右);if isempty(k), k=13; endz=input(每个信号取样点数2z, zk);if isempty(z), z=5; endaa=input(滚降系数=0.5);if aa=,aa=0.5;endN=2k;L=2z;M=N/L;Na=4; %示波器扫描宽度为4个码元Rb=2; %码速率是2Mb/sTs=1/Rb; %码元间隔dt=Ts/L; %时域采样间隔df=1/(N*dt); %频域采样间隔T=N*dt; %截短时间Bs=N*df/2; %系统带宽f=-Bs+df/2:df:Bs; %频域横坐标t=-T/2+dt/2:dt:T/2; %时域横坐标g=sin(pi*t/Ts).*cos(pi*t*aa/Ts)./pi*t/Ts.*(1-4*t.2*aa2/Ts2);%升余弦脉冲波形GG=t2f(g); EP=zeros(size(f)+eps;for ii=1:30b=sign(randn(1,M)+1; %发送码元取值为0、2s=zeros(1,N); %产生冲激序列s(L/2:L:N)=b;SS=t2f(s);S=SS.*GG; %升余弦信号的傅氏变换a=f2t(S); %升余弦信号的时域波形a=real(a);P=S.*conj(S)/T; %升余弦信号的功率谱EP=(EP*(ii-1)+P+eps)/ii;figure(1) %输出的功率画图plot(f,10*log10(EP),b);gridaxis(-3,+3,-100,20)xlabel(f (KHz);ylabel(Ps(f) (dBm/MHz) ;title(功率谱图)figure(2)hold ontt=0:dt:Na*L*dt; %系统的眼图画图 for jj=1:Na*L:N-Na*Lplot(tt,a(jj:jj+Na*L);title(升余弦滚降波形的眼图);endend8psk的源程序：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(temp(二),主程序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);g