function [g] = G_t( t,a,Ts)%UNTITLED1 Summary of this function goes here%  Detailed explanation goes here%  a為滾降系數%Ts =1; %receive signal symbol %ycx本函數輸入t是一段離散時間量，a為一個滾降系數，接收Ts采樣？？ 輸出RRC g為離散的根余弦 %******************************Rising cosine pulse************************************% t1=pi*t/T;% g=sin(t1)./t1.*cos(a*t1)./(1-4*a^2*t.^2/T^2);   %******************************Root Rising cosine pulse**RRC!******************************% for u =1:length(t)%     if t(u)==0%         g(u)=(1+(4/pi-1)*a)/Ts;              % ycx  與t無關！！%     else if t(u)==1/(4*a)|| t(u)==-1/(4*a)%             g(u)=sqrt(2)*a*((pi-2)*cos(pi/(4*a))+(pi+2)*sin(pi/(4*a)))/(2*pi);  %ycx  與t無關！！   %         else%             g(u) = 1/(pi*t(u)*(1-(4*a*t(u)/Ts)^2))*(4*a*t(u)/Ts*cos(pi*(1+a)*t(u)/Ts)+sin(pi*(1-a)*t(u)/Ts));%         end;%     end;% end;for i = 1:length(t)    if t(i) == 0        g(i) = ( 1 + (4/pi-1)*a )/Ts ;    else if(abs(abs(t(i))-Ts/(4*a))<1e-20)          % ycx ???            g(i) = sqrt(2)*a/(2*pi*Ts)*( (pi-2)*cos(pi/(4*a)) + (pi+2)*sin(pi/(4*a)) ) ;         else             g(i) = ( 4*a*t(i)/Ts*cos( pi*(1+a)*t(i)/Ts)+ sin(pi*(1-a)*t(i)/Ts) )/(pi*t(i)*(1-(4*a*t(i)/Ts)^2) ) ;         end     end end