hn.asv

来自「通信仿真： 利用Matlab实现的脉冲成型仿真」· ASV 代码 · 共 50 行

T=1;
a=0.3;
N=32;
for n=1:N
    t(n)=(n-N/2-1)*T/4;
    y(n)=t(n)/T;
    if(t(n)~=0)
    %temp_sin=sin(pi*(1-a)*y(n));
    %temp_cos=4*a*y(n)*cos(pi*(1-a)*y(n));
    %h(n)=(temp_sin+temp_cos)/(1-4*a*y(n))/(pi*y(n));
     temp_sin=sin(pi*t(n));
     temp_cos=cos(a*pi*t(n));
     temp1=pi*t(n);
     temp2=1-4*a*a*t(n)*t(n);
     h(n)=temp_sin*temp_cos/temp1/temp2;
    else 
       h(n)=1; 
    end
end
stem(t,h)

for i=0:255
    tt=256;
    ttt=i;
    for ii=1:8
        tt=tt/2;
        d(ii)=floor(ttt/tt);
        ttt=ttt-tt*d(ii);
    end   
    for j=0:3
        s(i*4+j+1)=d(8)*h(29+j)+d(7)*h(25+j)+d(6)*h(21+j)+d(5)*h(17+j)+d(4)*h(13+j)+d(3)*h(9+j)+d(2)*h(5+j)+d(1)*h(1+j);
        ss(i+1,j+1)=d(8)*h(29+j)+d(7)*h(25+j)+d(6)*h(21+j)+d(5)*h(17+j)+d(4)*h(13+j)+d(3)*h(9+j)+d(2)*h(5+j)+d(1)*h(1+j);
    end
end 
%yy=[ss(76,1:4) ss(38,1:4) ss(19,1:4) ss(137,1:4) ss(196,1:4) ss(98,1:4) ss(49,1:4) ss(152,1:4)];
yy=[ss(125,1:4) ss(38,1:4) ss(19,1:4) ss(137,1:4) ss(196,1:4) ss(98,1:4) ss(49,1:4) ss(152,1:4)];

plot(yy)