hn.m

通信仿真： 利用Matlab实现的脉冲成型仿真

T=1;
a=0.22;
N=32;
for n=1:N
    t(n)=(n-N/2)*T/4;
    y(n)=t(n)/T;
    if(t(n)~=0)
     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
figure
stem(t,h)
max=0;min=0;
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);
        if s(i*4+j+1)>max
            max=s(i*4+j+1);
        end
        
        if s(i*4+j+1)<min
            min=s(i*4+j+1);
        end 
        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 
sss=floor((s-min)*255/(max-min));
f=fopen('text.txt','wt');
fprintf(f,'    %d :  %d;\n',[0:1023;sss]);
fclose(f);

yy=[ss(77,1:4) ss(39,1:4) ss(20,1:4) ss(138,1:4) ss(197,1:4) ss(99,1:4) ss(50,1:4) ss(153,1:4) ss(77,1:4) ss(39,1:4) ss(20,1:4) ss(138,1:4) ss(197,1:4) ss(99,1:4) ss(50,1:4) ss(153,1:4)]
%yy=[ss(1,1:4) ss(3,1:4) ss(5,1:4) ss(9,1:4) ss(17,1:4) ss(33,1:4) ss(65,1:4) ss(129,1:4)  ss(1,1:4) ss(3,1:4) ss(5,1:4) ss(9,1:4) ss(17,1:4) ss(33,1:4) ss(65,1:4) ss(129,1:4)];
figure
plot(yy)