升余弦滾降系統的源程序

1、升余弦滾降系統的源程序：(一)，用到的函數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);

2、60;  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(時域取樣間隔) &

3、#160; H=fft(x);   X=H(N/2+1:N),H(1:N/2)*dt;   end(二),主程序   global dt t df N   close all   k=input('取樣點數2k, k取13左右');   if isempty(k), k=13; end   z=input('每個信號取樣點數2z, z<k'

4、);   if isempty(z), z=5; end   aa=input('滾降系數=0.5');   if aa=,aa=0.5;end   N=2k;   L=2z;M=N/L;   Na=4;                

5、0;         %示波器掃描寬度為4個碼元   Rb=2;                         %碼速率是2Mb/s   Ts=1/Rb;   

6、0;                   %碼元間隔   dt=Ts/L;                       %時域采樣間隔 &#

7、160; df=1/(N*dt);                   %頻域采樣間隔   T=N*dt;                      

8、  %截短時間   Bs=N*df/2;                     %系統帶寬   f=-Bs+df/2:df:Bs;            %頻域橫坐標  

9、; 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:30   b=sign(

10、randn(1,M)+1;          %發(fā)送碼元取值為0、2   s=zeros(1,N);                  %產生沖激序列   s(L/2:L:N)=b;   SS=t2f(s);  

11、; S=SS.*GG;                      %升余弦信號的傅氏變換   a=f2t(S);                   

12、   %升余弦信號的時域波形   a=real(a);   P=S.*conj(S)/T;                %升余弦信號的功率譜   EP=(EP*(ii-1)+P+eps)/ii;   figure(1)      

13、;                %輸出的功率畫圖   plot(f,10*log10(EP),'b');   grid   axis(-3,+3,-100,20)   xlabel('f (KHz)');ylabel('Ps(f) (dBm/MHz)') ;title('功率譜圖')   figure(2)   hold on   tt=0:dt:Na*L*dt;